matlab矩阵与方程组.pptVIP

下载本文档

5
0
约9.68千字
约 30页
2017-09-10 发布于四川
举报
版权申诉

matlab矩阵与方程组.ppt

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

matlab矩阵与方程组

第四讲矩阵与线性代数计算研究内容：线性方程组的求解与矩阵特征值问题一、矩阵定义二、矩阵生成生成规则生成方式：由元素列表直接生成矩阵。从外部数据文件读入矩阵。用户自编M文件生成矩阵。利用小阵生成大阵。利用系统内部函数生成矩阵。四、矩阵运算矩阵的行列式与转置阵行列式：det(A) A1=[1 2 3;4 5 6;7 8 0]; A10=det(A1) A10 = 27 A2=[1 2 3;4 5 6;7 8 9]; A20=det(A2) A20 = 0 syms c A3=[1 2 3;4 5 6;7 8 c]; A30=det(A3) A30 =-3*c+27 转置矩阵：A’ B1=[1 3 5;2 6 10]; BB=B1,AA=A1 BB = 1 2 3 6 5 10 AA = 1 4 7 2 5 8 3 6 0 矩阵的加、减、乘及乘方运算注：与数学中的运算符大致一样，可开任意有理数根。矩阵求逆与除法运算(满秩方阵) A-1 =inv(A) , A/B=A*inv(B) , A\B=inv(A)*B 矩阵秩与奇异值 rank(A) A1=[1 2 3;4 5 6;7 8 0];A2=[1 2 3;4 5 6;7 8 9]; B1=[1 3 5;2 6 10]; rank(A1) ans = 3 rank(A2) ans = 2 rank(B1) ans = 1 奇异值矩阵范数与条件数范数：norm(A,p) b=[1 2 3];A12=norm(A1,2);A22=norm(A2,2);b2=norm(b,2); A12,A22,b2 A12 = 13.2015 A22 = 16.8481 b2 = 3.7417 条件数：cond(A) K1=cond(A1) K1 = 35.1059 K2=cond(A2) K2 = 3.8131e+016 病态方程组五、矩阵特征值和特征向量 [P,R]=eig(A) 实对称矩阵：A1=[5 -2 0;-2 6 2;0 2 7]; [P,R]=eig(A1) P = 0.6667 0.6667 -0.3333 0.6667 -0.3333 0.6667 -0.3333 0.6667 0.6667 R = 3.0000 0 0 0 6.0000 0 0 0 9.0000 六、矩阵分解与广义逆阵广义逆阵：设A为m*n阶矩阵，如果存在n*m阶矩阵G,满足条件 (1)AGA=A, (2)GAG=G, (3)(AG)*=AG, (4)(GA)*=GA, *表示共轭后再转置，则称G为A的广义逆阵，记为A+ =G。 A+ =pinv(A) A1=[1 2 3;4 5 6;7 8 0]; G1=pinv(A1) G1 = -1.7778 0.8889 -0.1111 1.5556 -0.7778 0.2222 -0.1111 0.2222 -0.1111 LU分解 A=LU，A是方阵，U是上三角阵，L是一个经过行交换的下三角阵。 [L,U]=lu(A) A1=[5 -2 0;-2 6 2;0 2 7]; [L1,U1]=lu(A1) L1 = 1.0000 0 0 -0.4000 1.0000 0