文档介绍：《数值分析》
实验报告
学院：计算机科学与软件学院
姓名：XXX
班级：计算机XX班
学号：XXXXXX
实验一：舍入误差与数值稳定性
实验目的：
通过上机编程，复****巩固以前所学程序设计语言；
通过上机计算，了解舍入误at x[],float y[])
(
float l=0,w;
int i,j;
for (i=0;i<=n;i++)
(
w=1;
for(j=0;j<=n;j++)
(
if(i!=j)
w=w*((a-x[j])/(x[i]-x[j]));
)
l=l+w*y[i];
)
return l;
)
实验结果：
六晨"Cz \Progra^ Files\HicrDsof t Visual. S-tudio\ByPro j ect sYZlDebu^VS1. exe
超目二技下^数巡
：- -
：1-0
作二次插值，并来乂＞2 , 乂二瓦亡2 ,居时的函数近似值，
知劭蜘
函函函
及及及
翻点点
:::
数数数
函期期
5M
及及及
点点点
12 3
戮数时的的的
一*rHltlt
数数数U 个*丁丁值七一r_"Tth- Tth-x 的的的的
*要要要要酿要
\.\.\.\\t_5\
?????? 22
:■9: jnt-同>_k_同数数数 S函$ 女二二节节节 12 3 第第第5
T
Ms:2 十个个值 3 3 3制 rFrpferp 用匚工厂工厂X 的的的的 ..
(83入人入入科入 £国 “+
2、牛顿插值
实验内容：
按下列数据
xi

V\

3

作五次插值，并求x 1=0,46, X2=, X3 =。
实验源程序：
#include<>
#define M 6
float w(float X,int n,float x[]);
float quotient(int k,int i,float f[][M],float x[],float y[]);
float newton(float X,int n,float f[][M],float x[],float y[]);
void main()
float x[M],y[M],f[M][M];
x[0]=;y[0]=;
x[1]=;y[1]=;
x[2]=;y[2]=;
x[3]=;y[3]=;
x[4]=;y[4]=;
x[5]=;y[5]=;
printf("x:
'n"); printf("y:
\n");
printf(" 做五次插值，并求x=,x=,x= 时的函数近似
值.\n\n");
for(int i=0;i<M;i++)
f[0][i]=y[i];
float N[3];
N[0]=newton(,5,f,x,y);
N[1]=newton(,5,f,x,y);
N[2]=newton(,5,f,x,y);
printf("x= 时函数的近似值为 %f\n",N[0]);
printf("x= 时函数的近似值为 %f\n",N[1]);
printf("x= 时函数的近似值为 %f\n",N[2]);
)
float w(float X,int n,float x[])
(
float w=;
for (int i=0;i<n;i++) w=w*(X-x[i]);
return w;
)
float quotient(int k,int i,float f[][M],float x[],float
y[])
(
if(k==0)
f[0][i]=y[i];
else
f[k][i]=(quotient(k-1,i,f,x,y)-quotient(k-1,i-1,f,x,y))/ (x[i]-x[i-k]);
return f[k][i];