做fft图像问题。三角函数

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

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

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('频谱图');

效果

img

您好，我是有问必答小助手，您的问题已经有小伙伴帮您解答，感谢您对有问必答的支持与关注！
PS：问答VIP年卡 【限时加赠：IT技术图书免费领】，了解详情>>> https://vip.csdn.net/askvip?utm_source=1146287632