Interpolations

%---------------------------------------------------------
% interpolation linéaire, spline, cubique
% d'aprés 'Apprendre et maitriser Matlab'
%-----------------------------------------------------------
% remarque : en 3D utiliser interp2
%-----------------------------------------------------------

echo on;
%------------------------------------
% interpolations diverses 
% fonction interp1
%------------------------------------
echo off;

clear all;
x=0:4*pi;  		% domaine
z=0:0.25:4*pi; 		% pas inférieur à celui de x

n=length(x);
y=(cos(x)).^2;  		% courbe cos²x
	% interpolation linéaire
y1=interp1(x,y,z,'linear');  	% linear pas obligatoire car par défaut
	% tracé pts réels en diamant rouge , interp. linéaire en noir
plot(x,y,'rd',z,y1,'k'); 
hold on; 		% pour un autre tracé sur le même graphique
	% splines cubiques
y2=interp1(x,y,z,'spline'); 	% ou encore spline(x,y,z)
plot(z,y2,'b');		% spline en bleu
	% cubiques
y3=interp1(x,y,z,'cubic'); 
plot(z,y3,'m');		% cubic en magenta

text(2,1.6,'linéaire en noir');
text(2,1.4,'spline en bleu');
text(2,1.2,'cubique en magenta');
grid;