各位大神,麻烦帮我看看我的这个MATLAB为什么画不出图来,谢谢了!!!!
lx=51;ly=26;%把取b=pi/2,则m=50,n=25,h=l=pi/50.v1=0;v1(:,1)=0;v1(:,51)=0;forj=2:50v1(ly,j...
lx=51;ly=26;%把取b=pi/2,则m=50,n=25,h=l=pi/50.
v1=0;
v1(:,1)=0;v1(:,51)=0;
for j=2:50
v1(ly,j)=(5*(10^-2).*(j*pi/50).*(4-(j*pi/50)))./0.1249;
end
v2=v1;
maxt=1;t=0;
k=0;
while(maxt>1e-6)
k=k+1;
maxt=0;
for j=2:lx-1;
v1(1,j)=(v1(2,j)+2.*v1(3,j)-v1(4,j))./2;
v1(1,j)=sin((j*pi)/50);
v2(16,j)=(v1(16,j+1)+v1(17,j)+v1(15,j)+v1(16,j-1))/4;
t=abs(v2(16,j)-v1(16,j));
if(t>maxt)
maxt=t;
%end
end
end
v1=v2;
end
subplot(1,2,1),plot(j,v2)
axis([0,60,0,5]) 展开
v1=0;
v1(:,1)=0;v1(:,51)=0;
for j=2:50
v1(ly,j)=(5*(10^-2).*(j*pi/50).*(4-(j*pi/50)))./0.1249;
end
v2=v1;
maxt=1;t=0;
k=0;
while(maxt>1e-6)
k=k+1;
maxt=0;
for j=2:lx-1;
v1(1,j)=(v1(2,j)+2.*v1(3,j)-v1(4,j))./2;
v1(1,j)=sin((j*pi)/50);
v2(16,j)=(v1(16,j+1)+v1(17,j)+v1(15,j)+v1(16,j-1))/4;
t=abs(v2(16,j)-v1(16,j));
if(t>maxt)
maxt=t;
%end
end
end
v1=v2;
end
subplot(1,2,1),plot(j,v2)
axis([0,60,0,5]) 展开
1个回答
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询