东北大学matlab,东北大学matlab上机作业.doc

东北大学matlab上机作业.doc

1、安装MATLAB软件，应用demo命令了解主要功能，熟悉基本功能，会用help命令。2、用MATLAB语句输入矩阵和AB，1423A4j13j214jj前面给出的是矩阵，如果给出命令将得出什么结果？5)6,(A解：>>A=[1,2,3,4;4,3,2,1;2,3,4,1;3,2,4,1]A=1234432123413241>>B=[1+4j,2+3j,3+2j,4+1j;4+1j,3+2j,2+3j,1+4j;2+3j,3+2j,4+1j,1+4j;3+2j,2+3j,4+1j,1+4j;]B=1.0000+4.0000i2.0000+3.0000i3.0000+2.0000i4.0000+1.0000i4.0000+1.0000i3.0000+2.0000i2.0000+3.0000i1.0000+4.0000i2.0000+3.0000i3.0000+2.0000i4.0000+1.0000i1.0000+4.0000i3.0000+2.0000i2.0000+3.0000i4.0000+1.0000i1.0000+4.0000i>>A(5,6)=5A=1234004321002341003241000000053、假设已知矩阵，试给出相应的MATLAB命令，将其全部偶数行提取出来，A赋给矩阵，用命令生成矩阵，用上述命令检验一下结果是Bmagic(8)A不是正确。解：>>A=magic(8)A=64236160675795554121351501617474620214342244026273736303133323435292838392541232244451918484915145253111056858595462631>>B=A(2:2:end,:)B=955541213515016402627373630313341232244451918488585954626314、用数值方法可以求出，试不采用循环的630632842iS形式求出和式的数值解。由于数值方法是采用double形式进行计算的，难以保证有效位数字，所以结果不一定精确。试采用运算的方法求该和式的精确值。解：>>sum(sym(2).^[1:63])ans=184467440737095516145、选择合适的步距绘制出下面的图形。(1)，其中；(2)，其中。)/sin(t)1,(t)tan(si)si(t),(t解：(1)t=-1:0.03:1;y=sin(1./t);plot(t,y)-1-0.8-0.6-0.4-0.200.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.60.81(2)x=[-pi:0.05:-1.8,_1.799:.001:-1.2,-1.2:0.05:1.2,1.201:0.001:1.8,1.81:0.05:pi];y=sin(tan(x))-tan(sin(x));plot(x,y)-4-3-2-101234-3-2-101236、试绘制出二元函数的三维图和三22)1()1(),(yxyxyxfz视图。解：xx=[-2:.1:-1.2,-1.1:0.02:-0.9,-0.8:0.1:0.8,0.9:0.02:1.1,1.2:0.1:2];yy=[-1:0.1:-0.2,-0.1:0.02:0.1,0.2:.1:1];[x,y]=meshgrid(xx,yy);z=1./(sqrt((1-x).^2+y.^2))+1./(sqrt((1+x).^2+y.^2));surf(x,y,z),shadingflat;zlim([0,15])-2-1012-1-0.500.51051015[x,y]=meshgrid(-2:.1:2);z=1./(sqrt((1-x).^2+y.^2));subplot(221),surf(x,y,z),view(0,90);subplot(222),surf(x,y,z),view(90,0);subplot(223),surf(x,y,z),view(0,0);7、试求出如下极限。(1)；(2)；(3)。xx1)93(lim1lim0xy2)(cos1lim20yxyxe解：(1)symsx;f=(3^x+9^x)^(1/x);limit(f,x,inf)ans=9(2)symsxy;fb=x*y/(sqrt(x*y+1)-1);limit(limit(fb,x,0),y,0)ans=2(3)symsxy;fc=(1-cos(x^2+y^2))*exp(x^2+y^2)/(x^2+y^2);limit(limit(fc,x,0),y,0)ans=08、已知参数方程，试求出和。ttyxsincolxyd3/2t解：>>symst;x=log(cos(t));y=cos(t)-t*sin(t);diff(y,t)/diff(x,t)ans=-(-2*sin(t)-t*cos(t))/sin(t)*cos(t)>>f=diff(y,t,2)/diff(x,t,2);subs(f,t,sym(pi)/3)ans=3/8-1/24*pi*3^(1/2)9、假设，试求。xytef0d),(2222yfxfy解：symsxytf=int(exp(-t^2),t,0,x*y);x/y*diff(f,x,2)-2*diff(diff(f,x),y)+diff(f,y,2)simple(ans)ans=2*x^2*y^2*exp(-x^2*y^2)-2*exp(-x^2*y^2)-2*x^3*y*exp(-x^2*y^2)10、试求出下面的极限。(1)；1)2(16412lim22nn(2)。)3(li22nn解：(1)symskn;symsum(1/((2*k)^2-1),k,1,inf)ans=1/2(2)>>symsknlimit(n*symsum(1/(n^2+k*pi),k,1,n),n,inf)ans=111、试求出以下的曲线积分。(1)，为曲线，，lsyxd)(2l)sin(cottax)cos(inttay。0t(2)，其中为正向上半lyyxeex)d2()(33l22cybx椭圆。解：(1)>>symsat;x=a*(cos(t)+t*sin(t));y=a*(sin(t)-t*cos(t));f=x^2+y^2;I=int(f*sqrt(diff(x,t)^2+diff(y,t)^2),t,0,2*pi)I=2*a^2*pi^2*(a^2)^(1/2)+4*a^2*pi^4*(a^2)^(1/2)(2)>>symsxyabct;x=c*cos(t)/a;y=c*sin(t)/b;P=y*x^3+exp(y);Q=x*y^3+x*exp(y)-2*y;ds=[diff(x,t);diff(y,t)];I=int([PQ]*ds,t,0,pi)I=-2/15*c*(-2*c^4+15*b