getfft.m

来自「FM 的仿真程序」· M 代码 · 共 16 行

function [M,m_new,df_new]=GetFFT(m,t0,t1,ts,df)
%对信号求FFT，先得出FFT的点数
t=t0:ts:t1;
fs=1/ts;
n1=length(t);
    if nargin==4
        n2=0;
    else
        n2=fs/df;
    end
n=2^(max(nextpow2(n1),nextpow2(n2)));
M=fft(m,n); %m的n点FFT
m_new=[m,zeros(1,n-length(t))]; %新的m，有可能对m进行扩展
df_new=fs/n; %新的频率分辨率