Error using ss. How do i fix this?

Question

Anna Zito 2020-6-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/557371-error-using-ss-how-do-i-fix-this

评论： Star Strider 2020-7-2

My current code is below. This is my error code:

Error using ss (line 345)
The values of the "a" and "b" properties must be matrices with the same number of rows.
Error in Project43 (line 33)
sys=ss(A,BS,Chead,D);

>>

I have tried everything I know how to do to fix it, but it is not working. Any suggestions?

m1=5.3;
m2=28.5;
m3=8.6;
m4=12.8;
c1=400;
c2=4750;
c3=4585;
c4=4475;
k1=310000;
k2=183000;
k3=162800;
k4=90000;
t=0:0.01:5;
A=[0 0 0 0 1 0 0 0; 0 0 0 0 0 1 0 0; 0 0 0 0 0 0 1 0; 0 0 0 0 0 0 0 1; -k1/m2 k1/m2 0 0 -c1/m2 c1/m2 0 0; k1/m2 -k1/m2 k2/m2 0 c1/m2 (-(c1/m2)-(c2/m2)) c2/m2 0; 0 k2/m3 (-(k2/m3)-(k3/m3)) k3/m3 0 c2/m3 (-(c2/m3)-(c3/m3)) c3/m3; 0 0 k3/m4 (-(k3/m3)-k4/m4) 0 0 c3/m4 (-(c3/m4)-(c4/m4))];
D=[0];
EV=eig(A);
dt=0.01;
tf=5;
t=0:dt:tf;
%for head displacement x1
Chead=[1 0 0 0 0 0 0 0]; 
%For Us 
BS=[0 0 0 0 0 0 k4/m4 0]; 
sys=ss(A,BS,Chead,D); 
[yS,t,x]= step(sys,t) ; 
%For UI 
BI=[0 0 0 0 0 0 0 c4/m4 ]; 
sys=ss(A,B,Chead,D); 
[yI,t,x]= impulse(sys,t) ; 
Yhead=yS+yI; 
plot(t,Yhead) 
% For torso displacement
Ctorso=[ 0 1 0 0 0 0 0 0];
%For Us 
BS=[0 0 0 0 0 0 k4/m4 0]; 
sys=ss(A,BS,ctorso,D); 
[yS,t,x]= step(sys,t) ; 
%For UI 
BI=[0 0 0 0 0 0 0 c4/m4 ]; 
sys=ss(A,BS,Ctorso,D); 
[yI,t,x]= impulse(sys,t); 
Yhead=yS+yI; 
plot(t,Ytorso-Yhead)
%Repeat for head acceleration x1double dot using Cheadacc=[-k1/m1 k1/m1 0 0 –c1/m1 c1/m1 0 0] and plot(t,x1doubledot/9.81) 
Cheadacc=[-k1/m1 k1/m1 0 0 -c1/m1 c1/m1 0 0];
%For Us 
BS=[0 0 0 0 0 0 k4/m4 0]; 
sys=ss(A,BS,Cheadacc,D); 
[yS,t,x]= step(sys,t); 
%For UI 
BI=[0 0 0 0 0 0 0 c4/m4 ]; 
sys=ss(A,BS,Cheadacc,D); 
[yI,t,x]= impulse(sys,t) ; 
Yhead=yS+yI; 
plot(t,x1doubledot/9.81)