3次スプライン補間ですが、3次関数で補間するので既知の3点より、
補間したい値を求めると思っていたのですが、
下記の資料を見ると、既知の2点から値を求めています。
接線を使って、計算しているみたいですが、イマイチ分かりません。
分かりやすく教えていただけないでしょうか?
http://www.caero.mech.tohoku.ac.jp/publicData/Da …
No.1ベストアンサー
- 回答日時:
> 3次関数で補間するので
ここを誤解なさってるようです。そうはなくて、「2階微分まで連続な区分的3次関数」で補間するんです。
「n-1階微分まで連続な区分的n次関数」f(x)というのは、あらかじめ点列 k[j] (j=0,1,…, N) (これらを"knot" (「つなぎ目」という意味)と言います)が決めてあるものとし、また、f[j](x)をn次関数の列 (j=1,2,…,N)として、 k[0] ≦ x ≦ k[N]で定義される関数 f(x) が
(1) k[j-1] ≦ x ≦ k[j] のとき f(x) = f[j](x)
であって、しかも
(2) fのxによる0階微分f(x), 1階微分f'(x), 2階微分f''(x), …, n-1階微分(f^(n-1))(x)がどれもk[0] ≦ x ≦ k[N]の範囲で至る所連続
である、という意味。
もちろん(1)から、k[j-1] < x < k[j] ならf(x), f'(x), f''(x), …, (f^(n-1))(x)が連続なのは自明。だから、xが丁度knot k[1], …, k[N-1]に等しいときにf(x)が(2)を満たすことが、f(x)が(2)を満たすための必要十分条件。つまり、
(f[j]^(r))(k[j+1]) = (f[j+1]^(r))(k[j+1]) (r=0,1,2,…, n-1; j=1,…,N-1) …(★)
ということです。
n次関数ならn+1個の係数がある。(ご質問には「3次関数で補間するので既知の3点」とありますが、勝手な3次関数なら既知の4点が必要。これはラグランジュ補間ですね。)
しかし(2)の制約条件が付くため、f[j](x)は勝手なn次関数では駄目です。関数がN個あるんだから全部で(n+1)N個の係数がある。けれども、k[1], …, k[N-1]のN-1個のknotそれぞれについてn個の拘束条件(★)があるために、自由度は(n+1)N - n(N-1) = N+n だけしかない。この自由度を使って補間するんだから、結局「既知のN+n点 <x[i], y[i]> (i=1,2,…,N+n)を通るようにする」という補間をする。すなわち、
f(x[i]) = y[i] (i=1,2,…,N+n)
と上記の式★とを連立した連立方程式を解くことによってfが決まるんです。
なお問題によっては、一番端のknot k[0], k[N]におけるfのr階微分が予め指定されていたり(たとえば (f^(r))(k[0])=(f^(r))(k[N])=0 (r=2,3,…,n)とか)、あるいはfが周期的( (f^(r))(k[0])=(f^(r))(k[N]) (r=0,1,…,n) )であったり、といった条件が付く場合もあり、このときにはさらに自由度が減ることになります。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 ヒストスプライン平滑化をする際の節点の決め方ついて教えてください。 9 2022/08/08 16:17
- 統計学 確率統計の問題です。 3 2022/04/07 04:39
- 高校 三次関数のグラフにつきまして 3 2022/05/15 11:14
- 数学 2次関数y=ax^2のグラフは点A(4,2)を通っている。y軸上に点BをAB=OB(Oは原点)となる 1 2022/04/08 00:05
- 数学 三次関数 点A(2、a)を通って、曲線y=x^3-3xに三本の接線が引けるようなaの値を求めよ 少し 2 2023/07/02 16:09
- 統計学 統計学、エクセルがわかりません!解答と詳しい解説をお願いします! (1)それぞれの地域別に記述統計量 9 2022/08/21 16:30
- Excel(エクセル) 数列の数値補間 2 2022/10/27 16:38
- 数学 数2Bの数列の問題です。 自分は、 まず数列 an=ar^(n-1)と置き こちらの問題の、y= の 1 2022/07/07 16:26
- 数学 【 数I 放物線と直線の共有点 】 問題 放物線y=x²+ax+bが点(1,1)を通り, 直線y=2 4 2022/07/18 09:57
- 数学 数学1の問題がわかりません。 次の関数において、頂点の座標と、[]内のxの値に対するyの値を求めよ。 3 2023/02/13 00:36
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
これについて質問です。微分を...
-
2階微分d^2y/dx^2を詳しく教え...
-
高校数学が日常で役立つ場面を...
-
3階微分って何がわかるの??
-
授業で「yをxで微分する」とい...
-
二回微分して 上に凸下に凸 が...
-
微分について
-
y=log(logx)の微分について
-
これらの数式を声に出して読む...
-
なぜ微分したら円の面積が円周...
-
微分積分を理解できない人って...
-
位置を微分したら速度?
-
微分積分って何?
-
分母が文字の分数を微分する方...
-
二次関数 y=x^2 を微分すると---
-
dxやdyの本当の意味は?
-
微分積分を習わずに大学に入っ...
-
y^2をxについて微分してください
-
log(1+x)の微分
-
数学の微分の範囲で 増減を調べ...
おすすめ情報