南京邮电大学.docx

数值计算实践I、方程求根一、实验目的熟悉和掌握Newton法，割线法,抛物线法的方法思路，并能够在matlab上编程实现二、，分别用Newton法,割线法，抛物线法求解.
方程的构造方法：
根:方程的根为学号的,le^8)S
ans=，源程序代码为：
functionx=paowuxian(f,x0,x1,x2,e)ifnargin<4,e=1e-4;endx=x2;y=x1;z=x0;i=0;whileabs(x-y)>ei=i+1;h1=y-z;h2=x-y;c1=(feval(f,y)-feval(f,z))/h1;c2=(feval(f,x)-feval(f,y))/h2;d=(c1-c2)/(h1+h2);w=c2+h2*d;xi=x-(2*feval(f,x))/(w+(w/abs(w))*sqrt(wA2-4*feval(f,x)*d));z=y;y=x;x=xi;end在matla做件中执行下歹U语句并得到的最终结果截图:
＞〉■£皿=Lnlin§('耳*+')；
»i)aowuxian(furijCi,L5#1,le-8)
ans=：只改变初值由上述结果可知，，;割线法的初值设为0,;抛物线法的初值设为0,,;牛顿法根据问题1中牛顿法的程序，在matlab软件中执行下列语句并得到的最终结果截图：
»共[1,2,1,-]:
»E;
»
»epsl=l&-8:
»iui=100d；
»)(1割线法根据问题1中割线程序，在matlab软件中执行下列语句并得到的最终结果截图:
^un=ijiline(kJ+2+x*"2-Hk-0_335G):
)>gexianCfunj0^O*4S]e-8)
ans=，在matlab软件中执行下列语句并得到的最终结果截图：
?>fufL-irilin&(3a'3+2*k"24k-):
»paouracianCtuiij0,0*1b-S)
=-8改为1e-50和1e-100观察迭代次数有何变化牛顿法：
根据问题1中的牛顿法的程序，在matlab软件中执行下列语句并得到的最终结果截图：
精度为1e-50时
»w2,b-];
»村3：
»x0=0,5;
»rm=100(J；
»epsl=le-50;
y=(a,n,sOjiin,,eps1)L二7
-100时
»ep^l=le-100:
)>y=Itewto(n_l(a,xD,mijepsCL=7
y-，在matlab软件中执行下列语句并得到的最终结果截图：
精度为1e-50时
»fiwpinline(Jx<3+-2*it*2-kt—:
»gexianffunjOjljla-50)
ans二0,2240精度为1e-100时
,Eun^inlixte(*x*342^k*2+k-0_3356*):
»geKiantfunjOf1,le-100)
ana=,在matlab软件中执行下列语句并得到的最终结果截图：
精度为1e-50时
»fun^mlLnet7s*3+2*k2+x—0,3356T):
»paovuxianCfuTLj0,5,Lle-50)
10
ans--100时
paowjx]_an(£。事(L5』1,le—100)
JL二
10
ans=，当初值定得越靠近初值，迭代次数就越少
在只改变精度时，当精度越来越大时，迭代次数并几乎不变。
综上所述，初值对迭代次数的影响比较大，精度对迭代次数影响不大。
问题二问题描述根据所给的要求，可知待求的方程为：(1)中方法求解这一问题，并利用图解法找到初值，通过观察图像，将newton法初值设为：,割线法初值设为：0,。抛物线法初值设为：0,,。图像见下图：
-
Newton法问题一的牛顿法的求解只适用丁线性方程，所以在问题二中用其他方法来求解方程。
，源程序代码为：
fun