補間と近似.ppt

ラグランジュの補間法（ｎ次関数） y0 (x – x1)(x – x2)??? (x – xn) y = + ???? (x0 – x1)(x0 – x2)??? (x0 – xn) yn (x – x0)(x – x1)???(x – xn) + (xn – x0)(xn – x1)???(xn – xn-1) yi (x – x0)???(x – xi)???(x – xj)???(x – xn) + (xi – x0)???(xi – xi)???(xi – xj)???(xi – xn-1) i = j のときなし！ プログラミング化 n yi (x – x0)???(x – xi)???(x – xj)???(x – xn) y = S i = 0 (xi – x0)???(xi – xi)???(xi – xj)???(xi – xn-1) N-1 (x-x[0]) (x-x[j]) (x-x(N)) = S y[i] ?? ?? i = 0 (x[i]-x[0]) (x[i]-x[j]) (x[i]-x[j]) #include stdio.h #include math.h #define N 6 double x[N]={0.0, 1.0, 2.0, 3.0, 3.1, 5.0}; double y[N]={0.0, 1.1, 2.5, 4.0, 4.1, 5.0}; double lagrange( double ); int main( ){ double xx, yy; printf(XX\t\tYY\n); for(xx=0.0; xx=5.0; xx+=.2) { yy = lagrange( xx ); printf(%8.2lf\t%8.2lf\n,xx, yy); } return 0; } double lagrange( double xx ) { double z[N]; double yy = 0.0; int i, j; for(i=0;iN;i++) { z[i] = 1.0; for(j=0;jN;j++) if(i!=j) z[i] *= (xx - x[j])/(x[i] - x[j]); yy += z[i]*y[i]; } return yy; } N-1 (x-x[0]) (x-x[j]) (x-x(N)) y = S y[i] ?? ?? i = 0 (x[i]-x[0]) (x[i]-x[j]) (x[i]-x[j]) スプライン補間 各区間を３次関数で近似する。 　各区間の端点で、 １．値が関数と同じ ２．傾きが隣の区間と同じ ３．曲率が隣の区間と同じ x y 異なった近似式 スプライン補間 x y y = a10 + a11x + a12x2 + a13x3 y = a20 + a21x + a22x2 + a23x3 y = a30 + a31x + a32x2 + a33x3 ① ② ③ ① ② ③ 未知数 １２個 スプライン補間 値が関数と同じ 条件　６個 傾きが隣と同じ x y ① ② ③ 曲率が隣と同じ 条件　２個 条件　２個 両端は曲率ゼロ 条件　２個 計１２個 エクセルで多項式近似とスプラインを比較 x y ０ ２ １ １ ２ ３ ３ ４ 多項式近似とスプライン 多項式近似 スプライン 関数補間と近似式 数個、数十個のデータ点から その特徴をつかむ 関数の近似法 補間 　データ点上を必ず通過する近似 近似（最小２乗法） 　データ点に近い点を通る近似 x y x y エクセルで求める x y ０ ２ １ １ ２ ３ ３ ４ エクセルによる補間法 エクセルによる補間法 エクセルによる補間法 エクセルによる補間法 エクセルによる補間法 エクセルによる補間法 エクセルによる補間法 補間法 データ点上を必ず通過する近似 　例：多項式近似　データ数４点　→　３次式 y = a0 + a1x + a2x2 + a3x3