2025年南京邮电大学数值计算课程实践深度解析报告.doc

数值计算实践

=1\*ROMANI、方程求根

试验目的

熟悉和掌握Newton法,割线法,抛物线法的措施思绪，并可以在matlab上编程实現

问題描述

(1).給定一种三次方程,分别用Newton法,割线法,抛物线法求解.

方程的构造措施:

(a)根:方程的根為学号的后三位乘以倒数第二位加1再除以1000.

假设你的学号為则根為141*(4+1)/1000=0.564

(b)方程:以你的学号的后三位数分别作為方程的三次项,二次项,一次项的系数,根据所給的根以及三个系数确定常数项.

例如:

你的学号是则你的方程是x3+4x2+x+a0=0的形式.

方程的根為0.564,因此有

0.5643+4*0.5642+0.564+a0=0,于是a0=-2.

你的方程為x3+4x2+x-2.=0.

(2)假设方程是sinx+4x2+x+a0=0的形式（三个系数分别是学号中的数字），重新处理类似的问題

(3)构造一种五次方程完毕上面的工作.

四次方程的构造:将三次多项式再乘以(x-p*)2得到对应的五次多项式(p*為已经确定的方程的根,显然,得到的五次方程有重根).

(4)将（2）中的方程同样乘以(x-p*)得到一种新的方程来求解

注:(1)Newton法取0.5為初值,割线法以0,1為初值,抛物线法以0,0.5,1為初值,

(2)计算精度尽量地取高.

终止准则：根据来终止

(3)可供研究的问題：

（一）的取值不一样对收敛速度有多大的影响

（二）将注（1）中的初值该為其他的初值，对收敛性以及收敛速度有无影响

（三）?能否求出方程的所有的根

（4）试验汇报的撰写

试验汇报包括的内容：（一）试验目的（二）问題描述（三）算法简介（包括基本原理）（四）程序（五）计算成果（六）成果分析（七）心得体会

算法简介

在本问題中，我們用到了newton法，割线法，抛物线法。

1.Newton法迭代格式為：

当時值与真解足够靠近，newton迭代法收敛，对于单根，newton收敛速度很快，对于重根，收敛较慢。

2.割线法：為了回避导数值的计算，使用上的差商替代，得到割线法迭代公式：

割线法的收敛阶虽然低于newton法，但迭代以此只需计算一次函数值，不需计算其导数，因此效率高，实际问題中常常应用。

抛物线法：可以通过三点做一条抛物线，产生迭代序列的措施称為抛物线法。其迭代公式為：

其中是一阶均差和二阶均差。

收敛速度比割线法更靠近于newton法。

对于本问題的处理就以上述理论為根据。终止准则為：

本題中所有精度取1e-8。

程序计算成果

问題一

根据所給的规定，可知待求的方程為：

牛顿法

建立newton_1.m的源程序，源程序代码為：

functiony=newton_1(a,n,x0,nn,eps1)

x(1)=x0;

b=1;

i=1;

while(abs(b)eps1*x(i))

i=i+1;

x(i)=x(i-1)-n_f(a,n,x(i-1))/n_df(a,n,x(i-1));

b=x(i)-x(i-1);

if(inn)error

return;

end

end

y=x(i);

建立n_f.m的源程序来求待求根的实数代数方程的函数，源程序代码為：

functiony=n_f(a,n,x)%待求根的实数代数方程的函数

y=0.0;

fori=1:(n+1)

y=y+a(i)*x^(n+1-i);

end

建立n_df.m的源程序来方程的一阶导数的函数，源程序代码為：

functiony=n_df(a,n,x)%方程的一阶导数的函数

y=0.0;

fori=1:n

y=y+a(i)*(n+1-i)*x^(n-i);

end

在matlab软件中执行下列語句并得到的最终止果截图：

割线法

建立gexian.m的源程序，源程序代码為

functionx=gexian(f,x0,x1,e)

ifnargin4,e=1e-4;end

y=x0;x=x1;i=0;

whileabs(x-y)e

i=i+1;

z=x-(feval(f,x)*(