OFDM信号频谱图

OFDM的原理网上讲解的很多，但是真的画频谱图，连个基本的程序都没有找到。下面附上自己写的MATLAB

clc
clear all;


j=sqrt(-1);

B=100e6;%100Mhz
N=64;%子载波数number of subcurries
T=(N+1)/B;%码元持续时间
f=1/T;
fs=4*B;%采样率
Ts=1/fs;
sample=T/Ts;%采样点数

t=linspace(-T/2,T/2,sample);
st=zeros(1,length(t));
for ii=-N+1:N-1
    st=st+(1/sqrt(N))*exp(j*2*pi*f*ii*t);
end
 
%% 频谱
freq=linspace(-fs/2,fs/2,length(st));
fft_y=abs(fft(st));
plot(freq*1e-6,fftshift(fft_y/max(fft_y)));%归一化的频谱
% plot(freq*1e-6,fftshift(fft_y));
xlabel('Frequency in MHz');
title('Magnitude spectrum of chirp signal');
grid on;axis tight;

这只是基本的载波信号的频谱，当然要加相应的信号，但是加过随机信号之后的频谱怎么画带宽啥啥的都统统不对。再研究，之后附上。