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

数学物理建模与计算机辅助设计 数学物理建模与计算机辅助设计 上机练习专题5：数学物理方程的计算机求解和可视化 Page * 对代数表达式进行可视化 二维本征值问题 矩形区域的本征模与本征振动 边长为b和c的的四周固定的矩形膜振动的本征值问题为 采用分离变量法可以得到本征模和本征值为 Page * 对代数表达式进行可视化 二维本征值问题 绘制前4个本征函数的图形 %P70_1.m a=2; b=1; [m,n]=meshgrid(1:3); L=((m*pi./b).^2+(n*pi./b).^2); x=0:0.01:a; y=0:0.01:b; [X,Y]=meshgrid(x,y); w11=sin(pi*Y./b).*sin(pi*X./a); w12=sin(2*pi*Y./b).*sin(pi*X./a); w21=sin(pi*Y./b).*sin(2*pi*X./a);w22=sin(pi*Y./b).*sin(3*pi*X./a); figure subplot(2,2,1); mesh(X,Y,w11); subplot(2,2,2); mesh(X,Y,w12); subplot(2,2,3); mesh(X,Y,w21); subplot(2,2,4); mesh(X,Y,w22); Page * 对代数表达式进行可视化 二维本征值问题 动态可视化：设具有时间因子 %P71_1.m b=2; c=1; x=0:0.02:b;y=0:0.02:c; [X,Y]=meshgrid(x,y); Z=zeros(51,51); p=moviein(2*3*60); for m=1:2 for n=1:3 for i=1:60 a=sqrt((m*pi/c).^2+(n*pi/b).^2); Z=sin(a*i*.02*pi)*sin(m*pi*Y./c).*sin(n*pi*X./b); mesh(X,Y,Z); t=[本征振动：,m=,int2str(m), n=,int2str(n)]; title(t); axis([0 b 0 c -1 1]); p(:,((m-1)*3+(n-1))*60+i)=getframe; end end end MOVIE2AVI(p,D:\A.avi) Page * PDE工具箱求解数学物理方程 PDE：Partial Differential Equation PDE工具箱功能 设定二维定解区域、边界条件以及方程的形式和系数 用有限元法(FEM)生成网格、方程离散化并求出数值解 对所求解的可视化，图像和动画 PDE工具箱可以求解的4种偏微分方程 椭圆型方程 抛物型方程 双曲形方程 特征值方程 Page * PDE工具箱求解数学物理方程 PDE工具箱给定的两种边界条件 Dilichlet条件 Neumann条件 用PDE工具箱求解二维本征值问题 Page * PDE工具箱求解数学物理方程 求解二维本征值问题 PDE工具箱求解数学物理方程 1.求解定解问题 Page * 取 PDE工具箱求解数学物理方程 2.求矩形薄膜的横向振动。定解问题为 Page * 取 Page * 差分法求解数学物理方程 导数的差分公式 在x附近将函数f(x)展开成泰勒公式 前差公式 后差公式 二阶导数的差分公式(中心差分公式) Page * 差分法求解数学物理方程 用数值方法求解微分方程的时候步骤： (1)变量空间网格化 (2)微分方程变差分方程 (3)定解条件设置 (4)整理得到递推公式 Page * 差分法求解数学物理方程 一维波动问题 两端固定的弦的振动问题 设定解条件为 解的代数表达式为 Page * 差分法求解数学物理方程 两端固定的弦的振动问题的求解步骤 (1)变量空间网格化 (2)微分方程变差分方程 Page * 差分法求解数学物理方程 两端固定的弦的振动问题的求解步骤 (3)定解条件设置 思考： 1. 求时间分量的递推公式需要几个初始值？ 2. 对坐标分量呢? 3. 两个初始值时，另一个初始值如何处理？ Page * 差分法求解数学物理方程 两端固定的弦的振动问题的求解步骤 (4)整理得到递推公式 Page * 差分法求解数学物理方程 %P164_1.m N=4001; dx=0.0024; dt=0.0005; c=dt*dt/dx/dx; x=linspace(0,1,420); u(1:420,1)=0; u(181:240,1)=0.05*sin(pi*x(181:24