MATLAB数值域的绘图演算法.docVIP

第四章 数值域的绘图演算法 设A为 n阶复数矩阵, 其数值域(numerical range)为复数子集合 W(A)={x*Ax: x(Cn, x*x=1}. 为人所知, 此集合为包含A的固有值的凸集合. 当n=2时, W(A)为一椭圆盘, 其焦距为两个固有值(, (, 短轴长为 |A|22-|(|2- |(|2 之根号. 本章将讨论W(A)的绘图演算法, 重点是方法二与方法三的MATLAB演算法. 方法一 : Theorem 1 ([7]) Let A be an n-square complex matrix. Then W(A) = ( W(Axy), x,y are o.n. where Axy = [x*Ax x*Ay y*Ax y*Ay] ( M2 应用 Theorem 1, 只需用BASIC程式即可绘W(A), 参阅[7]的BASIC程式: step1 随机取两个orthonormal向量 x, y