[求助] 求lagrange插值编码

» 猜你喜欢

Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
2025冷门绝学什么时候出结果已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
论文投稿，期刊推荐已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复

1楼2011-07-02 01:33:19

libralibra

至尊木虫 (著名写手)

骠骑将军

程序强帖: 40
应助: 817 (博后)
金币: 12914.1
红花: 64
帖子: 2238
在线: 287.3小时
虫号: 696514
注册: 2009-02-05
专业: 计算机软件

【答案】应助回帖

★
余泽成(金币+1): 谢谢参与应助！ 2011-07-02 09:45:44
狼馨520(金币+2): 2011-07-03 14:18:59

看看lagrange插值介绍吧,这个很简单了

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com

2楼2011-07-02 02:41:02

文俊点点

木虫 (著名写手)

应助: 30 (小学生)
金币: 2529.5
散金: 1796
红花: 1
帖子: 1151
在线: 500.1小时
虫号: 1219128
注册: 2011-03-02
性别: GG
专业: 文艺美学

【答案】应助回帖

★
jjdg(金币+1): 感谢参与程序语言版讨论 2011-07-03 13:09:14
狼馨520(金币+3): 2011-07-03 14:19:08

c的写过········vb表示不会··········

这是你所拥有的时间·····这是你所能改变的生活·········

3楼2011-07-02 18:10:09

狼馨520

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1217.7
散金: 145
帖子: 57
在线: 63.3小时
虫号: 1238429
注册: 2011-03-19
性别: GG
专业: 能源化工

★
dubo(金币+1): 欢迎常来程序语言版讨论 2011-07-31 13:47:07

不会写帮忙搞定一下吧

4楼2011-07-03 01:04:19

jjdg

版主 (知名作家)

程序强帖: 9
应助: 223 (大学生)
贵宾: 3.926
金币: 51002.6
红花: 238
帖子: 7227
在线: 474.2小时
虫号: 27000
注册: 2003-10-21
专业: 药物代谢与药物动力学
管辖: 程序语言

【答案】应助回帖

★ ★ ★ ★ ★
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 (小学生)
金币: 20544.9
散金: 4751
红花: 29
帖子: 7413
在线: 635.8小时
虫号: 1314366
注册: 2011-06-03
性别: GG
专业: 无机材料化学

【答案】应助回帖

★ ★ ★ ★ ★
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

赞一下(2人)

我生活在一个经常爆发地震的年代

6楼2011-07-04 02:11:21