| 查看: 1290 | 回复: 5 | ||
| 本帖产生 2 个 程序强帖 ,点击这里进行查看 | ||
狼馨520金虫 (小有名气)
|
[求助]
求lagrange插值编码
|
|
谁能指教一下lagrange插值的vb 编写方法,我不会,急需!!! |
» 猜你喜欢
材料085601调剂
已经有28人回复
344 材料专业 求调剂211 无地域要求
已经有8人回复
279求调剂
已经有9人回复
求调剂,985材料与化工348分
已经有11人回复
一志愿中南大学 0855 机械 286 求调剂
已经有8人回复
085600材料与化工349分求调剂
已经有6人回复
考研求调剂
已经有5人回复
调剂求收留
已经有30人回复
297工科,求调剂?
已经有9人回复
0854调剂
已经有11人回复
libralibra
至尊木虫 (著名写手)
骠骑将军
- 程序强帖: 40
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件

2楼2011-07-02 02:41:02
文俊点点
木虫 (著名写手)
- 应助: 30 (小学生)
- 金币: 2529.5
- 散金: 1796
- 红花: 1
- 帖子: 1151
- 在线: 500.1小时
- 虫号: 1219128
- 注册: 2011-03-02
- 性别: GG
- 专业: 文艺美学

3楼2011-07-02 18:10:09
狼馨520
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1217.7
- 散金: 145
- 帖子: 57
- 在线: 63.3小时
- 虫号: 1238429
- 注册: 2011-03-19
- 性别: GG
- 专业: 能源化工
4楼2011-07-03 01:04:19
jjdg
版主 (知名作家)
- 程序强帖: 9
- 应助: 223 (大学生)
- 贵宾: 3.926
- 金币: 51005.2
- 红花: 238
- 帖子: 7227
- 在线: 474.2小时
- 虫号: 27000
- 注册: 2003-10-21
- 专业: 药物代谢与药物动力学
- 管辖: 程序语言
【答案】应助回帖
★ ★ ★ ★ ★
dubo(金币+5, 程序强帖+1): 感谢应助 2011-07-31 13:47:46
dubo(金币+5, 程序强帖+1): 感谢应助 2011-07-31 13:47:46
|
N 整型变量,输入参数,节点个数 XA() 输入参数,存放插值节点 YA() 输入参数,存放函数值 X 输入参数,插值自变量 Y 输出参数,所求值 DY 输出参数,误差估计 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 |

5楼2011-07-03 13:43:47
lilac_c
至尊木虫 (知名作家)
- 程序强帖: 1
- 应助: 26 (小学生)
- 金币: 20572.5
- 散金: 4751
- 红花: 29
- 帖子: 7420
- 在线: 636.4小时
- 虫号: 1314366
- 注册: 2011-06-03
- 性别: GG
- 专业: 无机材料化学
【答案】应助回帖
★ ★ ★ ★ ★
dubo(金币+5, 程序强帖+1): 感谢应助 2011-07-31 13:48:11
dubo: 增加了程序的注释。程序的注释其实比程序跟重要!!这是需要鼓励的!!希望大家多写一些程序的注释 2011-07-31 13:54:33
dubo(金币+5, 程序强帖+1): 感谢应助 2011-07-31 13:48:11
dubo: 增加了程序的注释。程序的注释其实比程序跟重要!!这是需要鼓励的!!希望大家多写一些程序的注释 2011-07-31 13:54:33
|
Sub POLINT(XA(), YA(), N, X, Y, DY) 子程序 Dim C(10), D(10) 定义两个数组 NS = 1 DIF = Abs(X - XA(1)) ! abstract 绝对值. For I = 1 To N 循环体开始,从1,到n, 这个n赋值是由外部传入,可以看子程序对应变量. 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!!这个地方看不很明白,以前自学的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楼2011-07-04 02:11:21













回复此楼