Difference between revisions of "Team:Tokyo Tech/Collaborations code"

 
(2 intermediate revisions by the same user not shown)
Line 5: Line 5:
 
   <meta name="generator" content="MediaWiki 1.16.5" />
 
   <meta name="generator" content="MediaWiki 1.16.5" />
 
   <style type="text/css">
 
   <style type="text/css">
 +
p{
 +
  margin-left: 50px;
 +
}
 
   </style type>
 
   </style type>
 
  </head>
 
  </head>
Line 19: Line 22:
 
　　<div class="textarea">
 
　　<div class="textarea">
 
         <h3 class="sub5">Matlab Code</h3>
 
         <h3 class="sub5">Matlab Code</h3>
 +
<p>
  
  
−
ふーせさんコードをよろしく。
+
%% code<br>
 +
X=52; %1=2mm<br>
 +
Y=52;<br>
 +
Z=21;<br>
 +
<br>
 +
T=60*12; %24h<br>
 +
%clear value of u<br>
 +
u=zeros(T,X,Y,Z);<br>
 +
%diffusion coefficient<br>
 +
D=0.027*4/5; <br>
 +
 
 +
<br>
 +
for i=6:1:46<br>
 +
for j=6:1:46<br>
 +
<br>
 +
    u(1,i,j,1)=8*182.5e-02;<br>
 +
end<br>
 +
end<br>
 +
<br>
 +
<br>
 +
 
 +
for i=1:X;<br>
 +
for k=1:Z;<br>
 +
    u(1,i,1,k)=0;<br>
 +
    u(1,i,Y,k)=0;<br>
 +
end<br>
 +
end<br>
 +
for j=1:Y;<br>
 +
for k=1:Z;<br>
 +
    u(1,1,j,k)=0;<br>
 +
    u(1,X,j,k)=0;<br>
 +
end<br>
 +
end<br>
 +
<br>
 +
for t=2:1:T<br>
 +
for i=1:X;<br>
 +
  for j=1:Y;<br>
 +
  for k=1:Z;<br>
 +
    u(t,1,j,k)=0;<br>
 +
    u(t,i,1,k)=0;<br>
 +
    u(t,X,j,k)=0;<br>
 +
    u(t,i,Y,k)=0;<br>
 +
  end<br>
 +
  end<br>
 +
end<br>
 +
<br>
 +
for i=2:1:X-1;<br>
 +
  for j=2:1:Y-1;<br>
 +
  for k=1:Z;<br>
 +
    if k==1;<br>
 +
    u(t,i,j,k)=u(t-1,i,j,k)+D*(u(t-1,i+1,j,k)-2*u(t-1,i,j,k)+u(t-1,i-1,j,k)+u(t-1,i,j+1,k)-2*u(t-1,i,j,k)+u(t-1,i,j-1,k)+u(t-1,i,j,k+1)-2*u(t-1,i,j,k)+0);<br>
 +
    else<br>
 +
    if k==Z;<br>
 +
        u(t,i,j,k)=u(t-1,i,j,k)+D*(u(t-1,i+1,j,k)-2*u(t-1,i,j,k)+u(t-1,i-1,j,k)+u(t-1,i,j+1,k)-2*u(t-1,i,j,k)+u(t-1,i,j-1,k)+0-2*u(t-1,i,j,k)+u(t-1,i,j,k-1));<br>
 +
    else<br>
 +
        u(t,i,j,k)=u(t-1,i,j,k)+D*(u(t-1,i+1,j,k)-2*u(t-1,i,j,k)+u(t-1,i-1,j,k)+u(t-1,i,j+1,k)-2*u(t-1,i,j,k)+u(t-1,i,j-1,k)+u(t-1,i,j,k+1)-2*u(t-1,i,j,k)+u(t-1,i,j,k-1));<br>
 +
    end<br>
 +
    end<br>
 +
  end<br>
 +
  end<br>
 +
end<br>
 +
