冪法求多項式方程的模大根matlab實現(xiàn)

1、冪法求多項式方程的模大根matlab實現(xiàn)作者:日期:幕法求多項式方程的模最大根matlab實現(xiàn)要求：利用 matlab編寫通用子程序，利用幕法求多項式方程的解：f (x) = xnan/Xn1 ax a0 = 0思想：i首先要將多項式轉(zhuǎn)化成矩陣形式。通過老師上課講的內(nèi)容。將上述多項式轉(zhuǎn)化成為如家格 式的矩陣：0 - a。1 - a?01 -an此矩陣的特征值，就是上述多項式的解。2幕法的思想就不多介紹了，書上講的很詳細，主要運用書上6.2.6的迭代公式:yk = Auk,A二m(k),yk的模最大分量Uk 二 y, %,實驗代碼：詳見附錄 i實驗結(jié)果：(代碼詳見附錄)(i) x3 x2 -5x

2、 3 = 08» jiij Xj spovermethodOj a, le-5) m =-3. 0000-O 40S20. 8165-0. 4082解： 其中m是模最大特征值，x是m對應的特征向量，s是迭代次數(shù)15。精度為1e-5(ii) x3 -3x -仁 0結(jié)果：>> 皿 s=povermthad(3,為le5)n =1. 87940. 24250.G5650.4557其中：m是模最大特征值（多項式模最大根），x是m對應的特征向量，s是迭代次數(shù)為57, 精度為1e-5./.、X8 101x7 208.01X6 10891.01X5 9802.08x4 (iii)+ 7

3、9108.9x-99902x2 +790x-1000 = 0結(jié)果：>>s=povermethod( S； a, 1-10)-100* 0000-0.00730.0062HL 77550.521S0.06990.08380>. 00030.00081212其中：m是模最大特征值（多項式模最大根），x是m對應的特征向量，s是迭代次數(shù)次，精度為1e-10.結(jié)論：幕法求多項式模最大根的效果還是很不錯的，迭代次數(shù)也不多，收斂比較快。附錄1冪法:fun ctio nm,x,s=powermethod( n,a,eps) %A轉(zhuǎn)化后的矩陣%x0迭代初向量 %l模最大特征值%n為最咼次幕A=

4、zeros (n);M = 500000;%v為主特征向量%迭代步數(shù)限制l = 0;for i=1: nA(i, n)=-a(i);endfor i=2: nfor j=1:n-1if i-j=1A(i,j)=1;endendends=0;n=max(size(A);u=ones(n ,1);y=o nes( n,1);%初始化，初始值是多少不重要beta 1=0;eta=no rm(u,2);y=u./eta;u=A*y;beta2=y'*u;while s<=Mif abs(beta2-beta1)/beta1)>epsbeta1=beta2;eta=no rm(u,2);y=u./eta;u=A*y;beta2=y'*u;ends=s+1;if (abs(beta2-beta1)/beta