原创力文档- 1、本文档共23页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验02MATLAB矩阵分析与处理(第3章)
实验02 MATLAB矩阵分析与处理
（第3章 MATLAB矩阵分析与处理）
一、实验目的
1. 掌握生成特殊矩阵的方法。
2. 掌握矩阵分析的方法。
3. 用矩阵求逆法解线性方程组。
4. 掌握MATLAB各种表达式的书写规则以及常用函数的使用。
二、实验内容
1. 分块矩阵
，其中E、R、O、S分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证。
命令窗口中的执行过程：
format compact
E=eye(3);
R=rand(3,2);
O=zeros(2,3);
S=diag([1,2]);
A=[E,R;O,S]
A =
Columns 1 through 4
1.0000 0 0 0.6787
0 1.0000 0 0.7577
0 0 1.0000 0.7431
0 0 0 1.0000
0 0 0 0
Column 5
0.3922
0.6555
0.1712
0
2.0000
A^2
ans =
Columns 1 through 4
1.0000 0 0 1.3575
0 1.0000 0 1.5155
0 0 1.0000 1.4863
0 0 0 1.0000
0 0 0 0
Column 5
1.1767
1.9664
0.5136
0
4.0000
C=[E,R+R*S;O,S^2]
C =
Columns 1 through 4
1.0000 0 0 1.3575
0 1.0000 0 1.5155
0 0 1.0000 1.4863
0 0 0 1.0000
0 0 0 0
Column 5
1.1767
1.9664
0.5136
0
4.0000 2. 希尔伯特矩阵、帕斯卡矩阵及其行列式的值和条件数
产生5阶希尔伯特矩阵H和5阶帕斯卡矩阵P，且求其行列式的值Hh和Hp以及它们的条件数Th和Tp，判断哪个矩阵性能更好（注：条件数越接近1，矩阵的性能越好）。
命令窗口中的执行过程：
H=hilb(5);P=pascal(5);
Hh=det(H)
Hh =
3.7493e-012
Hp=det(P)
Hp =
1
Th=cond(H)
Th =
4.7661e+005
Tp=cond(P)
Tp =
8.5175e+003 3. 求矩阵的行列式值、迹、秩和范数
建立一个5×5矩阵，求它的行列式值、迹、秩和范数。
命令窗口中的执行过程：
A=rand(5);
B=det(A),R=rank(A),T=trace(A),V=norm(A)
B =
0.0249
R =
5
T =
2.0045
V =
2.5684 4. 求A的特征值及特征向量
已知
命令窗口中的执行过程：
A=[-29,6,18;20,5,12;-8,8,5];
[V,D]=eig(A,nobalance)
V =
1.0000 0.3564 0.3132
-0.8533 -1.0000 1.0000
0.4890 0.6992 0.4642
D =
-25.3169 0 0
0 -10.5182 0
0 0 16.8351 5. 解线性方程组
下面是一个线性方程组：
(1) 求方程的解。
(2) 将方程右边向量元素b3改为0.53再求解，并比较b3的变化和解的相对变化。
(3) 计算系数矩阵A的条件数并分析结论。
命令窗口中
文档评论（0）