Discussion :

[karika]

Bonjour,voila l'article que je suis entrain de simuler,http://perso.ensem.inpl-nancy.fr/Jos...lita_ragot.pdf
et voila le programme que j'ai fais:

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
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
N=100;k=1:N; ndim=7;
v1(k)=rand(1,100);
v2(k)=rand(1,100);
v3(k)=rand(1,100);
%LES VARIABLES
x1=0.4*v1(k)+sin(k/N);
x2=v2(k)-2*cos(k/4);
x3=0.2*v3(k)-1;
x4=x1+x2;
x5=x2+x3;
x6=2*x1+x3;
x7=x4+x5;
x=[x1' x2' x3' x4' x5' x6' x7'];
%AJJOUT DES DEFAUTS
x1(:,8:15)=x1(:,8:15)+3;
x3(:,34:41)=x3(:,34:41)+3;
x5(:,58:65)=x5(:,58:65)+3;
%PLOTER LES VARIABLES
subplot(711);plot(x1);axis([0 100 -4 4]);
subplot(712);plot(x2);axis([0 100 -4 4]);
subplot(713);plot(x3);axis([0 100 -4 4]);
subplot(714);plot(x4);axis([0 100 -4 4]);
subplot(715);plot(x5);axis([0 100 -4 4]);
subplot(716);plot(x6);axis([0 100 -4 4]);
subplot(717);plot(x7);axis([0 100 -4 4]);
%MATRICE DE CORRELATION DES DONNEES NON BRUITES
a=x'*x;
%VALEURS ET VECTEURS PROPRES DE LA ATRICE DE CORELATION
[v,lamda]=eig(a);
%DALTA DE X,LE BRUIT
dalta=rand(100,7);
%AJJOUT DU BRUIT
 X=x+dalta;
%NOUVELLE MATRICE DE DONNEES
X1=X(1,:);
X2=X(2,:);
X3=X(3,:);
X4=X(4,:);
X5=X(5,:);
X6=X(6,:);
X7=X(7,:);
X=[X1' X2' X3' X4' X5' X6' X7'];
%SENSIBILITE DE LA MATRICE DE CORRELATION AU BRUIT(DALTA1=DALTA DE A
dalta1=dalta'*x+x'*dalta+dalta'*dalta;
%NOUVELLE MATRICE DE CORRELATION
A=a+dalta1;
i=1:7;
%DALTA3=DALTA DE lamda
dalta3(i,:)=v(i,:)'*dalta1*v(i,:);
%MINIMUM ET LE MAXIMUM DES VALEURS PROPRES
LAMDA=lamda+dalta3;
mLAMDA=max((LAMDA));
nLAMDA=min((LAMDA));
%DALTA4=DALTA DE V
w=inv([a-lamda*ones(7,7) -v;v' zeros(7,7)])*[-dalta1*v;zeros(7,7)];
dalta4=w(1:7,:);
%LE MINIMUM ET LE MAXIMUM DES VECTEURS PROPRES
V=v+dalta4;
for i=1:7;
     j=1:7;
nV(i,j)=min(V(i,j));
 
end 
 for i=1:7;
     j=1:7;
mV(i,j)=max(V(i,j));
 end

mais je sais pas pour quoi les colonnes du minimum et du maximum des vecteurs propres(de la ligne59 au 65) sont identique

[le fab]

salut

a priori v(i,j) est un scalaire donc

Code :

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

v(i,j) == min(v(i,j)) == max(v(i,j))

Fabien

[Invité]

Bonsoir,

Code :

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

1
2
3
4
5
6
7
8
X1=X(1,:);
X2=X(2,:);
X3=X(3,:);
X4=X(4,:);
X5=X(5,:);
X6=X(6,:);
X7=X(7,:);
X=[X1' X2' X3' X4' X5' X6' X7'];

ou plus simplement : X = X(1:7,:).';

[karika]

voici le résultat de nV:

Code :

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

1
2
3
4
5
6
7
8
9
10
 
nV =
 
   -0.5648   -0.5648   -0.5648   -0.5648   -0.5648   -0.5648   -0.5648
   -0.4168   -0.4168   -0.4168   -0.4168   -0.4168   -0.4168   -0.4168
   -0.6090   -0.6090   -0.6090   -0.6090   -0.6090   -0.6090   -0.6090
   -0.4808   -0.4808   -0.4808   -0.4808   -0.4808   -0.4808   -0.4808
   -0.4334   -0.4334   -0.4334   -0.4334   -0.4334   -0.4334   -0.4334
   -0.3274   -0.3274   -0.3274   -0.3274   -0.3274   -0.3274   -0.3274
   -0.0842   -0.0842   -0.0842   -0.0842   -0.0842   -0.0842   -0.0842

on remarque que la même colonne est entrain de ce répéter,et je sais pas pour quoi ???

[FLB]

Bonjour,
peux tu nous montrer la nouvelle version de ton code, qui prend en compte les remarques de le fab et de Winjerome?

[Invité]

As-tu pris en compte la remarque de le fab par rapport à ces lignes

Code :

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

1
2
3
4
for i=1:7;
    j=1:7;
    nV(i,j)=min(V(i,j));
end

?
@le fab : j=1:7; donc V(i,j) est bien un vecteur et non un simple scalaire.
@karika : Tu calcules avec ceci le maximum de chaque ligne de V, unique valeur sur chaque ligne donc que tu recopies sur chaque colonne de nV nV(i,j)=min(V(i,j)); avec j=1:7;

[karika]

merci beaucoup,voila ce que je souhaite faire exactement,et dites moi si c'est possible ou non.
Je veux que chaque composante du vecteur V auras un minimum et un maximum,c-à-d comme le vecteur V et de dimension 7*7 alors il ya 49 composantes,donc je souhaite avoir 49 valeur minimum et 49 valeur maximum,par exemple:
la composante V(1,1) aura une valeur minimum et une valeur maximum,et pour la composante V(1,2) aussi;et ainsi de suite pour toutes les composantes du vecteur V.
merci.

[le fab]

Salut
c'est possible si tu rajoutes une dimension à V
V(i,j,1) contient le minimum et V(i,j,2) le maximum par exemple
mais V(i,j) ne peux pas contenir deux valeurs

sinon tu peux aussi faire un tableau de structure comme ceci:
V(i,j).min et V(i,j).max

Fabien

[FLB]

Bonjour

il ya 49 composantes,donc je souhaite avoir 49 valeur minimum et 49 valeur maximum

J'ai du mal à te suivre, tu cherches le minimum et le maximum de quoi? Pour une unique composante (= 1 scalaire), tu auras maximum = minimum = composante.

[karika]

Je suis en train de simuler cet article http://perso.ensem.inpl-nancy.fr/Jos...lita_ragot.pdf, et comme je suis débutant en MATLAB j'ai besoin de quelqu'un qui peut m'orienter ; si vous le pouviez je vous serais très reconnaissant

[karika]

bonjour,
v(i,j) min = infk v(i,j) + δ(vi,j) (k)
v(i,j) max = supk v(i,j) + δ(vi,j) (k)
comment peut étre le code matlab de ces 2 équations???merci