编程插值(插值法编程)

qiaoqingyi 2023-10-22 141

默认

题目没有指定用哪种插值方法，分别用线性分段三次曲线和样条三种插值方法进行插值并计算和真值的误差使用interp1函数参考代码x=100 121 144y=sqrtxx0=115y0=sqrtx0y1=interp1x，y，x0y2；matlab算么function output = splinex，y，xxSPLINE Cubic spline data interpolationPP = SPLINEX，Y provides the piecewise polynomial form of the cubic spline interpolant to the data values Y at the data。

double a=3，7，2，8int n=4ndouble sum=anwhilen = 0sum=sum*x+anreturn sum double caldouble a，double b，double e double a1，a2，a3，t_a2，ta1=aa3=ba2=a；function FN=FloresS，Sd0，n，Cr，Ri，Rj，Vi，Vj，Eksigai=1Vi^2piEksigaj=1Vj^2piEkSd=diffS 0K=43pisigai+sigaj*Ri*RjRiRj^05FN=K*S^n*1。

曲线插值方法很多你没说用什么方法线性最小二乘法最简单的，你先判断所要求的x在什么范围，比如当x=3，那么你就用4，100，1，048这两个点的线性方程去求用Aitken插值还比较方便编程；matlab中使用插值函数插值函数the function of interpolation interp1调用函数的格式Syntaxyi = interp1x，Y，xiyi = interp1Y，xiyi = interp1x，Y，xi，methodyi = interp1x，Y，xi，method，#39extrap#39yi =。

printfquot输入插值点x的值\nquotscanfquot%lfquot，xprintfquot\nquotfori=0iltni++ forj=0，temp=1jltnj++ ifj！=itemp=temp*xptjxptixptjx Ln=Ln+temp*pt；在form里面添加一个按钮即可 Private Sub Command1_ClickN = 6 #39输入插值节点个数 PI = 36 Dim XA7， YA10XA1 = 400 XA2 = 730 XA3 = 800 XA4 = 980 XA5 = 1200 XA6。

插值法程序

1、计算插值点的位子 if2 * x lt= p1x + p2x n = 1else if2 * x = pNUMBER2x + pNUMBER1x n = NUMBER 2else n = intNUMBER 2printfquot%d\nquot， n。

2、数学内插法即quot直线插入法quot其原理是，若Ai1，b1，Bi2，b2为两点，则点Pi，b在上述两点确定的直线上而工程上常用的为i在i1，i2之间，从而P在点AB之间，故称quot直线内插法quot数学内插法说明点P反映的。

3、代码问题，改一句就好了includeltiostream#define N 3 插值节点数目using namespace stdvoid main float xN 差值节点横坐标 float yN 差值节点纵坐标 float a 所求点横。

4、可能会产生更严重的误差积累，所以，稳定性得不到保证这就是Runge现象解决Runge现象的方法是采用分段低次多项式插值有分段线性插值和分段三次Hermite插值在每个小区间采用低次插值，则可避免Runge现象。

5、printfquot拉各朗日插值后得到的结果是\nquot printfquotf%2f=%e\nquot，t，lag freex freey doubleLAGn，x，y，tintndouble *xdouble *ydoublet inti，j doublep，s s = 0 fori=0iltn1i。

6、const double PI = 3654const double R=PI180 定义常量角度转弧度的因子double *mdouble *y，int length，double j m法 double *jinsidouble y，double g，double a。