Discussion :

[bendesarts]

Bonjour,

Je dispose d'une fonction de deux variables p1 et p3 qui s'appelle Pb1.
Je cherche à déterminer les paramètres p1 et p3 qui minimise cette fonction. J'utilise pour cela la fonction fminsearch. Par contre, je n'arrive pas bien à m'en servir.

Code :

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

1
2
3
4
5
6
7
8
9
F=100;
omega=[0.001:0.2:10];
m=10;
k=100;
p1=omega/(sqrt(k/m));
p3 = [0.01:0.1:10];
function y = Pb1(p1,p3)
y = sqrt(1-2*p1.^2+p1.^4+4*p3.^2.*p1.^2)./(p3.*p1.^3);
[x,fval]=fminsearch(Pb1,[0.5, 0.5]);

J'obtiens l'erreur suivante :

Code :

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

1
2
??? Error: File: optimisation.m Line: 12 Column: 1
Function definitions are not permitted at the prompt or in scripts.

Pouvez-vous me dire ce qui ne va pas dans mon code ? Et comment je dois faire pour arriver à ce que la fonction fminsearch fonctionne ?

Merci pour votre aide

[shaiHulud]

Bonjour,

L'erreur n'a rien à voir avec l'optimisation. Tu définis une fonction ligne 7, mais tu es dans un script. 3 solutions possibles
- déplace la fonction dans un fichier à part
- transforme ton script en fonction
- Utilise une fonction en ligne

Code :

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

Pb1= @(p1,p3) sqrt(1-2*p1.^2+p1.^4+4*p3.^2.*p1.^2)./(p3.*p1.^3);

[bendesarts]

OK.

Comment fais-t-on pour faire une fonction en ligne ?
Je suis très novice en matlab. Je l'utilises ici ponctuellement pour faire juste une optimisation.

[FLB]

Tu n'as qu'à remplacer ton code

Code :

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

1
2
function y = Pb1(p1,p3)
y = sqrt(1-2*p1.^2+p1.^4+4*p3.^2.*p1.^2)./(p3.*p1.^3);

par la ligne que VV33D t'as donné.

[bendesarts]

J'ai mis maintenant le code suivant :

Code :

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

1
2
3
4
5
6
7
8
F=100;
omega=[0.001:0.2:10];
m=10;
k=100;
p1=omega/(sqrt(k/m));
p3 = [0.01:0.1:10];
Pb1= @(p1,p3) sqrt(1-2*p1.^2+p1.^4+4*p3.^2.*p1.^2)./(p3.*p1.^3);
[x,fval]=fminsearch(Pb1,[0.5, 0.5]);

J'obtiens alors le message d'erreur suivant :

Code :

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

1
2
3
4
5
6
7
8
9
10
??? Input argument "p3" is undefined.
 
Error in ==> @(p1,p3)sqrt(1-2*p1.^2+p1.^4+4*p3.^2.*p1.^2)./(p3.*p1.^3)
 
 
Error in ==> fminsearch at 205
fv(:,1) = funfcn(x,varargin{:});
 
Error in ==> optimisation at 13
[x,fval]=fminsearch(Pb1,[0.5, 0.5]);

Pour p3 est bien un vecteur ?

[Invité]

Bonjour,

Attention à Pb1= @(p1,p3) : la fonction fminsearch n'effectue son optimisation que par rapport au premier paramètre de la fonction fournie. Les autres éventuels paramètres n'étant que des constantes servant au calcul.
ici bendesarts

cherche à déterminer les paramètres p1 et p3 qui minimise cette fonction.

il faut donc rentrer un vecteur de deux valeurs tel que p = [p1 p3]. Soit p1 = p(1) et p3 = p(2).
On aura alors au final quelque chose comme :

Code :

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

Pb1= @(p) sqrt(1-2*p(1)^2+p(1)^4+4*p(2)^2*p(1)^2)/(p(2)*p(1)^3);

Pour en savoir plus sur les fonctions en ligne, je te suggère la FAQ : Que représente le symbole @ ?

[bendesarts]

OK. Une autre solution est alors peut être de définir ma fonction d'une autre manière que avec le @.

Comment peut-on définir alors une fonction de 2 variables sans le @?

Merci d'avance pour votre aide.