实验二MATLAB矩阵分析和处理.docxVIP

下载本文档

23
0
约2.46千字
约 10页
2019-07-17 发布于广东
举报
版权申诉

实验二MATLAB矩阵分析和处理.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

实验二MATLAB矩阵分析和处理实验二MATLAB矩阵分析和处理一、实验目的掌握生成特殊矩阵的方法。掌握矩阵分析的方法。用矩阵求逆法解线性方程组。二、实验内容 ?E3?31.设有分块矩阵A???O2?3R3?2?,其中E,R,O,S分别为单位矩阵、随机矩阵、零矩阵和S2?2?? 2对角矩阵，试通过数值计算验证A???E ?OR?RS?o ?2S? 解:gt;gt;E=eye(3)； gt;gt;R=rand(3,2)； gt;gt;O二zeros⑵3)； gt；gt;S=diag(l:2)； gt;gt;A=[ER;OS]； A= 1.0000 gt；gt;AA2 ans 二 1.0000001.90031.457901.0000000001.0000000.95010.23110.60681.000000.4 8600.89130.762102.0000 1.0000 000 01.0000 00 0.46231.21371.0000 2.67392.2863 04.0000 gt；gt;[ER+R*S;OSA2]ans= 1.0000 0000 所以A?? 2 01.0000 000 001.0000 00 1.90030.46231.21371.0000 1.45792.67392.2863 04.0000 ?E?OR?RS? ,S2?? 2.产牛5阶希尔伯特矩阵H和5阶帕斯卡矩阵巳求其行列式的值Hh 和Hp以及他们的条件数Th和Tp,判断哪个矩阵性能更好，为什么？解：H=hilb⑸H二 1.00000.50000.33330.25000.2000gt;gt;P 二 pascal(5)P 二 0.50000.33330.25000.20000.1667 0.33330.25000.20000.16670.1429 0.25000.20000.16670.14290.1250 0.20000.16670.14290.12500.1111 1 1 1 1 11234513610151410203515153570 gt;gt;Hh=det(H) Hh= 3.7493e-012 gt;gt;Hp=det(P) Hp= 1 gt;gt;Th=cond(H) Th= 4.7661e+005 gt;gt; Tp=cond(P) Tp= 8.5175e+003 由上式看出：帕斯卡矩阵性能更好，因为行列式为1,而希尔伯特矩阵条件数很差，使用一般方法求逆矩阵会因原始数据的微小变动产生不可靠的结果建立一个5x5矩阵，求它的行列式的值、迹、秩和范数解：a=magic(5)235 46 1012 1118 gt;gt;b=det(a) b= 5070000 gt;gt;c=trace(a) 65 gt;gt;d=rank(a) d= gt; gt;e=norm(a) e= 65.0000 gt;gt；71319251420212162239 ??29618???4.已知A?20512,求特征值和特征向量，并分析其数学意义?????885?? 解:A=[-29618;20512;-885]A= -29618 20512 ?885 gt； gt;b=eig(A) b= ? 25.3169 -10.5182 16.8351 gt；gt;[v,b]=eig(A) v= 0.71300.28030.2733-0.6084-0.78670.87250.34870.55010.4050b= ?25.316900 0-10.518200016.8351 下面是一个线性方程组 ?l/21/31/4??xl??0.95??l/31/41/5??x???0.67????2?????l/41/51/6???? x3????0.52?? 将方程右边向量第三个元素0.52改为0.53,并比较解的变化计算系数矩阵A的条件数并分析结论解:(l)A=[l/21/31/4; 1/31/41/5; 1/41/51/6] A= 0.50000.33330.2500 0.33330.25000.2000 0.25000.20000.1667 gt;gt;C=[0.95;0.67;0.52] C= 0.9500 0.6700 0.5200 gt； gt;B=A\C B= 1.2000 0.6000 0.6000 gt；gt; (2)A=[l/21/31/4; 1/31/41/5; 1/41/51/6] A= 0.50000.33330.2500 0.33330.25000.2000 0.25000.20000.1667 gt;gt;C=[0.95;0.67;0.53] C= 0.9500 0.6700 0.5