end<br>
 +
<br>
 +
figure<br>
 +
<br>
 +
w=zeros(X,Y);<br>
 +
v=zeros(X,Y);<br>
 +
 
 +
for x=1:X;<br>
 +
for y=1:Y;<br>
 +
    w(x,y)=u(T,x,y,1);<br>
 +
    <br>
 +
end<br>
 +
end<br>
 +
figure(1)<br>
 +
pcolor(w);<br>
 +
<br>
 +
<br>
 +
for z=1:Z;<br>
 +
for y=1:Y;<br>
 +
    v(y,z)=u(T,26,y,z);<br>
 +
    <br>
 +
end<br>
 +
end<br>
 +
figure(2)<br>
 +
pcolor(v)<br>
 +
<br>
 +
for x=1:X;<br>
 +
for y=1:Y;<br>
 +
    ww(x,y)=u(T,x,y,Z);<br>
 +
end<br>
 +
end<br>
 +
<br>
 +
figure(3)<br>
 +
pcolor(ww)<br>
 +
<br>
 +
<br>
 +
</p>
 
     </div>
 
     </div>
 
   　<div class="textbottom">
 
   　<div class="textbottom">

Latest revision as of 08:59, 6 October 2015

Collaborations

We were helped Nagahama!!

　　

Matlab Code

%% code
X=52; %1=2mm
Y=52;
Z=21;

T=60*12; %24h
%clear value of u
u=zeros(T,X,Y,Z);
%diffusion coefficient
D=0.027*4/5;

for i=6:1:46
for j=6:1:46

u(1,i,j,1)=8*182.5e-02;
end
end


for i=1:X;
for k=1:Z;
u(1,i,1,k)=0;
u(1,i,Y,k)=0;
end
end
for j=1:Y;
for k=1:Z;
u(1,1,j,k)=0;
u(1,X,j,k)=0;
end
end

for t=2:1:T
for i=1:X;
for j=1:Y;
for k=1:Z;
u(t,1,j,k)=0;
u(t,i,1,k)=0;
u(t,X,j,k)=0;
u(t,i,Y,k)=0;
end
end
end

for i=2:1:X-1;
for j=2:1:Y-1;
for k=1:Z;
if k==1;
u(t,i,j,k)=u(t-1,i,j,k)+D*(u(t-1,i+1,j,k)-2*u(t-1,i,j,k)+u(t-1,i-1,j,k)+u(t-1,i,j+1,k)-2*u(t-1,i,j,k)+u(t-1,i,j-1,k)+u(t-1,i,j,k+1)-2*u(t-1,i,j,k)+0);
else
if k==Z;
u(t,i,j,k)=u(t-1,i,j,k)+D*(u(t-1,i+1,j,k)-2*u(t-1,i,j,k)+u(t-1,i-1,j,k)+u(t-1,i,j+1,k)-2*u(t-1,i,j,k)+u(t-1,i,j-1,k)+0-2*u(t-1,i,j,k)+u(t-1,i,j,k-1));
else
u(t,i,j,k)=u(t-1,i,j,k)+D*(u(t-1,i+1,j,k)-2*u(t-1,i,j,k)+u(t-1,i-1,j,k)+u(t-1,i,j+1,k)-2*u(t-1,i,j,k)+u(t-1,i,j-1,k)+u(t-1,i,j,k+1)-2*u(t-1,i,j,k)+u(t-1,i,j,k-1));
end
end
end
end
end
end

figure

w=zeros(X,Y);
v=zeros(X,Y);
for x=1:X;
for y=1:Y;
w(x,y)=u(T,x,y,1);

end
end
figure(1)
pcolor(w);


for z=1:Z;
for y=1:Y;
v(y,z)=u(T,26,y,z);

end
end
figure(2)
pcolor(v)

for x=1:X;
for y=1:Y;
ww(x,y)=u(T,x,y,Z);
end
end

figure(3)
pcolor(ww)