Código fuente del centro de práctica

Dame un ejemplo, espero que te ayude. Ejemplo de código 1: Comparación de métodos de interpolación para funciones de interpolación unidimensionales.

clc borrar todo

x = 0:10;

y = sin(x); :10;

Defina el método de interpolación como una matriz de celdas.

strmod={ 'nearest', 'linear', 'spline', 'cubic' };

Defina el icono como una matriz de celdas.

str1b={'(a) método=más cercano','(b) método=lineal',...

(c) método=spline','(d ) método=cubic'};

Figura (1);

Porque i=1:4

yi=interp1(x, y,, strmod { I });

Figura (1);

Gráfico de rama (2, 2, I);

plot(x, y, 'ro', , yi, ' b '), xlabel(str 1b(I));

title(strmod{i}, 'FontWeight', 'Bold');

Fin

p>

Ejemplo 2: Interpolación spline cúbica

x0 = 0:10;

y0 = sin(x0);

x = 0: 25:10;

y=spline(x0, y0, x);

Figura (2); ', x, y, ' k ');

title('spline', 'FontWeight', 'Bold');

Igual que interp1.

Interpolación 2.2D

Se utiliza en campos como el procesamiento de gráficos e imágenes y el ajuste de curvas tridimensionales, y se implementa mediante interp2. El formato general es:

zi = interp2 (x, y, z,, yi, método) x, y - una matriz de variables independientes del mismo tamaño.

xi, yi: matriz de variables independientes de puntos de interpolación

Método: opciones del método de interpolación, 4 tipos.

*Interpolación del punto más cercano: método =' más cercano '

*Interpolación lineal: método =' lineal 'Este método es el método predeterminado de la función interp2.

*Interpolación spline cúbica: método = 'spline'

* Interpolación cúbica: método = 'pchip' o 'cubic'

Ejemplo: interpolación bidimensional Comparación de cuatro métodos.

Clc todo borrado;

[x, y, z] = picos(7);

Figura (3)

Neto Grid(x,y,z);

Título('Origen','FontWeight','Bold');

[xi,yi]=meshgrid(-3:0.2 :3,-3:0.2:3);

z1=interp2(x,y,z,,yi,'más cercano');

z2=interp2(x,y ,z,,yi,'lineal');

z3=interp2(x,y,z,,yi,'spline');

z4=interp2(x, y , z,, yi, 'cubic');

Figura (4)

Gráfico de rama (2, 2, 1)

Cuadrícula (, Fácil, z1)

Título('Reciente', 'FontWeight', 'Bold');

Subtrama(2, 2, 2)

Grid(xi, yi,z2)

Título('lineal','FontWeight','Bold');

Subtrama(2, 2, 3)

p>

Grid(xi,yi,z3)

title('spine', 'FontWeight', 'Bold');

Subplot(2, 2, 4)

Grid(,yi,z4)

Título('Cubo', 'FontWeight', 'Bold' Result

);