第1个回答 2011-08-24
你只要把图形的数据,带入fft的表达式就行。
例如:y=图形数据;
fy=fft(y(5001:7001));%5001:7001是你要进行傅里叶变换的数据,即第5001个点,到7001个点。因为傅里叶变换是最小二乘的,所以数据长度最好是2的n此方。
f=(0:length(fy)-1)*100000/length(fy); %100000为采样频率
Pyy=abs(fy);
figure(1)
plot(f(1:1001),Pyy(1:1001,'r') %傅里叶变换出来的图像是对称的,你只要取一半的点图形显示 就可以了。
grid on;
第2个回答 2011-08-22
x=0:2*pi/100:2*pi;
y=x.*sin(x);
yy=abs(fft(y));
subplot(211)
plot(x,y)
subplot(212)
stem(0:0.5/pi:50/pi,yy);
xlim([0,25/pi]);
本回答被提问者采纳