Discussion :

[Sammy44]

bonjour j'ai programme qui marchait parfaitement avec matlab 7 mais je n'ai plus acces qu'a la version 6 et je dois le faire tourner de facon assez urgente, j'ai essayé de rajouter des variables globales et tout mais rien n'y fait il y a toujours un pb, my_formquad n'est pas utilisée dans fminsearch...
Si vous pouvez m'aider.....

Code :

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

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
29
30
31
32
33
34
35
36
37
38
39
40
41
 
 function param= estimGMMtauxRestreint(R,T)
%R est le vecteur colonne des observation de R
%T est le nombre d'observation sur lequel on va travailler
%on initialise le vecteur de paramètres alpha, beta, gamma, sigma
x0=[0.1252;-0.9590;1.1189]
nbrTour=100 %nombre de tour de boucles
[n,m]=size(R) 
W=eye(3) %initialisée à la matrice identité (de la meme taille coté que longueur de g)
 
    my_formequad=useful_fcn(W);
 
    x=fminsearch(@my_formequad,x0)
 
 
for i=1:1:nbrTour
    W = matricepoids(T,x,R)
 
 
    my_formequad=useful_fcn(W);
 
    x=fminsearch(@my_formequad,x)
end
 
 
 
 
param=x
 
function fcn_hdl = useful_fcn(W_IN)
global W
        W = W_IN;
 
        fcn_hdl = @formequad;
 
        function f = formequad(x)
            f=grestreint(x,R)'*W*grestreint(x,R);
        end
    end
 
end

[Sammy44]

Pour faire marcher si vous voulez m'aider remplacer grestreint(x,R) par x et matrice de poids par W=eye(3)

la derniere partie marche separemment par exemple si je cree une fonction
useful_fct.m

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
 
function fcn_hdl = useful_fcn(W_IN)
 
        W = W_IN;
 
        fcn_hdl = @formequad;
 
        function f = formequad(x)
            f=grestreint(x,R)'*W*grestreint(x,R);
        end
    end
 
end

mais quand je l'appelle :

Code :

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

1
2
3
 
result=useful_fct(eye(3));
result([1;2;3])

retourne une erreur...

[Caro-Line]

Quel est le message d'erreur complet ?

[Jerome Briot]

Et quelle version exacte de MATLAB 6 ? 6.1 ? 6.5 ?

[Sammy44]

matlab 6.5