- 1、本文档共43页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
MATLAB矩阵分析与处理3.1特殊矩阵3.2矩阵的结构变换3.3矩阵求逆3.4矩阵求值3.5矩阵的特征值与特征向量3.6矩阵的超越函数3.1特殊矩阵3.1特殊矩阵常用的产生通用特殊矩阵的函数有:
zeros:产生全0矩阵(零矩阵)。
ones:产生全1矩阵(幺矩阵)。
eye:产生单位矩阵。
rand:产生0~1间均匀分布的随机矩阵。
randn:产生均值为0,方差为1的标准正态分布随机矩阵。3.1.1通用的特殊矩阵zeros(m)zeros(m,n)zeros(size(A))例:设A=[123;456]为2×3矩阵,分别建立3×3、3×2和与矩阵A同样大小的零矩阵。
3.1特殊矩阵解:zeros(3)
zeros(3,2)A=[123;456]zeros(size(A))例建立随机矩阵:
(1)在区间[20,50]内均匀分布的5阶随机矩阵。
(2)均值为0.6方差为0.1的5阶正态分布随机矩阵解:x=20+(50-20)*rand(5)
y=0.6+sqrt(0.1)*randn(5)
3.1特殊矩阵rand与randn的不同常用的函数还有reshape(A,m,n),它在矩阵总元素保持不变的前提下,将矩阵A重新排成m×n的二维矩阵。例:设A=[123;456]为2×3矩阵3.1特殊矩阵》reshape(A,1,6)ans=142536》reshape(A,3,2)ans=142536(1)魔方矩阵
魔方矩阵有一个有趣的性质,其每行、每列及两条对角线上的元素和都相等。对于n阶魔方阵,其元素由1,2,3,…,n2共n2个整数组成。MATLAB提供了求魔方矩阵的函数magic(n),其功能是生成一个n阶魔方阵。3.1.2用于专门学科的特殊矩阵3.1特殊矩阵例将101~125等25个数填入一个5行5列的表格中,使其每行每列及对角线的和均为565。解:M=100+magic(5)
M=
117124101108115
123105107114116
104106113120122
110112119121103
111118125102109
3.1特殊矩阵(2)范得蒙矩阵
范得蒙(Vandermonde)矩阵最后一列全为1,倒数第二列为一个指定的向量,其他各列是其后列与倒数第二列的点乘积。可以用一个指定向量生成一个范得蒙矩阵。在MATLAB中,函数vander(V)生成以向量V为基础向量的范得蒙矩阵。例如,A=vander([1;2;3;5])即可得3.1特殊矩阵A=111184212793112525513.1特殊矩阵(3)希尔伯特矩阵
在MATLAB中,生成希尔伯特矩阵的函数是hilb(n)。希尔伯特矩阵的每个元素为
希尔伯特矩阵是一个条件数很差的矩阵,使用一般方法求逆会因为原始数据的微小扰动而产生不可靠的计算结果。因此,在MATLAB中,有一个专门求希尔伯特矩阵的逆的函数invhilb(n),其功能是求n阶的希尔伯特矩阵的逆矩阵。3.1特殊矩阵例求4阶希尔伯特矩阵及其逆矩阵。解:
formatrat%以有理形式输出
H=hilb(4)H=11/21/31/41/21/31/41/51/31/41/51/61/41/51/61/7
H=invhilb(4)
3.1特殊矩阵H=16-120240
文档评论(0)