《数学物理建模与计算机辅助设计》课堂作业布置2.docVIP

课堂授课专题4：特殊函数的可视化 对复变量Gamma函数进行可视化。 答：源码：z=5*cplxgrid(30); f=mfun(gamma,z); cplxmap(z,f); view(60, 30) axis([-5 5 -5 5 -10 10]) 采用与教材例题不同的其它方法对球函数进行可视化。 答：l=6; m=2; R=4; A=3;theta0=0:pi/40:pi; phi=0:2*pi/40:2*pi; [phi,theta]=meshgrid(phi, theta0); Ylm=legendre(l,cos(theta0)); Ylm=Ylm(m+1,:); L=size(theta,1); yy=repmat(Ylm,1,L); Reyy=yy.*cos(m*phi); Imyy=yy.*sin(m*phi); ReM=max(max(abs(Reyy))); Rerho=R+A*Reyy/ReM; Rer=Rerho.*sin(theta); Rex=Rer.*cos(phi); Rey=Rer.*sin(phi); Rez=Rerho.*cos(theta); subplot(1,2,1); surf(Rex,Rey,Rez); light; lighting phong; axis(square); axis([-5 5 -5 5 -5 5]); view(40,30) title(实球谐函数); ImM=max(max(abs(Imyy))); Imrho=R+A*Imyy/(ImM+eps*(ImM==0)); Imr=Imrho.*sin(theta); Imx=Imr.*cos(phi); Imy=Imr.*sin(phi); Imz=Imrho.*cos(theta); subplot(1,2,2); surf(Imx,Imy,Imz); light; lighting phong; axis(square); axis([-5 5 -5 5 -5 5]); view(40,30) title(虚球谐函数); 拉盖尔—高斯光束(Laguerre-Gaussian Beam)是激光器的特征模式，其横截面标量场分布为： 其中，m为非负整数，n为整数，为常系数，是拉盖尔多项式（Laguerre Polynomials）也为一种特殊函数。自行查找资料，对拉盖尔—高斯光束对其横截面进行仿真，绘出不同的m和n值对应的光场的振幅分布、相位分布，并简单总结出其特征规律。 答： 课堂授课专题5：数学物理方程的仿真与可视化 1. 已知边长为b和c的的四周固定的矩形膜振动的本征值问题为 绘制前4个本征函数的图形，并利用PDE工具箱对该二维本征值问题进行求解与可视化。 【提示】本征值问题泛定方程写为：，采用分离变量法可以得到本征模和本征值为： 答： A=4; B=2;x=0:0.01:A; y=0:0.01:B; [m,n]=meshgrid(1:3); Lamda=((m*pi./B).^2+(n*pi./B).^2) %求本征值 [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 %绘画前4个本征函数的图像 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); 2.已知初始速度为零，初始位移只在区间上不为零，且在处达到最大值，求此振动过程中的位移的动态图像。 【提示】定解问题（初始位移如图1所示），： 根据达朗贝尔公式即得位移为： . 作图时，先画出的图形，然后再用或替代其中的，改变的值，画出不同时刻和的图形。 3. 设初始位移为零即，而且初速度也只在区间上不为零 的无界弦振动，求此振动过程的位移动态分布. 【提示】根据达朗贝尔公式得 于是有 这里指的是图2的曲线。 可作出和两个图形，让它们以速度分别向左、右两个方向移动，两者的和就描画出各个时刻的波形，由此即得出位移分布．