求帮忙求一下x1和h的幅频特性曲线，初始值随便取，谢谢啊...
ms=200; 
ks=30000;
mu=40;
Fm=200;
kt=181000;
c0=1000;
B=100;
h=1;
%FM=(Fm*(1-exp(-B*(x(2)-x(4))))/(1+exp(-B*(x(2)-x(4)))));
xr=0.01*sin(2*pi*h*t);

......................