上机练习专题6：数学物理方程的计算机求解和可视化.ppt

数学物理建模与计算机辅助设计;Page *;达朗贝尔解的物理意义;实例1：一维无界弦振动定解问题:;clear u0=2; u(1:560)=0; x1=-1;x2=1; x=linspace(-2,2,560); u(141:280)=2*u0*(x(141:280)-x1)/(x2-x1); u(281:420)=2*u0*(x2-x(281:420))/(x2-x1); uu=u; h=plot(x,u,linewidth,3); axis([2*x1,2*x2,-4,4]); set(h,EraseMode,xor) for at=2:140 lu(1:560)=0; ru(1:560)=0; lx=[141:420]-at; rx=[141:420]+at; lu(lx)=0.5*uu(141:420); ru(rx)=0.5*uu(141:420); u=lu+ru; set(h,Xdata,x,YData,u); drawnow; pause(0.1) end ;某时刻无界弦振动图型;实例2：一维无界弦振动定解问题:;其中: 　　　　;clear x1=0;x2=1;t=0:0.005:8;x=-10:0.1:10;a=1; [X,T]=meshgrid(x,t); xpat=X+a*T; xpat(find(xpat=x1))=0;xpat(find(xpat=x2))=1; xmat=X-a*T; xmat(find(xmat=x1))=0;xmat(find(xmat=x2))=1; jf1=1/2/a*(xpat-x1);jf2=1/2/a*(xmat-x1);jf3=1/2/a*((xpat-x1)-(xmat-x1)); subplot(3,1,1) h1=plot(x,jf1(1,:),linewidth,3); set(h1,erasemode,xor); axis([-10 10 -1 1]) hold on subplot(3,1,2) h2=plot(x,jf2(1,:),linewidth,3); set(h2,erasemode,xor); axis([-10 10 -1 1]) hold on subplot(3,1,3) h3=plot(x,jf3(1,:),linewidth,3); set(h3,erasemode,xor); axis([-10 10 -1 1]) hold on for j=2:length(t) pause(0.01) set(h1,ydata,jf1(j,:)); set(h2,ydata,jf2(j,:)); set(h3,ydata,jf3(j,:)); drawnow; end ;左行波: ;利用分离变量法求解数学物理定解问题;特解改写为： ;振幅: ;function zb t=0:0.005:2.0; x=0:0.001:1; ww1=wfun(1,0);ww2=wfun(2,0); ww3=wfun(3,0);ww4=wfun(4,0); ymax1=max(abs(ww1)); figure subplot(4,1,1) h1=plot(x,ww1,r,linewidth,5); axis([0,1,-ymax1,ymax1]) subplot(4,1,2) ymax2=max(abs(ww2)); h2=plot(x,ww2,r,linewidth,5); axis([0,1,-ymax2,ymax2]) subplot(4,1,3) ymax3=max(abs(ww3)); h3=plot(x,ww3,r,linewidth,5); axis([0,1,-ymax3,ymax3]) subplot(4,1,4) ymax4=max(abs(ww4)); h4=plot(x,ww4,r,linewidth,5); axis([0,1,-ymax4,ymax4]) ;n=1;Page *;Page *;Page *;Page *;Page *;比较下面两副图的区别： figure(4) subplot(2,2,1); [C,h]=contour(X,Y,w11,20); clabel(C,h,manual); subplot(2,2,2); [C,h]=contour(X,Y,w12,20); clabel(C,h,manual); subplot(2,2,3); [C,h]=contour(X,Y,w21,20); clabel(C,h,manual); subplot(2,2,4); [C,h]=contour(X,Y,