Dynamical System behaving strangely under high frequencies input
Show older comments
Dear all,
I have been simulating the following dynamical system in matlab:
w=0.04*(1-cos(8*pi*t)).*(heaviside(t-0.5)-heaviside(t-0.9));
g=9.81;
m1 = 290;
m2 = 59;
k1 = 16812;
k2 = 190000;
b1 = 1000;
alfa=4.515*(10^13);
beta=1;
gama=1.545*(10^9);
tau=(1/30);
Ps=10342500;
A=3.35*(10^(-4));
dxdt=zeros(6,1);
dxdt(1) = x(2);
dxdt(2) = -(k1/m1)*x(1)-(b1/m1)*x(2)+(k1/m1)*x(3)+(b1/m1)*x(4)+(A/m1)*x(5);
dxdt(3) = x(4);
dxdt(4) = (k1/m2)*x(1)+(b1/m2)*x(2)-((k1+k2)/m2)*x(3)-(b1/m2)*x(4)-(A/m2)*x(5)+(k2/m2)*w;
dxdt(5) = -alfa*A*x(2)+alfa*A*x(4)-beta*x(5)+gama*x(6)*sign(Ps-(x(5)*sign(x(6))))*sqrt(abs(Ps-(x(5)*sign(x(6)))));
dxdt(6) = -(1/tau)*x(6)-(1/tau)*((K(1,1))*x(1)+(K(1,2))*x(2)+(K(1,3))*x(3)+(K(1,4))*x(4)+(K(1,5))*x(5)+(K(1,6))*x(6));
This is a controller for a quarter car vehicle suspension. Its behaviour is fine as long as I keep the input frequency (inside the cos function) low. As I increase the frequency, states behaviour change and it turns out that oscilations never die. I have been thinking for days and I can not find a reason for this to happen. Could someone, please help me to discover what is going on?
Thank you all very much.
3 Comments
Walter Roberson
on 5 Jul 2016
Could you clarify whether this is a different question than your earlier http://www.mathworks.com/matlabcentral/answers/293568-problem-simulating-model-in-simulink-strange-results-in-high-frequency that looks like the same topic?
Allan Andre do Nascimento
on 5 Jul 2016
Allan Andre do Nascimento
on 5 Jul 2016
Answers (0)
Categories
Find more on Ordinary Differential Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!