移动通信系统OFDM系统仿真与实现基于MATLAB文档格式.docx

1、尽管峰值功率出现的概率较低，但为了不失真地传输这些高PAPR的OFDM信号，从而导致发送效率极低，接收端对前端放大器以及A/D变换器的线性度要求也很高。因此提出了基于信号畸变技术、信号扰码技术和基于信号空间扩展等降低OFDM系统PAPR的方法。（5） 自适应技术在OFDM系统中使用自适应技术，还应考虑频率分组、时间间隔、信道总延迟和信道估计误差等因素，其中信道估计误差对性能的影响较大。2.2 OFDM的优缺点OFDM主要有下列一些优点:（1） OFDM在对抗干扰及衰落的优势：把高速数据流通过串并变换，使得每个子载波上的数据符号持续长度相对增加，可有效对抗信号波形间的干扰ISI，适用于多径信道存

2、在频率选择性，而所有子载波都处于深衰落的概率极小，OFDM系统可通过动态比特分配和动态子信道分配的方法，充分利用信噪比较高的子信道，提高系统性能。（2） OFDM系统由于子载波之间存在正交性，允许子信道的频谱相互重叠，因此与常规的频分复用相比，OFDM可以最大限度地利用频谱资源。（3） 各子信道的正交调制和解调可通过离散傅利叶反变换（IDFT）和离散傅利叶变换（DFT）实现。对N很大（N32）的系统，可以通过快速傅立叶变换（FFT）来实现。基于FFT的系统在计算方面更有效，并且随着大规模集成电路技术和DSP的发展，IFFT和EFT都非常容易实现。（4） OFDM的开放灵活性：无论从无线数据业务

3、的使用需求,还是从移动通信系统自身要求，都希望物理层支持非对称高速数据传输,而OFDM系统可以很容易使用不同数量的子载波来实现上下链路中不同的传输速率。OFDM技术的不足: 由于OFDM系统内存在多个正交子载波，而且其输出信号是多个子信道的叠加，因此对子信道的正交性有严格要求。而由于无线信道的时变性，还有发射机载波和本地振荡器的频率偏差，所以OFDM易受频率偏差的影响。如果多个子信号的相位一致时，所得到的叠加信号的瞬时功率会远大于信号的功率，出现较大的峰值与均值功率比（PAR），这个比值的增大会降低射频放大器的功率效率，使系统性能恶化。本文主要研究了OFDM调制解调系统，对噪声信道抗干扰能力的

4、研究以及在实际中的基本应用。利用Matlab进行仿真和性能分析。2.3 原理及数学描述2.3.1 OFDM的载波调制正交频分复用OFDM（OrthogonalFrequencyDivisionMultiplex）是一种多载波调制方式，通过减小和消除码间串扰的影响来克服信道的频率选择性衰落。它的基本原理是将信号分割为N个子信号，然后用N个子信号分别调制N个相互正交的子载波。由于子载波的频谱相互重叠，因而可以得到较高的频谱效率。下图是OFDM基带信号处理原理图。其中，（a）是发射机工作原理，（b）是接收机工作原理。图2.1 OFDM基带信号处理原理图2.3.2 OFDM的调制解调原理OFDM通过把

5、需要发射的数据流分解为若干个并行的数据子流，这样每个数据子流在速率上就会降低很多，然后再进行相关调制，将它们调制到一组总数为N，频率之间的间隔相等，且又两两正交的子载波上。OFDM的调制可以用离散傅里叶反变化（IDFT）来实现，相应的有在解调端可以用离散傅里叶变换（DFT）来实现。而这两种傅里叶变换都有相应的快速算法，在系统的效率和相应时间上会有所提高。系统的实际应用中一般采用IFFT和FFT技术。2.3.3串并转换数据传输的典型形式是串行数据流，符号被连续传输，每一个数据符号的频谱可占据整个可利用的带宽。但在并行守护据传输系统中，许多符号被同时传输，减少了那些在串行系统中出现的问题。在OFD

6、M系统中，每个传输符号速率的大小大约在几十bps到几十Kbps之间，所以必须进行串并变换，将输入串行比特流转换为可以传输的OFDM符号。由于调制模式可以自适应调节，所以每个子载波的调制模式是可变化的，因而每个子载波可传输的比特数也是可以变化的，所以串并转换需要分配给每个子载波数据段的长度是不一样的。在接收端执行相反的过程，从各个子载波处来的数据被转换回原始的串行数据。2.3.4保护间隔和循环前缀应用OFDM的一个重要原因在于它可以有效的对抗多径时延扩展。通过把输入数据流串并变换到N个并行的子信道中，使得每一个调制子载波的数据周期可以扩大为原始数据符号周期的N倍。为了最大限度的消除符号间干扰，还

