MATLAB习题及答案.docx

1、MATLAB习题及答案习题:1， 计算与的数组乘积， 对于,如果，,求解X。3， 已知:,分别计算的数组平方和矩阵平方，并观察其结果。 4，角度,求的正弦、余弦、正切和余切。(应用sin,o,tan。)5, 将矩阵、和组合成两个新矩阵：(1)组合成一个43的矩阵,第一列为按列顺序排列的a矩阵元素,第二列为按列顺序排列的b矩阵元素,第三列为按列顺序排列的c矩阵元素，即 ()按照、b、的列顺序组合成一个行矢量,即 6， 将(x-6)(x-3)(x-8）展开为系数多项式的形式。（应用oly，olyva）,求解多项式3x2240的根。(应用roots) 8， 求解在x=8时多项式（x-1)(x-2)

2、（x-3)（x-4)的值。（应用pl，olyval） 9，计算多项式的微分和积分。(应用polyr,polyint，poysm) 0, 解方程组。(应用x=ab), 求欠定方程组的最小范数解。(应用pinv)， 矩阵,计算a的行列式和逆矩阵.（应用dt，inv)3, y=si（x）,x从0到2，x=0。0，求的最大值、最小值、均值和标准差(应用ma,in，mn,std), 参照课件中例题的方法，计算表达式的梯度并绘图(应用meshgrid， raient，contur, old n， iver）1,用符号函数法求解方程at2b*+c。（应用solve)6, 用符号计算验证三角等式:（应用sym

3、s，iple)17， 求矩阵的行列式值、逆和特征根.（应用syms，d，iv,ei）18, 因式分解: (应用ys, fato)19, ,用符号微分求d/x.（应用sym,dff）20， 符号函数绘图法绘制函数x=si()cos（t)，sin(3t）in()的图形,t的变化范围为0,。(应用sm,zpot)2， 绘制曲线,x的取值范围为5,5。（应用pt)22, 有一组测量数据满足,t的变化范围为010，用不同的线型和标记点画出a0。1、a=。和a=5三种情况下的曲线,在图中添加标题,并用箭头线标识出各曲线a的取值,并添加标题和图例框。(应用lo，title,xt，legnd)3,表中列出了个

4、观测点的6次测量数据,将数据绘制成为分组形式和堆叠形式的条形图。第次第2次第3次第4次第5次第次观测点1742观测点267247观测点3725观测点46427424,x=66 49 71 56 8，绘制饼图，并将第五个切块分离出来。 ， 用phere函数产生球表面坐标,绘制不通明网线图、透明网线图、表面图和带剪孔的表面图。(应用phere，ms, hdde of, urf, NaN)26， 编制一个解数论问题的函数文件:取任意整数，若是偶数,则用除，否则乘加1，重复此过程,直到整数变为1。2,有传递函数如下的控制系统,用ulink建立系统模型，并对系统的阶跃响应进行仿真。 27, 建立一个简单

5、模型，用信号发生器产生一个幅度为2V、频率为0。5Hz的正弦波，并叠加一个0.1V的噪声信号,将叠加后的信号显示在示波器上并传送到工作空间。28 建立一个模拟系统，将摄氏温度转换为华氏温度（T = 9/5Tc+32)。答案:1， 计算与的数组乘积。 a 9 ; 75;=2 4;46 8; a。*ban= 2 36 8 4 402，对于,如果,求解X。 A= 2;7 6 ;3 57；B=37 2 28; X=ABX 0.18 4.047 13183， 已知：,分别计算a的数组平方和矩阵平方,并观察其结果. a=1 2;4 6；78 ;.2ans 1 4 9 16 25 3 49 64 1as 3

