C++_数值分析_三次样条插值_自动选取步长梯形法_Romberg求积法_列主元高斯消去法_列主元LU分解法_Jacobi迭.docxVIP

PAGE  PAGE 16 数值分析上机实验报告之样条插值 1.三次样条插值（初值条件1）： P52.9、给定函数y=fx的函数表和边界条件s75=0，s80=0，求三次样条插值函数s(x)，并求f(78.3)的近似值。 函数表 x757677787980y=fx2.7682.8332.9032.9793.0623.153 源代码：yangtiao.cpp #includeiostream.h #includemath.h void main() { int choice = 0; int n = 2; double xx,*x, *y,*a,*b,*a1,*b1,*h,*m; cout请输入插值节点个数n：endl; cin n; x = new double[n]; y = new double[n]; a = new double[n]; b = new double[n]; a1 = new double[n]; b1 = new double[n]; h = new double[n-1]; m = new double[n+1]; cout请输入n个插值的节点(xi，yi)：endl; for (int i = 0 ; i n ; i++) { cinx[i]y[i]; } for (int j = 0 ; j n-1 ; j++) { h[j] = x[j+1] - x[j]; } cout请输入待估点xx：endl; cinxx; cout请选择边界条件：endl; cinchoice; switch(choice) { case 1: { double temp1,temp2; a[0] = 0; a[n-1] = 1; cout请输入边界条件的两个一阶微商值s(x1)与s(xn)：endl; cintemp1temp2; b[0] = 2*temp1; b[n-1] = 2*temp2; break; } case 2: { a[0] = 1; a[n-1] = 0; b[0] = 3/h[0]*(y[1]-y[0]); b[n-1] = 3/h[n-2]*(y[n-1]-y[n-2]); break; } } for (int k = 1 ; kn-1 ; k++) { a[k] = h[k-1]/(h[k-1]+h[k]); b[k] = 3*( (1-a[k])/h[k-1]*(y[k]-y[k-1]) + a[k]/h[k]*(y[k+1]-y[k]) ); } a1[0] = -a[0]/2; b1[0] = b[0]/2; for (int l = 1; ln ;l++) { a1[l] = -a[l]/(2+(1-a[l])*a1[l-1]); b1[l] = ( b[l] - (1-a[l])*b1[l-1] )/( 2+ (1-a[l]) * a1[l-1] ); } m[n] = 0; for (j = n-1 ; j=0; j--) { m[j] = a1[j] * m[j+1] + b1[j]; } //判别xx所在区间并输出结果 cout\n插值结果为：; for(k = 0 ;k n-1 ;k++) { if (x[k] = xx x[k+1] xx ) { double output = 0; output = (1+2*(xx-x[k])/(x[k+1]-x[k]))* pow(((xx-x[k+1])/(x[k]-x[k+1])),2) *y[k] +(1+2*(xx-x[k+1])/(x[k]-x[k+1]))* pow(((xx-x[k])/(x[k+1]-x[k])),2) *y[k+1] +(xx-x[k])* pow(((xx-x[k+1])/(x[k]-x[k+1])),2) *m[k] +(xx-x[k+1])* pow(((xx-x[k])/(x[k+1]-x[k])),2) *m[k+1]; cout output endl; break; } } delete x; delete y; delete a; delete b; delete a1; delete b1; delete h; delete m; } 运行结果截图： 2.三次样条插值（初值条件2）： P52.10、给定函数y=fx的函数表和边