做fft图像问题。三角函数

Question

y=x^2*cos(x)对y做fft图。设定fft的点数为1024

Answer 1

代码供参考，有帮助望采纳哟

x = linspace(0,2*pi,1024);
Fs = 1/(x(2)-x(1));
y = x.^2.*cos(x);
N = length(y);               % 求取抽样点数
yy = fft(y);                  % 对信号进行傅里叶变换
f = Fs/N*(0:round(N/2)-1);   % 显示实际频点的一半，频域映射，转化为Hz
plot(f,abs(yy(1:round(N/2))),'b');% 画频谱图
xlabel('频率/Hz');ylabel('幅值');
title('频谱图');

效果