| 查看: 1291 | 回复: 5 | ||
| 本帖产生 2 个 程序强帖 ,点击这里进行查看 | ||
狼馨520金虫 (小有名气)
|
[求助]
求lagrange插值编码
|
|
谁能指教一下lagrange插值的vb 编写方法,我不会,急需!!! |
» 猜你喜欢
322求调剂,08工科
已经有4人回复
材料工程281还有调剂机会吗
已经有30人回复
化学070300 求调剂
已经有16人回复
化工学硕294分,求导师收留
已经有12人回复
一志愿哈工大 085600 277 12材科基求调剂
已经有32人回复
268分085602化学工程调剂
已经有31人回复
一志愿安大生物学07初试322、本科二本、调剂求助
已经有11人回复
食品与营养(0955)271求调剂
已经有7人回复
085600材料与化工349分求调剂
已经有9人回复
一志愿中南大学 0855 机械 286 求调剂
已经有10人回复
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













回复此楼
5