三次樣条C语言的实现.docVIP

实验二 .三次样条C语言的实现 1.实验目的： 进一步熟悉三次样条插值法。 掌握编程语言字符处理程序的设计和调试技术。 2.实验要求： 已知：某些点的坐标以及初值条件。 输入：条件点数以及这些点的坐标，边界条件 输出：根据三次样条公式求出多项式公式。 3.程序流程： （1）输入已知点的个数； （2）分别输入已知点的X，Y坐标； （3）输入边界条件； （4）通过调用函数，求三次样条函数。 三次样条函数如下： 常见边界条件有三种： 第1种边界条件： 第2种边界条件： 周期边界条件 对于第1，2种条件可写成矩阵形式： 对于第3种条件可写成矩阵形式： 通过以上方程求出Mi 其中： 将求出的解带入： 即得所求函数 程序流程图： ↓ ↓ ↓ ↓ 程序如下： #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); int main(){ int n,i; char ch; do{ cout请输入节点个数:; cinn; n=n-1; for(i = 0; i = n; i++){ cout请输入 Xi:; cinx[i]; //