7、可以在每个OFDM符号间插入保护间隔（GI），而且该保护间隔长度一般要大于无线信道中的最大时延扩展，这样一个符号的多径分量就不会对下一个符号造成干扰。在这段保护间隔内，可以不插入任何信号，即是一段空闲的传输时段。2.3.5 OFDM的基本参数各种OFDM参数的选择就是需要在多项要求冲突中进行折中考虑。通常来讲，如前所述，首先要确定3个参数：带宽、比特率以及保护间隔，按照惯例，保护间隔的时间长度应该为应用移动环境信道的时延均方根值的24倍。一旦确定了保护间隔，则OFDM符号周期长度就可以确定。为了最大限度的减少由于插入保护间隔所带来的信噪比的损失，希望OFDM符号周期长度要远远大于保护间隔长度。

8、但是符号周期长度又不可能任意大，否则OFDM系统中包括更多的子载波数，从而导致子载波间隔相应减少，系统的实现复杂度增加，而且还加大了系统的峰值平均功率比，同时使系统对频率偏差更加敏感。因此在实际应用中，一般选择符号周期是保护间隔长度的5倍，这样由于插入保护比特所造成的信噪比损耗只有1dB左右。2.3.6调制方式可以通过改变发射的射频信号的幅度、相位和频率来调制信号。对于OFDM系统来说，只能采用前两种调制方法，而不能采用频率调制的方法，这是因为子载波是频率正交，而且携带独立的信息，调制子载波频率会破坏这些子载波的正交特性，这是频率调制不能在OFDM系统中采用的原因。2.3.7 16QAM调制1

9、6QAM 是用两路独立的正交 4ASK 信号叠加而成，4ASK 是用多电平信号去键控载波而得到的信号。它是 2ASK 体制的推广，和 2ASK 相比，这种体制的优点在于信息传输速率高。16QAM 信号采取正交相干解调的方法解调，解调器首先对收到的 16QAM 信号进行正交相干解调，一路与 cos c t 相乘，一路与 sin c t 相乘。然后经过低通滤波器，低通滤波器 LPF 滤除乘法器产生的高频分量，获得有用信号，低通滤波器LPF 输出经抽样判决可恢复出电平信号。16 进制的正交振幅调制是一种振幅相位联合键控信号。16QAM 的产生有 2 种方法：（1）正交调幅法，它是有 2 路正交的四电

10、平振幅键控信号叠加而成；（2）复合相移法：它是用 2 路独立的四相位移相键控信号叠加而成。2.3.8 QPSK调制QPSK是在2PSK（二相调相）的基础上发展而来的一种多进制相位调制。QPSK是利用载波的四种不同相位来表征传送的数字信息。在QPSK调制中，首先对输入的二进制数据按二位数字编成一组，以此构成双比特码元。其组合共有4种，即有4种不同状态。故可以用M=4种相位或相位差来表示。这里M=4，故称为四相调相。同样，若采用八相调制方式，在一个码元时间内可传送3位码，其信息传送速率是二相调制方式的3倍。由此可见，采用多相调制的级数愈多，系统的传输速率愈高，但相邻载波之间的相位差愈小，接收时要区

11、分它们的困难程度就愈大，将使误码率增加。3. 基于MATLAB的OFDM系统仿真实验3.1 MATLAB的特点与功能MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计

12、算软件的先进水平。3.2 MATLAB仿真系统3.2.1 仿真结构图图3.1 系统仿真原理图3.2.2 主要模块（1）串并转换 OFDM系统是多载波传输的一种特殊方式，而发送端产生的数据为串行的，因此需要进行串并转换。OFDM将高速输入的串行数据比特流转换成并行传输的低速数据流，而且串并转换之后会提高系统的抗干扰能力。串并转换的示意图如图3.2。图3.2 串并转换示意图（2）QPSK调制QPSK调制的方法有两种，一种是数字方法，另一种是相位选择法。本程序中选择的是利用数字的方法来调制信号。这里通过调用一个子函数qpskmod，程序如下：function iout,qout=qpskmod（pa

