解：
程序设计：
clear;
a=floor(rand(1)*90+10)
b=floor(rand(1)*90+10)
c=input('请输入运算符号：','s');
ifc=='+'
s=a+b
elseifc=='-'
s=a-b
elseifc=='*'
s=a*b
elseifc=='\'
s=a\b
else
disp('error')
Columns 31 through 35
+ + + + +
Columns 36 through 40
+ + + + +
Columns 41 through 45
+ + + + +
Columns 46 through 50
+ + + + +
Columns 51 through 55
+ + + + +
Matlab
实验一 MATLAB运算基础
1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。
(1)
>> z1=2*sin(85*pi/180)/(1+exp(2))
z1 =
(2) ，其中
>> x=[2 1+2i; 5];
>> z2=1/2*log(x+sqrt(1+x^2))
>> D=C(3:5,2:3)
D =
520 397
705 557
890 717
>> whos
Name Size Bytes Class Attributes
A 5x5 200 double
B 5x3 120 double
C 5x3 120 double
D 3x2 48 double
4. 完成下列操作：
else
y=x.^2-x-1;
end
y
运行结果：
y =
2. 输入一个百分制成绩，要求输出成绩等级A、B、C、D、E。其中90分~100分为A，80分~89分为B，79分~79分为C，60分~69分为D，60分以下为E。
要求：
(1) 分别用if语句和switch语句实现。
程序设计：
（一）用if语句
clear;
5. 下面是一个线性方程组：
(1) 求方程的解。
(2) 将方程右边向量元素b3改为再求解，并比较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/2 1/3 1/4
score=input('请输入百分制成绩：');
a=fix(score/10);
ifa==9 | a==10
grade='A';
elseifa==8
grade='B';
elseifa==7
grade='C';
elseifa==6
grade='D';
else
grade='E';
end
grade
（二）用if语句
>> ch(wz)=[]
ch =
aegb0je23nc
实验
1. 设有分块矩阵 ，其中E、R、O、S分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证 。
>> E=eye(3);
>> R=rand(3,2);
>> O=zeros(2,3);
>> S=diag(1:2);
>> A=[E,R;O,S]
1 3 6 10 15
1 4 10 20 35
1 5 15 35 70
>> Hh=det(H)
Hh =ห้องสมุดไป่ตู้
>> Hp=det(P)
Hp =
1
>> Th=cond(H)
Th =
+05
>> Tp=cond(P)
Tp =
+03
答：5阶帕斯卡矩阵P的性能好。矩阵的性能是由条件数决定的，条件数越接近于1其性能就越好。由上机操作求得Th=+005，Tp=+003。Tp的值更接近于1则其性能要好。所以5阶帕斯卡矩阵P的性能好。
6. 建立A矩阵，试比较sqrtm(A)和sqrt(A)，分析它们的区别。
实验
1. 求分段函数的值。
用if语句实现，分别输出x=,,,,,,时的y值。
程序设计：
clear;
x=[ ];
ifx<0&x~=-3
y=x.^2+x-6;
elseifx>=0&x<5&x~=2&x~=3
y=x.^2-5*x+6;
(3)
>> A^3
ans =
37226 233824 48604
247370 149188 600766
78688 454142 118820
>> A.^3
ans =
1728 39304 -64
39304 343 658503
27 274625 343
(4)
>> A/B
ans =
>> B\A
ans =
3. 建立一个5×5矩阵，求它的行列式值、迹、秩和范数。
>> A=rand(5)
A =
>> det(A)
ans =
>> trace(A)
ans =
>> rank(A)
ans =
5
>> norm(A)
ans =
4. 已知
求A的特征值及特征向量，并分析其数学意义。
>> A=[-29,6,18;20,5,12;-8,8,5]
wages=84*120+(time-120)*84*(1+;
else
wages=time*84;
end
disp=(['员工工号：'num2str(num)'应发工资：'num2str(wages)])
运行结果：
请输入员工工号：4
请输入员工工时数：35
disp =
员工工号：4 应发工资：2240
4. 设计程序，完成两位数的加、减、乘、除四则运算，即产生两个两位随机整数，再输入一个运算符号，做相应的运算，并显示相应的结果。
clear;
score=input('请输入百分制成绩：');
ifscore>=90&score<=100
grade='A';
elseifscore>=80&score<90
grade='B';
elseifscore>=70&score<80
grade='C';
elseifscore>=60&score<70
(2) 工作时数低于60小时者，扣发700元。
(3) 其余按每小时84元计发。
试编程按输入的工号和该号员工的工时数，计算应发工资。
程序设计：
clear;
num=input('请输入员工工号：');
time=input('请输入员工工时数：');
iftime<60
wages=time*84-700;
elseiftime>120
grade='A';
case8
grade='B';
case7
grade='C';
case6
grade='D';
otherwise
grade='E';
end
grade
else
disp('error')%成绩不合理时输出出错信息
end
3. 硅谷公司员工的工资计算方法如下：
(1) 工作时数超过120小时者，超过部分加发15%。
1/3 1/4 1/5
1/4 1/5 1/6
>> format %恢复默认格式
>> b=[;;]
b =
>> X=A\b
X =
>> b2=[;;]
b2 =
>> X2=A\b2
X2 =
>> D=cond(A)
D =
+03
矩阵的条件数决定矩阵的性能，条件数越接近于1其性能越好，系数矩阵A的条件数为+003，和1相差很大，则其性能不好。因此b矩阵个别元素的微小变动，对方程的解（X值）影响很大。
A =
-29 6 18
20 5 12
-8 8 5
>> [V,D]=eig(A)
V =
D =