题目

求助：求三次样条插值函数的C++程序
已知：
α 30° 40° 50° 60° 70° 80° 90°
m 2.73 2.14 1.75 1.49 1.28 1.13 1.00
n 0.49 0.57 0.64 0.72 0.80 0.89 1.00
用三次样条差值函数求对于不同的α所对应的m、n的值,求此三次样条差值函数的C++程序.

提问时间：2021-03-21

答案

#include
#include
using namespace std;
const int MAX = 50;
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 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 = 0; i--)
fxym[i] = fxym[i] - B[i]*fxym[i+1];
}
void printout(int n);
int main(){
int n,i; char ch;
do{
coutn;
for(i = 0; i