关于脉冲响应函数与频响函数的关系

eric415

书上说脉冲响应函数的fft就是频响函数，我用matlab的impluse函数做的再用fft，发现怎么都得不到，请老师们帮我看看。
我想求出对应于一号点激励一号点响应的一个五峰值系统位移频响函数曲线：
M=diag([50 50 50 50 50]);
K=1e7*[5.8 -2.9 0 0 0; -2.9 5.8 -2.9 0 0;0 -2.9 5.8 -2.9 0 ;0 0 -2.9 5.8 -2.9;0 0 0 -2.9 2.9];
C=1e3*[2 -1 0 0 0;-1 2 -1 0 0;0 -1 2 -1 0;0 0 -1 2 -1;0 0 0 -1 1];

A=[zeros(5,5) ones(5,5);-inv(M)*K -inv(M)*C];
B=[zeros(5,5);inv(M)];
C=[1 zeros(1,9)];
D=zeros(1,5);
sys=ss(A,B,C,D);
[y,t]=impulse(sys);
n=length(y);
Y=fft(y(:,:,1));
plot(1:n,abs(Y))