| 查看: 1271 | 回复: 5 | ||
| 本帖产生 2 个 程序强帖 ，点击这里进行查看 | ||
狼馨520金虫 (小有名气)
|
[求助]
求lagrange插值编码
|
|
谁能指教一下lagrange插值的vb 编写方法，我不会，急需！！！ |
回复此楼» 猜你喜欢
一志愿北京化工大学070300 学硕336求调剂
已经有8人回复
-
一志愿东华大学化学070300，求调剂
已经有6人回复
-
318求调剂
已经有7人回复
-
求调剂一志愿海大，0703化学学硕304分，有大创项目，四级已过
已经有6人回复
-
一志愿武理材料工程348求调剂
已经有5人回复
-
求调剂
已经有7人回复
-
275求调剂
已经有6人回复
-
招08考数学
已经有7人回复
-
303求调剂
已经有6人回复
-
306求调剂
已经有4人回复
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 (小学生)
- 金币: 20567.5
- 散金: 4751
- 红花: 29
- 帖子: 7419
- 在线: 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 循环体开始，从１，到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！！这个地方看不很明白，以前自学的ＢＡＳＩＣ都忘记差不多拆． 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
