2018/8/27
5
Fc1 - F ╳ - Fc2
1/(m2
c2 s
s2 )
x2
╳ -
c1 s
Fc1
╳ - Fk
1/(m1s2) k
x1
clc sys7=tf(d,n) sys0=tf(1,[3.5 0 0]) sys8=series(sys7,sys3) sys1=feedback(sys0,7,-1) sys=feedback(sys5,sys8,-1) sys05=tf([0.5 0],1) sys2=series(sys1,sys05) sys3=feedback(sys2,1,-1) sys4=tf(1,[5.6 0 0]) sys02=tf([0.2 0],1) sys5=feedback(sys4,sys02,-1) [n,d]=tfdata(sys1,'v')
2
function[sysm1, sysm2]=modelm4(sysp) sysp=[12 38 1000 0.1]; m1=sysp(1); m2=sysp(2); [X1_F,X2_F]=modelm4(sysp) k=sysp(3); c=sysp(4); A=[0 0 1 0;0 0 0 1;-k/m1 k/m1 -c/m1 c/m1;k/m2 -k/m2 c/m2 -c/m2 ]; B=[0 0 0 1/m2]'; C=[0 1 0 0 ]; M=[1 0 0 0]; D=0; sys2=ss(A,B,C,D); sys1=ss(A,B,M,D); sysm2=tf(sys2); sysm1=tf(sys1);
2018/8/27 1
function y=f(x) sym x if x<5 y=0 elseif x>10 y=2*x else y=(1-x)/4 end