| 6 | 1/1 | 返回列表 |
| 查看: 1265 | 回復(fù): 5 | ||
| 本帖產(chǎn)生 2 個(gè) 程序強(qiáng)帖 ,點(diǎn)擊這里進(jìn)行查看 | ||
狼馨520金蟲 (小有名氣)
|
[求助]
求lagrange插值編碼
|
|
誰能指教一下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 ! end 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ā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑(0805) +6 | 煙汐憶海 2026-03-16 | 16/800 |
|
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +4 | Ymlll 2026-03-18 | 6/300 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +10 | yangfz 2026-03-17 | 10/500 |
|
|
[考研] 材料專碩英一數(shù)二306 +4 | z1z2z3879 2026-03-18 | 4/200 |
|
|
[考研] 化工學(xué)碩306求調(diào)劑 +10 | 42838695 2026-03-12 | 10/500 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[基金申請(qǐng)] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動(dòng)化所實(shí)習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|