MATLAB程序设计与应用刘卫国编课后实验答案.docx

1、MATLAB程序设计与应用刘卫国编课后实验答案实验一 LB运算基础、 先求下列表达式得值,然后显示MTA工作空间得使用情况并保存全部变量.(1） (2） ，其中（3） （4） ,其中t=0：0、：2、5解:M文件：=*in（85p/0)/（1+exp（）)x +2；-、5 5；z212log（xsqrt(1x） a=3、：0、1：、0；3(ep（、3、*a）-exp(、3、a）、/2、sn(a+0、3）og(（0、3+a）、/)t=0：0、5:、5；z=(0t1）、*（t、2）+（t=1&2)、*（、21）+(=t3） 、*(t、-*t+1)运算结果:z1=2*sin(pi1)/（xp（）x=

2、2 1+2*i；、5 5；z2=12lg(x+sqrt(+x2）) 3、0：、：3、；3（exp(0、3、）-ex(0、3、*a)、2、sin(a0、3)log(（0、3+a）、/)=0:0、5:2、；z4=(t0t1）、（、2)+(t=1&=2&t3） 、(t、2+)z 0、275z2 = 0、7114 0、0253 、898 、36i 0、39 + 0、93i 、541 0、0044iz3 Co hrog 4 0、7388 +、116i 0、7696 3、1416i 、781 +3、146 0、913 3、141 Columns5 trough 8 、7822 + 、1416 0、7602

3、 3、16 0、254+ 3、141i 、6784+3、416i olumns 9 through 12 0、616 +3、1416 、5496+3、146i 0、468 + 3、1416i 0、780 3、416i Clums 13 troug 6 0、275 + 3、1416i 0、1680 +、1416i 0、49 + 3、6i 0、0771 +3、416 Colums 17 rough 20 0、124 + 3、416i 0、366+ 3、1416i 0、104+3、16i0、752 +3、1416i olumns 21 rough24 0、86+ 、1416i 1、097 + 3、1

4、416 -、70+ 3、146i -1、271 + 3、141i Column 25 hrouh2 、843+ 3、116 -2、227 + 3、116 2、983 +3、116 37、0245 Colmns 29 thrugh 2 、0017 2、3085 、871 -、5978 Columns 3 trough 36 1、355 1、151 -0、93 0、8083 Cumns 37 hogh 40 -、6567 0、5151 0、3819 -0、2561 Columns 41 trog44 -、134 0、0255 0、0792 0、6 Columns 45 through 48 、26

5、63 0、478 0、06 0、4841 Coumns9 rogh 2 0、537 、515 、6145 0、6366 Comns5 trugh 0、6474 0、6470 0、3 、61 Columns throu 0 0、777 0、5327 、47 0、4126 Column 0、3388 = 0 0、200 0 、2500 1、000 2、25002、 已知：求下列表达式得值：（) A6*与A-B+I（其中I为单位矩阵）(） A与A、*B（） A3与、3（4）/及BA(） A，B与A(1,3，：）；B2解: 文件：A=12 34 -4;34 787；3 65 ;B=1 1；20 3；3

6、 7;A+6、*B -Beye(3) A*B A、*B 3 A、3 A/B A，B （1，,:）;B2运算结果:A=12 3 -4；4 7 87；3 ;=1 3 -；2 03; 2 ;A+6、B A+e(）B A、B A3 A、3 B BAA， A（,3,:）；Bns = 18 2 10 46 7 105 21 53 49ans = 31 32 8 8 0 7 1an = 68 44 62 309 2 596 15 5 21as = 12 10 6 0 26 9 -1 49ns 326 3324 480 2730 14918 6006 78688 44142 1188ns = 178 3930

7、4 64 3904 33 5503 274625 3an = 16、400 -、600 7、6000 35、8000 -7、2000 0、2000 67、0 -34、0000 68、00ans = 109、400 13、200 22、800 53、00 8、00 171、000 61、6000 9、8000 186、2000ans 2 34 -4 1 34 0 3 3 65 7 7an 12 4 - 3 65 7 5 1 11 0 19 20 5 403、 设有矩阵A与B（1) 求它们得乘积C.（2） 将矩阵C得右下角32子矩阵赋给。（3） 查瞧MATLAB工作空间得使用情况.解：、 运算结果

8、:E=(reshae(1：1：,，)）；F=3 16；7 9；0 2 -4;9 7 0; 3 11;C E*H=C（3：5，：3)C 9 150 77 258 335 27 423 520 397 58 705 557 75 890 17H = 97 7 55 89 7174、 完成下列操作：(） 求100，9之间能被21整除得数得个数。（2） 建立一个字符串向量,删除其中得大写字母.解：(1) 结果:m=10：99；n=find(mod（m,21）=0)；lenh（n）n 43（）、 建立一个字符串向量 例如：h=AB12345F；则要求结果就是：ch=ABC23456Fg9；k=fid(c

9、hA&hTp，所以pascal矩阵性能更好.3、 建立一个55矩阵，求它得行列式值、迹、秩与范数.解: M文件如下：输出结果为:A= 17 24 8 5 23 5 7 14 1 4 6 13 20 22 1 12 19 21 1 5 2 d = 5070t = 6c1 = 6、500c2 = 5、4618cinf= 6、504、已知求得特征值及特征向量,并分析其数学意义。解：M文件如图:输出结果为:V = 0、0 0、803 0、2733 -0、6084 、7867 0、8725 0、3487 、5501 0、40 -5、36 0 10、18 0 0 0 16、831数学意义:得个列向量就是A

10、得特征向量，得主对角线上个就是A得特征值,特别得，V得个列向量分别就是D得3个特征值得特征向量.、 下面就是一个线性方程组：（1)求方程得解.（2） 将方程右边向量元素3改为0、53再求解，并比较b得变化与解得相对变化。（） 计算系数矩阵A得条件数并分析结论。解： M文件如下：输出结果:X = 1、200 0、00 0、00X2 1、200 、600 0、6001、533e0由结果,X与得值一样,这表示b得微小变化对方程解也影响较小，而A得条件数算得较小，所以数值稳定性较好，就是较好得矩阵。、 建立矩阵，试比较sqrt（A）与rt（）,分析它们得区别。解：文件如下：运行结果有：A = 16 1

11、8 20 5 2 9 8 1 = 3、8891 -0、1102 3、3 3、17 2、136 0、698 0、385 、70 、735b2 = 4、0002、4495 4、2426 4、4721 、236 3、641 、0000 2、284 2、361 = 1、000 、0000 、000 20、00 5、0000 、00 9、00 8、000 、000分析结果知：qrm（）就是类似得数值平方根（这可由b1bA得结果瞧出）,而sqrt（)则就是对中得每个元素开根号,两则区别就在于此。实验三 选择结构程序设计一、实验目得1、 掌握建立与执行M文件得方法。2、 掌握利用i语句实现选择结构得方法.3

12、、 掌握利用wich语句实现多分支选择结构得方法.、 掌握try语句得使用。二、实验内容1、 求分段函数得值。用if语句实现，分别输出=-5、0，-3、0，、0,2、0,2、5,、0,、0时得y值.解:M文件如下：运算结果有:f(5）y= 14 （3） 11 （）y = 2 f（2）y 1 f（2、5）y 0、2500 f（3)y = 5 (5） 192、输入一个百分制成绩，要求输出成绩等级A、B、C、D、。其中90分10分为A,80分分为B，分7分为C,60分69分为，6分以下为E.要求：（1） 分别用if语句与with语句实现。(） 输入百分制成绩后要判断该成绩得合理性，对不合理得成绩应输

13、出出错信息.解:M文件如下试算结果：score=8gae Bscore=23错误:输入得成绩不就是百分制成绩3、 硅谷公司员工得工资计算方法如下：（1)工作时数超过120小时者，超过部分加发15。(2）工作时数低于60小时者，扣发700元。（3） 其余按每小时84元计发。试编程按输入得工号与该号员工得工时数，计算应发工资。解：M文件下4、 设计程序,完成两位数得加、减、乘、除四则运算,即产生两个两位随机整数，再输入一个运算符号,做相应得运算，并显示相应得结果。解:M文件如下；运算结果例：a 3b = 3输入一个运算符：c =falsea 2 4输入一个运算符：+c = 15、 建立56矩阵,要

14、求输出矩阵第n行元素。当n值超过矩阵得行数时，自动转为输出矩阵最后一行元素,并给出出错信息。解：M文件如下:运算结果如下：输入一个5行6列矩阵A=2 345 5;2 5 7 6;2 22 2 2 3；11 2 7 3； 4 5 6 7输入一正整数n=4 9 7 输入一个5行6列矩阵A=12 3 5 ;23 4576；2 2 22;11 2 39 7 3; 4 6 7输入一正整数n=6 2 3 4 5 7as=Error sng =disTomany t arguets、实验四 循环结构程序设计一、实验目得1、 掌握利用fr语句实现循环结构得方法.2、掌握利用whle语句实现循环结构得方法。、

15、熟悉利用向量运算来代替循环操作得方法。二、实验内容1、 根据，求得近似值。当n分别取10、1000、1000时,结果就是多少?要求:分别用循环结构与向量运算（使用um函数)来实现。解：M文件如下：运行结果如下：K %循环结构计算p值y0;n=npu（n=);fr i=1：n y+1/i/i；endp=sqt（6*y)n100 3、132n100i 3、106n1000p=、115向量方法计算i值ninput(n=）;1、/（1：n）、2；s=sum（i)；pi=qrt（6s）n=1p = 、121n=100i = 3、140n=1000pi =3、14152、 根据,求:（） y n=n1;n

16、dny = 3、003n= 57n = 53、 考虑以下迭代公式：其中a、b为正得学数。（) 编写程序求迭代得结果,迭代得终止条件为xnx15,迭代初值x0=1、0,迭代次数不超过50次.(2） 如果迭代过程收敛于,那么r得准确值就是，当(a，b）得值取(1，1)、(8，3）、(1,、1）时，分别对迭代结果与准确值进行比较。解:M文件如下:运算结果如下；请输入正数a=1请输入正数b=1 、6180 = 0、6180 、016r= 0、80 -1、6180 -、00 2、26请输入正数a8请输入正数b=3 = 1、06r 、7016 1、610r = 1、01 4、716=1.0 -6、01请输

17、入正数1请输入正数b=0、1x 3、127r 3、127 4、16r = 3、112 、22 = -0、0000 6、32544、 已知求1f00中:(1） 最大值、最小值、各数之与。(2） 正数、零、负数得个数。解:M-文件以下就是运算结果:ax（f)=43738635min（)=8994121152u（f）-427560151c149c223495、若两个连续自然数得乘积减就是素数，则称这两个边疆自然数就是亲密数对,该素数就是亲密素数.例如，23-15，由于5就是素数，所以2与3就是亲密数,5就是亲密素数。求2,0区间内:() 亲密数对得对数.（2）与上述亲密数对对应得所有亲密素数之与。解

18、:M文件：运算结果为:j= 29s = 2361实验五 函数文件一、实验目得1、理解函数文件得概念.2、 掌握定义与调用MATAB函数得方法。二、实验内容1、 定义一个函数文件,求给定复数得指数、对数、正弦与余弦，并在命令文件中调用该函数文件。解：M文件如下：函数shu、M文件:function e,l,s,c=fu()fush 复数得指数，对数,正弦，余弦得计算%e 复数得指数函数值 复数得对数函数值s 复数得正弦函数值c 复数得余弦函数值e=p(z）;l=og（z）;s=in(z)；c=cos(z）；命令文件M:=input（请输入一个复数z)；a，b,c,dfushu（z）运算结果如下:

19、z=npu(请输入一个复数z=);a,c,d=ushu（z)请输入一个复数z1+i = 1、467 2、87ib= 0、3466 + 0、7854ic= 、2985 +、635id 0、8337 0、98i2、 一物理系统可用下列方程组来表示:从键盘输入m1、m2与得值,求1、a2、N1与N2得值.其中g取9、8,输入时以角度为单位。要求：定义一个求解线性方程组A=B得函数文件,然后在命令文件中调用该函数文件。解：文件函数fc、文件:funcio X= f（A,）fc fc就是求解线性方程得函数 A就是未知矩阵得系数矩阵X=AB； 命令M文件：cc;m1input(输入1=)；m2=inut(输入2=)；thetainpu(输入the=)；x=het*i/18；g=9、8;A1*os（） m1 -in（x）0 msin（x） s（x） 0 0 m -sin（x） 0 0 cos（) ;B=0;m1*g；0；m2g