Discussion :

[florelle]

Can you please help me to correct this code. I have an error in the boucle for/end and the function plot. Because xtemp and xem is vector matlab can't plot this. Thank you
--------
Pouvez-vous s'il vous plait m'aider à corriger ce code. J'ai une erreur au niveau de la boucle for/end et la fonction plot. xtemp et xem étant des vecteurs, MATLAB ne peut pas plotter cela. Merci.

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
%Initial param
yield;
randn('state',24)
theta=0.37;
sigma_r=0.01;
sigma_s=0.6;
sigma_l=0.4;
sigma_e=0.21;
lambda_s=0.06;
lambda_l=0.05;
lambda_e=0.12;
rho=0.001;
T=5;
dt= 1/250;
t=1:dt:T;
N=length(t);
%Calcul pi_s and pi_l
B=zeros(1,N);
gamma=  0.5;
B =(gamma/theta).*(exp(theta*(T-t))-1);
k =(1-exp(-theta*(T-t)))/theta;
 
pi_s= 1/(1-gamma)*(lambda_s/sigma_s)-(gamma/(1-gamma))* sigma_r/sigma_l .*k;
 
pi_l= 1/(1-gamma)*(lambda_l/sigma_l)-(gamma/(1-gamma))* sigma_r/sigma_l .*k;
 
pi_b= 1-pi_s-pi_l;
%calcul capital adequacy
x = pi_s(1:end-864);
y = pi_l(1:end-864);
alpha1= r' + (x + y).*lambda_e;
alpha2=(x+y).*sigma_e;
alpha3=rho;
beta1=0.2.*x*(r+lambda_s)+0.5.*y*(r+lambda_l)-(0.2.*x*sigma_s).^2-(0.5.*y*sigma_l).^2;
beta2=0.5.*y*sigma_l;
beta3=0.2.*x*sigma_s;
 
xzero=0.92;
T=5;
N=2^7;
dt=T/N;
dw1=sqrt(dt)*randn(1,N);
dw2=sqrt(dt)*randn(1,N);
dw3=sqrt(dt)*randn(1,N);
%Simulation
R=4; dt=R*dt; L=N/R;
xtemp=zeros(1,L);
xtemp=xzero;
for j=1:L
    winc1=sum(dw1(R*(j-1)+1:R*j));
    winc2=sum(dw2(R*(j-1)+1:R*j));
    winc3=sum(dw3(R*(j-1)+1:R*j));
    xtemp=xtemp + (xtemp.*(alpha1 - beta1)- alpha3)*dt - xtemp.*(beta2*winc1 + beta3*winc2 - alpha2*winc3);
    xem(j)=xtemp';
end
plot([0:dt:T],[xzero,xem],'k')
y=0.8;
hold on
hline=refline([0 y]);
set (hline, 'color', 'r')
xlabel('time','Fontsize',12)
ylabel('k(t)', 'Fontsize', 16)
title('Simulation de la dynamique CAR')

[Jerome Briot]

Pourrais-tu nous montrer le message d'erreur complet ?