MatLab语言及应用--习题及参考答案

1、MatLab语言及应用-习题1.产生一个1x10的随机矩阵，大小位于(-5 5),并且按照从大到小的顺序排列 好！ clear all close all clc %A=rand(l, 10)sort (A, 2) A = 0.81470. 90580. 12700.91340.63240.09750. 27850.54690. 95750.S649arts =0.09750. 12700.27850.516g0.63240.8U70. 90580.91310. 95750.96491. A = 21 2 3; 4 35 6; 7 8 49;B = 求A2、A的2次数事，并求矩阵A的秩、A =

2、 21 2 3： 4 35 6: 7 8 49;B = 3; 2; 1;% 求 AB、ETA、AT 但)与 ATE (b)AB=ABBTA=B，*A a=inv(A) b=a*B%求A2(c)、ACd)的2次数曷c=A4d=Aft2%求矩随A的秩(e)、转送和行列式(g) e=rank(A) 餐A g=det (A)3; 2; 1;设求AB、BTA、A-l与ATB 转置和行列式。AB =708S86ETA =7884700.0489-0. 0022-0. 0027-0.00450.0296-0. 0033-0.0063-0. 00450. 02130. 139830. 0423-0, 0065

3、 f =21472358372960390428722856364973511820311131994568171715626441786570760g =d =3406047013622226612S151652268624702.请产生一个100*5的矩阵，矩阵的每一行都是1 2 3 4 5repeat (1:5, 100, 1)4545454512345d512345445512345451234545123454455I234545123454.请修改下面的程序，让他们没有for循环语句!A=l 2 3; 4 5 6; 7 8 9;r c=size(A);for i=l:1:rfor

4、j=l:1:cif (A(i,j)>8 A(i, j)<2) A(i, j)=0;end%A(L 9) = 0 0end end5. a=l 2 3 4 5, b=a (1) *a (5) +a (2) *a (4) +a (3) *a (3) +a (4) *a (2) +a (5) *a (1). 试用MATLAB中最简单的方法计算b,注意最简单。XXa=l 2 3 4 5; p=flipud(a")： b=a*pb =356.求下列联立方程的解 3x+4y-7z-12w=4 5x-7y+4z+ 2w=一3 x +8z- 5w=9 -6x+5y-2z+l0w=-8 求

5、系数矩阵的秩；求出方程组的解。a=3 4 -7 -12:5 -7 4 2:1 0 8 -5：-6 5 -2 10 ; %系数矩临Ra=.rarxk (a) %系数矩阵的秩b=4：-3：9;-8; %常数项矩阵X=ab %方程组的解-1.4841-0.68160.5337-1.24297.设 y=cos 0. 5+ (3sinx) / (l+x*2) 把x=02 n间分为101点，画出以x为横坐 标，y为纵坐标的曲线；XXx=linspace (0, 2的pi, 101)：y=cos(0. 5+3*sin(x)./(1+x. A2);figureplot k y)8 .设£&）

6、内-5-4*-4+3乂-2+2乂+6取乂=-2,8之间函数的值（取100个点），画出 曲线，看它有几个零点。XXx=linspace (-2, 8, 100):y=x. a5-4*x.m4+3*x. a2+2*x+6;figureplot (xy)axis(-2, 8,-100, 300)grid on%由图知共有三个零目的三维图形。9 .在-10, 10； TO, 10范围内画出函数X, Y =m.eshgrid(-10: 0. 5: 10): a=sqrt (X. "2十Y/2)十eps ： Z=sin(a) ./a;meshCX, Y, Z)0.5x,x<2 /(x) =

7、 51.5-0.25%, 2 < x < 610.编写一个函数，使其能够产生如下的分段函数:056 < x并调用此函数,绘制在X = 9 + 2范围的f（x） f（J + 2）曲线。（10分）fun.m* 勰 + function f=fun<x)- f=0. 5*x.*(x<=2) + ( 1.5-0. 25*x).*(Q2). *(x<=6)+0.5*&>6)：- endUntitled* X +XXx=0:. 1:2;y=fun(x).*fun(x+2)：plot (x, y)11.将一个屏幕分4幅，选择合适的步长在右.上幅与左下幅绘制出

