MATLAB题目和答案

2008-2009-2MATLAB初步入门习题1.已知矩阵5432141097539108627810715675A，矩阵60151403514436136349624B，解线性方程X，使得A.X=B，并将矩阵A的右下角2×3子矩阵赋给矩阵C。A=[5,7,6,5,1;7,10,8,7,2;6,8,10,9,3;5,7,9,10,4;1,2,3,4,5];B=[24,96;34,136;36,144;35,140;15,60];X=A\BX=1.00004.00001.00004.00001.00004.00001.00004.00001.00004.0000C=[9,10,4;3,4,5]C=91043452.根据下面两个矩阵987654321,136782078451220124BA执行下列的矩阵运算命令，并回答有关的问题（1）A+5*B和A-B+I分别是多少（其中I为单位矩阵）？（2）A.*B和A*B将分别给出什么结果，它们是否相同，为什么？得出A.^B、A/B及A\B的结果，并分别解释它们的物理意义。A=[4,12,20;12,45,78;20,78,136];B=[1,2,3;4,5,6;7,8,9];I=[1,0,0;0,1,0;0,0,1];A+5*Bans=92235327010855118181A-B+Ians=41017841721370128A.*Bans=42460482254681406241224A*Bans=1922282647388731008128415181752A.^Bans=1.0e+019*0.00000.00000.00000.00000.00000.00000.00000.00011.5917X=A/BWarning:Matrixissingulartoworkingprecision.X=NaN-InfInfNaN-InfInfNaN-InfInfA\BWarning:Matrixisclosetosingularorbadlyscaled.Resultsmaybeinaccurate.RCOND=1.328379e-018.ans=0.0194-0.2737-0.5894-0.09432.10314.34540.1027-1.1071-2.33943.分别用for和while循环结构编写程序，求出6362326302222212iiK并考虑一种避免循环的简洁方法来求和，并比较各种算法的运行时间。（注：一共是三种方法）n=0;s=1;whilen64n=n+1;s=1+2^n;endss=1.8447e+019n=0;s=0;forn=0:63s=s+2^n;endss=1.8447e+019n=64;a1=1;q=2;s=a1*(1-q^n)/(1-q);ss=1.8447e+0194.画图题要求：1、两根曲线画在一幅图中，不采用子图的方式。2、横坐标（0~5π），纵坐标为（-1.5~+1.5）3、要求给横坐标轴加说明“t(deg)”，纵坐标轴加“magnitude”，并且加网格4、在（π/2，1.2）处加创建说明性的文字“这是我的程序运行结果！”5、图中，sint用红色*线绘制，cost用篮色实线绘制。6、给图形加标题“sinewavefromzeroto5/pi”x=0:0.1:5*pi;f=sin(x);g=cos(x);plot(x,f,'r*',x,g,'b');axis([0,5*pi,-1.5,1.5]);xlabel('t(deg)');ylabel('magnitude');text(pi/2,1.2,'这是我的程序运行结果!');title('sinewavefromzeroto5/pi');5.体会各种绘图命令及效果。假设用户有下面的绘图数据向量t=0:0.1:10;y=sin(tan(t))-tan(sin(t));且绘图函数采用下面各个函数，如polar()，bar()，stem()，stairs()，那么试得出并解释所得出的结果。t=0:0.1:10;y=sin(tan(t))-tan(sin(t));subplot(2,2,1);polar(t,y);subplot(2,2,2);bar(t,y);subplot(2,2,3);stem(t,y);subplot(2,2,4);stairs(t,y);6.证明：函数z=xy的图形是双曲抛物面。（提示：在区域-2≤x≤2,-2≤y≤2上作出它的图形。）x=-2:0.1:2;y=-2:0.1:2;[X,Y]=meshgrid(x,y);Z=X.*Y;surf(X,Y,Z);7.求解微分方程2/5)1(12xxydxdyy=dsolve('Dy-2*y/(x+1)=(x+1)^(5/2)')y=-1/2*(x+1)^(5/2)*x-1/2*(x+1)^(5/2)+exp(2/(x+1)*t)*C18.用subplot分别在不同的坐标系下作出下列四条曲线，为每幅图形加上标题，①概率曲线2xey（-5≤x≤5）；②四叶玫瑰线sin2r（0≤x≤2pi）；③叶形线;13,13323ttyttx（-2≤x≤2）；④曳物线22111lnyyyx（-1≤x≤1）。x=1:0.01:5;y=exp(-x.*x);subplot(2,2,1)plot(x,y);title('概率曲线');subplot(2,2,2)a=linspace(0,2*pi);b=sin(2*a);polar(a,b);title('四叶玫瑰线');subplot(2,2,3)ezplot('3*t/(1+t.^3)','3*(t.^2)/(1+t.^3)',[020]);title('叶形线');subplot(2,2,4)f1=inline('log((1+sqrt(1-y.^2))./y)-sqrt(1-y.^2)-x');f2=inline('log((1-sqrt(1-y.^2))./y)+sqrt(1-y.^2)-x');ezplot(f1);holdon;ezplot(f2);holdon;title('曳物线');9.求11111111111122222222ddddccccbbbbaaaaD，并简化结果。symsabcdD=[a^2+1/a^2a1/a1;b^2+1/b^2b1/b1;c^2+1/c^2c1/c1;d^2+1/d^2d1/d1];D=simple(det(D))D=-(-d+c)*(-d+b)*(-c+b)*(-d+a)*(-c+a)*(-b+a)*(c*d*b*a-1)/a^2/c^2/d^2/b^210.11.simulink仿真题（选做）要求如下：1、绘制出系统结构图2、输入信号为step信号3、输出接scope和workspace4、要使能够看出系统响应的趋势，仿真时间如何处理？（给出具体措施）5、输出到workspace的保存类型（svaeformat）设定为何种？并使用plot语句绘制出系统响应曲线（用命令）21s92ss+-