重根的二阶迭代法matlab,计算方程重根的一个高阶迭代程序

亏l我们知道,用通常的迭代程序(例如Newton一RaPhson程序)去求方程的重根,或者由于程序收敛甚慢而浪费机器的宝贵时间,或者导致程序发散.对代数方程而言,重根使方程具有“病态”特性,这时可能使求得的根值不可信,甚至可能改变根的性质(例如,实根变为复根).因此,人们关注着方程重根的计算.在所有计算方程重根的迭代程序中,基本上可以分为两类:一类是在使用时与根的重数无关而仍有较高的敛速阶数;另一类是在上机前或在机器运行中需要确定所求根的重数才具有较高的敛速阶数.从已知的程序来看,前一类程序的基本思想是将方程的重根化为单根来处理,常用的办法是将求方程f(二)一。的根化为求F(劝~t(劝/t’(幻~o的根.作为这一类迭代程序的例子是xn+,~x,一F(x。)/F’(x。)一x,一了(x。)f‘(x。)f‘(x。)12一f(x,;)f“(x,)(,=O,l,…),(l)其中二。为所求根x*的某一初始近似(下同),程序(l)具有平方敛速‘,,.作为这一类迭代程序的另一个例子是X刀+!~X”x。一丸一,f(x。)卫匕2一_f(x,一;)f‘(x。)f,(x,)f,(x。一l)(。~z,2,…),其中x。,x,是所求根二*的两个相异的初始近似.这个程序具有敛速阶数为1.618+12].作为第二类迭代程序的例子是t‘,3,”,”丸二;~x),一。l(乱)/j,(介)(n~o,1,…),__1_./。_‘、f(心)_~;厂(礼)门(x,)、,几J了牛1孟”I移、J刀多JIHl.12一f‘(xn)2!f‘(x。)f‘(xn)7(2)X左+1~X叮一(n一0,l,…),,(,。+一)于(x,,)f{川)(,n+l)f‘(x,:)f{砚)一f(x。)f{跳‘,(,~o,l,…),(2’)(2,‘)其中程序(2)具有平方敛速,程序(2’),(2“)都具有立方敛速,。为所求根的重数.*19夕9年l月23日收到.3期陈永昌:计算方程重根的一个高阶迭代程序52我们知道,当求方程f(幻~。的一个单根a*时,下列程序[4l盆”+i=戈.-卫业-万f‘(x。)L每步只需计算两个矛(幻f(x。一了(x。)/l,(x。))一f(x。)Zf(x,一f(x。)/声,(x。))一f(x。)值和一个f’(幻值,但它却具有(n~0,1,2,…)(*)4阶敛速。现设广是方程矛(幻~O的一个m级重根,即f(二)~(x一x*)”g(二)~0,(3)m为一非负整数,g(x*)铃0.我们通过F(幻~f(幻/厂(幻将x*单根化,然后将程序(*)使用于F(幻,即得迭代中序且、程X月十1~X雌一~旦远2.J业应二卫兰兰立F‘(x,)ZF(z。)一F(x。)(4)劣.+1~X.一f‘(x。)f(x。)[f,(x。)l“一了(x。)f,,(x。)f(。,)f’(x。)一f(x。)f,(z,)Zf(二。)f’(x。)一f(x,夕f’(:,)’(。~o,z,2,…).其中F(二)~f(,)/f(二),。,~另,一F(x,)/F‘(x。).程序(的的优点是显然的:它是第一类迭代程序;敛速阶数是4;计算量只比程序(1)、多计算一次F(介),因此,具有较高的计算效能.事实上,程序(钓的计算效能E;~4就1+口1)+口:故对任何函数了(幻,其计算效能总大于程序(1)的计算效能凡表示计算尸(x)(i~1,2)所花的代价._2‘十纷气,其中。,分别程序(4)具有四阶敛速是显然的.为此,我们考虑更为一般的迭代程序:劣”+i~尤.一口F(x。)孟(l+几,)F(:。)一孟,F(z。)F‘(x,)几(l+又,)F(x。)一几:(l+又