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

3 vues (au cours des 30 derniers jours)
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?

Réponses (0)

Catégories

En savoir plus sur Modeling dans Help Center et File Exchange

Produits


Version

R2020b

Community Treasure Hunt

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

Start Hunting!

Translated by