Empty matrix output

Version imprimable

21/08/2009, 16h37
george216

Empty matrix output

Salut les amis,
j'ai défini une fonction intitulée "matrix_exp" qui retourne une matrice x de dimension (T-m,2*p+8) avec T, m et p sont connus. Lorsque j'ai appelé cette fonction et en choisissant p=4, j'ai obtenu une matrice x:"empty 0 by 16".
Cependant, lorsque j'ai copié les lignes définissant cette fonction dans un autre fichier, j'ai obtenu une matrice x bien définie ( càd je l'ai obtenue sans aucun appel de la fonction matrix_exp). J'ai conclu que ce que j'ai écrit pour définir cette fonction est juste mais je ne sais pas où réside le problème.
Esq vous pouvez m'aider svp. Merci d'avance.
George.
21/08/2009, 17h50
phryte

Bonjour.
Fais nous voir ta function.

Voilà Phryte, je vous donne la fonction ainsi que le fichier à partir duquel je l'ai appelée

Merci.

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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
function x=matrix_exp(y,p,d_F,s_F,d_G,s_G,gamma_1_F,gamma_2_F,c_1_F,  c_2_F,gamma_1_G,gamma_2_G, c_1_G,  c_2_G)
 
 
 
 
global T m
 
 
stdS=inv(std(s_F));
stdG=inv(std(s_G));
f=(exp(-gamma_1_F*stdS*(s_F-c_1_F))+exp(gamma_2_F*stdS*(s_F-c_2_F)))./(1+exp(-gamma_1_F*stdS*(s_F-c_1_F))+exp(gamma_2_F*stdS*(s_F-c_2_F)));
g=(exp(-gamma_1_G*stdG*(s_G-c_1_G))+exp(gamma_2_G*stdG*(s_G-c_2_G)))./(1+exp(-gamma_1_G*stdG*(s_G-c_1_G))+exp(gamma_2_G*stdG*(s_G-c_2_G)));
 br=zeros(T-m,4);
for k=1:T-m
     for j=1:4
        if k==j
            br(k,j)=1;
        elseif rem(abs(k-j),4)==0
            br(k,j)=1;
        else
            br(k,j)=0;
        end
     end
end
brf=zeros(length(br),4);
for i=1:length(br)
for k=1:4
    brf(i,k)=br(i,k)*f(i,1);
end
end
x=zeros(T-m,2*p+8);
for i=1:T-m
  for j=1:4
 x(i,j)=br(i,j);
  end
end
 
for i=1:T-m
  for j=5:8
x(i,j)=brf(i,j);
  end
end
 
for i=1:T-m
  for j=1:p
    x(i,j+8)=y(m+1-j,1);
    x(i,j+p+8)=x(i,j+1)*g(i,1);
  end
  end
return

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
yy=[10 11 10 11 12 13 14 15 16 12 11 12 14 15 16 17 13 14 15 16 17 18 17 20 21 22 20 21 23 24 25 26 20 16 17 18 19 30 31 32 30 33 34 35 31 29 28 27 31 33 36 30 ]';
zz=zeros(length(yy),1);
for i=5:length(yy)
    zz(i)=yy(i)-yy(i-4);
end
y=yy(5:length(yy));
s_F=zz(5:length(zz));
s_G=y;
d_F=1;d_G=2;p=4;
gamma_1_F=0.1;gamma_2_F=0.1;c_1_F=0;  c_2_F=0;gamma_1_G=0.1;gamma_2_G=0.1; c_1_G=0;  c_2_G=0;
T=length(y);
m=max([p d_F d_G]);
x=matrix_exp(y,p,d_F,s_F,d_G,s_G,gamma_1_F,gamma_2_F,c_1_F,  c_2_F,gamma_1_G,gamma_2_G, c_1_G,  c_2_G)

C'est l'utilisation de GLOBAL qui est mauvaise.
Quand tu arrives dans ta fonction T et m sont vides.

Il faut éviter d'utiliser les variables globales, on se retrouve dans ce genre de situation après.
Pourquoi ne passes-tu pas T et m en paramètres de ta fonction : cela résoudra ton problème dans un 1er temps.

Mais pas dans un 2nd puisque ça donne une autre erreur.
A mon avis ton code doit pouvoir être amélioré.
Si j'ai du temps je regarde.

[EDIT] quelques améliorations :

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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
function truc()
yy=[10 11 10 11 12 13 14 15 16 12 11 12 14 15 16 17 13 14 15 16 17 ...
    18 17 20 21 22 20 21 23 24 25 26 20 16 17 18 19 30 31 32 30 33 ...
    34 35 31 29 28 27 31 33 36 30 ]';
y=yy(5:end);
s_F=y-yy(1:end-4);
s_G=y;
d_F=1;
d_G=2;
p=4;
gamma_1_F=0.1;
gamma_2_F=0.1;
c_1_F=0;
c_2_F=0;
gamma_1_G=0.1;
gamma_2_G=0.1;
c_1_G=0;
c_2_G=0;
x=matrix_exp(y,p,d_F,s_F,d_G,s_G,gamma_1_F,gamma_2_F,c_1_F,  c_2_F,gamma_1_G,gamma_2_G, c_1_G,  c_2_G);
 
 
function x=matrix_exp(y,p,d_F,s_F,d_G,s_G,gamma_1_F,gamma_2_F,c_1_F,  c_2_F,gamma_1_G,gamma_2_G, c_1_G,  c_2_G)
T=length(y);
m=max([p d_F d_G]);
%global T m
 
stdS=inv(std(s_F));
stdG=inv(std(s_G));
f=(exp(-gamma_1_F*stdS*(s_F-c_1_F))+exp(gamma_2_F*stdS*(s_F-c_2_F)))./...
    (1+exp(-gamma_1_F*stdS*(s_F-c_1_F))+exp(gamma_2_F*stdS*(s_F-c_2_F)));
g=(exp(-gamma_1_G*stdG*(s_G-c_1_G))+exp(gamma_2_G*stdG*(s_G-c_2_G)))./...
    (1+exp(-gamma_1_G*stdG*(s_G-c_1_G))+exp(gamma_2_G*stdG*(s_G-c_2_G)));
% br=zeros(T-m,4);
% for k=1:T-m
%      for j=1:4
%         if k==j
%             br(k,j)=1;
%         elseif rem(abs(k-j),4)==0
%             br(k,j)=1;
%         else
%             br(k,j)=0;
%         end
%      end
% end
%**
%"diagonales" à 1
%Hypothèse FORTE : T-m est un multiple de 4
minS = min(size(br));
maxS = max(size(br));
br = repmat(eye(minS),maxS/minS,1);
%**
brf=zeros(length(br),4);
for i=1:length(br)
    for k=1:4
        brf(i,k)=br(i,k)*f(i,1);
    end
end
 
x=zeros(T-m,2*p+8);
x(1:T-m,1:4) = br(1:T-m,:);
% for i=1:T-m
%     for j=1:4
%         x(i,j)=br(i,j);
%     end
% end
 
x(1:T-m,5:8) = brf(1:T-m,:);
% for i=1:T-m
%     for j=5:8
%         x(i,j)=brf(i,j-4);
%     end
% end
 
for i=1:T-m
    for j=1:p
        x(i,j+8)=y(m+1-j,1);
        x(i,j+p+8)=x(i,j+1)*g(i,1);
    end
end
return

22/08/2009, 12h56
george216

Merci bien Caroline,
j'ai utilisé l'expression " global T m", parce que j'ai d'autres fonctions qui utilisent en commun ces paramètres. Je vais essayer de supprimer cette expression.