(完整版)MATLAB)课后实验[1]

实验一 MATLAB 运算基础

1. 先求以下表达式的值，尔后显示 MATLAB 工作空间的使用情况并保存全部变量。

(1)

2sin 85 z

1 2

1 e

(2) 1

2

z ln( x 1 x ) ，其中

2

2 x

2 1 2i

5

(3)

a

e e a

z sin( a 0.3) ln , a 3.0, 2.9, L , 2.9, 3

2 2

2

t 0 t 1

(4) 2

z t 1 1 t 2

4

2

t 2t 1 2 t 3

，其中

解：

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>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(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>='A'&ch<='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

x

1

1 1 1

x

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

1 其他

x x

用 if 语句实现，分别输出时的 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) y<3 时的最大 n 值。

(2) 与(1)的 n 值对应的 y 值。解：M—文件以下：

3. 考虑以下迭代公式：

x n 1

a

b x

n

其中 a、b 为正的学数。

(1) 编写程序求迭代的结果，迭代的停止条件为 |x n+1 -x n|≤10 -5，迭代初值 x0，迭代次数不高出 500 次。

(2) 若是迭代过程收敛于 r，那么 r 的正确值是

2 4

b b a

2

，当(a,b) 的值取

(1,1) 、

(8,3) 、(10,0.1) 时，分别对迭代结果和正确值进行比较。

解：

M 文件以下：