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)