2011041020007_宫克_第二次作业.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值对应的光场的振幅分布、相位分布，并简单总结出其特征规律。 答： function main z=cplxgrid(30); r=abs(z);phi=angle(z); r2=r.^2;count=1; for m=2:2:6 for n=2:2:6 Cmn=sqrt(2/pi*gamma(m+1)./gamma(m+n+1));%常系数Cmn的建立 L=Laguerrel(m,n,2*r2); z1=Cmn*(sqrt(2)*r).^2.*L.*exp(i*n*phi).*exp(-r2); figure(count) cplxmap(z,z1,5*pi),colorbar(vert) %绘图 str=[m=,num2str(m),，n=,num2str(n)]; title(str) count=count+1; end end end function L=Laguerrel(n, alpha, x) %拉盖尔多项式的建立 L= 1; bin= 1; for i= n:-1:1 bin=bin.*(alpha+ i)./(n+ 1- i); L= bin-x.* L./i; end end 由图像可以知道m跟振幅、相位有关；n表示有n个极大值点和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*