Free Hosting

Free Web Hosting with PHP, MySQL, Apache, FTP and more.
Get your Free SubDOMAIN you.6te.net or you.eu5.org or...
Create your account NOW at http://www.freewebhostingarea.com.

Cheap Domains

Cheap Domains
starting at $2.99/year

check

四點拉格朗日插值多項式(I)

最新更新日期: 2011年4月29日

程式需要在 REG Lin 模式下執行,因此在選擇新程式位置後,按 5 1 選用REG Lin模式。

注意: 藍色的英文字為統計模式中的變數(Σx 按 Shift 1 2,Σy 按 Shift 1 → 2,Σxy 按 Shift 1 → 3)。若果不需要記存方程的答案,程式中綠色部份可以不輸入。

第一個程式 (148 / 144 bytes)

ClrStat: ?→A: ?→B: ?→C: ?→D: ?→X: ?→Y:

(Y - D) ÷ (X - C , D - B DT: ?→D: ?→M:

(M - Y) ÷ (D - X→M: (M - Σx) ÷ (D - C→B:

(D - A)-1(B - (X - A)-1(Σx - Σy ÷ (C - A→A◢

B - A(D +X + C→B◢ D + X: M - BAns - A(DAns + X2→C

Y - AX3 - BX2 - AnsX→D

 

第二個程式 (169 / 163 bytes,程式亦可以計算三點拉朗日格插值多項式)

ClrStat: ?→X: ?→Y: ?→D: ?→M: ?→B: ?→A:

(M - Y) ÷ (D - X , B DT: (X - B)-1(Σx - (A- M) ÷ (B - D→M:

?→B: If B≠Σy: Then ?→C:  (D - B )-1(M - (X - B)-1( (A - Y)

÷ (Σy - X) - (C - A) ÷ (B - Σy→A◢ Else 0→A: IfEnd:

M - A(X + D + Σy→B◢ D + X: Σx - MAns + A(ΣyAns + DX→C

Y + X(MD - AΣyD - Σx→D

 

例題1: 求 (0 , - 4), (1 , - 4)、(2 , 2) 及 (3 , 26)四點的拉格朗日插值多項式。

按 Prog 1  再按 0 EXE - 4 EXE 1 EXE - 4 EXE 2 EXE 2 EXE 3 EXE 26 EXE

(顯示x3的係數為2) EXE (顯示x2的係數為 -3)

EXE (顯示x的係數為 1) EXE (顯示常數項為 - 4)

所以拉格朗日插值多項式為 2x3 - 3x2 + x - 4

計算完結後按 MODE 1 返回正常計算模式。

 

例題2: 求過點 (0 , 5), (1 , 4) 及 (-1 , 10)三點的拉格朗日插值多項式。

注意: 只適用於第二個程式

按 Prog 1 再按 0 EXE 5 EXE 1 EXE 4 EXE

- 1 EXE 10 EXE EXE (不輸入數值代表計算三點拉格朗日插值多項式 ,顯示x2的係數為2)

EXE (顯示x的係數為 -3) EXE (顯示常數項為5)

所以拉格朗日插值多項式為  2x2 - 3x + 5

計算完結後按 MODE 1 返回正常計算模式。

 

註1: 若程式包含綠色程式碼,程式執行完成後,按 RCL A 、RCL B、RCL C及RCL D顯示多項式的係數。

註2: 基於計算方法本身限制,輸入的四點中有任何兩點的x座標相同,程式會出現 Math ERROR,不過若果第三點的x值與第四點的x值相等時,第二個程式會自動假定是計算三點拉格朗日插值多項式,因此亦不能輸入第四點的y值。

 

返回 CASIO fx-50FH、fx-3650P II、fx-50FH II及fx-50F PLUS 程式集