| 查看: 1261 | 回复: 5 | |||
| 本帖产生 2 个 程序强帖 ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
狼馨520金虫 (小有名气)
|
[求助]
求lagrange插值编码
|
||
谁能指教一下lagrange插值的vb 编写方法，我不会，急需！！！ |
回复此楼» 猜你喜欢
311求调剂
已经有6人回复
-
0856调剂
已经有4人回复
-
材料学硕318求调剂
已经有8人回复
-
307求调剂
已经有5人回复
-
281求调剂
已经有4人回复
-
304求调剂
已经有7人回复
-
材料工程274求调剂
已经有3人回复
-
303求调剂
已经有4人回复
-
求调剂
已经有6人回复
-
课题组接收材料类调剂研究生
已经有4人回复
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
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
10