Discussion :

[karika]

Bonjour tout le monde,

J'ai 3 variable v1(k),v2(k) et v3(k) qui suivent une loi normale de moyenne nulle et de variance &^2,et k=1:100;
Comment je plot ces variables avec MATLAB ?

Merci pour votre aide.

[le fab]

Code :

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

1
2
3
4
5
figure;
hold on
plot(v1);
plot(v2);
plot(v3);

[karika]

merci,et comment j'identifie ces 3 variable avant plot

[FLB]

Bonjour,
pour générer tes vecteurs v, utilise la fonction randn.

Code :

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

doc randn;

[karika]

dans mon exemple les trois variables suivent une lois normale de moyenne=0 et de variance inconnu,comment je fait dans ce cas...et merci beaucoup pour votre aide.

[FLB]

de variance inconnu

Tu souhaites générer ces trois variables, ou tu les as déjà?

[karika]

Code :

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

1
2
3
4
5
6
7
x1(k) = 0.4*v1(k)+sin(k/N),        v1(k) ∼ N(0, σ2)
x2(k) = v2(k)−2*cos(k/4),           v2(k) ∼ N(0, σ2)
x3(k) = 0.2*v3(k) − 1,                 v3(k) ∼ N(0, σ2)
x4(k) = x1(k) + x2(k)
x5(k) = x2(k) + x3(k)
x6(k) = 2*x1(k) + x3(k)
x7(k) = x4(k) + x5(k)

voila mon exemple, je veux simuler les 7 variables, mais je ne sais pas comment identifier dans MATLAB v1(k),v2(k) et v3(k) qui suivent une lois normale de moyenne =0 et de variance inconnue.....

[karika]

je veut générer ces 3 variable v1(k),v2(k) et v3(k) de moyenne=0 et de variance inconnu

[shaiHulud]

Bonjour,

Je ne comprends pas trop ce que tu cherches à faire
- Estimer la variance de x1: var(x1)
- Tester la Gaussianité : kstest(), jbtest()
- afficher l'histogramme : hist(x1)

[karika]

en principe c'est bon ça marche merci beaucoup....maintenant j'ai un autre problème
voila ce que j'ai fait comme programme:

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
N=100;k=1:N;
v1(k)=0+var(v1(k)).*randn(1,k);
v2(k)=0+var(v2(k)).*randn(1,k);
v3(k)=0+var(v3(k)).*randn(1,k);
 
x1(k)=0.4*v1(k)+sin(k/N);
x2(k)=v2(k)-2*cos(k/4);
x3(k)=0.2*v3(k)-1;
x4(k)=x1(k)+x2(k);
x5(k)=x2(k)+x3(k);
x6(k)=2*x1(k)+x3(k);
x7(k)=x4(k)+x5(k);
 
x=[x1(k)' x2(k)' x3(k)' x4(k)' x5(k)' x6(k)' x7(k)'];
 
x(8:15,1)=x(8:15,1)+3;
 
subplot(711);plot(x1(k));
subplot(712);plot(x2(k));
subplot(713);plot(x3(k));
subplot(714);plot(x4(k));
subplot(715);plot(x5(k));
subplot(716);plot(x6(k));
subplot(717);plot(x7(k));

dan cette ligne

Code :

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

x(8:15,1)=x(8:15,1)+3;

j'ai voulu ajouter une constante qui est 3 a l'instant de 8:15 pour x1(k) mais ça na pas fonctionné,pouvez vous m'indiquer ou est mon erreur ??

[le fab]

bah si ca marche très bien
ce qui ne marche pas c'est

Code :

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

v1(k)=0+var(v1(k)).*randn(1,kN);

v1(k) n'est pas défini et tu le passe en argument à var et tu passes un vecteur à randn

[karika]

dacord merci...et ce qui concerne l'ajout de la constante 3 à la variable x1(k) à l'instant de 8:15 comment ça ce fait??

[le fab]

et ce qui concerne l'ajout de la constante 3 à la variable x1(k) à l'instant de 8:15 comment ça ce fait??

ce que tu as écrit marche ... mais ne fait peut etre pas ce que tu attend
après je ne comprend pas ce que tu veux dire par l'instant de 8:15

[karika]

je veux ajouter une constante "3" au vecteurs: de la 1 ligne et de colonnes de 8 jusqu’à 15.comment je fait ?

1 ligne et de 8:15 colonnes de la matrice x

[Invité]

Bonjour,

Comme te l'a dit le fab la ligne :

Code :

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

x(8:15,1)=x(8:15,1)+3;

marche très bien

Un exemple :

Code :

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

1
2
x = (1:20).'
x(8:15,1)=x(8:15,1)+3

[le fab]

je veux ajouter une constante "3" au vecteurs: de la 1 ligne et de colonnes de 8 jusqu’à 15.comment je fait ?

1 ligne et de 8:15 colonnes de la matrice x

dans ce cas

Code :

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

x(1,8:15)=x(1,8:15)+3;