| 查看: 1526 | 回复: 5 | |||
| 本帖产生 2 个 程序强帖 ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
狼馨520金虫 (小有名气)
|
[求助]
求lagrange插值编码
|
||
谁能指教一下lagrange插值的vb 编写方法，我不会，急需！！！ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
2027博士申请-临床检验诊断学、病原生物学、微生物方向
已经有5人回复
-
【长期招聘】复旦大学 AI 制药实验室 博士后招聘
已经有10人回复
-
各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧
已经有4人回复
-
越来越担心硕士毕业后没有竞争力了
已经有15人回复
-
申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请
已经有3人回复
【答案】应助回帖
★ ★ ★ ★ ★
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