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

1、MATLAB程序设计与应用刘卫国编课后实验答案81653实验一 MATLAB运算基础1、先求下列表达式的值，然后显示MATLAB工作空间的使用情2sin85(1) zi况并保存全部变量Z21ln(x1x2),其中x 21 2i20.455Z3e0.3a e0.3asin (a0.3 a0.3) In ,2a 3.0, 2.9,川,2.9, 3.02t20 t1Z4t2 11 t2，其中 t=0：0、5:2、5t2 2t12 t3解：1e2M文件：z1=2*si n(85*pi/180”(1+exp(2)x=2 1+2*i;-、45 5;z2=1/2*log(x+sqrt(1+xA2)a=-3、

2、0:0、1:3、0;z3=(exp(0、3、*a)-exp(-0、3、*a)、/2、*sin(a+0、3)+log(0、3+a)、/2) t=0:0、5:2、5;z4=(t=0&t=1 &t=2&t=0&t=1 &t=2&t26630、34780 42060 4841Columns 49 through 520、53790、58150、61450 6366Columns 53 through 560、64740、64700 63510、6119Columns 57 through 600、57770、53270 47740 4126Colum n 610、 3388z4 =0 0、 25000

3、1、 2500 1、 00002、 25002、已知:12344131A 34787 ,B 2033657327求下列表达式的值：(1)A+6*B与A-B+l(其中I为单位矩阵)(2)A*B 与 A、*B(3)AA3 与 A、八3A/B及BA(5) A,B与A(1,3,:);BA2解：M文件：A=12 34 -4;34 7 87;3 65 7;B=1 3 -1;2 0 3;3 -2 7;A+6 、 *BA-B+eye (3)A*BA、*BAA3A、A3A/BBAA,BA(1,3,:);BA2运算结果：A=12 34 -4;34 7 87;3 65 7;B=1 3 -1;2 0 3;3 -2 7

4、;A+6、*BA-B+eye(3)A*BA、*BAA3A、A3A/BBA A,BA(1,3,:);BA2ans =1852-10467105215349ans =1231-3328840671ans =684462309-72596154-5241ans =1210246802619-13049ans =372262338244860424737014918860076678688454142118820ans =172839304-643930434365850327274625343ans =16、4000-13、60007、 600035、8000-76、200050、 200067、0

5、000 -134、000068、 0000ans =109、4000 -131、2000322、 8000-53、000085、0000 -171、0000-61、600089、8000 -186、2000ans =1234-413-13478720336573-27ans =1234-436574511101920-5403、设有矩阵A与B1234530166789101769A 1112131415 ,B 02341617181920970212223242541311(1)求它们的乘积C。(2)将矩阵C的右下角3X2子矩阵赋给D(3)查瞧MATLAB工作空间的使用情况。解:、运算结果：E

6、=(reshape(1:1:25,5,5);F=3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11; C= E*FH=C(3:5,2:3)C =9315077258335237423520397588705557753890717H5203977055578907174、完成下列操作：(1)求100,999 之间能被21整除的数的个数。(2)建立一个字符串向量，删除其中的大写字母。解:(1) 结果：m=100:999;n=fin d(mod(m,21)=0);len gth( n)ans =43(2)、建立一个字符串向量 例如：ch二ABC123d4e56Fg9：则要求结

7、果就是：ch=ABC123d4e56Fg9; k=fi nd(ch=A&ch f(-3) y =11 f(i) y = fy =1 f(2、5)y =-0、 2500 f(3)y =5 f(5)y =192、 输入一个百分制成绩，要求输出成绩等级 A、B、C、D、E。其中90分100分为A,80 分89分为B,79分79分为C,60分69分为D,60分以下为 E。要求：(1)分别用if语句与switch语句实现。(2)输入百分制成绩后要判断该成绩的合理性 ，对不合理的成绩应输出出错信息。解:M文件如下试算结果：score=88grade =B score=123错误:输入的成绩不就是百分制成绩

8、3、硅谷公司员工的工资计算方法如下 ：(1)工作时数超过120小时者,超过部分加发15%。(2)工作时数低于60小时者，扣发700元。(3)其余按每小时84元计发。试编程按输入的工号与该号员工的工时数 ，计算应发工资。解:M文件下4、设计程序，完成两位数的加、减、乘、除四则运算 ，即产生两个两位随机整数，再输入一个 运算符号，做相应的运算，并显示相应的结果。解：M文件如下；运算结果例：a =38b =33输入一个运算符：人c =falsea =92b =40输入一个运算符:+c =5、建立5 x 6矩阵,要求输出矩阵第n行元素。当n值超过矩阵的行数时，自动转为输出矩 阵最后一行元素，并给出出错

9、信息。解：M文件如下:运算结果如下输入一个 5 行 6 列矩阵 A=1 2 3 4 5 5;2 3 4 5 7 6;2 2 2 2 2 3;11 2 3 9 7 3;2 3 4 5 6 7 输入一正整数 n=411 2 3 9 7 3输入一个 5 行 6 列矩阵 A=1 2 3 4 5 5;2 3 4 5 7 6;2 2 2 2 2 3;11 2 3 9 7 3;2 3 4 5 6 7 输入一正整数n=62 3 4 5 6 7 ans =Error using = dispToo many in put argume nts、实验四循环结构程序设计、实验目的1、 掌握利用for语句实现循环结构

10、的方法。2、 掌握利用while语句实现循环结构的方法。3、 熟悉利用向量运算来代替循环操作的方法。、实验内容21、根据一632III12n，求n的近似值。n 分别取 100、1000、10000时,结果就是多少?要求:分别用循环结构与向量运算 （使用sum函数）来实现。解:M文件如下:运行结果如下K %循环结构计算pi值y=o；n=in put( n=);for i=1: ny=y+1/i/i;endpi=sqrt(6*y)n=100pi =3、 1321n=1000pi =3、 1406n=10000pi =3、1415%向量方法计算Pi值n=in put( n=);i=1、/(1: n)

11、、A2; s=sum(i);pi=sqrt(6*s)n=100pi =3、 1321n=1000pi =3、 1406n=10000pi =3、141511 12、根据y 1 ，求:3 5 I 2 n 1y y=0; n=0; while y3n=n-1;endn3、 0033n =57 n =563、考虑以下迭代公式：aXn 1b Xn其中a、b为正的学数。(1)编写程序求迭代的结果，迭代的终止条件为|Xn+1-Xn|W 10-5,迭代初值X0=1、0,迭代次 数不超过500次。(2)如果迭代过程收敛于 r,那么r的准确值就是 一 ，当 (a,b)的值取(1,1)、2(8,3)、(10,0、

12、1)时,分别对迭代结果与准确值进行比较。解：M文件如下：运算结果如下请输入正数a=1请输入正数b=1x =0、 6180r =0、 6180-4、7016r =0、 6180-1、6180s =-0、 0000-2、2361请输入正数a=8请输入正数b=3x =1、 7016r =1、 7016-1、6180r =1、 7016-4、7016s =0.0 -6、4031请输入正数a=10请输入正数b=0、1x =3、 1127r =3、1127-4、7016r =3、 1127-3、2127s =-0、 0000-6、32544、 已知f11n1f20n2f31n3fnfn 1 2 fn 2fn 3 n3求 f1 f 100 中：(1)最大值、最小值、各数之与。(2)正数、零、负数的个数。解:M文件以下就是运算结果max(f)=437763282635 mi n(f)=-899412113528 sum(f)=-7427456019