二维驻波的波动方程仿真

jainxing

关于图片上的公式，我用下面的代码做对吗？

1. m=2; %设置m，n值
   
2. n=3;
   
3. t=0;
   
4. a=4;  %设置板大小
   
5. b=4;
   
6. x=0:.1:a;  %节点
   
7. y=0:.1:b;
   
8. [X,Y]=meshgrid(x,y);  %X，Y对应x，y产生矩阵。
   
9. Modal_shape=4*sin(m*pi*X/a).*sin(n*pi*Y/b).*cos(sqrt((m*pi/a)^2+(n*pi/b)^2)*t);  %模态形状的表示
   
10. for ii=1:2  %循环次数
    
11.     for i=-100:2:100
    
12.         surf(X,Y,0.01*i*Modal_shape)  %作图
    
13.         axis([0 a 0 b -4 4])  %坐标轴大小
    
14.         xlabel('x','FontSize',11);  %设置x轴标签
    
15.         ylabel('y');  %设置y轴标签
    
16.         zlabel('z');  %设置z轴标签
    
17.         title('驻波图(2,3)');  %设置标题
    
18.         pause(0.001)  %暂停程序执行的秒数
    
19.     end
    
20. end
    

复制代码

Ryan2017

wawawawawawaw

jainxing

Ryan2017 发表于 2018-2-4 09:29
wawawawawawaw

你也是要做这个吗？