东南大学数值分析上机实验题(下).doc

《东南大学数值分析上机实验题(下).doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、数值分析上机报告姓名：学号：2013年12月22日第四章38.（上机题）3次样条插值函数（1）编制求第一型3次样条插值函数的通用程序;（2）已知汽车曲线型值点的数据如下：0123456789102.513.304.044.705.225.545.785.405.575.705.80端点条件为=0.8，=0.2。用所编制程序求车门的3次样条插值函数S(x),并打印出S(i+0.5)(i=0,1,…9)。解：（1）#include#includefloatx1[100],f1[100],f2[99],f3[98],m[100],a[100][101],x

2、,d[100];floatc[99],e[99],h[99],u[99],w[99],y_0,y_n,arr,s;inti,j,k,n,q;voidselectprint(floaty){if((y>0)&&(y!=1))cout<<"+"<

3、1]/2.0*(x-x1[j-1])*(x-x1[j-1])+e[j-1]*(x-x1[j-1])*(x-x1[j-1])*(x-x1[j-1]));break;}returns;}voidmain(){do{cout<<"请输入n值:";cin>>n;if((n>100)

4、

5、(n<1))cout<<"请重新输入整数(1..100):"<100)

6、

7、(n<1));cout<<"请输入Xi(i=0,1,...,"<>x1[i];cout<

8、n<<"n):";for(i=0;i<=n;i++)cin>>f1[i];cout<>y_0>>y_n;cout<

9、(f2[i+1]-f2[i])/(x1[i+2]-x1[i]);//二阶差商for(i=1;i

10、];for(i=1;i<=n;i++)//用追赶法解方程,得M[i]{arr=a[i][i-1];for(j=0;j<=n+1;j++)a[i][j]=a[i][j]-a[i-1][j]*arr/a[i-1][i-1];}m[n]=a[n][n+1]/a[n][n];for(i=n-1;i>=0;i--)m[i]=(a[i][n+1]-a[i][i+1]*m[i+1])/a[i][i];for(i=0;i

11、=(m[i+1]-m[i])/(6*h[i]);for(i=0;i