matlab题目和答案

《matlab题目和答案》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、2008-2009-2MATLAB初步入门习题1.已知矩阵，矩阵，解线性方程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=ABX=1.00004.00001.00004.00001.00004.00001.00004.00001.00004.0000>>C=[9,10,4;3,4,5]C=9104345

2、>>2.根据下面两个矩阵执行下列的矩阵运算命令，并回答有关的问题（1）A+5*B和A-B+I分别是多少（其中I为单位矩阵）？（2）A.*B和A*B将分别给出什么结果，它们是否相同，为什么？得出A.^B、A/B及AB的结果，并分别解释它们的物理意义。>>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=92235327010855118181>>A-B+Ians=410178417213701

3、28>>A.*Bans=42460482254681406241224>>A*Bans=1922282647388731008128415181752>>A.^Bans=1.0e+019*0.00000.00000.00000.00000.00000.00000.00000.00011.5917>>X=A/BWarning:Matrixissingulartoworkingprecision.X=NaN-InfInfNaN-InfInfNaN-InfInf>>ABWarning:Matrixisclosetosingu

4、larorbadlyscaled.Resultsmaybeinaccurate.RCOND=1.328379e-018.ans=0.0194-0.2737-0.5894-0.09432.10314.34540.1027-1.1071-2.3394>>3.分别用 for和while循环结构编写程序，求出并考虑一种避免循环的简洁方法来求和，并比较各种算法的运行时间。（注：一共是三种方法）>>n=0;>>s=1;>>whilen<64n=n+1;s=1+2^n;end>>ss=1.8447e+019>>>>n=0;>>s=0;

5、>>forn=0:63s=s+2^n;end>>ss=1.8447e+019>>>>n=64;>>a1=1;>>q=2;>>s=a1*(1-q^n)/(1-q);>>ss=1.8447e+019>>4.画图题要求：1、两根曲线画在一幅图中，不采用子图的方式。2、横坐标（0~5π），纵坐标为（-1.5~+1.5）3、要求给横坐标轴加说明“t(deg)”，纵坐标轴加“magnitude”，并且加网格4、在（π/2，1.2）处加创建说明性的文字“这是我的程序运行结果！”5、图中，sint用红色*线绘制，cost用篮色实线绘制。

6、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:1

7、0;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的图形是双

8、曲抛物面。（提示：在区域-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.求解微分方程>>y=dsolve('Dy-2*y/(x+1)=(x+1)^(5/2)')y=-1/2