matlab 比较矩阵差异,Matlab矩阵操作

1求特征值

[v,d]=eig(A);

A为你的矩阵，V为特征向量矩阵，D为特征值矩阵

2除法运算

Matlab提供了两种除法运算：左除(\)和右除(/)。一般情况下，x=a\b是方程a*x

=b的解，而x=b/a是方程x*a=b的解。

例：a=[1 2 3; 4 2 6; 7 4 9]

b=[4; 1; 2];

x=a\b

则显示：x=

-1.5000

2.0000

0.5000

如果a为非奇异矩阵，则a\b和b/a可通过a的逆矩阵与b阵得到：

a\b = inv(a)*b

b/a = b*inv(a)

数组除法：

A./B表示A中元素与B中元素对应相除。

3矩阵乘方

运算符：^

运算规则：

(1)当A为方阵，P为大于0的整数时，A^P表示A的P次方，即A自乘P次；P为小于0的整数时，A^P表示A-1的P次方。

(2)当A为方阵，p为非整数时，则 其中V为A的特征向量， 为特征值对角矩阵。如果有重根，以上指令不成立。

(3)标量的矩阵乘方PA，标量的矩阵乘方定义为 式中V，D取自特征值分解AV=AD。

(4)标量的数组乘方P.^A，标量的数组乘方定义为 数组乘方：A.^P：表示A的每个元素的P次乘方。

4矩阵函数

命令 方阵指数

函数 expm

格式 Y = expm(A) %使用Pade近似算法计算eA，这是一个内部函数，A为方阵。

Y=expm1(A) %使用一个M文件和内部函数相同的算法计算eA

Y=expm2(A) %使用泰勒级数计算eA

Y=expm3(A) %使用特征值和特征向量计算eA

命令 矩阵的对数

函数 logm

格式 Y = logm(X) %计算矩阵X的对数，它是expm(X)的反函数。

[Y,esterr] = logm(X)

%esterr为相对残差的估计值：norm(expm(Y)-X)/norm(X)