matlab

[afef1]

bonjour tout le monde

Xn+1=A*(Sn+1)*Xn + b * (Sn+1) ................. (1)
et on peut la écrire a Xn+1=A(i)Xn + b(i)
les variables aléatoires Sn sont iid a valeur dans (1,2,.......,K)

sont programme qui permet de simuler une réalisation x1;.........;xn du modèle (1) pour une taille donnée n .

function [X]=Exo401(n,x0,A,B,p)
X(:,1)=x0;
P=cumsum(p);
l=size(A);
m=length(p);
for i=2:n+1
    y=rand;
    j=1;
    while y>P(j)
        j=j+1;
    end
    a=A(:,(l(1)*(j-1))+1:l(1)*j);
    b=B(:,j);
    X(:,i)=a*X(:,i-1)+b;
end
X=X(:,2:n+1);
plot(X)

s'il vous plaît est qu'il y'a quelqu'un qui peut me l'explique
-----