2014年数学物理建模期末考试试题.docxVIP

数学物理建模与计算机仿真考试试题 （2014——2015第一学期） （n为自然数）1、计算机仿真计算积分： （n为自然数） 积分方向为正方向（10分）。 解: clear al 1 弔利用留数定理求解有方向的闭合路径的积分计算 m=zeros（L101）； m(l) = l； m(101) = l； [R,P,K]=residue([l 0],m)； int=2*pi*sum（R） （不要投影图）2、在同一幅图中绘制函数于⑵二 2 和平面x = O的图像（10分）? + 3^ （不要投影图） 解: 2畑3?z) 15- 程序： function text2 clear al 1 z=2*cplxgrid(30)； cplxmapl (z,2. /((5+3*z) +eps? (z=~3/5)))； title(r2./(5-3*z)*)； hold on [y, z]=meshgrid(-3：0. 5：3.-3：0? 5:3): x二zeros(size (y))； mesh(x.y.z)； end function cplxmapl(z,w.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([-l 1 -1 1 m M])； caxis([-l 1])； s = ones(size(z))； surf(xty.u.v)； colormap(hsv(64)) End 3、在同一幅图中绘制绘制勒让德函数 环心妁 的图形(10分)。 解：  程序: clear al 1 t=0：0. 01：l； y1=Legendre （1,t）； y2=legendre（2t t）； y3=legendre （3rt）； plot（tf yl （1. ：）, . t,yl （2. ：） J y」J linewidth t3）; hold on plot（t,y2（l, ：） f t,y2（2, ：） t rk-, tty2（3, ：） , m-, * linewidth* ,3）； hold on plot（t,y3（l, ：）, t,y3（2f ：） / g-, tty3（3, ,t,y3（4, ：） /- 11 linewidth* ? 3）； legend（?p_「0「P_「「，°P_2V.?P 21「P :f T , 32\*P_3^ 3J; title （勒让德’）； 3、 求任意两个圆相交的区域，并求出相交区域的面积。（注: 如果两个圆不相交，输出“两个圆不相交”）（20分）。 程序： function text4 clear al 1 char ch； 弔画两个圆 心 100； rl=100； xO二0；y0=0；xl=40；yl=60； ch=*g*； circle(rO.xO,yO.ch)； ch=，r*； circle(rltxl,yl,ch)； 弔判断两个圆是否相交并求面枳 x3=[]； y3=[J； for x=-rl：rl for y=-rl：rl dl=sqrt((x-xO) 2+ (y-yO) 2)； d2=sqrt((x-xl)*2+(y-yl) 2)； if((dl=rO)(d2=rl))； x3=[x3.x]； y3=[y3.y]； end end end if (size(x3) ==0) fprintf (r可’、囲不想交’)； else 弔计算面积 syms a b [al,bl]=solve((a-xO)2+(b-y0厂2-rO 2.(a~xl) 2+(b-yl)*2-rl 2)； al=double(al)； bl=double(bl)； plot(x3,y3. k. *)； juli=sqrt((al(l)-al(2))^2+(bl(l)-bl(2))2)； hO=sqrt(r0*2-(juli/2) 2)； hl=sqrt(rl 2-(juli/2)2)； S=rO*rO*acos(hO/rO)/2+rl*r1*acos(h1/rl)/2-(hO+hl)*juli/2； fprintf (*输出面积的大小为:\n%f\n .S)； end end function circle(r.xO