ss狀態空間方程matlab計算

k1=1000; k2=2000; k3=3000; m1=1;m2=2;alpha=1; beta=0.005; M=[m1,0;0,m2]; K=[k1+k2,-k2;-k2,k2+k3]; n=size(M,1); C=alpha*M+beta*K; ssA=[zeros(n),eye(n);-M\K,-M\C]; ssB=[C,M;M,zeros(n)]\[eye(n);zeros(n)]

>>阅读原文<<