幂法求多项式方程的模最大根matlab实现.doc

幂法求多项式方程的模最大根matlab实现 要求:利用maｔｌab编写通用子程序,利用幂法求多项式方程的解: 　 思想： 首先要将多项式转化成矩阵形式。通过老师上课讲的内容。将上述多项式转化成为如家格式的矩阵： 　 此矩阵的特征值，就是上述多项式的解。 幂法的思想就不多介绍了,书上讲的很详细，主要运用书上６.2.6的迭代公式： 　　 实验代码：详见附录１ 实验结果：(代码详见附录) （i） 解： 其中m是模最大特征值，x是m对应的特征向量，s是迭代次数１5。精度为１ｅ-5 （ｉｉ) 结果： 其中：m是模最大特征值（多项式模最大根）,x是m对应的特征向量，s是迭代次数为57，精度为１e-5. (iｉi) 结果: 其中:m是模最大特征值（多项式模最大根),x是ｍ对应的特征向量，s是迭代次数12次，精度为1ｅ－１0. 结论:幂法求多项式模最大根的效果还是很不错的，迭代次数也不多,收敛比较快。 附录1 幂法: ｆuncｔｉｏｎ ［ｍ,ｘ,s]=powermethoｄ(ｎ,ａ,epｓ) %Ａ转化后的矩阵 ％x0迭代初向量 %ｌ模最大特征值 %ｎ为最高次幂 A=zｅrｏs(n);　　 　 %v为主特征向量 M ＝ 500００0；　 　 　 ％迭代步数限制 　 l =　0； ｆoｒ i＝1:n A(i,ｎ)=-a(i）； end for i=2：n ｆｏr ｊ=1:n-1 　 　 　 if　ｉ-j=＝1 　 　 A(i,j)＝1； 　 　　 enｄ end eｎd s=0； ｎ＝ｍax（sｉze（A))； 　 　u＝ｏnes(n，１)； 　 　y=oneｓ(ｎ,１）；%初始化,初始值是多少不重要 　　　 beta1＝0； eta=ｎoｒm(u,２); 　y＝u./ｅｔａ; ｕ＝A*y; 　 beta2=y*u; while s=Ｍ 　 　if abs（(bｅta2－beta１）／beta1)eps 　　 　 　 bｅｔa1=beta2； 　　　　　eta＝ｎorｍ(u,２); 　 y=ｕ./eta; 　　 u=A＊ｙ; 　　 beta２=y＊ｕ; 　 enｄ 　s＝s+1; 　 iｆ(abｓ（（beta2-beｔａ1)/betａ1)＜=eps) 　 　 breａk； 　 　enｄ end if s=Ｍ 　 　ｍ＝beta2; ｘ＝ｙ; 　 ｅlse 　 m=ｂetａ2; 　 x=y; end