Cómo utilizar el lenguaje R para escribir la fórmula de interpolación de Newton para interpolar valores faltantes
Polinomio de Lagrange?lt;-?function(x,y)?{ len?=?length(x) if(len?!=?length(y))
stop( "¿longitud? ¿no? igual!") if(len?lt;?2)
stop("dim?size?must?more?than?1") #pretreat?data?abd?alloc? memery xx?lt;-?paste("(","a?-",x,")") m?lt;-?c(rep(0,len)) #combin?express for(i?in? 1:len)?{
td?lt;-?1
tm?lt;-?"1"
for(j?in?1 :len)?{ if(i?=?j)?{
td?lt;-?td*(x[i]?-?x[j])
tm?lt;-?paste(tm, "*", xx[j]) }
}
tm?lt;-?paste(tm, "/", td )
m[i]lt;-tm?#m[i]?lt;-?parse(text=tm) } #combin?the?exrpession m?lt;-?paste(m , "*",y) r?lt;-?paste(m,collapse=" ") #combin?la?función fbody?lt;-?paste("{?return(",r,")}") f ?lt;-?function(a)?{} #llenar?el?cuerpo de la?función body(f)?lt;-?parse(text=fbody) return(f)
} p >
¿Este es el algoritmo de interpolación polinómica de Lagrang? Consúltelo
.