本人最近搞懂了OFDM的一些知识，便给本章取名为第7章——OFDM信道估计、均衡。本文所有代码下载地址是：123kevin456/OFDM-

一、OFDM相关书籍

OFDM书籍随便一查有非常多，缺少均衡、同步相关的代码实现，便难体会一些通信过程。我最近正在看的两本OFDM书籍，都是理论+代码结合来讲解，推荐给你。

1、《MIMO-OFDM无线通信技术及MATLAB实现》

此书最大特点便是有MATLAB源代码下载！下载！！下载！！！但由于代码版本可能老旧，因此需要偶尔改改代码适应新的MATLAB版本才可以正常运行。我试过OFDM的几个代码，结果还是不错的。

2、《Simulation and Software Radio for Mobile Communications》

此书我目前还没有发现中文版本，因此直接看的英文版本。这本书介绍了PSK调制、OFDM、CDMA、多址接入协议等内容的仿真，均有MATLAB源代码。

就OFDM来说，第1本书比第2本书，难度更大，且在均衡、同步，比如载波同步、位定时同步等内容，讲得更加详细。由于难度太大，我已琢磨许久，偶有所得，便在此更新。

因此，本次文章后面分享的MATLAB可运行代码主要参考是第2本书。

二、OFDM基本原理

关于OFDM原理的基本介绍，知乎上的 @子木在其《OFDM完整仿真过程及解释（MATLAB）》 已经介绍非常清楚了，本系列文章主要是作为补充，尤其是将OFDM经过瑞利信道时的误码率性能仿真，以及如何通过信道估计来对信道进行补偿（compensation）。

图1 OFDM仿真系统框图

在这次的OFDM仿真中，为方便调试程序，首先做如下假设：
（1）调制：采用QPSK调制，可自己编写，也可以使用系统自带的pskmod函数；

（2）信道编码：一般采用（2，1，7）卷积码，本次仿真无编码；

（3）没有上载波与下载波过程，同时也没有发送端的上采样，接收端的匹配滤波，寻找最佳采样点过程；

待最基本的程序仿真正确后，再对其加入相应模块，以做完善。

三、OFDM经过高斯白噪声信道的误码率分析

本次仿真主要集中在将OFDM信号经过瑞利信道，是否有均衡，来观察误码率情况。

当然，在经过瑞利信道之前，不妨先测试一下OFDM在高斯白噪声信道下的传输。代码无误后，再将信道改成瑞利衰落信道。

先给出主要代码，各个函数的代码附在最后：

