数值分析实验误差分析.docVIP

实 验 报 告 课程名称 数值分析实验 实验项目 误差分析 专业班级 姓 名 学 号 指导教师 成 绩 日 期 一、C、C++或MATLAB的计算机。 三、实验原理 根据不同的算法，得到的结果的精度是不一样的。 四、实验内容步骤 求方程ax2+bx+c=0的根，其中a=1，b= -(5×108+1)，c=5×108 采用如下两种计算方案，在计算机上编程计算，将计算结果记录下来，并分析产生误差的原因。 方案一： 方案二： 要求：编写程序实现该算法；调试程序，检查输出结果。 五、实验结果及分析 比起方案一，方案二更为精确。 六、#includestdio.h #includemath.h void main() { float a=1,b=-(500000000+1),c=500000000,q; float x1,x2,x3,x4; q=b*b-4*a*c; x1=(-b+sqrt(q))/2; x2=(-b-sqrt(q))/2; x3=-(b-sqrt(q))/2; x4=c/x1; printf(%f\n,x1); printf(%f\n,x2); printf(%f\n,x3); printf(%f\n,x4); } 佛山科学技术学院 实 验 报 告 课程名称 数值分析实验 实验项目 Lagrange插值 专业班级 09计算机科学与技术2班 姓 名 张国烽 学 号 2009314215 指导教师 成 绩 日 期 一、），最后输入要求的自变量x的值，输出对应的函数值。 二、实验设备 安装有C、C++或MATLAB的计算机。 三、实验原理 插值的基本原理（求解插值问题的基本思路） 构造一个函数y=f(x)通过全部节点，即 (i=0、1、… n) 再用f(x)计算插值，即 拉格朗日（Lagrange）多项式插值 Lagrange插值多项式： 3．牛顿（Newton）插值公式 四、实验内容步骤 1．给定，构造Lagrange插值函数计算。 2．已知4个点的函数值如下表，用Newton插值法求x=0.596时的函数值。 i 0 1 2 3 xi 0.40 0.55 0.65 0.80 yi 0.41075 0.57815 0.69675 0.88811 五、实验结果及分析 1.拉格朗日插值 2.牛顿插值 六、 #includeiostream #includestring #includecmath #define Max 100 using namespace std; float xi[Max],yi[Max],value_x; int n; float fenzi[Max],fenmu[Max][Max]; void getdata() { cinn; for(int i=0;in+1;i++) { cinxi[i]yi[i]; } cinvalue_x; } void init() ， { int i,j; for(i=0;in+1;i++) { fenzi[i]=value_x-xi[i]; } for(i=0;in+1;i++) { for(j=0;jn+1;j++) { fenmu[i][j]=xi[i]-xi[j]; } } } void result() { float sum=1,result=0; for(int k=0;kn+1;k++) { sum=1; for(int j=0;jn+1;j++) { if(j==k) continue; sum*=fenzi[j]/fenmu[k][j]; } sum*=yi