求矩阵特征值算法及程序.docxVIP

— PAGE \* Arabic 1 — 求矩阵特征值算法及程序 求矩阵特征值算法及程序简介  1.幂法  1、幂法规范化算法  （1）输入矩阵A 、初始向量)0(μ  ，误差eps ； （2）1?k ；  （3）计算)1()(-?k k A V μ；  （4）)max (,)  max ()1(1)(--??k k k k V m V m ； （5）k k k m V /)()(?μ；  （6）如果eps m m k k m,m2=x[[k]];m=m1],  {k,1,Length[x]}];  m2]  v=a.u;  m0=fmax[u];  m1=fmax[v];  t=Abs[m1-m0]//N;  k=0;  While[teps  v=a.u;  k=k+1;  m0=m1;  m1=fmax[v];  t=Abs[m1-m0]//N;  Print[k=,k, 特征值=,N[m1,10], 误差=,N[t,10]];  Print[ 特征向量=,N[u,10]]];  If[k ≥nmax,Print[迭代超限]]  说明：本程序用于求矩阵A 按模最大的特征值及其相应特征向量。程序执行后，先通过键盘输入矩阵A 、迭代初值向量)0(μ、精度控制eps 和迭代允许最大次数max n ，程序即可给出每次迭代的次数和对应的迭代特征值、特征向量及误差序列，它们都按10位有效数输出。其中最后输出的结果即为所求的特征值和特征向量序列。如果迭代超出max n 次还没有求出满足精度的根则输出迭代超限提示，此时可以根据输出序列判别收敛情况。  程序中变量说明  a:存放矩阵A ；  u:初始向量)0(μ和迭代过程中的向量)(k μ及所求特征向量；  v:存放迭代过程中的向量)(k V ；  m1:存放所求特征值和迭代过程中的近似特征值；  nmax:存放迭代允许的最大次数；  eps:存放误差精度；  fmax[x]: 给出向量x 中绝对值最大的分量；  k:记录迭代次数；  t1:临时变量；  注：迭代最大次数可以修改为其他数字。  3、例题与实验  例1. 用幂法求矩阵????  ? ??=9068846544  1356133A 的按模最大的特征值及其相应特征向量，要求误差410-m,m2=x[[k]];m=m1],{k,1,Length[x]}];  m2];  v=a.u;  a1=Inverse[a];  m0=fmax[u];  m1=fmax[v];  t=Abs[m1-m0]//N;  k=0;  While[teps  v=a1.u;  k=k+1;  m0=m1;  m1=fmax[v];  t=Abs[m1-m0]//N;  t1=Abs[1/m1-1/m0]//N;  Print[k=,k, 特征值=,N[1/m1,10], 误差=,N[t1,10]];  Print[ 特征向量=,N[u,10]]];  If[k ≥nmax,Print[迭代超限]]  说明：本程序用于求矩阵A 按模最小的特征值及其相应特征向量。程序执行后，先通过键盘输入矩阵A 、迭代初值向量)0(μ、精度控制eps 和迭代允许最大次数max n ，程序即可给出每次迭代的次数和对应的迭代特征值、特征向量及误差序列，它们都按10位有效数输出。其中最后输出的结果即为所求的特征值和特征向量序。如果迭代超出max n 次还没有求出满足精度的根则输出迭代超限提示，此时可以根据输出序列判别收敛情况。  程序中变量说明  a ：存放矩阵A  u ：初始向量)0(μ和迭代过程中的向量)(k μ及所求特征向量  v:存放迭代过程中的向量)(k V  a1：存放逆矩阵1-A  m1: 存放所求特征值和迭代过程中的近似特征值  nmax:存放迭代允许的最大次数  eps:存放误差精度  fmax[x]: 给出向量x 中绝对值最大的分量  k:记录迭代次数  t1:临时变量  注：迭代最大次数可以修改为其他数字。   3、例题与实验  例3. 用反幂法求矩阵????