%%%%%%%%%%%%%%%%%%%%% OFDM仿真 %%%%%%%%%
%%%%%%%%%%%%%%%%%%%%% ofdm.m  %%%%%%%%%
%%%%%%%%%  date:2020年10月11日  author:飞蓬大将军 %%%%%%%%%%%%%%%程序说明
%%%高斯白噪声信道下OFDM传输%%%%    仿真环境
%软件版本：MATLAB R2019a%********************** 程序主体 ************%%%%%%%%%%%%%%%%%%%%%%  参数设置   %%%%%%%%%%%%%%%%%%%
para = 128;   %Number of parallel channel to transmit
fftlen = 128;  %FFT length
noc = 128;    %Number of carrier 
nd = 6;  %Number of information OFDM symbol for one loop
ml = 2;   %Modulation：QPSK
sr = 250000;  %Symbol rate 符号速率
br = sr.*ml;  %Bit rate per carrier
gilen = 32; %length of guard interval ebn0_temp = 3:1:10;
for kkk = 1:length(ebn0_temp)ebn0 = ebn0_temp(kkk);  %Eb/No 
% ebn0 = 8;  %Eb/Nonloop = 10000; %Number of sumulation loops 
noe = 0;   %Number of error data  
nod = 0;   %Number of transmitted data 
eop = 0;   %Number of error packet  
nop = 0;   %Number of transmitted packetfor iii = 1:nloop%%%%%%%%%%%%%%%%%  发射机  %%%%%%%%%%%%%%%%%%%seldata = rand(1,para*nd*ml)>0.5;  %串行数据paradata = reshape(seldata,para,nd*ml); %串并转换[ich,qch] = qpskmod(paradata,para,nd,ml); %调制kmod = 1/sqrt(2);ich1 = ich.*kmod;qch1 = qch.*kmod;%%%%%%%%%%%% IFFT %%%%%%%%%%x = ich1 + qch1 *1j;spow1 = sum(sum(ich1.^2+ qch1.^2))/nd./para;y = ifft(x);ich2 = real(y);qch2 = imag(y);spow2 = sum(sum(ich2.^2+ qch2.^2))/nd./para;%%%%%%%% 添加保护间隔 %%%%%%%%[ich3,qch3] = giins(ich2,qch2,fftlen,gilen,nd);fftlen2 = fftlen + gilen;%     figure(1);
%     plot(ich3,'-+');
%     hold on;%****************  Attenuation Calculation **************spow = sum(ich3.^2+ qch3.^2)/nd./para;attn = 0.5*spow*sr/br*10.^(-ebn0/10);attn = sqrt(attn);%%%%%%%%%%%%%%%%%%%%%%%  接收机  *************%%%%%%%%%%%%%%%%% AWGN addition *******************[ich4,qch4] = comb(ich3,qch3,attn);%     plot(ich4,'-*');%%%%%%%% 去掉保护间隔 %%%%%%%%[ich5,qch5] = girem(ich4,qch4,fftlen2,gilen,nd);%%%%%%%%%%%%%%  FFT %%%%%%%%rx = ich5 + qch5.*1i;ry = fft(rx);ich6 = real(ry);qch6 = imag(ry);%%%%%%%%%%%%%% demoluation %%%%%%%%%%%%%%ich7 = ich6./kmod;qch7 = qch6./kmod;demodata = qpskdemod(ich7,qch7,para,nd,ml);%%%%%%%%%%%%%% 并串转换 %%%%%%%%%demodata1 = reshape(demodata,1,para*nd*ml);%%%%%%%%%%%%%%% Bit Error Rate %%%%%%%%%%%noe2 = sum(abs(demodata1-seldata));nod2 = length(seldata);%%%%cumulative the number of error and data in noe and nodnoe = noe + noe2;nod = nod + nod2;%%%计算PERif noe2~=0eop = eop +1;
%     else
%         eop = eop;endnop = nop + 1;
%     fprintf('%f\t%e\t%d\n',iii,noe2/nod2,eop);
end%**************   output result *************
per = eop/nop;
ber = noe/nod;fprintf('%f\t%e\t%e\t%d\t\n',ebn0,ber,per,nloop);
fid = fopen('BERofdm.dat','a');
fprintf(fid,'%f\t%e\t%e\t%d\t\n',ebn0,ber,per,nloop);
fclose(fid);end%************************* 画误码率曲线进行对比 ******************%
% ebn0 = 3:1:10;
% awgn_theory = [0.0228784075610853,0.0125008180407376,0.00595386714777866,0.00238829078093281,0.000772674815378444,0.000190907774075993,3.36272284196176e-05,3.87210821552205e-06];
% awgn_no_compensation = [3.698496e-02,2.254329e-02,1.226654e-02,5.823633e-03,2.305339e-03,7.492187e-04,1.757812e-04,3.170573e-05];
% rayleign_one_path_theory = [0.125000000000000,0.100000000000000,0.0833333333333333,0.0714285714285715,0.0625000000000000,0.0555555555555556,0.0500000000000000,0.0454545454545455]; 
% 
% semilogy(ebn0,awgn_theory,'-*',ebn0,awgn_no_compensation,'-+');
% xlabel('比特信噪比');
% ylabel('误码率');
% title('不同信噪比下误码率仿真曲线');
% legend('理论曲线','实验曲线');
% grid on;%***********   高斯白噪声下的OFDM采用QPSK调制的误码率值 **********%%%%%%%%%%%%%%%%%%%%%%%%      理论值          **************%
%%%%%%%%%%%%%     EbN0(dB)      误码率        
%%%%%%%%%%%%%       3        0.0228784075610853
%%%%%%%%%%%%%       4        0.0125008180407376
%%%%%%%%%%%%%       5        0.00595386714777866
%%%%%%%%%%%%%       6        0.00238829078093281
%%%%%%%%%%%%%       7        0.000772674815378444
%%%%%%%%%%%%%       8        0.000190907774075993
%%%%%%%%%%%%%       9        3.36272284196176e-05
%%%%%%%%%%%%%      10        3.87210821552205e-06%%%%%%%%%%%%%%%%%%%%%        实验值        *******%%%%%%%%%%%%
%%%%%%%%%%%%%      EbN0(dB)      误码率         误包率        循环次数
% % % % % % % % % 3.000000	3.698496e-02	1.000000e+00	10000	
% % % % % % % % % 4.000000	2.254329e-02	1.000000e+00	10000	
% % % % % % % % % 5.000000	1.226654e-02	1.000000e+00	10000	
% % % % % % % % % 6.000000	5.823633e-03	9.999000e-01	10000	
% % % % % % % % % 7.000000	2.305339e-03	9.709000e-01	10000	
% % % % % % % % % 8.000000	7.492187e-04	6.837000e-01	10000	
% % % % % % % % % 9.000000	1.757812e-04	2.384000e-01	10000	
% % % % % % % % % 10.000000	3.170573e-05	4.780000e-02	10000	    %%%%%%%%%%%%%%%%%   结论    %%%%%%%%%%%%%%%%%%%%
%完成了OFDM经过高斯白噪声信道的误码率仿真
%由于OFDM中引入保护间隔，是一种冗余信息，因此相比于理论误码率曲线有10*log（160/128）=0.969dB的损失
%对应到误码率曲线上，即本实验误码率曲线相比理论误码率曲线向右平移0.969dB
%2020年11月7日   

图2 OFDM采用QPSK在AWGN下的误码率曲线图

在上面的代码中，注意两点：

（1）当信号经过FFT之后，或者IFFT之后，功率会发生变化，这是为什么呢？功率的数量关系如何？

