Discussion :

[cpalperou]

Bonjour,
je rencontre un problème pour fitter des données à un modèle.
Il s'agit de valeurs N dépendant de 2 variables, la latitude (l) et le temps (t).
Le modèle est le suivant:

Nmod(l,t)=N0+N1cos[(2*pi*t/12)+t1]+N2*(l+pi/2)(l-pi/2)cos[(2*pi*t/6)+t2]

Les paramètres à ajuster étant N0, N1, N2, t1 et t2.
Je n'ai aucune idée de comment faire ni aucun point de départ.
Pouvez vous m'aider?
Merci

[le fab]

salut

commence par :

Code :

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

1
2
3
4
 
x0 = [0,0,0,0,0]; % améliorer si tu as une idée de point de départ meilleur
fun = @(x,xdata) x(1)+x(2)*cos[(2*pi*xdata(2)/12)+x(4)]+x(3)*(xdata(1)+pi/2)(xdata(1)-pi/2)cos[(2*pi*xdata(2)/6)+x(5)]
x = lsqcurvefit(fun,x0,xdata,ydata);

ou x represente les coefficients [N0, N1, N2, t1, t2], x0 les coefficients de départ pour l'optimisation
xdata est une matrice m*2 [l t] de tes entrées
ydata une matrice m*1 [N] du résultat à fitter

bon pour ça faut avoir l'optimisation toolbox
si tu as la curvefitting toolbox, ça doit être encore plus simple

[cpalperou]

Bonjour et merci leFab.
Je ne comprends pas ce que doivent être xdata et ydata.
Ce que j'ai: un vecteur l contenant les latitudes de dimensions 1*36, un vecteur time de dimension 1*124 et une matrice N de dimensions 36*126.

D'après la façon dont tu écris fun, xdata(2) contient le vecteur temps (dimension 1*36) et xdata(1) les latitudes (dim 1*124). Comment écrire alors xdata, si ces 2 composantes n'ont pas les même dimensions?
Ce que tu appelles ydata correspond, je pense à ma matrice N. Est-ce correct?

Merci de m'aider

[le fab]

c'est ca ...

solution un peu bourrin pour construire xdata et ydata :

Code :

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

1
2
3
4
 
[L,T] = meshgrid(vLatitudes,vTemps)
xdata = [L(:) T(:)];
ydata = N(:);

avec peut être des transposés à mettre, peut être inverser L T dans meshgrid....
le but est de créer 2 vecteur colonnes de 36*124 éléments pour la latitude et le temps et un vecteur N de même taille avec correspondance

y a peut etre mieux et plus propre, mais c'est comme ca que j'ai l'habitude de faire

[cpalperou]

Bonjour,
le vecteur L contient 36 éléments.
Le vecteur T contient 124 éléments.

Logiquement, il est impossible de les concaténer:
Par exemplz:

Code :

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

1
2
3
4
5
6
7
 
L=ones(36,1);
T=ones(124,1);
xdata = [L(:) T(:)];
??? Error using ==> horzcat
All matrices on a row in the bracketed expression must have the 
 same number of rows.

Je ne comprends pas ta méthode!

[le fab]

la première ligne est indispensable (meshgrid)
L et T sont alors des matrices de 36*124 ou 124*36 suivant l'ordre d'appel dans la fonction meshgrid
ensuite L( : ) et T( : ) sont des vecteurs colonnes de 36*124 éléments, tout comme N( : ) et tu as alors une correspondance élément par élément (L(12) & T(12) -> N(12))
et donc tu obtient des xdata et ydata cohérents que tu peux utiliser dans lsqcurvefit