矩阵分析与处理 (2).ppt

3.4矩阵求值3.4.1方阵的行列式

把一个方阵看作一个行列式，并对其按行列式的规则求值，这个值就称为所对应的行列式的值。在MATLAB中，求方阵A所对应的行列式的值的函数是det(A)。A=[1-11;5-43;211]det(A)第31页，共43页，星期日，2025年，2月5日关于矩阵分析与处理(2)第1页，共43页，星期日，2025年，2月5日第3章MATLAB矩阵分析与处理3.1特殊矩阵3.2矩阵结构变换3.3矩阵求逆与线性方程组求解3.4矩阵求值3.5矩阵的特征值与特征向量3.6矩阵的超越函数第2页，共43页，星期日，2025年，2月5日3.1特殊矩阵

1．通用的特殊矩阵

常用的产生通用特殊矩阵的函数有：

zeros：产生全0矩阵(零矩阵)。

ones：产生全1矩阵(幺矩阵)。

eye：产生单位矩阵。

rand：产生0～1间均匀分布的随机矩阵。

randn：产生均值为0，方差为1的标准正态分布随机矩阵。第3页，共43页，星期日，2025年，2月5日例3.1分别建立3×3、3×2和与矩阵A同样大小的零矩阵。

(1)建立一个3×3零矩阵。

zeros(3)

(2)建立一个3×2零矩阵。

zeros(3,2)

(3)设A为2×3矩阵，则可以用zeros(size(A))建立一个与矩阵A同样大小零矩阵。

A=[123;456]%产生一个2×3阶矩阵A

zeros(size(A))%产生一个与矩阵A同样大小的零矩阵第4页，共43页，星期日，2025年，2月5日例3.2建立随机矩阵：

(1)在区间[20,50]内均匀分布的5阶随机矩阵。

(2)均值为0.6、方差为0.1的5阶正态分布随机矩阵。

解：x=20+(50-20)*rand(5)

y=0.6+sqrt(0.1)*randn(5)

第5页，共43页，星期日，2025年，2月5日此外，常用的函数还有reshape(A,m,n)，它在矩阵总元素保持不变的前提下，将矩阵A重新排成m×n的二维矩阵。如，前面的A,reshape(A,1,6)ans=142536第6页，共43页，星期日，2025年，2月5日3.12．用于专门学科的特殊矩阵

(1)魔方矩阵

魔方矩阵有一个有趣的性质，其每行、每列及两条对角线上的元素和都相等。对于n阶魔方阵，其元素由1,2,3,…,n2共n2个整数组成。MATLAB提供了求魔方矩阵的函数magic(n)，其功能是生成一个n阶魔方阵。第7页，共43页，星期日，2025年，2月5日例3.3将101~125等25个数填入一个5行5列的表格中，使其每行每列及对角线的和均为565。解：M=100+magic(5)

M=

117124101108115

123105107114116

104106113120122

110112119121103

111118125102109

第8页，共43页，星期日，2025年，2月5日(2)范得蒙矩阵

范得蒙(Vandermonde)矩阵最后一列全为1，倒数第二列为一个指定的向量，其他各列是其后列与倒数第二列的点乘积。可以用一个指定向量生成一个范得蒙矩阵。在MATLAB中，函数vander(V)生成以向量V为基础向量的范得蒙矩阵。例如，A=vander([1;2;3;5])即可得第9页，共43页，星期日，2025年，2月5日A=11118421279311252551第10页，共43页，星期日，2025年，2月5日(3)希尔伯特矩阵

在MATLAB中，生成希尔伯特矩阵的函数是hilb(n)。

使用一般方法求逆会因为原始数据的微小扰动而产生不可靠的计算结果。MATLAB中，有一个专门求希尔伯特矩阵的逆的函数invhilb(n)，其功能是求n阶的希尔伯特矩阵的逆矩阵。第11页，共43页，星期日，2025年，2月5日例3.4求4阶希尔伯特矩阵及其逆矩阵。

命令如下：

formatrat%以有理形式输出

H=hilb(4)

H=invhilb(4)

第12页，