文档介绍：Vzjgqq数值分析实验报告(1)
2
———————————————————————————————— 作者：
———————————————————————————————— 日期：
个人收集整理 勿做商业用途
个人收集整理 勿做商业用途
个人收集整理 勿做商业用途
生命是永恒不断的创造，因为在它内部蕴含着过剩的精力，它不断流溢，越出时间和空间的界限,它不停地追求，以形形色色的自我表现的形式表现出来。
－－泰戈尔
《数值分析》
实 验 报 告 册
姓名:
学号：
专业:
年级:
计算机科学学院
计算机应用教研室
个人收集整理 勿做商业用途
个人收集整理 勿做商业用途
个人收集整理 勿做商业用途
2008 年 春季 学期
目 录
实验一 3
实验二 5
实验三 7
实验四 10
实验五 12
实验六 15
实验七……………………………………………………………………18
个人收集整理 勿做商业用途
个人收集整理 勿做商业用途
个人收集整理 勿做商业用途
实验一
一、课题名称
非线性方程数值解法
二、目的和意义
学会常用的插值方法，求函数的近似表达式,以解决其它实际问题；明确插值多项式和分段插值多项式各自的优缺点;熟悉插值方法的程序编制；如果绘出插值函数的曲线，观察其光滑性。
三、计算公式
Lagrange插值公式：
牛顿插值公式:
四、结构程序设计
程序设计:
#include""
float f(float x）
{
return（（x*x＊x—1)/3); /*牛顿迭代函数*/
}
main()
｛
float x1,x2，eps,d;int k=0;
clrscr()；
个人收集整理 勿做商业用途
个人收集整理 勿做商业用途
个人收集整理 勿做商业用途
printf（”\n input x1="); /＊输入迭代初值*/
scanf（”%f”,＆x1)；
printf（”\n input eps="）; /＊输入求解精度eps*/
scanf(”%f”，&eps）;
do{
k++;
x1=x2；
x2=f(x1）;
printf（”\n ％d %f\n”,k,x2)；
}while(fabs(x2—x1)〉=eps）；
printf（"the root of f(x)=0 is x=%f,k=%d\n"，x2,k)； /＊输出x和迭代次数k＊/
getch(）；
}
五、结果讨论和分析
计算结果分析：
将六种迭代格式分别代入程序试验：
第一种格式：无论何值都无法求出，即发散
第二种格式：初值为任意的x（x2<=1），精度为0。00001
X=-，k=6
其他值为发散.
第三种格式：初值为任意的x（x〉0),精度为0。0001
X=，k=10
其他值为发散。
第四种格式：初值为任意值,
X=-0。347296，k=5
（5)第五种格式：初值为任意值，精度为0。00001
X=—,k=4
第六种格式：初值为任意值，
X=—0。347296，k=4
由此可知不同的初值对公式的计算有影响，当初值不满足函数的收敛条件时，无法计算结果，函数发散。
精度的大小不同也使迭代函数迭代的次数不同,从而影响xn的近似程度