6、 6 4 6 8 6 12 1261504， 角度,求x的正弦、余弦、正切和余切。 =30 45 6;1=x180p; si(x1) = 0。000 0。771 .8cos（1)ans = 0。8660 。70 0。0 an(x)ns 5774 1。0000 1。321cot(x）as = .731 .0000 05745, 将矩阵、和组合成两个新矩阵:()组合成一个43的矩阵,第一列为按列顺序排列的a矩阵元素,第二列为按列顺序排列的b矩阵元素，第三列为按列顺序排列的c矩阵元素,即 （2)按照a、的列顺序组合成一个行矢量,即 =4 2; 7; =7 1; 3; c59;6 2;% (） d=(

7、:) b（:) c(:) d = 4 5 5 8 2 1 9 7 3 2 (） =a(:）;b(:)；c(:) e = 5 2 7 8 1 3 5 6 9 2 或利用(1）中产生的 e=reshape(,1，12) as= 4 5 2 7 7 8 5 9 6, 将(6）(x-3）(x-8)展开为系数多项式的形式。 a=6 ; pa=ply(a)； 也可以用pa=oly（ 3 8)来替换1,2两行 ppa=poy2sym(pa) ppa x317*x2+9*x-4,求解多项式x-7x2+2x+40的根。 r=170; rots（r） = 5。0000 .000 2。0008, 求解在8时多项式(

8、x1)(x2) (x-3）(x4)的值。 =poly( 2 3 4)； oyvalm（p,8) ns = 84 9， 计算多项式的微分和积分。 clearf=y(*4-12*x14*x2+5*x+9）dif（f)in(f） ans 16*x336*x228*x+5 s 4/*5-*x4-1/3+52x2+x10，解方程组. 2 ;3 11；2 2 ; b=3 6 6; =a x 7。400 -0.2000 -1。400011,求欠定方程组的最小范数解。 a=4 7 4；9 5 6; b=85; x=nv(a）b x -02151 0。49 0。799 0。270712，矩阵，计算的行列式和逆矩

9、阵. a=4 6；7 5 4 ；3 4 9; ad=det(a) ai=nv(a) a = -64 ai = -.4531 0.6562 0。593 0.69-0。87 0。9062 0。201 0。562 0。03713 ysin（x）,x从0到2，x=0。2，求的最大值、最小值、均值和标准差。 =0:0.02*i:2*i； =sin(x); ymax=max（y） ymin=mi() ymeanmean() ytd=std（y) ax 1 ymn= -1 ymen = 2。295e17 ystd = 0。77114, 参照课件中例题的方法，计算表达式的梯度并绘图。 v 2：2:; x，y

10、=meshrid(v）; z=10（.5）exp（x。2y。）; px,py =adet(z，。，。）; conou(x,y,z） holon quiver(,px,py） hol ff1, 下面三种表示方法有什么不同的含义?(1)x+5x2(2)f=3x25*x+2(3)x=ym() f=3*2+5*x+2(1)f=*x2+*+表示在给定x时,将3*x2+5*x的数值运算结果赋值给变量,如果没有给定x则指示错误信息。（2）3x25*x表示将字符串3*x2+*x+赋值给字符变量f，没有任何计算含义,因此也不对字符串中的内容做任何分析.(3）xsy(x） f=3*+5*x+2表示x是一个符号变量

11、，因此算式=3*x5*x+2就具有了符号函数的意义,f也自然成为符号变量了。16， 用符号函数法求解方程at+*c=0。 r=solv（+b*tc0，t) r 1/a（b+(b2*a)（/)） 1/2a(-b-（4*a*c)(/2)）, 用符号计算验证三角等式：(应用syms,siml） sin(1)co(2）-cos(1)sin(2)si(-2) syms pi1 h2； ysimple（sn(hi1)*cs(ph）cos（pi1）*n(hi2)） y = sin(ph-ph)18, 求矩阵的行列式值、逆和特征根。 syms a11a121 a2; Aa1,a2;a21,a22 D=de(A

12、） %行列式 AI=inv(A) 逆 E=eg（A) % 特征值 A = 11， 2 a21, a2 = a11*a22a1*a1 A -a22/（a1*22a1*a1), a12/(a1*22+a12*a2) a21/(a1*a2a12*a21), -a11/(-a11a2+12*a21) AE = 12*1+/2*a21/*(a112a1a22+a22+4*12*21）（1/2） 1/*1/2a22/2*(a11-2a1*a22+22+4*2*a2）(1/2）19，因式分解： symsx; =x4x35*56; actor（f） as = ()(x2)（）*（x+)20，,用符号微分求f/

13、dx。（应用yms,iff） smsa x; f=a，x2, 1/x; exp(a*x), lo(x）， in（x）; df=iff(f) d 0， *, 1/x2 aep(ax), 1/x， cos（x)21, 符号函数绘图法绘制函数x=sn（3t）cos(t）,sin(3t)si(t)的图形，的变化范围为0,。 sm ezplt（sin(*t）cs(t）,sin（3*t）*sin(t)，0，p) 2, 绘制曲线,x的取值范围为5,。 ：0。2:5； =x。3+x1; plot(x,）3,有一组测量数据满足,t的变化范围为010，用不同的线型和标记点画出a=0。1、a0.2和a0.三种情况下

14、的曲线,在图中添加标题，并用箭头线标识出各曲线的取值,并添加标题和图例框. 0:0。:10; y=exp（。1*)； y=xp（0。*t); y3=ep(05t); plot（t,y，b，t,y2，:*,y3,.g） ttle(tyrm=e-itat) title(tyrmetat，otSze，2) tt(t(6)，y(6)，leftarroim=0.1,FontSze,11) tex(t(6)，y(6）,earowiarm=0。2,FnSze,11) tet（t(6),y3(6）,leftrrowitam=0。,Fntize,1) ttle（itrm=e-ta，FontSze,12） leg

15、end（=0.1，a02，a=.5）5，表中列出了4个观测点的次测量数据,将数据绘制成为分组形式和堆叠形式的条形图。第1次第2次第3次第4次第5次第6次观测点1367428观测点2673247观测点397254观测点4643274 y=3 6 9 6；67 7 4;732 3;4 ;2 8 7;8 74 4；bar(）26， = 9 1 56 3，绘制饼图，并将第五个切块分离出来. x=66 49 1 56 3； 000 0 1； pi(x，L)7, 用spher函数产生球表面坐标,绘制不通明网线图、透明网线图、表面图和带剪孔的表面图。 x,y,sphre(0); m(，y,z) ms（,y,

16、),hiddoff sf（，y，z） z(18:30，1:)=NN*ones（1,)； surf(x,y，z)8，有一周期为的正弦波上叠加了方差为01的正态分布的随机噪声的信号，用循环结构编制一个三点线性滑动平均的程序。（提示：用0.1*an（1，）产生方差为0。1的正态分布的随机噪声;三点线性滑动平均就是依次取每三个相邻数的平均值作为新的数据,如x1()=(x（）x(2)+x()3，1(3)=(x(2)x(3)+x（)3） =:pi/5:4*p； lengt(t); y=s(t)+0。1*rndn（1,）； ya(1)=（1); for i:n y(i)=sum（y（1:i+1)）/3； e

17、d ya（)=（n); plot(t,y，c,t,ya,r,linidth,)9, 编制一个解数论问题的函数文件:取任意整数，若是偶数,则用2除,否则乘3加1,重复此过程，直到整数变为1 nconollaz(n) % cllatz % Clsic “3+1 Ploblemfom nuber heo =n； whle n1 m(n,2)=0 n=n/； ese n=3*+1; ed c=n； end3,有传递函数如下的控制系统，用Simulink建立系统模型,并对系统的阶跃响应进行仿真。 31, 建立一个简单模型,用信号发生器产生一个幅度为2V、频率为.5z的正弦波,并叠加一个。1V的噪声信号，将叠加后的信号显示在示波器上并传送到工作空间。 32,建立一个模拟系统,将摄氏温度转换为华氏温度(Tf = 95Tc+3)。