数学物理方法matlab建模.docxVIP

数学物理建模与计算机仿真考试试题 （2011——2012第二学期） 1、计算机仿真计算积分:向为正方向 1、 计算机仿真计算积分: 向为正方向（10分〉。 （n为自然数〉 同=2 解：利用留数定理计算闭合路径积分，^c/(z)dz = 2^Res[/(z),z/ 设 g=1+zn,仁z八(2*n)/(1+zn) 求出所有n个奇点并放在矩阵A中，A二solve(g) 再求出每个奇点kEA对应的留数，Res (k) = l imit (f*(z~A(k)), z, A(k)) 故所求的积分结果为：INT=2niZRes(k) 对应程序为： n=input ( 5= T ； syms z g=l+zAn; f=zA(2*n)/(1+zAn); A=solve (g); sum=0; for k=l: n Res=limit(f*(z-A(k))A z,A (k)); sum=sum+Res; end int=2*pi*i*sum INT=simplify(int) 2、在同一幅图中绘制函数 2、 在同一幅图中绘制函数 和平面x = 0的图像 0 分〉。 解：参照教材第二页，用cplxgrid指令绘图，去掉投影，程序如下: function cplxmap2(z,wz B) blue = 0.2; x = real (z); y = imag(z); u = real(w); v = imag(w); if nargin 2 k = find((abs(w) B) | isnan (abs (w))); if length(k) 0 u(k) = B*sign(u (k)); v(k) = zeros (size (k)); v = v/max(max(abs (v))); v(k) = NaN*ones(size (k)); end end M = max(max(u)); m = min(min(u)); axis ( [-1 1 -1 1 m M]); caxis ( [-1 1]); s = ones (size (z)); surf(xz yz v); colormap(hsv (64)) 再参照第8页绘出2/ (5-3z)的图像及x=0的图像，程序如下: z=cplxgrid(30); cplxmap2(z, 2?/((5-3 ?*z)+eps* (z==5/3))); colorbar(1 vert1) title ( * 2/ (5-3z) *) hold on [y,z]=meshgrid(-1:0.01:1,-1:0.01:1); a=0; x2=a*ones (size (y)); mesh(x2z yr z); hold off 输出图像为： 2/(5-3z)-1 -10.80.60.40.2002-0.4?0.608 2/(5-3z) -1 -1 0.8 0.6 0.4 0.2 0 02 -0.4 ?0.6 08 3、 绘制贝塞尔函数丿0，丿1，丿2，，3的图形并标注出零点(10分》。 解：参考书上29页绘图，程序如下： y=besselj (0:3, (0:0.2:10) 1 ); figure (1); plot ( (0:0.2:10) *,y); legend(1 JO 1, 1 JI1, 1J21, 1J31) hold on x=0:0.05:10; LD=[]; for n=0:3 y=besselj(nA x); for k=l:length (y)-1 if y (k)*y (k+1)0 h = interpl (y(k:k+l)zx(k:k+l),0); LD=[LDZ h]; end end end plot(LDZ 0, * o *) hold off 输出图像为： 果两个圆不相交，输出“两个圆不相交”〉（20分〉。 解：先写一个绘制任意圆的函数，程序如下： function circle(x0,yO,R) fai=0:2*pi/100:2*pi; y=yO+R*sin (fai); x=xO+R*cos (fai); plot (xz y, bj 1 LineWid th , 5); %fill(x,y,c1) hold on; axis equal end 根据输入绘制两个圆并计算面积： function twocircle xl=input(Xl= *); yl=input(Yl= *); rl=input(1Rl= *); x2=input(1X2=1); y2=input(* Y2=); r2=input(R2= *); circle(xlA ylz rl); circle(x2r y2z r2); Num=0; for x=xl-r1:0 ?1:xl + r1;