[Débutant] Condition initiale fonction de transfert

Bonjour,

S'il vous plait j'ai pas réussi à imposer les condition initiale pour mon système. Pouvez vous m'aider , mon code est le suivant, Je veux que la réponse commence par la valeur 2

Code:

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
 
clear all;
close all;
 
%__________nombre d'itérations____________%
N =10000; 
%__________Allocation des variables _______%
Rdr=zeros(1,N); 
Cdr=zeros(1,N);
C0r=zeros(1,N);
R0r=zeros(1,N);
SOC=zeros(1,N);
%______________________Choix du modèle de système________________________%
 
 
 
        r0 = 6e-3;
        rd = 2e-4;
        cd = 2000;
        c0 = 7e5;
 
b2=r0;
b1=(r0*c0+rd*c0+rd*cd)/(rd*cd*c0);
b0=1/(rd*cd*c0);
a1=1/(rd*cd);
a0=0;
 
 
%______________________génération du Signal d'entrée_______________________%
I=idinput(N,'PRBS',[0 1/120],[0 100]); %Pseudo Random Binary Signal
 
%__________Période d'échantillonnage
T=1;
%________Fonction continue, Fonction discrète, Géneration du sortie________%
 
%Fonction de transfert continue
syst=tf([b2 b1 b0],[1 a1 a0]); 
 
t=0:1:N-1;
[A,B,C,D] = tf2ss([b2 b1 b0],[1 a1 a0]);
sys = ss(A,B,C,D);     
Y=lsim(sys, I, t,[2 2]) ;
plot(Y,'G');