原创力文档- 1、本文档共73页，可阅读全部内容。
- 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数值计算1
第5章 MATLAB数值计算;5.1 特殊矩阵
5.1.1对角阵与三角阵
1. 矩阵的对角元素
(1)提取矩阵的对角线元素
设A为m×n矩阵，diag(A)函数用于提取矩阵A主对角线元素产生一个具有min(m,n)个元素的列向量。
diag(A)函数还有更进一步的形式diag(A,k)，其功能是提取第k条对角线的元素。
(2)构造对角矩阵
设V为具有m个元素的向量，diag(V)将产生一个m×m对角矩阵，其主对角线元素即为向量V的元素。
diag(V)函数也有更进一步的形式diag(V,k)，其功能是产生一个n×n(n=m+|k|)对角阵，其第k条对角线的元素即为向量V的元素。;A= 1 2 3 4
5 6 7 8
5 9 2 3
diag(A)
ans =
1
6
2;diag([1 2 -1 4])
ans =
1 0 0 0
0 2 0 0
0 0 -1 0
0 0 0 4;例5.1 先建立5×5矩阵A，然后将A的第1行元素乘以1，第2行乘以2，…，第5行乘以5。
命令如下：
A=[17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3;11,18,25,2,19];
D=diag([1,2,3,4,5])
C=D*A ; 2. 矩阵的三角阵
(1)下三角矩阵
求矩阵A的下三角阵的MATLAB函数是tril(A)。
tril(A)函数也有更进一步的一种形式tril(A,k)，其功能是求矩阵A的第k条对角线以下的元素。
(2)上三角矩阵
在MATLAB中，提取矩阵A的上三角矩阵的函数是triu(A)和triu(A,k)，其用法与提取下三角矩阵的函数tril(A)和tril(A,k)完全相同。;C =
17 0 1 0 15
46 10 14 28 32
12 0 39 0 66
40 48 76 84 12
55 90 125 10 95
tril(C)
M =
17 0 0 0 0
46 10 0 0 0
12 0 39 0 0
40 48 76 84 0
55 90 125 10 95; ; A=magic(5)
A= 17 24 1 8 15
23 5 7 14 16
4 6 13 20 22
10 12 19 21 3
11 18 25 2 9
B=100+magic(5)
B = 117 124 101 108 115
123 105 107 114 116
104 106 113 120 122
110 112 119 121 103
111 118 125 102 109; ;5.2 矩阵分析;A =
1 1 1 1
8 4 2 1
27 9 3 1
64 16 4 1
C=A
C =
1 8 27 64
1 4 9 16
1 2 3 4
1 1 1 1; ; ; 如果矩阵A不是一个方阵，或者A不是一个满秩矩阵，矩阵A没有逆矩阵，但可以找到一个与A的转置矩阵A’同型的矩阵B,使得：
ABA=A
BAB=B
此时称矩阵B为矩阵A的伪逆，记作A+;A=[0,0,0;0,1,0;0,0,1];
pinv(A)
ans =
0
文档评论（0）