| 6 | 1/1 | 返回列表 |
| 查看: 1266 | 回復(fù): 5 | ||
| 本帖產(chǎn)生 2 個(gè) 程序強(qiáng)帖 ,點(diǎn)擊這里進(jìn)行查看 | ||
狼馨520金蟲 (小有名氣)
|
[求助]
求lagrange插值編碼
|
|
誰(shuí)能指教一下lagrange插值的vb 編寫方法,我不會(huì),急需。。 |
至尊木蟲 (著名寫手)
驃騎將軍

木蟲 (著名寫手)

金蟲 (小有名氣)
版主 (知名作家)
|
N 整型變量,輸入?yún)?shù),節(jié)點(diǎn)個(gè)數(shù) XA() 輸入?yún)?shù),存放插值節(jié)點(diǎn) YA() 輸入?yún)?shù),存放函數(shù)值 X 輸入?yún)?shù),插值自變量 Y 輸出參數(shù),所求值 DY 輸出參數(shù),誤差估計(jì) Sub POLINT(XA(), YA(), N, X, Y, DY) Dim C(10), D(10) NS = 1 DIF = Abs(X - XA(1)) For I = 1 To N DIFT = Abs(X - XA(I)) If DIFT < DIF Then NS = I DIF = DIFT End If C(I) = YA(I) D(I) = YA(I) Next I Y = YA(NS) NS = NS - 1 For M = 1 To N - 1 For I = 1 To N - M HO = XA(I) - X HP = XA(I + M) - X W = C(I + 1) - D(I) DEN = HO - HP If DEN = 0# Then Print "PAUSE" Exit Sub End If DEN = W / DEN D(I) = HP * DEN C(I) = HO * DEN Next I If 2 * NS < N - M Then DY = C(NS + 1) Else DY = D(NS) NS = NS - 1 End If Y = Y + DY Next M End Sub |

至尊木蟲 (知名作家)
|
Sub POLINT(XA(), YA(), N, X, Y, DY) 子程序 Dim C(10), D(10) 定義兩個(gè)數(shù)組 NS = 1 DIF = Abs(X - XA(1)) ! abstract 絕對(duì)值. For I = 1 To N 循環(huán)體開始,從1,到n, 這個(gè)n賦值是由外部傳入,可以看子程序?qū)?yīng)變量. DIFT = Abs(X - XA(I)) If DIFT < DIF Then 如果dift小于dif NS = I DIF = DIFT End If 。 nd if 是 if ...endif 不能缺少. C(I) = YA(I) D(I) = YA(I) Next I Y = YA(NS) NS = NS - 1 For M = 1 To N - 1 For I = 1 To N - M HO = XA(I) - X HP = XA(I + M) - X W = C(I + 1) - D(I) DEN = HO - HP If DEN = 0# Then!!這個(gè)地方看不很明白,以前自學(xué)的BASIC都忘記差不多拆. Print "PAUSE" Exit Sub End If DEN = W / DEN D(I) = HP * DEN C(I) = HO * DEN Next I If 2 * NS < N - M Then DY = C(NS + 1) Else DY = D(NS) NS = NS - 1 End If Y = Y + DY Next M End Sub |

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +6 | 小材化本科 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +10 | SUICHILD 2026-03-12 | 10/500 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 10/500 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|