做fft图像问题。三角函数

y=x^2*cos(x) 做y在fft后的振幅谱。设定fft的点数为1024。

fs = 50;%采样频率
Ts = 1/fs;%采样周期
n = 1024;
x = 0:Ts:(n*Ts);
y = power(x,2).*cos(x);
plot(y);
title('原始信号');
figure;

%傅里叶变换处理
y = fft(y);
f = (0:length(y)-1)*fs/length(y);

plot(f,abs(y))
xlabel('频率(Hz)')
ylabel('振幅')
title('傅里叶变换后振幅谱')