【数值计算】五次样条插值推导

函数形式

$p_k(x) = a_kx^5 + b_kx^4 + c_kx^3 + d_kx^2 + e_kx + l_k$

其中 $k = 1, 2, 3, ... , n-1$ 。

已知条件

经过的 $n$ 个点： $(x_j, y_j)$ ，满足 $p_{k}(x_{k}) = y_{k}$ ， $p_{k+1}(x_{k+1}) = y_{k+1}$ ，其中 $j = 1,2,3,...,n$ 。

四阶及以下导数连续： $p^{(p)}_{k}(x_{k+1}) = p^{(p)}_{k+1}(x_{k+1})$ ，其中 $p = 0,1,2,3,4$ 。

边界条件： $p^{'}_{1}(x_1) = y^{'}_1$ ， $p^{'}_{n-1}(x_n) = y^{'}_n$ ， $p^{''}_{1}(x_1) = y^{''}_1$ ， $p^{''}_{n-1}(x_n) = y^{''}_n$ 。

推导过程

改写一下 $p_k(x)$ 的形式：

$p_k(x) = y_{k}f_{k}(x) + y_{k+1}f_{k+1}(x) + y^{'}_{k}g_{k}(x) + y^{'}_{k+1}g_{k+1}(x) + y^{''}_{k}h_{k}(x) + y^{''}_{k+1}h_{k+1}(x)$

满足

$\begin{cases} p_{k}(x_{k}) = y_{k} \\ p_{k}(x_{k+1}) = y_{k+1} \\ p^{'}_{k}(x_{k}) = y^{'}_{k} \\ p^{'}_{k}(x_{k+1}) = y^{'}_{k+1} \\ p^{''}_{k}(x_{k}) = y^{''}_{k} \\ p^{''}_{k}(x_{k+1}) = y^{''}_{k+1} \\ \end{cases}$

则有

$\begin{cases} f_{k}(x_{k})=1,&f_{k+1}(x_{k})=0,&g_{k}(x_{k})=0,&g_{k+1}(x_{k})=0,&h_{k}(x_{k})=0,&h_{k+1}(x_{k})=0\\ f_{k}(x_{k+1})=0,&f_{k+1}(x_{k+1})=1,&g_{k}(x_{k+1})=0,&g_{k+1}(x_{k+1})=0,&h_{k}(x_{k+1})=0,&h_{k+1}(x_{k+1})=0\\ f^{'}_{k}(x_{k})=0,&f^{'}_{k+1}(x_{k})=0,&g^{'}_{k}(x_{k})=1,&g^{'}_{k+1}(x_{k})=0,&h^{'}_{k}(x_{k})=0,&h^{'}_{k+1}(x_{k})=0\\ f^{'}_{k}(x_{k+1})=0,&f^{'}_{k+1}(x_{k+1})=0,&g^{'}_{k}(x_{k+1})=0,&g^{'}_{k+1}(x_{k+1})=1,&h^{'}_{k}(x_{k+1})=0,&h^{'}_{k+1}(x_{k+1})=0\\ f^{''}_{k}(x_{k})=0,&f^{''}_{k+1}(x_{k})=0,&g^{''}_{k}(x_{k})=0,&g^{''}_{k+1}(x_{k})=0,&h^{''}_{k}(x_{k})=1,&h^{''}_{k+1}(x_{k})=0\\ f^{''}_{k}(x_{k+1})=0,&f^{''}_{k+1}(x_{k+1})=0,&g^{''}_{k}(x_{k+1})=0,&g^{''}_{k+1}(x_{k+1})=0,&h^{''}_{k}(x_{k+1})=0,&h^{''}_{k+1}(x_{k+1})=1 \end{cases}$

解得

$\begin{cases} f_{k}(x) = -6x^5+15x^4-10x^3+1\\ f_{k+1}(x) = 6x^5-15x^4+10x^3\\ g_{k}(x) = -3x^5+8x^4-6x^3+x\\ g_{k+1}(x) = -3x^5+7x^4-4x^3\\ h_{k}(x) = -0.5x^5+1.5x^4-1.5x^3+0.5x^2\\ h_{k+1}(x) = 0.5x^5-x^4+0.5x^3 \end{cases}$

