Matlab第2讲新.pptVIP

3、矩阵的特征值运算 利用eig和eigs两个函数来进行矩阵的特征值运算。 在MATLAB中，计算矩阵A的特征值和特征向量 的函数是eig和eigs，常用的调用格式有3种： 1) E=eig(A)：求矩阵A的全部特征值，构成列向量E。 2) [V,D]=eig(A)：生成两个矩阵V和D。矩阵A的全部特征值，构成对角阵D，A的特征向量构成V的列向量。 3) eigs函数使用迭代法求解矩阵六个最大的特征值及相应的特征向量。 例18 求矩阵X和A的特征值和特征向量。 X=magic(3), A=[1 0 0;0 0 3;0 9 0] X = 8 1 6 3 5 7 4 9 2 A = 1 0 0 0 0 3 0 9 0 E=eig(X) E = 15.0000 4.8990 -4.8990 [V,D]=eig(X) V = -0.5774 -0.8131 -0.3416 -0.5774 0.4714 -0.4714 -0.5774 0.3416 0.8131 D = 15.0000 0 0 0 4.8990 0 0 0 -4.8990 D=eigs(A) D = -5.1962 5.1962 1.0000 [V,D]=eigs(A) V = 0 0 1.0000 0.5000 0.5000 0 -0.8660 0.8660 0 D = -5.1962 0 0 0 5.1962 0 0 0 1.0000 5 线性系统与矩阵的特征值 4、矩阵的逆运算与伪逆运算 求矩阵逆运算的函数为 inv 如： B = inv(A) 如果矩阵A不是一个方阵，或者A是一个非满 秩的方阵时，矩阵A没有逆矩阵，但可以找到 一个与A的转置矩阵A’同型的矩阵B，使得：A·B·A=A ，B·A·B=B 此时称矩阵B为矩阵A的伪逆，也称为广义逆 矩阵。求一个矩阵伪逆的函数是pinv(A)。 5 线性系统与矩阵的特征值 5、矩阵的行列式运算 求矩阵行列式运算的函数为 det 如： x = det(A) 5 线性系统与矩阵的特征值 6、矩阵的秩 矩阵的秩通过函数rank来求得。 如：r=rank(T) 7、矩阵的迹-用函数trace求得。 矩阵的迹是指矩阵主对角线上所有元素的和。 4 符号微积分初步 3、符号方程的生成： 例8 使用sym函数生成符号方程。 equation1=sym(sin(x)+cos(x)=1) equation1 = sin(x)+cos(x)=1 4 符号微积分初步 符号微积分中常用的函数： 1、findsym函数用于寻找符号变量： 该函数用于找出一个表达式中存在哪些符号变量。注：使用findsym(S,n)可以找出表达式S中n个与x接近的变量。 2、limit(F,x,a)命令用来计算符号表达式当x→a时F=F(x)的极限值。 4 符号微积分初步 例9 使用findsym函数寻找符号表达式中的符 号变量。 syms a alpha b x1 y findsym(alpha+a+b) ans = a, alpha, b 4 符号微积分初步 findsym(cos(alpha)*b*x1+14*y,2) ans = x1,y findsym(cos(alpha)*b*x1+14*y,3) ans = x1,y,b 4 符号微积分初步 findsym(cos(alpha)*b*x1+14*y,4) ans = x1,y,b,alpha findsym(y*(4+3*i)+6*j) ans = y 4 符号微积分初步 3、 limit(F)命令，使用命令findsym(x)确定F中的自变量，设为变量x，再计算当x→0时F的极限值。 4、limit(F,x,a, ‘right’)或limit(F,x,a, ‘left’)命令用来计算符号函数F的单侧极限：左极限x→a-或右极限x→a+。 4 符号微积分初步 sy