直接使用fft()绘制% 两个频率分别为15HZ 和 20HZ 的正弦信号 Fs=50;%采样频率50Hz f1=15; f2=20; t = 0:1/Fs:10-1/Fs % 0-9.98s 一共500个点 x = sin(2*pi*f1*t) + sin(2*pi*f2*t);%原始信号 N=length(x)% N=500 f…
FFT 与matlab 1) fft(x(n),N) (对于输入序列长度为N的序列,不填写N会默认N点DFT),出来的结果是【1:N】个点,就刚好是序列的DTFT的第一个主值序列,所以,为了观测方便,只需要显示范围为【1 :N/2】个点就可以(我就很喜欢这么干!);或者想要转换成双边谱去看,也是可以的,直接整体向左平移N个点,亦或者,...
% 当画出的频谱图在0频率时出现峰值,可将该选项置1,可以有效抑制0频率的峰值,默认情况不去趋势% options.analysisType :频谱分析类型,当options.analysisType = 'fft'时,进行快速傅里叶变化FFT% 当options.analysisType = 'psd'时,计算功率谱% options.psdUnit :功率谱单位,当options.psdUnit = 1时,单位为...
AI代码助手复制代码 显然,该信号有三个频率分量,直流(0频),5Hz和10Hz,对应的幅度分别为3、1、2,其时域波形图如下: 转换为频域的关键函数是X = fftshift(fft(x./(n)));一定注意需要除以总样本数n,然后用fftshift将曲线挪一下位置。 做完这步操作后,得出的是双边频谱,频率范围从-fs/2到fs/2,这是因为...
直接使用fft()绘制 % 两个频率分别为15HZ 和 20HZ 的正弦信号 Fs=50;%采样频率50Hz f1=15; f2=20; t = 0:1/Fs:10-1/Fs % 0-9.98s 一共500个点 x = sin(2*pi*f1*t) + sin(2*pi*f2*t);%原始信号 N=length(x)% N=500 figure(1); plot(t,x); title('Original Signal'); xlabel...