数值计算实验-特征值.docxVIP

实验七特征值 姓名：方小开 班级：计科2班 学号：20060810202 【实验目的】 用MATLAB求特征值。 【实验原理】 幕法?主特征值： 是曰的最大分黄 当k充分大时Ubmaxlx^lJ ： 堆程实现： 输入：矩阵…、的初值,最大次数knw,误差epsilon 找出、的最大分量maxvaLk=0; 对y=x/max、?l,计算新一轮的?Ay i+Mx的新的最大分量maxvalnew 假设(maxvalnew-maxvalepsilon 那么结束循环 否那么 m?xval=maxvalnew k-k,1 回到第(2)步 掉法?主特征值?加速 如果能选取适当的冬，使得冬二那么用-4-V代替A来进行迭 代求得模最大的特征值万和特征向量V,那么相对A的模最大的特征值为片+舄, 特征向量为V。 察法?最小特征值 如果A的逆矩阵品存在得到、有.4* = = A的逆矩阵A *的特征值是A矩阵的特征值的倒数 对应的特征向量相同? A-1的模最大的特征值对应于AJ的模最小的特征华 Ax 【实验过程】 籍法?主特征值： function [lambdaz v)=powerchar1(a,x,epsilon,kmax) k=0; maxval=max(abs(x)); while(kkmax) y=x/maxval; ■对x进行归 化得到y I %10.5f %10.5f %10.5f\n,/krx(l)zx(2)/x(3)rmaxval, y(1) zy(2) zy(3) ; x-a*y;二计。卜?轮的x maxvalnew-max (abs (x) ; *沂?轮的IrlXffi if abs(maxvalnew-maxval)epsilon break; end k=k+l; maxval=maxvalnew; end lambda=maxvalnew; v=x/lambda; 运行结果： ? Cl y]-povoroliar 1 x, !? Oo-e, OOOO) 1 ■ 3. O. fAtm O. 640/ 1. oooo 皋法?主特征值?加速： 编程时选择适当值p，对A进行调整，计算完后加上P，可以不在代码中表达，只要运算时换消整一下系数炬阵 a=[l 1 0.S; 1 1 0.25; 0.S 0.25 2] x=[l 1 If循环22次后得结果 p=0.75; b=a-p*diag((l 11]) ll,v]=powercharl(brxz1.0e-6,6000);循环 15 次出结果 ?l=l+p= 2.536S 特征向量也，样 帮法?她小特征值 编程实现1：求A的逆B, b=inv(A)调用powercharl() 如 a=[l 10.5; 11 0.25; 0.5 0.25 2] x=[l 11]*； b=inv(a) (l,v]=power-charl(b#x/1.0e-6,6000); 编程实现方法2：由于求逆问题较多，故采用求方程的形式。而A是对称的，正交的。修改原来的powercharl() 中计算x=a，y为x=a\y或者调用上章的代码x=gau$sl(a,n),此处a是增广矩阵x=gauss2(a,n),x=gau$s」ordnl(a,n)与 gaussjordn function [lambda,v]=powerchar3a,x,epsilon,kmax) k=O;maxval=max(abs(x));n=length(x);b=zeros(n,n+1); for i=l:n for j=l:n b(i,j)=a(izj); end end y - zeros(nz1); while kkmax y=x/maxval; for i=l:n b(irn+l)=y(i); end x=gauss2(br n); maxva1new=max(abs(x)); if abs(maxvalnew-maxval)epsilon break; end k-k+1； maxva1-maxva1new; end lambda-maxvalnew; v=x/lambda; 运行结果： ? a=( 2 1 0; 1 3 l;0 1 4) ?x=(l 11), /100); 得到A的懿小特征值 1.00000 -0.73205 0.26795 归一化后的特征向 lit 幕法求特征ffi-matlab提供的eig(A) 函数格式 eig d = eiu(A) (l = eig(A.B) [\\D] = eig(A) [\：D] = eig(A/i)obalaiKer) %求矩阵A的特征值d,以向量形式存放d. Bd,以向量形式存放d.。。计算A的特征值对角阵D和特征向量V,使A\=\D成立. 。?当矩阵A中有与戡断误差数量级相差