牛顿法解方程组和三次样条插值代码.docVIP

已知： 0 0.15 0.30 0.45 0.60 1 0.97800 0.91743 0.83160 0.73529 求满足边界条件： ，的三次样条插值函数。 解：算法步骤为： 步1：始数据，， 和。 步2：从到计算及。 步3：从1到由及式计算。 步4：用追赶法公式解方程，求出。 步5：计算的系数或计算在若干点上的值，并输出结果 下图为输入输出结果： 输入数据 输出数据 二.用Newton法求解非线性方程组： 给定初值，误差极限为。 解：算法步骤为： 方程组 其中均为多元函数。若用向量记号记， 就可写成 当,且中至少有一个是自变量 的非线性函数时，则称方程组为非线性方程组。非线性方程组求根问题是前面介绍的方程即求根的直接推广，实际上只要把单变量函数看成向量函数则可将单变量方程求根方法推广到方程组。若已给出方程组的一个近似根 将函数的分量在用多元函数泰勒展开，并取其线性部分，则可表示为 令上式右端为零，得到线性方程组 其中 称为为雅可比（Jacobi）矩阵。求解线性方程组并记解为 这就是解非线性方程组的牛顿法。 程序执行结果为 附录（程序代码） 程序1：三次样条插值 #includeiomanip #includeiostream using namespace std; const int MAX = 250; float x[MAX], y[MAX], h[MAX]; float c[MAX], a[MAX], fxym[MAX]; float f(int x1, int x2, int x3){ float a = (y[x3] - y[x2]) / (x[x3] - x[x2]); float b = (y[x2] - y[x1]) / (x[x2] - x[x1]); return (a - b)/(x[x3] - x[x1]); } //求差分 void printout(int n){ coutsetprecision(6); for(int i = 0; i n; i++){ couti+1: [x[i] , x[i+1]]\n\t; float t = fxym[i]/(6*h[i]); if(t 0)coutt*(x[i+1] - x)^3; else cout-t*(x - x[i+1])^3; t = fxym[i+1]/(6*h[i]); if(t 0)cout + t*(x - x[i])^3; else cout - -t*(x - x[i])^3; cout\n\t; t = (y[i] - fxym[i]*h[i]*h[i]/6)/h[i]; if(t 0)cout+ t*(x[i+1] - x); else cout- -t*(x[i+1] - x); t = (y[i+1] - fxym[i+1]*h[i]*h[i]/6)/h[i]; if(t 0)cout + t*(x - x[i]); else cout - -t*(x - x[i]); coutendlendl; } coutendl; } void cal_m(int n){ //用追赶法求解出弯矩向量M…… float B[MAX]; B[0] = c[0] / 2; for(int i = 1; i n; i++) B[i] = c[i] / (2 - a[i]*B[i-1]); fxym[0] = fxym[0] / 2; for(i = 1; i = n; i++) fxym[i] = (fxym[i] - a[i]*fxym[i-1]) / (2 - a[i]*B[i-1]); for(i = n-1; i = 0; i--) fxym[i] = fxym[i] - B[i]*fxym[i+1]; } void printout(int n