MATLAB fft变换，横坐标为频率

显示全部楼层 · 2009-3-12 08:52:40

如何写呢？
clc;clear;close all;
t = 0:0.001:0.6;
x = sin(2*pi*50*t)+2*sin(2*pi*120*t);
Y = fft(x,512);
py1=Y/512;
py2=py1*2;
% 为什么乘以2
py3=abs(py2);
i=(0:256)*500/256;
% 为什么乘以500
plot(i,py3(1:257))
%py3=abs(py2)*2^(1/2);

千问 · 2009-3-12 08:52:40

可以用自功率密度函数或互功率谱密度函数，给你一个自功率谱：>> fni=input('随机信号谱分析：','s')>> fid=fopen(fni,'r')>> sf=200>> nfft=1024>>fno= ‘自功率谱密度 ‘>> a=fscanf(fid,'%f',[2,inf])>> status=fclose(fid)>> x=a(1,:)>> y=a(2,:)-a(1,:)>> f=0:sf/nfft:sf/2-sf/nfft>> w=hanning(nfft)>> z=psd(y,nfft,sf,w,nfft/2)>> nn=1:nfft/4>> subplot(2,1,1)>>...