电气系统仿真课件3 matlab数值计算教学讲义.pptVIP

6 矩阵的逆 例题： 设矩阵A和B满足关系式AB=A+2B 其中，A=[4 2 3;1 1 0;-1 2 3] 求矩阵B。 分析： 1 AB=A+2B =〉（A-2E）B=A 2 所以，B=（A-2E）-1A A=[4 2 3;1 1 0;-1 2 3] B=inv(A-2*eye(3))*A 7. 方阵的行列式 在MATLAB中，求方阵的行列式的函数是det，用法： det (a)。 a=magic(3) a = 8 1 6 3 5 7 4 9 2 det(a) ans = -360 8. 矩阵的特征值 如果AX=λX存在非零解向量，则称λ为A的一个特征值，相应的非零解向量X称为与λ相对应的特征向量。 在MATLAB中，以函数eig计算矩阵的特征值。 用法一： eig(a)，把它的特征值以列向量的形式显 示出来。 用法二：[v,d]= eig(a)，v为特征向量， 求如下矩阵的特征值和特征向量 G=[1 2 0; 2 5 -1; 4 10 -1]; eig(G) ans = 3.7321 0.2679 1.0000 [v,d]=eig(G) v = -0.2440 -0.9107 0.4472 -0.3333 0.3333 -0.0000 -0.9107 -0.2440 0.8944 d = 3.7321 0 0 0 0.2679 0 0 0 1.0000 [v,d]=eig(A) 9. 矩阵的特征多项式 设A为n阶方阵，如果λ和n维非零列向量X使关系式AX= λX成立，那么（A-λE）X=0称为矩阵A的特征方程，| A- λE |称为方阵A的特征多项式。 在MATLAB中，以函数poly计算矩阵的特征多项式。 a=[-2 1 1; 0 2 0;-4 1 3]; poly(a) ans = 1 -3 0 4 poly2sym(poly(a)) %把系数矩阵转为符号表达式 ans = x^3-3*x^2+4 10. 矩阵的秩 设在矩阵A中存在一个不等于0的r阶子式D，且所有（r+1）阶子式全为0，那么称D为方阵A的最高阶非零子式，r为矩阵A的秩。 在MATLAB中，以函数rank计算矩阵的秩。 a=magic(4) a = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1 rank(a) ans = 3 det(a) ans = 0 3.3.2 数组的运算 数组运算符与矩阵运算符的区别便是多一个小黑点。 .+ 加法 .- 减法 .* 乘法 .^ 幂 .\ 左除 ./ 右除 .’ 转置 数组的运算操作都是对元素逐个进行的。 a=ones(3); b=magic(3) b = 8 1 6 3 5 7 4 9 2 a*b ans = 15 15 15 15 15 15 15 15 15 a.*b ans = 8 1 6 3 5 7 4 9 2 b^2 ans = 91 67 67 67 91 67 67 67 91 b.^2 ans = 64 1 36 9 25 49 16