13、radata,para,nd,ml）m2=ml./2;paradata2=paradata.*2-1;count2=0;for jj=1:nd isi = zeros（para,1）; isq = zeros（para,1）; for ii = 1 : m2 isi = isi + 2. （ m2 - ii ） .* paradata2（1:para）,ii+count2）; isq = isq + 2. （ m2 - ii ） .* paradata2（1:para）,m2+ii+count2）; end iout（1:para）,jj）=isi; qout（1:para）,jj）=isq;

14、 count2=count2+ml;end（3）IFFT/FFT运算通过N点的IDFT运算，把频域数据符号变换为时域数据信号，经调制后发送到信道中去。在实际应用中，一般采用FFT/IFFT运算，因为它可以显著的降低运算的复杂度。在Matlab软件里可以使用函数fft（）和ifft（）来对数据进行FFT/IFFT运算，可以省去很多复杂的运算。（4）保护间隔和循环前缀 在OFDM系统中，保护间隔是一种循环复制增加了符号的波形长度，在符号的数据部分，每一个子载波内有一个整数倍的循环，此种符号的复制产生了一个循环的信号，即将每个OFDM符号的后一段时间的样点复制到OFDM符号的前面，形成前缀，在交接点

15、没有任何的间断。当信道的最大多径时延扩展小于保护间隔的时间时，OFDM系统可以完全克服ISI的影响。具体的m语言实现为：ich3=ich2（fl-gl+1:fl,:）;ich2;qch3=qch2（fl-gl+1:qch2;（5）并串转换 这一过程是串并转换的逆过程，将N个子载波的数据传送到一个载波信道中去，将并行数据转换为串行数据序列进行传输。示意图如图3.3。图3.3 并串转换示意图（6）加入高斯噪声 白噪声是根据噪声的功率谱密度是否均匀来定义的，而高斯噪声则是根据它的概率密度函数呈正态分布来定义的，高斯型白噪声称高斯白噪声，在通信系统的理论分析中，特别是在分析、计算系统抗噪声性能时，经常

16、假定系统中信道噪声（即前述的起伏噪声）为高斯型白噪声。在matlab软件里产生高斯噪声的函数一般有两个WGN和AWGN。Awgn的用法为：y = awgn（x,SNR） 在信号x中加入高斯白噪声。信噪比SNR以dB为单位。x的强度假定为0dBW。如果x是复数，就加入复噪声。y = awgn（x,SNR,SIGPOWER） 如果SIGPOWER是数值，则其代表以dBW为单位的信号强度；如果SIGPOWER为measured，则函数将在加入噪声之前测定信号强度。本程序正是采用的ReData=awgn（TrData，SNR，measured）来给复数发射数据TrData加入高斯噪声。（7）QPSK解

17、调在进行解调之前，信号要完成去掉保护间隔和进行FFT运算的任务。去掉保护间隔也就可以去掉符号间的干扰。这里通过调用一个子函数qpskdemod，具体的m语言实现为：function demodata=qpskdemod（idata,qdata,para,nd,ml）demodata=zeros（ para,ml*nd ）;demodata（1:para）,（1:ml:ml*nd-1）=idata（1:nd）=0;para）,（2:ml*nd）=qdata（1:（8）接收信号 完成解调之后下面就是进行解调信号的判决得到接收信号，程序中将发送段的原始信号和接收端的信号的波形输了出来。（9）MATL

18、AB仿真源程序：clear;clc;SNR=10;%信噪比为10fl=128;%FFT长度为128Ns=6;%OFDM信号的个数为6para=128;%子载波个数为128sr=250000;%符号速率为250000br=sr.*2;%子载波的比特率gl=32;%保护时隙的长度nloop=1;%循环次数noe=0;nod=0;eop=0;nop=0;for iii=1:nloop%发送和产生数据Signal=rand（1,para*Ns*2）0.5;%串并转换paradata=reshape（Signal,para,Ns*2）;%QPSK调制ich,qch=qpskmod（paradata,pa

19、ra,Ns,2）;kmod=1./sqrt（2）;ich1=ich.*kmod;qch1=qch.*kmod;qpsk_x=ich1+qch1.*sqrt（-1）;%IFFTfy=ifft（qpsk_x）;ich2=real（fy）;qch2=imag（fy）;%保护间隔设定%衰减的计算spow=sum（ich3.2+qch3.2）/Ns./para;attn=0.5*spow*sr/br*10.（-SNR/10）;attn=sqrt（attn）;%并串转换ich4=reshape（ich3,1,（fl+gl）*Ns）;qch4=reshape（qch3,1,（fl+gl）*Ns）;%通过形成

