电力系统数字仿真软件应用第3课-2MATLAB数值分析与可视化2013.pdfVIP

电力系统数字仿真软件应用 王西田副教授 x.t.wang@126.com 电信大楼1-322 2013/3/4 上海交通大学电力系统及其自动化教研室 第3章 MATLAB电力系统仿真与分析 —— 3.2 MATLAB数值分析与可视化 介绍MATLAB在数值计算分析与数据结果可视化 方面强大而独特的功能。 参考书目： 入门—— 随便选一本参考书 进阶—— MATLAB用户手册 2013/3/4 《电力系统数字仿真软件应用》 3.2.1 数的相关操作  求余运算 M = mod(X, Y) R = rem(X, Y) 当被除数X为负、除数Y为正时：mod(X, Y) = rem(X, Y) + Y  舍入法 Y = fix(X) 向0舍入 Y = floor(X) 向负无穷舍入 Y = ceil(X) 向正无穷舍入 Y = round(X) 向最接近的数4舍5入 2013/3/4 《电力系统数字仿真软件应用》 3.2.2 线性代数 矩阵属性 rank 矩阵的秩 cond 矩阵的条件数 norm 矩阵的范数 trace 矩阵的迹 sum(diag(A)) 2013/3/4 《电力系统数字仿真软件应用》 矩阵变换  Jordan标准型 [V, J] = jordan(A) -1 J —— Jordan标准型，J=V AV  Hessenberg矩阵 [P, H] = hess(A) H —— Hessenberg矩阵，A=PHP’，P’P=I  零空间 Z = null(A) Z —— 零空间的正交矩阵，AZ=0且Z’Z=I  正交空间 B = orth(A) B —— 正交空间的正交矩阵，B’B=eye(rank(A)) 2013/3/4 《电力系统数字仿真软件应用》 矩阵分解  LU分解 [L,U] = lu(A) 把矩阵分解为上三角矩阵和下三角矩阵：A=LU The LU factorization of A allows the linear system A*x = b to be solved quickly with x = U\(L\b)  奇异值分解 [U,S,V] = svd(A) A = U*S*V’ 其中U矩阵和V矩阵是正交矩阵，S矩阵是 对角矩阵，它的对角元素是A矩阵的奇异值。 2013/3/4 《电力系统数字仿真软件应用》 特征值和特征向量 [V,D]=eig(A) 求方阵A的特征向量V及特征值D，满足A*V=V*D；其 中D的对角线元素为特征值，V的列为对应的特征向量。 D=eig(A) 只返回特征值 [V,D]=eig(A,’nobalance’) A阵元素大小悬殊时适用 [V,D]