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

数学物理建模与计算机仿真考试试题 (2014——2015第一学期) 1、计翼机仿真计算积分: 1、计翼机仿真计算积分: (n为自然数) 积分方向为正方向(10分)。 解： c1 ear al 1 %利用留数定理求解有方向的闭合路径的积分计箪 m=zeros(l. 101)； m(l)=l； m(101) = l； [R.P.K]=residue([l 0].m)； inl=2*pi*sum(R) 2、在同一幅图中绘制函数y(,)= 2 (不要投影图) 和平面x=0的图像(10分)? + 3々 程序： function text2 clear al 1 z=2*cplxgrid(30)； cplxmap1(z,2./((5+3章z)+eps*(z==-3/5)))； title(2./(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): i f nargin 2 k = find((abs(w) B) | i snan(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 IT 1 m M])： caxis([-l 1]): s = ones(size (z))； surf(x.y,u.v)； colormap(hsv(64)) End 3、在同一幅图中绘制绘制勒让德函数的图形（10分）。 解:  程序: c1ear a11 t=0：0.01：l； yl=legendre（l,t）； y2=legendre（2tt）； y3=legendre（3.t）； plot（t,yl（L：）, 9r t t.y! （2, ：） / y-*, * i inewidth*. 3）: hold on plot （t,y2（l.：）, *b- . t.y2（2, ：）. k-\ t.y2（3. ：）, *m-. linewidth* .3）； hold on plot（t,y3（L：）, c-* tt.y3（2 ：） , t.y3（3. ：）/ — . t.y3（4? 1 inewidth9 .3）; legendCP「OLP 1 1, T 20*.*P 2T,「2 ? ?P 30.P 3°「.P 32,P 3 3）； title（勒让德.）； IM3、 求任意两个圆相交的区域，并求出相交区域的面积。（注: IM 如果两个圆不相交，输出“两个圆不相交”）（20分）。 程序： function text4 clear ali char ch： 先画两个圓 r0=100； rl=100； x0=0；y0=0；xl=40；yl=60； chf ; circle(rO.xO,yOtch)； ch=,r#: circle(rLxL yLch): %判断两个圆是否相交并求面枳 x3=[]； y3=[]； for x=-rl:rl for y=-rl：rl dl=sqrt ((x-xO)涉+ (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) fprintfC两个圆不想交 else 、计算面积 syms a b [altbl]=solve((a-xO)2+(b-y0)M2-r0*2,(a-xl) 2*(b-yl)*2-rl 2)； al二double(al): bl=double(bl): plot(x3ty3,k?’)； juli=sqrt ((al (1) -al (2)广 2+ (bl (1) -bl ⑵)”2); hO=sqrt(r0*2-(juli/2)2)； hl=sqrt(rl*2-(juli/2)2)； S=rO*rO*acos(hO/rO)/2+r1*rl*acos(h1/rl)/2-(hO+h1)*juli/2； fprintf(f输出面积的大小*：\n?f\n\S)； end end function circle(rtxO,yO.ch) theta=0：0. 01*pi：2?