FFT画图的方法

plot(-Fs/2:Fs/2-1, fftshift(abs(fft(dpsk,Fs))),'k');

Fs为采样频率

 

假设符号速率为40sps，采样频率为Fs，则每个符号的采样数为：

fd  = Fs/RB;

对于，经过采样后，每个bit扩展fd倍，总数量也扩展了fd倍，代码如下：

for i = 1:symlen
    sbas_rect((i-1)*fd+1:i*fd) = drc(i)*ones(1,fd);
end

再经过载波调制，经过实数载波调制为实信号，代码如下：

dpsk = sbas_rect.*cos(2*pi*fc*t); % 矩形成形

最后画出FFT图如第一行代码所示

若载波为复信号，则

dpsk_c = sbas_rect.*exp(1j*2*pi*fc*t);

复信号FFT频谱画图也和第一行一样。