矩阵本征值并行计算与工程计算的核心问题.docVIP

矩阵本征值并行计算与工程计算的核心问题 一、矩阵本征值问题 矩阵本征值问题是高等数学中的一个古典问题,又是当前迅速发展的计算机科学和数值代数中的一个活跃的研究课题。矩阵本征值问题的提法十分简单,最常见的由矩阵方程Ax=λBx定义。A、B是给定的n行n列方阵。方阵的元素通常是实数或复数。问题是要求解未知量本征值λ和未知的非零本征向量x,并称问题为矩阵偶A、B的广义本征值问题。假若B的对角线上的元素全是1,其他元素是零,方程退化为Ax=λx。 这个方程定义的被称为标准的矩阵本征值问题,简称为矩阵本征值问题。 上述两个定义方程可等价地写为第二种形式的定义方程 (A-λI)x=0,(A-λB)x=0 自然地可写出一般形式的方程: (A0-λA1-λ2A2-,…,-λkAk)x=0 这个方程定义的被称为高次矩阵本征值问题,其中A0,A1,A2,…,Ak是给定的n行n列方阵,求解λ与x。这个一般形式的矩阵本征值问题不单纯是出于逻辑推理和抽象的偏爱,而常常是更能精细地刻画客观的研究对象。 根据高等代数的理论,本征值问题可以等价地用矩阵行列式写成第三种形式的定义方程,常称为本征方程 det(A-λI)=0 det(A-λB)=0 det(A0-λA1-λ2A2-,…,-λkAk)=0 由行列式理论可知以上三个方程可以理解为λ的多项式方程,表面上它只定义了本征值,但是结合第二种形式的方程可知,它也隐含着定义了本征向量,它表现为λ矩阵的零空间是确定的。 二、矩阵本征值问题的应用 矩阵本征值(也称为矩阵特征值)。顾名思义,它是矩阵的一个最重要的本质特征。对于一切用矩阵描述的客观运动状态和因素之间关系的问题,矩阵本征值都将扮演重要角色。因此,矩阵本征值问题除了直接用于解决数学中诸如非线性规划、优化、常微分方程、几何,以及各类数学推算和计算方法中的问题之外,更多地用于科学研究与工程设计的计算工程之中。 本征值常常表现为频率,在建筑工程中它是至关重要的,设计高层建筑、桥梁、水坝时都要使其固有频率在限值之内,以免在承受复杂的动荷载时因发生共振而遭到破坏。在航空航天工程中,本征值要符合严格的要求,否则飞机、飞船、导弹等飞行物会因颤震而毁于瞬间。国民经济系统的本征值起着保持社会经济稳定发展的作用,而在股票、证券投资金融业中,若运作不能把握住市场系统的固有频率本征值,整个经营将惨败无疑。对于生物群体,若本征值超标,导致的后果不堪设想,不是群种灭绝就是胡生乱长,为害四方。 本征值几乎无处不在,以至于苦心探索微观世界分子、原子、电子运动规律的专家们,特别重视矩阵本征值问题。它被用来解决人体基因的密码、高级晶体材料的生长、蛋白质的结构,还有各种新型器材、药品的研制,等等。下面我们举几个最简单的矩阵本征值问题的构成例子。 例1.弦振动 假设受力F作用的弦上等距h分布的质点m1、m2、m3、m4作上下运动。由牛顿定律容易推导出运动方程 d2x D———=Tx dt2 其中,x=(x1,x2,x3,x4)T,di=mih/F,D=diag(d1,d2,d3,d4)T是三对角线矩阵,主对角线元素为-2,次对角线元素为1。 再设运动为标准的振动,即 d2x ———=-ω2x dt2 于是建立了矩阵偶T,D的广义本征值问题 Tx=-ω2Dx 例2.信息系统研究 假设一个信息系统模拟的一组部件(子系统)s1,s2,…,sm要完成一组操作O1,O2,…,Om相应的容量为V1,V2,…,Vm。 效率e=ct/n,其中c表示单位时间、单位重量的费用,t表示时间,n表示操作规模。 系统的效率矩阵可以用E=(eij)表示。eij=cijtij/nij,其中i表示部件号,j表示操作号,系统完成任务的费用是x=Ev,x=(x1,x2,…,xn)T,V=(v1,v2,…,vm)T,于是系统的性能由E,V,x确定。 若以Rayleigh商为度量性能标准,即 a2=(Ev)T(Ev)/vTv 系统的极大费用表示为矩阵ETE的本征值。达到极大值的本征向量v在系统上是一个临界容量。 例3.经济发展研究 假设工厂i对工厂j的销量为bij,对消费者的销售量为yi,总销售量为xi=yi+ibij,若按实际情况定义了输入系统模拟值aij,使得bij=aijxj,就可以建立起简单的Leontief矩阵模型,写出经济动态模拟方程 D[(A-I)x(t)+y(t)]=x(t+1)-X(t) 若中矩阵A=(aij),x=(x1,x2,…)T,y=(y1,y2,…)T,矩阵D是对角阵,对角元称为各工业的反应系数,这个方程中的矩阵(A—I)的本征值的实部大于零,则系统是不稳定的。那是表示输出量按时间的指数增长。 例4.动物群体研究 假设某种动物的最大年龄是n年,按年龄分群为n+1群,假定雌、雄动物数量成固定比例,因此只