Azzera filtri
Azzera filtri

How to change continuous state space model to discrete state space model with disturbance and display the discrete state space?

1 visualizzazione (ultimi 30 giorni)
This is parameter that i used in my state space model
% Parameter Massa
m1 = 8095; % massa train set 1 dalam kg
m2 = 8500; % massa train set 2 dalam kg
m3 = 8457; % massa train set 3 dalam kg
% Parameter Gaya
f1 = 205.*10^3; % Gaya train set 1 dalam N
f2 = 302.*10^3; % Gaya train set 2 dalam N
f3 = 302.*10^3; % Gaya train set 3 dalam N
% Parameter Resistansi
c_0_1 = 0.01176;
c_1_1 = 0.00077616;
c_2_1 = 4.48 ;
c_0_2 = 0.01176 ;
c_1_2 = 0.00077616;
c_2_2 = 4.48;
c_0_3 = 0.01176 ;
c_1_3 = 0.00077616;
c_2_3 = 4.48;
% Desired Speed
v_0 = 300;
% Desired Train Following Headway Time
hstar = 120;
a_1 = -1./m1.*(c_1_1 + 2.*c_2_1.*v_0);
a_2 = -1./m2.*(c_1_2 + 2.*c_2_2.*v_0);
a_3 = -1./m3.*(c_1_3 + 2.*c_2_3.*v_0);
a_1_head = 1-(a_1.*hstar);
a_2_head = 1-(a_2.*hstar);
a_3_head = 1-(a_3.*hstar);
b = 1;
p_1 = -1./m1.*(c_0_1 - c_2_1.*(v_0).^2);
p_2 = -1./m2.*(c_0_2 - c_2_2.*(v_0).^2);
p_3 = -1./m3.*(c_0_3 - c_2_3.*(v_0).^2);
A = [0 a_1_head 0 0 0 0;
0 0 0 0 0 0;
0 (a_2_head-1) 0 a_2_head 0 0;
0 0 0 0 0 0;
0 (a_3_head-1) 0 (a_3_head-1) 0 a_3_head;
0 0 0 0 0 0;
];
B = [-b*hstar 0 0;
b 0 0;
0 -b*hstar 0;
-b b 0;
0 0 0;
0 0 -b;
];
C = [1 0 0 0 0 0;
0 0 1 0 0 0;
0 0 0 0 1 0;
];
% W is bounded disturbance
W = [((a_1 - 1).*v_0) - (p_1.*hstar);
0;
((a_2 - 1).*v_0) - (p_2.*hstar);
((a_1 - 1).*v_0) - (p_1*hstar);
0;
((a_3 - 1).*v_0) - (p_3*hstar);
];
D = 0;
Ts = 100;
x0 = [0; 0; 100; 0;200;0];
and then i use idss to change my continuous state space model to discrete state space model, but without my disturbance
sys = idss(A,B,C,D);
display(sys);
How do i change my continuous state space model with disturbance to discrete state space moodel?
and
How i display my discrete time state space model like A_discrete, B_discrete, W_discrete?

Risposte (0)

Categorie

Scopri di più su Sequence and Numeric Feature Data Workflows in Help Center e File Exchange

Prodotti


Release

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by