1求特征量
[v，d]=EIG(a )；
a是你的矩阵，v是特征向量矩阵，d是特征值矩阵
2除法
Matlab有两种除法：左除法 (和右除法) /。 一般来说，x=ab是方程a*x
=b的解，而x=b/a是方程x*a=b的解。
例如： a=[1(2) 3； 4 ) 6； 7(9) )。
b=[4； 1； 2 )；
x=ab
显示x=
-1.5000
2.0000
0.5000
在a是非奇异矩阵的情况下，ab和b/a是由a的逆矩阵和b阵列得到的。
ab=inv(a ) b
b/a=b*inv(a ) ) )
数组除法：
A./B表示a中元素和b中的元素相对应。
三矩阵幂
运算符： ^
运算规则：
)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的每个元素的p次方。
四矩阵函数
命令方阵指数
函数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 ) x )