求助!看看这个程序有什么问题?运行不出结果!!

bhy521

1. clear
   
2. h=1;
   
3. v0=100;
   
4. hx=16;
   
5. hy=10;
   
6. v1=zeros(hy,hx);
   
7. v1(hy,:)=ones(1,hx)*v0;
   
8. v1(2:hy-1,2:hx-1)=ones(hy-2,hx-2)*v0/2;
   
9. v2=zeros(hy,hx);
   
10. for i=1:hy
    
11.     for j=1:hx
    
12.         if i==hy
    
13.            v2(i,j)=v0;
    
14.         elseif i==1|j==1|j==hx
    
15.             v2(i,j)=0;
    
16.         else
    
17.             v2(i,j)=[v1(i,j-1)+v1(i,j+1)+v1(i-1,j)+v1(i+1,j)]/4;
    
18.         end
    
19.     end
    
20. end
    
21. for k=1:500
    
22.     v1=v2;
    
23.   for i=1:hy
    
24.     for j=1:hx
    
25.         if i==hy
    
26.             v2(i,j)=0;
    
27.         elseif i==1|j==1|j==hx
    
28.             v2(i,j)=0;
    
29.         else
    
30.             v2(i,j)=[v1(i,j-1)+v1(i,j+1)+v1(i-1,j)+v1(i+1,j)]/4;
    
31.         end
    
32.    
    
33.   end
    
34. end
    
35. subplot(1,2,1),mesh(v2)
    
36. h1=mesh(v2);
    
37. set(h1,'linewidth',1,'edgecolor','b')
    
38. text(-6,10,120,'v','fontsize',10);
    
39. text(-4,12,0,'h','fontsize',10);
    
40. text(18,0,0,'b','fontsize',10);
    
41. hold on
    
42. axis([0 17 0 10 0 100]);
    
43. subplot(1,2,2),
    
44. contour(v2,'k-')
    
45. text(-2,10,'h','fontsize',10);
    
46. text(17,1,'b','fontsize',10);
    
47. hold on

复制代码

[ 本帖最后由 lxq 于 2006-11-22 21:31 编辑 ]

suffer

少了一个end，给改如下

1. clear
   
2. h=1;
   
3. v0=100;
   
4. hx=16;
   
5. hy=10;
   
6. v1=zeros(hy,hx);
   
7. v1(hy,:)=ones(1,hx)*v0;
   
8. v1(2:hy-1,2:hx-1)=ones(hy-2,hx-2)*v0/2;
   
9. v2=zeros(hy,hx);
   
10. for i=1:hy
    
11.     for j=1:hx
    
12.         if i==hy
    
13.            v2(i,j)=v0;
    
14.         elseif i==1|j==1|j==hx
    
15.             v2(i,j)=0;
    
16.         else
    
17.             v2(i,j)=[v1(i,j-1)+v1(i,j+1)+v1(i-1,j)+v1(i+1,j)]/4;
    
18.         end
    
19.     end
    
20. end
    
21. for k=1:500
    
22.     v1=v2;
    
23.   for i=1:hy
    
24.     for j=1:hx
    
25.         if i==hy
    
26.             v2(i,j)=0;
    
27.         elseif i==1|j==1|j==hx
    
28.             v2(i,j)=0;
    
29.         else
    
30.             v2(i,j)=[v1(i,j-1)+v1(i,j+1)+v1(i-1,j)+v1(i+1,j)]/4;
    
31.         end
    
32.     end
    
33.   end
    
34. end
    
35. subplot(1,2,1),mesh(v2)
    
36. h1=mesh(v2);
    
37. set(h1,'linewidth',1,'edgecolor','b')
    
38. text(-6,10,120,'v','fontsize',10);
    
39. text(-4,12,0,'h','fontsize',10);
    
40. text(18,0,0,'b','fontsize',10);
    
41. hold on
    
42. axis([0 17 0 10 0 100]);
    
43. subplot(1,2,2),
    
44. contour(v2,'k-')
    
45. text(-2,10,'h','fontsize',10);
    
46. text(17,1,'b','fontsize',10);
    
47. hold on

复制代码

yangji

这都能看出来，厉害