Discussion :

[Marcusss]

bonjour,

J'ai une fonction à plusieurs variables f(x1,x2,x3....) et je veux chercher le minimum par rapport à x1 seulement,avec comme contrainte 0<x1<1

J'ai essayé des fonctions tel que fminsearch,fminbnd,min et tout ça mais ça ne marche pas.
Quelle fonctions devrais-je utiliser s'il vous plait?

Merci

[rostomus]

Pourais tu nous donner un plus d'explication sur ton probleme??

un bout de code par exemple?

[Marcusss]

Bonjour,
voila une partie du code

Code Matlab :

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

1
2
3
4
5
6
7
8
9
h(1)=(1-lambda)*r(1)^2; 
for i=2:N-1
    h(i) = lambda*h(i-1)+(1-lambda)*r(i)^2; % The expression of volatility
end
 
RMSE=0;
for i=1:N-1
    RMSE=RMSE+sqrt((1/N)*(r(i+1)^2)-(h(i)^2))
end

Puisque h dépend de lambda donc RMSE dépend également de lambda,,,
et je cherche le lambda qui minimise le RMSE
PS: r est un vecteur donné

merci

[rostomus]

Salut,

Je ne suis pas sûr si ça va marcher ou non, essaie de faire ça:

Crée une fonction :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
function RMSE=essaie(lambda)
r=........
h(1)=(1-lambda)*r(1)^2; 
for i=2:N-1
    h(i) = lambda*h(i-1)+(1-lambda)*r(i)^2; % The expression of volatility
end
 
RMSE=0;
for i=1:N-1
    RMSE=RMSE+sqrt((1/N)*(r(i+1)^2)-(h(i)^2))
end

et tapes cette ligne:

Code :

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

1
2
 
lambda0=fminsearch(@essaie,.5)

[Marcusss]

En tapant la ligne de commande,j'ai eu ce message d'erreur:

Code :

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

1
2
3
4
5
6
7
??? Attempted to access r(1); index out of bounds because numel(r)=0.

Error in ==> dec at 11
h(1)=(1-lambda)*r(1)^2; 

Error in ==> fminsearch at 175
fv(:,1) = funfcn(x,varargin{:});

C'est du à quoi?

[rostomus]

C'est du à quoi?

regarde:

Code :

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

1
2
3
4
5
6
7
??? Attempted to access r(1); index out of bounds because numel(r)=0.

Error in ==> dec at 11
h(1)=(1-lambda)*r(1)^2; 

Error in ==> fminsearch at 175
fv(:,1) = funfcn(x,varargin{:});

Avec quoi t'as remplacé "r=........" ?