Matlab上机实验答案文档格式.docx
《Matlab上机实验答案文档格式.docx》由会员分享，可在线阅读，更多相关《Matlab上机实验答案文档格式.docx（36页珍藏版）》请在冰豆网上搜索。
Columns26through30
Columns31through35
Columns36through40
Columns41through45
Columns46through50
Columns51through55
Columns56through60
Column61
+
（4），其中t=0:
t=0:
z4=（t>
=0&
t<
1）.*（t.^2）+（t>
=1&
2）.*（t.^2-1）+（t>
=2&
3）.*（t.^2-2.*t+1）
z4=
00
2.已知：
求下列表达式的值：
（1）A+6*B和A-B+I（其中I为单位矩阵）
（2）A*B和A.*B
（3）A^3和A.^3
（4）A/B及B\A
（5）[A,B]和[A（[1,3],:
）;
B^2]
A=[1234-4;
34787;
3657];
B=[13-1;
203;
3-27];
A+6*B
ans=
1852-10
467105
215349
I=eye（3）;
A-B+I
1231-3
32884
0671
（2）
A*B
684462
309-72596
154-5241
A.*B
121024
680261
9-13049
（3）
A^3
3722623382448604
247370149188600766
78688454142118820
A.^3
172839304-64
39304343658503
27274625343
（4）
A/B
B\A
（5）
[A,B]
1234-413-1
34787203
36573-27
[A（[1,3],:
1234-4
3657
451
11019
20-540
3.设有矩阵A和B
（1）求它们的乘积C。
（2）将矩阵C的右下角3×
2子矩阵赋给D。
（3）查看MATLAB工作空间的使用情况。
A=（reshape（1:
25,5,5））'
B=[3016;
17-69;
023-4;
970;
41311];
C=A*B
C=
9315077
258335237
423520397
588705557
753890717
D=C（3:
5,2:
3）
D=
520397
705557
890717
whos
NameSizeBytesClassAttributes
A5x5200double
B5x3120double
C5x3120double
D3x248double
4.完成下列操作：
（1）求[100,999]之间能被21整除的数的个数。
（2）建立一个字符串向量，删除其中的大写字母。
n=100:
999;
l=find（rem（n,21）==0）;
length（l）
43
ch='
aegbBOIEG0je23RGnc'
wz=find（ch>
='
A'
&
ch<
Z'
ch（wz）=[]
ch=
aegb0je23nc
实验二MATLAB矩阵分析与处理
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]
A=
0000
0000
A^2
[E,R+R*S;
O,S^2]
A^2==[E,R+R*S;
11111
通过验证，矩阵成立。
2.产生5阶希尔伯特矩阵H和5阶帕斯卡矩阵P，且求其行列式的值Hh和Hp以及它们的条件数Th和Tp，判断哪个矩阵性能更好。
为什么
H=hilb（5）
H=
P=pascal（5）
P=
12345
1361015
14102035
15153570
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的性能好。
3.建立一个5×
5矩阵，求它的行列式值、迹、秩和范数。
A=rand（5）
det（A）
trace（A）
rank（A）
5
norm（A）
4.已知
求A的特征值及特征向量，并分析其数学意义。
A=[-29,6,18;
20,5,12;
-8,8,5]
-29618
20512
-885
[V,D]=eig（A）
V=
00
在数学上，特别是线性代数中，对于一个给定的线性变换，它的特征向量（本征向量或称正规正交向量）是这样一个非零的向量v：
当v经过这个线性变换的作用之后，得到的新向量（长度也许改变）仍然与原来的v保持在同一条线上。
一个特征向量的长度在该线性变换下缩放的比例称为其特征值（本征值）。
如果特征值为正，则表示v在经过线性变换的作用后方向也不变；
如果特征值为负，说明方向会反转；
如果特征值为0，则是表示缩回零点。
但无论怎样，仍在同一条直线上。
5.下面是一个线性方程组：
（1）求方程的解。
（2）将方程右边向量元素b3改为再求解，并比较b3的变化和解的相对变化。
（3）计算系数矩阵A的条件数并分析结论。
formatrat%用分数格式显示
A=[1/2,1/3,1/4;
1/3,1/4,1/5;
1/4,1/5,1/6]
1/21/31/4
1/31/41/5
1/41/51/6
format%恢复默认格式
b=[;
]
b=
X=A\b
X=
b2=[;
b2=
X2=A\b2
X2=
D=cond（A）
矩阵的条件数决定矩阵的性能，条件数越接近于1其性能越好，系数矩阵A的条件数为+003，和1相差很大，则其性能不好。
因此b矩阵个别元素的微小变动，对方程的解（X值）影响很大。
6.建立A矩阵，试比较sqrtm（A）和sqrt（A），分析它们的区别。
实验三选择结构程序设计
1.求分段函数的值。
用if语句实现，分别输出x=,,,,,,时的y值。
程序设计：
clear;
x=[];
ifx<
0&
x~=-3
y=x.^2+x-6;
elseifx>
x<
5&
x~=2&
x~=3
y=x.^2-5*x+6;
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语句
score=input（'
请输入百分制成绩：
'
a=fix（score/10）;
ifa==9|a==10
grade='
elseifa==8
B'
elseifa==7
C'
elseifa==6
D'
E'
grade
（二）用if语句
ifscore>
=90&
score<
=100
elseifscore>
=80&
90
=70&
80
=60&
70
（三）用switch语句
switchfix（score/10）
case{9,10}
g
升级会员 