Matlab上机实验答案Word格式.docx

1、0.5:2.5;z4=（t=0&t=1&2）.*（t.2-1）+（t=2&3）.*（t.2-2.*t+1）z4 =0 0 2.已知：12344A 34787, B65求以下表达式的值：（1）A+6*B 和 A-B+I （此中 I 为单位矩阵）（2）A*B 和 A.*B（3）A3 和 A.3（4）A/B 及 BA（5） A,B 和A（1,3,:）;B2A=12 34 -4;34 7 87;3 65 7;B=1 3 -1; 2 0 3;3 -2 7;A+6*Bans =1852-1046105215349I=eye（3）;A-B+I31-33288467（2） A*B684462309-72596

2、154-5241A.*B1022619-130（3） A3372262338244860424737014918860076678688454142118820A.3172839304-6434365850327274625（4） A/BBA（5）A,B-4-187-2A（1,3,:B2 ans =111920403.设有矩阵 A和B1661017A 11131415, B23222425（1）求它们的乘积 C。（2）将矩阵 C 的右下角 32 子矩阵赋给 D。（3）查察 MATLAB 工作空间的使用状况。A=（reshape（1:25,5,5）;B=3 0 16;17 -6 9;0 23 -

3、4;9 7 0;4 13 11;C=A*BC =93 150 77258 335 237423 520 397588 705 557753 890 717 D=C（3:5,2:3）D =520397705557890717whosName Size Bytes Class AttributesA5x5200doubleB5x3120CD3x2484.完成以下操作：（1）求100,999 之间能被 21 整除的数的个数。（2）建立一个字符串向量，删除此中的大写字母。n=100:999;l=find（rem（n,21）=0）;length（l）43 ch=aegbBOIEG0je23RGncwz=

4、find（ch=A&chZch（wz）=ch =aegb0je23nc实验二 MATLAB 矩阵解析与办理1. 设有分块矩阵 AE3R32 ，此中 E、R、 O、 S 分别为单位矩阵、随机矩阵、零矩O2S2阵和对角阵，试经过数值计算考据A2ER RSO。 E=eye（3）; R=rand（3,2）; O=zeros（2,3）; S=diag（1:2）; A=E,R;O,SA =0 0 0 0 A2E,R+R*S;O,S2A2=E,R+R*S;经过考据，矩阵 A2建立。2. 产生 5 阶希尔伯特矩阵 H 和 5 阶帕斯卡矩阵 P，且求其行列式的值 Hh 和 Hp 以及它们的条件数 Th 和 Tp

5、 ，判断哪个矩阵性能更好。为何？H=hilb（5） H =P=pascal（5） P =153570Hh=det（H） Hh =Hp=det（P） Hp =Th=cond（H） Th =4.7661e+05Tp=cond（P） Tp =8.5175e+03答： 5 阶帕斯卡矩阵 P 的性能好。矩阵的性能是由条件数决定的，条件数越凑近于 1 其性能就越好。由上机操作求得 Th=4.7661e+005 ， Tp=8.5175e+003 。 Tp 的值更凑近于 1 则其性能要好。所以 5 阶帕斯卡矩阵 P 的性能好。3.建立一个 5 5 矩阵，求它的行列式值、迹、秩和范数。A=rand（5） det

6、（A） trace（A） rank（A）norm（A） ans =4.已知29 6 18A205128 8 5求 A 的特色值及特色向量，并解析其数学意义。A=-29,6,18;20,5,12;-8,8,5-29-8V,D=eig（A）V = 00 0在数学上，特别是线性代数中，对于一个给定的线性变换，它的 特色向量 （本征向量或称正规正交向量 ）是这样一个非零的向量 v：当 v 经过这个线性变换的作用以后，获取的新向量（长度也许改变）依旧与本来的 v 保持在同一条线上。一个特色向量的长度在该线性变换下缩放的比率称为其 特色值 （本征值）。假如特色值为正，则表示 v 在经过线性变换的作用后方向

7、也不变；假如特色值为负，说明方向会反转；假如特色值为0，则是表示缩回零点。但无论如何，仍在同一条直线上。5.下边是一个线性方程组：x1x2x3（1）求方程的解。将方程右侧向量元素b3 改为 0.53 再求解，并比较b3 的变化和解的相对变化。（3）计算系数矩阵 A 的条件数并解析结论。 format rat % 用分数格式显示A=1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6 A =1/21/31/41/51/6 format% 恢复默认格式b=0.95;0.67;0.52b =X=Ab X = b2=0.95;0.53b2 = X2=Ab2X2 = D=cond（A）

8、1.3533e+03矩阵的条件数决定矩阵的性能，条件数越凑近于1 其性能越好，系数矩阵A 的条件数为1.3533e+003，和 1 相差很大，则其性能不好。所以解（ X 值）影响很大。b 矩阵个别元素的渺小变动，对方程的6.建立 A 矩阵，试比较 sqrtm（A） 和 sqrt（A） ，解析它们的差别。实验三选择结构程序设计1.求分段函数的值。x 6x0且 xyx25x 6 0x 5且 x2及 x 3x 1其余用 if 语句实现，分别输出时的 y 值。程序设计：clear;x=-5.0 -3.0 1.0 2.0 2.5 3.0 5.0;if xx=90&score=80&90=70&80=60

9、&（三） 用 switch 语句score=input（ switch fix（score/10）case9,10otherwise（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。=100%判断成绩合理性switchfix（score/10）disp（ error ） % 成绩不合理时输出出错信息3.硅谷公司员工的薪水计算方法以下：（1） 工作时数超出 120 小时者，超出部分加发 15% 。（2）工作时数低于 60 小时者，扣发 700 元。（3）其余按每小时 84 元计发。试编程按输入的工号和该号员工的工时数，计算应发薪水。num=input（ 请输入员工工号：t

10、ime=input（ 请输入员工工时数：iftimewages=84*120+（time-120）*84*（1+0.15）;wages=time*84;disp=（ 员工工号： num2str（num） 应发薪水： num2str（wages）请输入员工工号： 4请输入员工工时数： 35disp =员工工号： 4 应发薪水： 22404.设计程序，完成两位数的加、减、乘、除四则运算，即产生两个两位随机整数，再输入一个运算符号，做相应的运算，并显示相应的结果。解：a=floor（rand（1）*90+10）b=floor（rand（1）*90+10）c=input（ 请输入运算符号： , sc=

11、+s=a+bc= -s=a-b*s=a*bs=abdisp（a =3077请输入运算符号： -s =-475.建立 5 6 矩阵，要求输出矩阵第 n 行元素。当 n 值超出矩阵的行数时，自动转为输出矩阵最后一行元素，并给出出错信息。A=rand（5,6）n=input（ n=? ）if nn=5B=A（n,1:6）;erro!B=A（5,1:n =B = 实验四循环结构程序设计1. 依据L2 ，求的近似值。当 n 分别取 100 、 1000 、10000n时，结果是多少？分别用循环结构和向量运算（使用sum 函数）来实现。程序设计（一） ：for n=100,1000,10000sum=0;for x=1:sum=sum+1/（x.2）;result=sqrt（6*sum）result =程序设计（二） ：x=1:n;result=sqrt（6*sum（1./x./x）2. 依据 y 1 1，求：2n（1）y3 时的最大 n 值。（2）与 （1） 的 n 值对应的 y 值。y=0;n=1;while y10（-5）&=499x0=x1;y=abs（x1-x0）