《数学物理方法》补充教材一-弦振动方程的MATLAB仿真结果.docVIP

计算机仿真在数理方程中的应用 一 齐次波动问题 1 作出有界弦振动问题解析解的图形分布 ， 解析解：（取系数为1） 1）n=1本振解随时间变化的图形分布 clear a=1;l=1;x=0:0.05:1;t=0:0.001:3;u=0; n=1; [X,T]=meshgrid(x,t); u=(cos(n*pi*a*T/l)+sin(n*pi*a*T/l)).*sin(n*pi*X/l) figure(1) axis([0 1 -0.05 0.05]) mesh(X,T,u) title(n=1本振随时分布的透视图) xlabel(x) ylabel(t) zlabel(u) （图形输出） 2）n=2本振解随时间变化的图形分布 clear a=1;l=1;x=0:0.05:1;t=0:0.001:3;u=0; n=2; [X,T]=meshgrid(x,t); u=(cos(n*pi*a*T/l)+sin(n*pi*a*T/l)).*sin(n*pi*X/l) figure(1) axis([0 1 -0.05 0.05]) mesh(X,T,u) title(n=2本振随时分布的透视图) xlabel(x) ylabel(t) zlabel(u) （图形输出） 3）前10个本征振动迭加解随时分布图 clear a=1;l=1;x=0:0.05:1;t=0:0.001:3;u=0; [X,T]=meshgrid(x,t); for n=1:10; u0=(cos(n*pi*a*T/l)+sin(n*pi*a*T/l)).*sin(n*pi*X/l) u=u+u0 end figure(1) axis([0 1 -0.05 0.05]) mesh(X,T,u) title(前10个本振迭加解随时分布的透视图) xlabel(x) ylabel(t) zlabel(u) （仿真结果） 2 作出磁致伸缩换能器作纵振动的定解问题解的图形分布（P185例1） ， 解析解：（取系数为1） 1）n=1本振解随时间变化的图形分布 clear a=1;c=1;d=2;l=1;x=0:0.05:1;t=0:0.001:3;u=0; n=1; [X,T]=meshgrid(x,t); u0=(cos(n*pi*a*T/l)+sin(n*pi*a*T/l)).*cos(n*pi*X/l) u= u0+c+d*T figure(1) axis([0 1 -0.05 0.05]) mesh(X,T,u) title(n=1本振随时分布的透视图) xlabel(x) ylabel(t) zlabel(u) （图形输出） 2）n=2本振解随时间变化的图形分布 clear a=1;c=1;d=2;l=1;x=0:0.05:1;t=0:0.001:3;u=0; n=2; [X,T]=meshgrid(x,t); u0=(cos(n*pi*a*T/l)+sin(n*pi*a*T/l)).*cos(n*pi*X/l) u= u0+c+d*T figure(1) axis([0 1 -0.05 0.05]) mesh(X,T,u) title(n=2本振随时分布的透视图) xlabel(x) ylabel(t) zlabel(u) （图形输出） 3）前10个本征振动迭加解随时分布图 clear a=1;c=1;d=2;l=1;x=0:0.05:1;t=0:0.001:3;u1=0; [X,T]=meshgrid(x,t); for n=1:9; u0=(cos(n*pi*a*T/l)+sin(n*pi*a*T/l)).*sin(n*pi*X/l) u1=u1+u0 end u=c+d*T+u1 figure(1) axis([0 1 -0.05 0.05]) mesh(X,T,u) title(前10个本振迭加解随时分布的透视图) xlabel(x) ylabel(t) zlabel(u) （仿真结果） 二 齐次输运问题 1 作出热传导方程解的图形（P188，例2） ， 解析解： 1）k=0本征解随t变化的图形 a=1;l=40;u0=10;k=0; x=0:0.5:40; t=0.0001:0.00001:0.005; [X,T]=meshgrid(x,t) u=(2*u0/(pi^2))*((-1)^k)/((k+1/2)^2)*exp(-((k+1/2)^2*pi^2*a^2/1600)*T).*sin((k+1/2)*pi*X/l) mesh(X,T,u) xlabel(x) ylabel(t) zlabel(u) 2）k=2本征解随t变化的图形 a=1;l=40;u0=10;k=2; x=0:0.5:40; t=0.0001