WPF 牛頓多項式插值和三次樣條插值法的簡單實現

2022-01-12 05:12:31 字數 1122 閱讀 4241

折騰不少時間

牛頓插值法主要是在於差商,三次樣條主要矩陣,發現自己不會,又回頭整矩陣的,什麼線性方程的解法..簡直頭大

牛頓的差商,其實就是遞迴,雖然計算一次是可以用陣列存起來,但是寫出來的時候是遞迴感覺很爽,就滅有儲存各階段的結果,導致計算量大增,應該有個n方的時間把...

剩下就是三次樣條,對於公式來看起來也是很簡單,但是坑不少,到計算的時候發現自己不會矩陣,回頭整了矩陣,發現解法不是求逆或者伴隨矩陣這樣的計算量極大的方法,而是三角分解,lup分解或者追趕法等。

所以三次樣條除去自身對邊界條件的處理,剩下就是矩陣計算了。特別的,三次樣條在用的時候不是按段輸入值,而是直接輸入值,所以做了下對輸入值的進行查詢段的方法,就是判斷輸入的值屬於哪兩個相鄰的x點之間。

簡單弄了demo。**上沒什麼好說的,有不對請指出

示例

控制**

public

icommand drawtypecommand

//依次刪除插值方法產生的點

for (int i = 0; i < maindata.count; i++)

}switch

(a)

break

;

case

interpolationfunction.三次樣條自然條件插值法:

for (int i = 0; i < copyarry.length - 1; i++)

}break

;

default

:

break

; }

isdrawing = false

; redraw = true

; }, (a) => !isdrawing && maindata.count > 3

); }

}

Matlab 三次樣條插值多項式表示

原文自 matlab spline pp.coefs 如何運用matlab 三次樣條插值的問題,今天做作業,突然想用matlab搞搞。題目如下 清華大學出版社的 數值分析 第5版 p49,20題。x 0.25 0.3 0.39 0.45 0.53 y 0.5 0.5477 0.6245 0.6708...

實驗三 牛頓插值多項式

實驗三 牛頓插值多項式 實驗內容 天安門廣場公升旗的時間是日出的時刻,而降旗的時間是日落時分,根據天安門廣場管理委員會的公告,某年10月份公升降旗的時間如下 日期 1 15 22 公升旗 6 09 6 23 6 31 降旗 17 58 17 36 17 26 根據上述資料構造newton插值多項式,...

三次樣條插值

條件 1 輸入 x y f x 0 leq i leq n 2 要求擬合的曲線 s x 滿足 對於任意的 1 leq i leq n 1 在 x 處一階二階導數連續,s x 也連續,且 s x f x s x f x 求解過程 設 s m 對於區間 x x s x 是 x x 上的線性函式,所以設 ...