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 |
Partager