Discussion :

[mounabs]

Bonjour,

j'ai un process de diffusion dont la moyenne est nonlineaire de la forme:
mu(dt) = (alpha1+alpha2*x+alpha3*(x2)+alpha4/x) dt

je veux simuler cette moyenne j'ai alors fait le programme suivant
mais j'obtient un graphique vide

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
clf
load TRB.txt				% Load Data
 
X = TRB(1:T,1)/100;
T = rows(TRB);
 N = 2087; Delta = 1/T;           
% simulation du processus de CIR-SR
alpha1 = -4.0224; alpha2 = 0.08893; alpha3 = -0.0608; alpha4 = 7.1945;  
 
for X = 1:T
 
    mu = (alpha1 +alpha2*X + alpha3*(X^2) + alpha4/X)* Delta;
plot(X,mu);
 
end

est ce que vous pouvez me dire où est le problème.

Merci

[Jean Dumoncel]

Bonjour,

en fait ton graphique n'est pas totalement vide : il contient les points que tu traces a chaque itération.
Pour tracer une courbe, il faut stocker les points et afficher la courbe après la boucle :

Code :

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

1
2
3
4
for X = 1:T
    mu(X) = (alpha1 +alpha2*X + alpha3*(X^2) + alpha4/X)* Delta;
end
plot(1:T,mu);

Par contre c'est normal que tu définisses X puis que tu l'utilises aussi comme variable d'itération? J'ai l'impression que tu souhaites faire :

Code :

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

1
2
3
4
for k = 1:T
    mu(k) = (alpha1 +alpha2*X(k) + alpha3*(X(k)^2) + alpha4/X(k))* Delta;
end
plot(X,mu);

non?

[mounabs]

Merci pour la rectification mais bizarre le graphique qui apparait reste toujours vide. Comment faire apparaitre la courbe?

[Jean Dumoncel]

As-tu regardé ce que tu obtiens dans le vecteur mu?

[mounabs]

non comment faire, rien ne m'apparait

[Jean Dumoncel]

Déjà essaie de faire :
clear all
close all
puis relances ton programme.

Si il n'y a toujours rien qui s'affiche :
pour afficher mu, si il ne contient pas beaucoup de valeurs, tu peux faire dans la fenêtre de commande :

Code :

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

>> mu

ou bien double-clic sur son nom dans la fenêtre nommée Workspace.

[mounabs]

j'ai changer le programme comme suit:

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
clear all
clf
load TB3S.txt				% Load Data
T = rows(TB3S);
		
x = TB3S(1:T,1)/100;
x0 = 0.081;
x = zeros(T,1);
 N = 2087; Delta = 1/N;         
%simate the ait sahalia drift
alpha1 = 0.006; alpha2 = 0.0721; alpha3 = 0.141; alpha4 = 0.0003;  
x0 = 0.081;
x = zeros(T,1);
x(1) = x0;
for k = 1:T
    mo(k) = (alpha1 +alpha2*x(k) + alpha3*(x(k)^2) + alpha4/x(k))* Delta;
end
plot(1:T,mo);
close all
    

mais le problème persiste le graphique apparait vide et le vecteur mo apart la première valeur , toutes les valeurs sont inf. je ne sais pas où est le problème

[Jean Dumoncel]

As-tu regardé les valeurs contenues dans x? Elles te paraissent correctes? Une division par 0 donne Inf par exemple...

[membreComplexe12]

je ne sais pas si c'est ça le probleme mais remplace ceci:

Code :

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

plot(1:T,mo);

par:

Code :

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

1
2
3
4
 
figure (1)
hold on
plot(1:T,mo);

[Jean Dumoncel]

je ne sais pas si c'est ça le probleme mais remplace ceci:

Code :

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

plot(1:T,mo);

par:

Code :

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

1
2
3
4
 
figure (1)
hold on
plot(1:T,mo);

Je ne vois pas bien le rapport avec le problème de mounabs?

[membreComplexe12]

je me suis dit que ses points n'apparaiseent pas car il a peut etre ecrasé un précédent graphique ou un truc dans ce genre. Je voulias juste verifier que ça soit pas un probleme dans ce genre...

[Jean Dumoncel]

Non, les point n'apparaissent pas car :

e vecteur mo apart la première valeur , toutes les valeurs sont inf.

[mounabs]

bonjour:
le code est devenu comme suit:

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
clf
 
 load TB3S.txt % read your data here: create a vector x where x(1) = first observations, x(n) = last observations
T = rows(TB3S);			% nb. obs.
x = TB3S(1:T,1)/100;
 
 Delta = 1/T;
 % simulation du processus de CIR-SR
 alpha1 = 0.0143; alpha2 = -0.0006; alpha3 = -0.0136; alpha4 = 0.0003;  
 
for i = 1:T
 
drift(i) = (alpha1 +alpha2*x(i) + alpha3*(x(i)^2) + alpha4/x(i))* Delta;
 plot(x(i),drift(i));
 hold on;
 
end
 
hold off;

la figure apparait mais pas assez lisible pourquoi?

[Jean Dumoncel]

As-tu compris ce que j'ai écris dans ce précédent message :

Pour tracer une courbe, il faut stocker les points et afficher la courbe après la boucle :

Code :

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

1
2
3
4
for X = 1:T
    mu(X) = (alpha1 +alpha2*X + alpha3*(X^2) + alpha4/X)* Delta;
end
plot(1:T,mu);

?

[mounabs]

Code :

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

1
2
3
4
clf
 
 load TB3S.txt 
 x = TB3S(1:T,1)/100;

avec x(1) = première observation, x(n) = dernière observation

Code :

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

1
2
3
T = rows(TB3S);			% nb. obs.
 
 Delta = 1/T;

valeurs initiales

Code :

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

alpha1 = 0.0143; alpha2 = -0.0006; alpha3 = -0.0136; alpha4 = 0.0003;

Code :

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

1
2
3
4
for X = 1:T
    mu(X) = (alpha1 +alpha2*X + alpha3*(X^2) + alpha4/X)* Delta;
end
plot(1:T,mu);

ici le X doit être les valeur dans le vecteur TB3S.txt de la première valeur jusqu'à la dernière.
comment faire

[Jean Dumoncel]

Par rapport à ton code du post #13, il suffisait juste de sortir le plot :

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
clf
 
load TB3S.txt % read your data here: create a vector x where x(1) = first observations, x(n) = last observations
T = rows(TB3S);			% nb. obs.
x = TB3S(1:T,1)/100;
 
 Delta = 1/T;
 % simulation du processus de CIR-SR
 alpha1 = 0.0143; alpha2 = -0.0006; alpha3 = -0.0136; alpha4 = 0.0003;  
 
for i = 1:T
 
drift(i) = (alpha1 +alpha2*x(i) + alpha3*(x(i)^2) + alpha4/x(i))* Delta;
 
end
plot(x,drift);