8、下列函数的图 形。Jcos(x), r 九冗、X , 一 |2 2 （曲线图）；X2 V2f(x，y) = K + K；(-2 X < 2,-4 < y < 4)2- 4一(曲面图)。xl=(l/(-2)：2)*pi；f=sqrt (cos (xl)： x2= (-2:2)：y=(-2:2):fx.y=x2. A2/2A2+y. A2/4*2：subplot (2, 2n 2), plotxlj f )jaxis(pi/ (-2),pi/2, -L )»grid on图，)subplot (2, 2,3)plot (x2, y), axis ( -2, 2-4,4)

9、, grid ontitle(2)")12 .编写一个判断闰年的函数来判断从2010年到2100年之间有哪些年份是闰年。(闰年可以被4整除但不能被100整除，或者可以被400整除)XXa = input ('请输入一个年份用来判断：'' )if (100) 二 0 &&= 0 ) 11 (rem(a 100) = 0 也也 rem (a, 400) = 0)disp C这年是闰年)elsedisp C这年不是闰年D请输入一个年份用来判断；18991899这年不是闰年13 .给出某地一月内温度的测定值(内含平均温度、最低、最高温度)，并求其平 均

10、温度、最低、最高温度的平均值、标准差等。(参考下表)周口四月份每天的温度报表/111L度 日期平均温度 ()最高温度 ()最低温度 ()110.2198214. 31910313.0199415.81912L 09. 715768.0164715.2208814.8207910. 31771012. 11971114. 72291215. 327121322.629171420. 727161522. 331191613.015121716.925121813.51610198. 7116207. 31142111.81582210.91482315.621122417.825132518.6

11、24142622. 729172723.528182824.529182918.627143018.92715A=10.2 19 8; 14. 3 19 10:13, 0 19 9:15. 8 19 12:9.7 15 7;8. 0 16 4:15.2 20 8： Id. 8 20 7; 10. 3 17 7; 12. 1 19 7：14, 7 22 9:15.3 27 12:22,6 29 17:20.7 27 16:22, 3 31 19;13,0 15 12:16. 9 25 12:13.5 16 10:8, 7 11 6:7.3 11 4:11. 8 15 8; 10. 9 14 S;

12、15. 6 21 12;17. 8 25 13; 18. 6 24 14;22. 7 29 17:23. 5 28 18:24,5 29 18:18.6 27 14:18.9 27 15;%平均温度、最低、最高温度的平均值标准差m.e an (A) |st d(A) ans =15.376721.200011. 10004.85065.78584.293914.用二分法求方程f (x)十3+4女-2-10=0在区间1,2内根的近似值，为使误 差不超过10-5时所需要的二分次数。%f=iiLliiie Cx. 3+-10.蚁-20);x=l. 4 ： 0. 05： 1. 9 ；a=l;b=2：c

13、=(a+b)/2：三 while abs (f (c) > 1 e-4 if f (a) *f (c) >0a=c：elseb=c;end|c=(a+b)/2;end1.594615 .取 = L5,用牛顿迭代法求/*) = 1+4/-10 = 0的跟的近似值。Ffunction new (tolerance) xO=l.5k=0a=xO - 3+4*x CT 2Tob=3*x0*2+8*x0xl=x0-a/b while (abs (xO-xl)>tolerance) x0=xlk=k+la=x。0 3+4 0 0 2Tob=3*:<0 A2+8*x0>> new(L 5)L 3733xl=x0- a/benddisp (xl)disp (k)|16 .求人) = (sin、)”-。.5|，|=。的根。KKsynis tft=sm(t) A2*exp(-0. Ut)-O. 5abs(t) ; S = ftS =S=s