一个碰撞振动的时域频域响应问题

guleiss

我在分析一个碰撞振动系统冲击响应，直接用ode45求解，时域的曲线没什么问题，但FFT之后，频域的很怪，不知道怎么解释，跟实验结果不同，按照实验结果，除了谐振频率外，振幅都很应该小，但fft分析的频域响应却在0Hz附近的低频时时，振幅也很大，不知道问题出在哪里？想请教大家，谢谢

matlab程序如下
function y=impact1(t,x)
gap=0;
k1=6.49;
k2=591;
f=1.7;
a=0;
c1=2
y=zeros(2,1)
y(1)=x(2)
y(2)=-2*c1*x(2)-(k1*x(1).*(x(1)<gap)+((k1+k2)*x(1)-k2*gap).*(x(1)>=gap))/0.008


[t,x]=ode45(@impact1, [0:0.01:10], [0 0.5])
x1=x(:,1);
x2=x(:,2);
fs=100
N=1024
n=0:N-1

subplot(231)
plot(t,x1)

y1=fft(x1,N);


mag=abs(y1);
f=(0:512)*fs/N;% one side
figure(1);
subplot(232);
plot(f, mag(1:513));

hcharlie

1）这是一个典型的非线性问题，不能用传统的线性问题的结论来思维。
2）碰撞力是单向的，所以会得到单向的0频的响应。