¿Qué función debería usarse para llamar a la interpolación de Newton en MATLAB?

El siguiente es el código del archivo de función. Puede crear un archivo de función usted mismo para guardarlo.

Función s = Newton(x, y, x0, nn)

Interpolación de Newton, donde x es el punto de interpolación e y es el valor de su función.

x0 es el valor de la coordenada x del punto que se va a interpolar y nn es el orden del polinomio de interpolación de Newton.

nx = longitud (x);

ny = longitud (y);

Si nx~=ny

Advertencia (" Los vectores x e y deben tener la misma longitud")

Retorno

Fin

m = longitud (x0);

Según la fórmula Calcule cada elemento del punto de interpolación deseado x0.

Porque i=1: m

t = 0.0

j = 1;

YY = y

kk = j;

Y (kk lt=nn)

kk = kk 1;

Para k=kk: nx

YY(k)=(YY(k)-YY(kk-1))/(x(k)-x(kk-1));

Fin

fin

t = YY(1);

Para k=2: nn

u = 1,0;

jj = 1;

mientras(jj lt;k)

u = u *(x0(I)-x(jj));

jj = jj 1;

Fin

t = t YY(k)* u;

Fin

s(I)= t;

Fin