关于用Ｒ实现拉格朗日插值的一些问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我在https://bbs.pinggu.org/thread-3076192-1-1.html贴中找到一个用Ｒ实现拉格朗日插值的代码，但把程序拷进Ｒ执行例子并没有得到-11而是得到-inf。感觉这个代码并没有什么错误，求大神指教

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：拉格朗日 拉格朗 pinggu thread Ping 程序

俩地方，一个是这里：
td <- td*(x - x[j])
tm <- paste(tm,"*",xx[j])  
还有就是：

tm <- paste(tm,"/",td)
m<-tm #m <- parse(text=tm)
  
最后就是这样的啦——
   
LagrangePolynomial <- function(x,y) {
  len = length(x)
  if(len != length(y))
    stop("length not equal!")
  
  if(len < 2)
    stop("dim size must more than 1")
  
  #pretreat data abd alloc memery
  xx <- paste("(","a -",x,")")
  m <- c(rep(0,len))
  
  #combin express
  for(i in 1:len) {
    td <- 1
    tm <- "1"
    for(j in 1:len) {
      if(i != j) {
       td <- td*(x - x[j])
        tm <- paste(tm,"*",xx[j])
      }
    }
    tm <- paste(tm,"/",td)
    m<-tm #m <- parse(text=tm)
  }
  
  #combin the exrpession
  m <- paste(m,"*",y)
  r <- paste(m,collapse="+")
  
  #combin the function
  fbody <- paste("{ return(",r,")}")
  f <- function(a) {}
  
  #fill the function's body
  body(f) <- parse(text=fbody)
  
  return(f)
}


a = 4:6
b = c(10, 5.25, 1)
f <- LagrangePolynomial(a,b)
f(18)

   
要等于-11，就是这个样子啦
看瞎了快...