由N点DFT的表达式，知： $X\left(k\right)=\sum _{n=0}^{N-1}x\left(n\right)e^{-j\frac{2\pi }{N}kn},k=0,1,\cdots ,N-1$

为理解信号经过FFT前后功率发生的变化，联想到概率论中的公式： $D\left(ax+by\right)=a^{2}D\left(x\right)+b^{2}D\left(y\right)$

因为 $E\left\{x{\left(0\right)}^2\right\}=E\left\{x{\left(1\right)}^2\right\}=\cdots =E\left\{x{\left(N-1\right)}^2\right\}$

所以有： $E\left\{X{\left(k\right)}^2\right\}=N\bullet E\left\{x{\left(0\right)}^2\right\}$

即信号经过FFT之后，功率变为原来的N倍。同理可以得到，信号经过IFFT后，功率变为原来的 $\frac{1}{N}$ 。

这样你便能理解ofdm.m代码中，发送端在经过IFFT后，才计算信号功率spow，再根据信噪比和spow，去计算相应的噪声功率。

（2）仿真曲线相比于理论曲线看起来是向右平移的关系，所以是平移的关系码？如果是平移，那到底是平移了多少dB呢？

由于在本次仿真中，所有子载波采用相同的调制方式。当然，没有规定就说所有的子载波一定要采用相同的调制方式。

因此，既然子载波是采用相同的QPSK调制方式，那么理论线就是QPSK的误码率曲线，QPSK和2PSK的误码率曲线是一样的。

如果对QPSK和2PSK的误码率曲线是否一致，还未想清楚的话，请点击《陈老湿：问答交流1：采用过采样的原因》，看问题2的解释。

结合《第二章：线性分组码》中讲过信息比特与信道比特的关系。在上面的代码中，由于加入了保护间隔，保护间隔是采用循环前缀（Cyclic Prelix,CP）方式，即将OFDM符号中的后一部分进行复制，加入到OFDM前面，构成一个完整的OFDM符号。

复制的这一段信号，占用了功率资源，没有用来发送新的信息比特，因此便有了功率损失。注意到子载波数目是128，CP长度是32，即CP长度占整个OFDM符号长度的1/5。

举个数值例子方便你理解，本来你有1w的功率用来发送1个比特，现在你只能有0.8w用来发这个比特了，因为还有0.2w要用去给CP了。

所以相比下来会损失：10log10(1)-10log10(4/5)=0.9691dB。对应到本次实验中，即实验曲线会相比理论曲线有向右平移0.9691dB。

你也可以根据《第二章：线性分组码》中Eb/N0与SNR的换算关系，来将实验曲线进行向左平移0.9691dB，便可发现二者相近了。

那么问题来了，CP有什么用呢？CP的长度如何取呢？（思考几秒钟？）

在多径信道中，当信号带宽大于无线信道的相干带宽时，链路会受到多径衰落的影响，从而造成ISI（符号间干扰）。插入CP，可以避免ISI。当然，这要求CP的长度要大于最大多径时延扩展。

上面这段话，我不知道你看着是什么感觉？书上都是这么写的，以及上课老师也是这么讲的？但都没有说清楚为什么CP长度大于最大多径时延扩展，便可以避免ISI。

最大多径时延扩展，可以理解对一个符号最晚到达路径相对于最早到达路径的时间差。

下面从数学表达式来帮助你理解，这需要你有循环卷积和线性卷积的概念。

如图4.14所示，如果保护间隔长度小于最大多径时延扩展，对于前一个OFDM符号的尾部已经影响到下一个OFDM符号的前端，这样就造成ISI。在实际中，可能会出现符号定时偏差（STO），STO会使OFDM信号的前端和FFT窗的起始点不一致。

如图4.15所示，若CP的长度大于最大多径时延扩展，也可能会出现ISI和ICI（子载波干扰），这取决于FFT的起始点的定时。如果FFT窗的起始点早于时延的前一个符号末端，则会出现ISI，如果FFT窗的起始点晚于符号的起始点，既会出现ISI，又会出现ICI。

因此，在本次的AWGN信道仿真中，如果不加入保护间隔，也是可以的，加入CP对AWGN信道来说没起到实际的用处。如果你把保护间隔去掉，便就可以得到误码率的实验曲线和理论曲线非常接近了。

因此，当保护间隔采用循环前缀（Cyclic Prelix,CP）方式时，即将OFDM符号中的后一部分进行复制，加入到OFDM前面，构成一个完整的OFDM符号，保证了发射采样和信道采样的循环卷积性质。

注意到循环卷积的理解，下一节将根据循环卷积原理来对信道进行补偿。

至此，已经讲完了CP的作用，以及怎么取。

三、总结

CP在高斯白噪声信道下起不到作用，那么如何看到OFDM在多径瑞利信道传输，CP发挥作用呢？

我们便首先要知道瑞利信道如何进行建模与仿真，之后便可以看到均衡对于信道补偿的作用。我们下次再讲。

也欢迎读者朋友就相关技术问题与我交流，一起学习，共同进步。请你也别忘了把这篇文章分享给你身边正在学习通信专业的同学们，也许能够帮到Ta。这是《陈老湿·通信MATLAB》仿真的第7章，期待下次更新见！