Matlab上机实验答案.docx
- 文档编号:2314812
- 上传时间:2022-10-28
- 格式:DOCX
- 页数:36
- 大小:348.98KB
Matlab上机实验答案.docx
《Matlab上机实验答案.docx》由会员分享,可在线阅读,更多相关《Matlab上机实验答案.docx(36页珍藏版)》请在冰豆网上搜索。
Matlab上机实验答案
Matlab上机实验答案
Matlab上机实验答案
实验一MATLAB运算基础
1.先求下列表达式的值,然后显示MATLAB工作空间的使用情况并保存全部变量。
(1)
>>z1=2*sin(85*pi/180)/(1+exp
(2))
z1=
(2),其中
>>x=[21+2i;5];
>>z2=1/2*log(x+sqrt(1+x^2))
z2=
-+
+-
(3)
>>a=:
:
;
>>z3=(exp.*a)-exp.*a))./2.*sin(a++log(+a)./2)
(>>z33=(exp*a)-exp*a))/2.*sin(a++log(+a)/2)
可以验证z3==z33,是否都为1)
z3=
Columns1through5
+++++
Columns6through10
+++++
Columns11through15
+++++
Columns16through20
+++++
Columns21through25
+++++
Columns26through30
+++++
Columns31through35
+++++
Columns36through40
+++++
Columns41through45
+++++
Columns46through50
+++++
Columns51through55
+++++
Columns56through60
+++++
Column61
+
(4),其中t=0:
:
>>t=0:
:
;
>>z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<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
ans=
1231-3
32884
0671
(2)
>>A*B
ans=
684462
309-72596
154-5241
>>A.*B
ans=
121024
680261
9-13049
(3)
>>A^3
ans=
3722623382448604
247370149188600766
78688454142118820
>>A.^3
ans=
172839304-64
39304343658503
27274625343
(4)
>>A/B
ans=
>>B\A
ans=
(5)
>>[A,B]
ans=
1234-413-1
34787203
36573-27
>>[A([1,3],:
);B^2]
ans=
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)
ans=
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=
00
00
00
0000
0000
>>A^2
ans=
00
00
00
0000
0000
>>[E,R+R*S;O,S^2]
ans=
00
00
00
0000
0000
>>A^2==[E,R+R*S;O,S^2]
ans=
11111
11111
11111
11111
11111
通过验证,矩阵成立。
2.产生5阶希尔伯特矩阵H和5阶帕斯卡矩阵P,且求其行列式的值Hh和Hp以及它们的条件数Th和Tp,判断哪个矩阵性能更好。
为什么
>>H=hilb(5)
H=
>>P=pascal(5)
P=
11111
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)
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]
A=
-29618
20512
-885
>>[V,D]=eig(A)
V=
D=
00
00
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]
A=
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)
D=
+03
矩阵的条件数决定矩阵的性能,条件数越接近于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>=0&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语句
clear;
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语句
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
grade='D';
else
grade='E';
end
grade
(三)用switch语句
clear;
score=input('请输入百分制成绩:
');
switchfix(score/10)
case{9,10}
g
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Matlab 上机 实验 答案
![提示](https://static.bdocx.com/images/bang_tan.gif)