Matlab上机练习二答案

1、求i44i93i49i67i23i57i41i72i53i84x的共轭转置。x=[4+8i3+5i2-7i1+4i7-5i;3+2i7-6i9+4i3-9i4+4i];x’ans=4.0000-8.0000i3.0000-2.0000i3.0000-5.0000i7.0000+6.0000i2.0000+7.0000i9.0000-4.0000i1.0000-4.0000i3.0000+9.0000i7.0000+5.0000i4.0000-4.0000i2、计算572396a与864142b的数组乘积。a=[693;275];b=[241;468];a.*bans=12363842403、对于BAX，如果753467294A，282637B，求解X。A=[492;764;357];B=[372628]’;X=A\BX=-0.51184.04271.33184、463521a，263478b，观察a与b之间的六种关系运算的结果。a=[123;456];b=[8–74;362];abans=010101a=bans=010101abans=101010a=bans=101010a==bans=000000a~=bans=1111115、7.0802.05a，在进行逻辑运算时，a相当于什么样的逻辑量。相当于a=[11011]。6、角度604530x，求x的正弦、余弦、正切和余切。x=[304560];x1=x/180*pi;sin(x1)ans=0.50000.70710.8660cos(x1)ans=0.86600.70710.5000tan(x1)ans=0.57741.00001.7321cot(x1)ans=1.73211.00000.57747、用四舍五入的方法将数组[2.45686.39823.93758.5042]取整。b=[2.45686.39823.93758.5042];round(b)ans=26498、设81272956313841A，793183262345B，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。A=[14813;-36-5-9;2-7-12-8];B=[543-2;6-23-8;-13-97];C1=A*B'C1=19-823012273-385429C2=A'*BC2=-1516-243663-1793-105226117-60194684-10C3=A.*BC3=51624-26-18-12-1572-2-21108-56inv(C1)ans=0.00620.0400-0.0106-0.00460.01690.00300.01680.02090.0150inv(C2)Warning:Matrixisclosetosingularorbadlyscaled.Resultsmaybeinaccurate.RCOND=8.997019e-019.ans=1.0e+015*-0.9553-0.2391-0.19970.27000.96670.24200.2021-0.2732-0.4473-0.1120-0.09350.1264-1.1259-0.2818-0.23530.3182inv(C3)???Errorusing==invMatrixmustbesquare.9、设x=rcost+3t,y=rsint+3,分别令r=2,3,4，画出参数t=0~10区间生成的x~y曲线。t=linspace(0,10);r1=2;x1=(r1*cos(t)+3*t);y1=r1*sin(t)+3;r2=3;x2=(r2*cos(t)+3*t);y2=r2*sin(t)+3;r3=4;x3=(r3*cos(t)+3*t);y3=r3*sin(t)+3;plot(x1,y1,'r',x2,y2,'b',x3,y3,'m')10、设f(x)=x5-4x4+3x2-2x+6(1)在x=[-2,8]之间取100个点，画出曲线，看它有几个过零点。（提示：用polyval函数）x=linspace(2,8,100);y=polyval([10-43-26],x);plot(x,y,'b',x,0,'y')(2)用roots函数求此多项式的根。t=[10-43-26]p=roots(t)12、求解多项式x3-7x2+2x+40的根。r=[1-7240];p=roots(r);-0.21510.44590.79490.270713、符号函数绘图法绘制函数x=sin(3t)cos(t)，y=sin(3t)sin(t)的图形，t的变化范围为[0,2]。symstezplot(sin(3*t)*cos(t),sin(3*t)*sin(t),[0,pi])14、设,求x=sym('x');y=(sin(x))^4+(cos(x))^4;diff(y,10)15、dxxx632)9(x=sym('x');int((sqrt((9-x^2)^3)/x^6),x)16、213xxdxx=sym('x');int(x+x^3,x,1,2)17、求级数的和:122nnnsym(‘n’);symsum(n+2/2^n,n,1,inf)18、利用函数int计算二重不定积分x=sym('x');y=sym('y');z=(x+y)*exp(-x*y);a=int(z,x);int(a,y)19、试求出如下极限。（1）xxxx1)93(lim；（2）11lim00xyxyyx；（1）symsx;f=(3^x+9^x)^(1/x);l=limit(f,x,inf)ans=9（2）symsxy;f=x*y/(sqrt(x*y+1)-1);limit(limit(f,x,0),y,0)ans=2（3）symsxy;f=(1-cos(x^2+y^2))*exp(x^2+y^2)/(x^2+y^2);limit(limit(f,x,0),y,0)ans=020、已知参数方程tttytxsincoscosln，试求出xydd和3/22ddtxysymst;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)21、假设xytteyxf0d),(2，试求222222yfyxfxfyxsymsxyts=int(exp(-t^2),t,0,x*y);x/y*diff(f,x,2)-2*diff(diff(f,x),y)+diff(f,y,2)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)22试求出下面的极限。dxdyeyxxy)(（1）1)2(1161141121lim2222nn；（2）)131211(lim2222nnnnnnnsymskn;symsum(1/((2*k)^2-1),k,1,inf)ans=1/2limit(symsum(1/((2*k)^2-1),k,1,n),n,inf)ans=1/2（2）limit(n*symsum(1/(n^2+k*pi),k,1,n),n,inf)ans=123、假设一曲线数据点为x=0:2:4*piy=sin(x).*exp(-x/5)试将x的间距调成0.1，并用下列方法进行内插：（1）线性内插法（method=’linear’）（2）样条内插法（method=’spline’）（3）三次多项式内插法（method=’cubic’）（4）多项式拟合法：直接利用6次多项式去通过7个数据点请将这些内插法的结果及原先的数据点画在同一个图上25、阅读以下程序并在指定位置上写注释（6分）t=(0:0.02:2)*pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z,'r-',x,y,z,'bd')%___(1)_绘制一条红色连续曲线和蓝色菱形离散点view([-82,58])%___(2)__以方位角-82度，俯角58度观察视图boxon%显示坐标轴的矩形框legend('图形','视角')%__(3)_标注图例______holdon%__(4)_保持图形_____x=-4:4;y=x;[X,Y]=meshgrid(x,y);%__(5)_生成平面网格坐标矩阵__Z=X.^2+Y.^2;surf(X,Y,Z);%__(6)_绘制三维曲面_26、编写一个M函数文件fun_es(x)，计算如下函数：)sin(5.023/xxeyx，其中参数可以为标量，也可以为向量。在Matlab指令窗口键入如下命令检验该函数：fun_es(3)functionfun_es(x)y=0.5.*exp(x./3)-x.^2.*sin(x)27、编写m函数文件求半径为r的圆的面积和周长。functionf=circle(r)S=pi*r^2L=2*pi*rx,yMatlab上机实验四28、编制m文件，等待键盘输入，输入密码123，密码正确，显示输入密码正确，程序结束；否则提示，重新输入。functionf=zy2()t=input('Inputkey(keyis123)')ift==123'good'else'Youarefool,inputagain'zy2;//从新输入；end4.3已知y=22221312111n。当n=100时，求y的值步骤：建立函数，命名为Cbb.m求平方倒数的和。其程序如下functionF=Cbb(n)F=0;k=1;whilek=nF=F+(1/k)/k;k=k+1;endend在工作空间中求y值y=Cbb(100)y=1.63504.5编制程序，计算1+2++n2000时的最大n值步骤：建立函数文件Dbb.m,求累加和小于n的最大数，使n=2000程序如下functionk=Dbb(n)k=1;f=0;whilefnf=f+k;k=k+1;endk=k-2;end将2000作为输入值，求n.在工作空间中输入如下：Dbb(2000)ans=624.6分别用for和while循环结构编写程序，求出K=6302ii分别用for和while循环结构编写命令文件Ebb,Fbb.并在工作空间中求KEbb程序如下f=1;J(1)=1;fori=1:63J(i+1)=J(i)*2;f=f+J(i+1);End求K程序：Ebb.mff=1.8447e+019Fbb程序如下I=1;H=1;i=1;whilei=63H=H*2;I=I+H;i=i+1;end在工作空间中求KFbb.mII=1.8447e+0191.倘若要是生成一个4行5列的矩阵a=rand(4,5)2.倘若要是对x进行赋值，从5到25，间隔是0.1x=5:0.1:25;倘若要是对x进行赋值，从2到20，中间间隔100个点x=linspace(2,20,100);4.用3种方法将图画在同一个窗口中x=5:0.1:25;y1=sin(x).*cos(2*x);y2=sin(2*x).*cos(x);plot(x,y1,x