自抗扰matlab程序代码

m1=100;
k1=490000;
yipu1=0.1;
c1=yipu1*2*sqrt（m1*k1）;

m2=100;
k2=490000;
yipu2=0.1;
c2=yipu2*2*sqrt（m2*k2）;

num1=[c1 k1];
den1=[m1 c1 k1];

num21=[c1k1];
num22=