(完整版)MATLAB)课后实验答案[1].docVIP

实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示 MATLAB 工作空间的使用情 况并保存全部变量。 (1) 0 2sin 85 z 1 2 1 e (2) 1 2 z ln( x 1 x ) ，其中 2 2 x 2 1 2i 0.45 5 (3) 0.3a 0.3 a e e 0.3 a z sin( a 0.3) ln , a 3.0, 2.9, L , 2.9, 3.0 3 2 2 2 t 0 t 1 (4) 2 z t 1 1 t 2 4 2 t 2t 1 2 t 3 ，其中 t=0:0.5:2.5 解： M 文件: z1=2*sin(85*pi/180)/(1+exp(2)) x=[2 1+2*i;-.45 5]; z2=1/2*log(x+sqrt(1+x^2)) a=-3.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=0t1).*(t.^2)+(t=1t2).*(t.^2-1)+(t=2t3) .*(t.^2-2*t+1) 4. 完成下列操作： (1) 求[100,999] 之间能被 21 整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。 解：(1) 结果： m=100:999; n=find(mod(m,21)==0); length(n) ans = 43 (2). 建立一个字符串向量 例如： ch=ABC123d4e56Fg9; 则要求结果是： ch=ABC123d4e56Fg9; k=find(ch=Ach=Z); ch(k)=[] ch = 123d4e56g9 实验二 MATLAB矩阵分析与处理 1. 设有分块矩阵 A E R 3 3 3 2 O S 2 3 2 2 ，其中 E、R、O、S 分别为单位矩阵、随机矩阵、零矩 阵和对角阵，试通过数值计算验证 2 A E R RS 2 O S 。 解: M文件如下； 5. 下面是一个线性方程组： 1 1 1 2 3 4 0.95 x 1 1 1 1 x 0.67 2 3 4 5 x 0.52 1 1 1 3 4 5 6 (1) 求方程的解。 (2) 将方程右边向量元素 b3 改为 0.53 再求解，并比较 b3 的变化和解的相对变化。 (3) 计算系数矩阵 A 的条件数并分析结论。 解： M 文件如下： 实验三 选择结构程序设计 1. 求分段函数的值。 2 x x 6 x 0且x 3 2 y x x x x x 5 6 0 5且 2及 3 2 x x 1 其他 用 if 语句实现，分别输出 x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0 时的 y 值。 解：M 文件如下： 2. 输入一个百分制成绩，要求输出成绩等级 A、B、C、D、E。其中 90 分~100 分为 A， 80 分~89 分为 B，79 分~79 分为 C，60 分~69 分为 D，60 分以下为 E。 要求： (1) 分别用 if 语句和 switch 语句实现。 (2) 输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 解：M 文件如下 3. 硅谷公司员工的工资计算方法如下： (1) 工作时数超过 120 小时者，超过部分加发 15% 。 (2) 工作时数低于 60 小时者，扣发 700 元。 (3) 其余按每小时 84 元计发。 试编程按输入的工号和该号员工的工时数，计算应发工资。 解：M 文件下 实验四 循环结构程序设计 1. 根据 2 1 1 1 1 L ，求π的近似值。当 n 分别取 100 、1000 、10000 2 2 2 2 6 1 2 3 n 时，结果是多少？ 要求：分别用循环结构和向量运算（使用 sum 函数）来实现。 解：M 文件如下： 运行结果如下： 2. 根据 y 1 1 1 1 L ，求： 3 5 2n 1 (1) y3 时的最大 n 值。 (2) 与(1)的 n 值对应的 y 值。 解：M—文件如下： 3. 考虑以下迭代公式： x n 1 a b x n 其中 a、b 为正的学数。 (1) 编写程序求迭代的结果，迭代的终止条件为 |xn+1 -xn|≤10 -5，迭代初值 x0=1.0 ，迭代 次数不超过 500 次。 (2) 如果迭代过程收敛于 r，那么 r 的准确值是 2 4 b b a 2 ，当(a,b) 的值取 (1,1) 、 (8,3) 、(10,0.1) 时，分别对迭代结果和准确值进行比较。 解： M 文件如下： 运算结果如下； 5. 若两个连续自然数的乘积减 1 是素数，则称这两个边疆自然数是亲密数对，该素数 是亲密素数。 例如，2× 3-