Discussion :

[membreComplexe12]

Salut tous,

j'aimerai faire différents tests d'interpolations sous matlab. J'aimerai tester :

1°) l'interpolation polynomiale classique
2°) l'interpolation par moindre carrés
3°) les splines
4°) l'interpolation par polynomes de Lagrange
5°) l'interpolation par polynome d'euler

sauriez vous comment faire tous ceci sous matlab ?

par exemple si je veux interpoler ces données:

x=1:1:10
y=5,4,10,4,5,8,9,3,50,4

pour les splines j'ai trouvé dans l'aide mais pour le reste je n'ai pas trouvé comment faire

merci d'avance pour toutes les pistes que vous pourrez me donner

A+

[duf42]

Bonjour,

Je pense qu'une partie de ton problème doit pouvoir être résolu avec les fonctions POLYFIT/POLYVAL et INTERP1.

Qu'utilises-tu pour les splines?

Duf

[membreComplexe12]

salut DUF et merci d'avoir pris le temps de répondre

Je pense qu'une partie de ton problème doit pouvoir être résolu avec les fonctions POLYFIT/POLYVAL et INTERP1.

=> j'ai regardé et effectivement POLYFIT permet bien de faire l'identification d'un polynome au sens des moindre carré.

=> par contre pour le reste je ne pense pas que c'est ce que je cherche. En fait j'aimerai une fonction qui fonctionne un peu comme polyfit mais avec des polynomes de Lagrange ou Euler...

Qu'utilises-tu pour les splines?

=> la fonction "spline" dans l'aide de matlab tout à la fin il y a un exemple qui explique comment utiliser les splines, c'est assez facile.

[Gakusei]

Bonjour,

Si tu regarde doc interp2, tu vois différents type d'interpolation 'spline', 'cubic',... Cela te résous certain points (les splines que tu as trouvé).

1) interpolation polynomiale : doc polyval
2) interpolation au sens des moindres carrés: doc polyfit
3) Les splines c'est fait
4) Lagrange:
On pose la température T et masse volumique R(T)
Et voici un script issue de la formule:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
function Ri=lagrange(T,R,Ti)
Ri=zeros(size(Ti)); % Initialisation des Ti 
n=length(R); % Nombre de points
for i=1:n 
    y=ones(size(Ti)); 
    for j=1:n 
        if i~=j 
            y=y.*(Ti-T(j))/(T(i)-T(j)); 
        end 
        Ri=Ri+y*R(i) 
    end 
end 
return

Désolé le temps que je tape le script il y a eu des postes donc du déjà dit. Sauf pour Lagrange

[membreComplexe12]

merci pour votre aide