20、复数发射数据TrData=ich4+qch4.*sqrt（-1）;%在此加入一个高斯白噪声ReData=awgn（TrData,SNR,idata=real（ReData）;qdata=imag（ReData）;idata1=reshape（idata,fl+gl,Ns）;qdata1=reshape（qdata,fl+gl,Ns）;idata2=idata1（gl+1:gl+fl,:qdata2=qdata1（gl+1:%FFTRex=idata2+qdata2.*sqrt（-1）;ry=fft（Rex）;ReIChan=real（ry）;ReQChan=imag（ry）;ReIchan=R

21、eIChan/kmod;ReQchan=ReQChan/kmod;%QPSK解调RePara=qpskdemod（ReIchan,ReQchan,para,Ns,2）;ReSig=reshape（RePara,1,para*Ns*2）;%误码比特率和及时的的误码和数据noe2=sum（abs（ReSig-Signal）;nod2=length（Signal）;noe=noe+noe2;nod=nod+nod2;if noe2=0 eop=eop+1;else eop=eop;eop;nop=nop+1;fprintf（%dt%et%dn,iii,noe2/nod2,eop）%输出结果per=e

22、op/nop;ber=noe/nod;figure（1）subplot（2,1,1）,stem（Signal）,grid minor;title（signalxlabel（x）,ylabel（ysubplot（2,1,2）,stem（ReSig）,grid minor;Resignal）figure（2）subplot（2,1,1）,stem（ich2）,grid minor;subplot（2,1,2）,stem（qch2）,grid minor;figure（3）subplot（2,1,1）,stem（idata）,grid minor;时间幅度加噪声后I路波形subplot（2,1,2）

23、,stem（qdata）,grid minor;加噪声后Q路波形figure（4）subplot（2,1,1）,stem（ReIChan）,grid minor;subplot（2,1,2）,stem（ReQChan）,grid minor;%QPSK星座图figure（5）for alfa=0:0.001*pi:2*pi plot（cos（alfa）,sin（alfa）,b hold ongrid;I路ylabel（Q路调制后的星座图hold off4. 仿真实验结果分析及总结4.1 实验分析通过图3.4可以看出，通过将图中发送数据和接收到的数据进行对比发现，经过OFDM系统的传输后，信号的

24、误码率为0，对抗码间干扰和时延扩展有很好的效果。实际OFDM系统中，子载波的数目较大时，系统的误码率也是非常低的。为了能够很好的观察到QPSK的调制，同时也得到了调制后的星座图。图3.4 发送短信号和接收短信号图3.5 进行IFFT后的I路和Q路波形图3.6 加入噪声后的I路和Q路波形图3.7 进行FFT变换后的I路和Q路波形图3.8 经过调制后的星座图通过仿真可以看出，OFDM在高斯信道具有比较良好的性能，信噪比在比较大的时候，误码率比较低。图3.4为发送信号和接收信号的波形图，从中可以看出，通过两者之间的对比，OFDM系统在对抗码间干扰中起到了很好的效果。同时通过图3.5和图3.6以及图3

25、.7进行观察I路和Q路的波形，并通过星座图作为参考，对进行IFFT变换、加入噪声、FFT变换过程中，OFDM系统对抗干扰的性能都有所体现。本次仿真利用了MATLAB强大的符号计算功能，由于在数值计算的同时需要获取问题的解析解，因此就需要涉及到符号运算，而MATLAB软件提供了便捷的符号运算，在利用MATLAB软件的过程中，关于矩阵的计算、多项式的运算以及数据统计分析等功能都起到了重要的作用。有优点的同时也存在很多不足，主要体现在一下几个方面：（1）编写代码过程比较复杂；（2）在仿真观察波形的时候修改参数比较复杂；（3）可视化程度不高。因此针对这些问题，下面通过SystemView进行一个更便捷的系统仿真，同时再次验证OFDM系统的抗干扰性能。4.2 实验总结主要工作有：（1）学习并总结OFDM的原理及其技术应用。（2）通过MATLAB平台编写m程序对QPSK的调制的OFDM进行仿真并对其性能方面的分析。体现了MATLAB软件强大的符号运算和数值方面的运算功能，借助了MATALAB运算的灵活性等特点。（3） 整个过程中主要基于操作方面的相对较多，锻炼了动手能力和对MATLAB软件的使用和了解，以及其他仿真的平台的使用，同时操作的过程中了解和掌握了OFDM的基本原理和相关应用。本文所做工作的