fftseq.m

来自「双边带调制解调信号代码 用matlab实现（通信原理课程设计）」· M 代码 · 共 13 行

%求傅里叶变换的子函数
function [M,m,df]=fftseq(m,ts,df)
fs=1/ts;
if nargin==2  %nargin为输入参量的个数
    n1=0;
else
    n1=fs/df;
end
n2=length(m);
n=2^(max(nextpow2(n1),nextpow2(n2)));%nextpow2(n)取n最接近的较大2次幂
M=fft(m,n);%M为信号m的傅里叶变换，n为快速傅里叶变换的点数，及为基n-FFT变换
m=[m,zeros(1,n-n2)];%构建新的m信号
df=fs/n;  %重新定义频率分辨率