Minimisation d'une fonction

Version imprimable

28/04/2011, 16h38
foufi86

Minimisation d'une fonction

Bonjour
j'espére bien que vous pouvez m'aider car vraiment je suis débutante en matlab et j'en ai besoin d'aide mon probleme en faite consiste a une minimisation d'erruer entre une courbe experimentale et autre analytique envue de determiner trois parametres avec la fonction fminsearch.Mes questions sont les suivantes :
1- c'est ou que je doit stoker mes .m fonctions
2-comment faire cette minimisation et commententrer les valeus expérimentales
merci.
28/04/2011, 17h40
Jean Dumoncel

Citation:

Envoyé par foufi86

1- c'est ou que je doit stoker mes .m fonctions

Peu importe ou tu les stockes du moments qu'ils sont "visibles" par MATLAB. Le plus simple, c'est de définir le répertoire courant de MATLAB ou tu veux et de placer tes fichiers dedans.

Citation:

Envoyé par foufi86

2-comment faire cette minimisation et commententrer les valeus expérimentales
merci.

Un exemple d'utilisation de fminsearch :
Nonlinear Curve Fitting
29/04/2011, 11h08
arbas

Bonjour:
fminsearch est à base de la méthode de moindres carrés non linéaire.
pour l'utiliser:
help fminsearch
Inconvenient de cette méthode:
Vous pouvez pas encadrez les paramètres cherchés dans un intervalle donnée , vous pouvez obtenir une meilleur identification mais les paramètres ne sont pas dans un intervalle qui satisfait des conditions sur les paramètres.
29/04/2011, 12h37
tonyjub

simplex

Bonjour,

c'est la méthode du simplex qui est utilisée par fminsearch, c'est une méthode très robuste par rapport au problème par contre elle n'est pas contrainte par défault, un user US a programmer cette fonction avec contrainte, il faut rechercher "boundary value fminsearch" sur le site de mathworks.

Le défaut de cette méthode, c'est que comparé à fmincon on a pas accès à des données qui peuvent s'avérer importante dans un problème comme le gradiant, le hessien, ect...

réponse

bonjour
je vous remercie tous pour vos réponses.
voila c'est le code que j'ai utilisé pour résoudre mon problème :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
function main
global ti Fi
fid =fopen('données simulat.m','r');
M = fscanf(fid,'%e',[1,1]);
xy=fscanf(fid,'%e%e',[2,M]);
xy1=xy';
ti=xy1(:,1);
Fi=xy1(:,2);
x0 = [0.08 0.05 0.05]; %initialisation aléatoire des paramètres dont je ne dispose d'aucune idée concernant leur estimation
xsol = fminsearch ( @F, x0)
plot ( f(xsol, ti))
 
 
function y = f(x,t)
y=(5480351341935199.*x(1).^2.*heaviside(t - 20).*(37./(975.*x(1) + 925.*x(2)) - (5476.*x(2).*x(3))./ ...
    (exp(((925.*x(2).^2 + 975.*x(1).*x(2)).*(t - 20))./(2886.*x(1).*x(3) + 2886.*x(2).*x(3))).* ...
    (2886.*x(1).*x(3) + 2886.*x(2).*x(3)).*(975.*x(1) + 925.*x(2)))))./3010462836850688 + ...
    (5.*1590.^(1./2).*x(1).*x(2).*heaviside(t - 20).*(37./(975.*x(1) + 925.*x(2)) - ...
    (5476.*x(2).*x(3))./(exp(((925.*x(2).^2 + 975.*x(1).*x(2)).*(t - 20))./(2886.*x(1).*x(3) + ...
    2886.*x(2).*x(3))).*(2886.*x(1).*x(3) + 2886.*x(2).*x(3)).*(975.*x(1) + 925.*x(2)))))./111 + ...
    (1622183997212819.*x(2).^2.*x(3).*heaviside(t - 20))./(8246337208320.*exp(((925.*x(2).^2 + ...
    975.*x(1).*x(2)).*(t - 20))./(2886.*x(1).*x(3) + 2886.*x(2).*x(3))).*(2886.*x(1).*x(3) + ...
    2886.*x(2).*x(3))) + (2.*1590.^(1./2).*x(2).^2.*x(3).^2.*heaviside(t - 20).* ...
    (25./(296.*x(2).*x(3).^2.*exp((25.*x(2).*(t - 20))./(74.*x(3)))) - (975.*x(1) + 925.*x(2))./ ...
    (4.*x(2).*x(3).*exp(((925.*x(2).^2 + 975.*x(1).*x(2)).*(t - 20))./(2886.*x(1).*x(3) + 2886.*x(2).*x(3))).* ...
    (2886.*x(1).*x(3) + 2886.*x(2).*x(3)))))./375 + (3332053615896601.*x(1).*x(2).*x(3).*heaviside(t - 20))./ ...
    (16492674416640.*exp(((925.*x(2).^2 + 975.*x(1).*x(2)).*(t - 20))./(2886.*x(1).*x(3) + 2886.*x(2).*x(3))).* ...
    (2886.*x(1).*x(3) + 2886.*x(2).*x(3)));

Code:

1
2
3
function out = F(x)
global ti Fi
out = sum ( (Fi - f(x,ti)).^2 );

le probleme c'est qu'il me donne un paramètre qui est négatif :cry:
Y-a-t'il quelqu'un qui peut me dire pourquoi
merci beaucoup