则

$p_{k}(x) = (-6y_{k}+6y_{k+1}-3s_{k}y^{'}_{k}-3s_{k}y^{'}_{k+1}-0.5s^{2}_{k}y^{''}_{k}+0.5s^{2}_{k}y^{''}_{k+1})t^{5}_{k}+(15y_{k}-15y_{k+1}+8s_{k}y^{'}_{k}+7s_{k}y^{'}_{k+1}+1.5s^{2}_{k}y^{''}_{k}-s^{2}_{k}y^{''}_{k+1})t^{4}_{k}\\ +(-10y_{k}+10y_{k+1}-6s_{k}y^{'}_{k}-4s_{k}y^{'}_{k+1}-1.5s^{2}_{k}y^{''}_{k}+1.5s^{2}_{k}y^{''}_{k+1})t^{3}_{k}+0.5s^{2}_{k}y^{''}_{k}t^{2}_{k}+s_{k}y^{''}_{k}t_{k}+y_{k}$

其中： $s_{k}=x_{k+1}-x_{k}$ ， $t_{k}=(x-x_{k})/s_{k}$ 。

由 $p^{(4)}_{k-1}(x_{k})=p^{(4)}_{k}(x_{k})$ 和 $p^{(3)}_{k-1}(x_{k})=p^{(3)}_{k}(x_{k})$ ，分别可推出以下两条等式

$-\frac{168}{s^{3}_{k-1}}y^{'}_{k-1} -(\frac{192}{s^{3}_{k-1}}+\frac{192}{s^{3}_{k}})y^{'}_{k} -\frac{168}{s^{3}_{k}}y^{'}_{k+1} -\frac{24}{s^{2}_{k-1}}y^{''}_{k-1} +(\frac{36}{s^{2}_{k-1}}-\frac{36}{s^{2}_{k}})y^{''}_{k}+\frac{24}{s^{2}_{k}}y^{''}_{k+1} =\frac{360}{s^{4}_{k-1}}y_{k-1} +(\frac{360}{s^{4}_{k}}-\frac{360}{s^{4}_{k-1}})y_{k} -\frac{360}{s^{4}_{s_k}}y_{k+1}$

$-\frac{24}{s^{2}_{k-1}}y^{'}_{k-1} +(\frac{36}{s^{2}_{k}}-\frac{36}{s^{2}_{k}})y^{'}_{k} +\frac{24}{s^{2}_{k}}y^{'}_{k+1} -\frac{3}{s^{2}_{k}}y^{''}_{k-1} +(\frac{9}{s_{k-1}}+\frac{9}{s_{k}})y^{''}_{k} -\frac{3}{s_{k}}y^{''}_{k+1} =\frac{60}{s^{3}_{k-1}}y_{k-1} -(\frac{60}{s^{3}_{k-1}}+\frac{60}{s^{3}_{k}})y_{k} +\frac{60}{s^{3}_{k}}y_{k+1}$

求解出 $y^{'}_{k-1}$ 、 $y^{'}_{k}$ 、 $y^{'}_{k+1}$ 、 $y^{''}_{k-1}$ 、 $y^{''}_{k}$ 、 $y^{''}_{k+1}$ 代入 $p_k(x)$ 即为所求的五次样条解析式。

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 212,657评论 6赞 492
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 90,662评论 3赞 385
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 158,143评论 0赞 348
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 56,732评论 1赞 284
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 65,837评论 6赞 386
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 50,036评论 1赞 291
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 39,126评论 3赞 410
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 37,868评论 0赞 268
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 44,315评论 1赞 303
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 36,641评论 2赞 327
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 38,773评论 1赞 341
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 34,470评论 4赞 333
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 40,126评论 3赞 317
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 30,859评论 0赞 21
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 32,095评论 1赞 267
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 46,584评论 2赞 362
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 43,676评论 2赞 351

【数值计算】五次样条插值推导

函数形式

已知条件

推导过程

推荐阅读更多精彩内容