移动通信系统OFDM系统仿真与实现(基于MATLAB)

1、OFDM系统仿真与实现1. OFDM的应用意义在近几年以内，无线通信技术正在以前所未有的速度向前发展。由于用户对各种实时多媒体业务需求的增加和互联网技术的迅猛发展，未来的无线通信及技术将会有更高的信息传输速率，为用户提供更大的便利，其网络结构也将发生根本的变化。随着人们对通信数据化、个人化和移动化的需求，OFDM技术在无线接入领域得到了广泛的应用。OFDM是一种特殊的多载波传输方案，它将数字调制、数字信号处理、多载波传输技术结合在一起，是目前已知的频谱利用率最高的一种通信系统，具有传输速率快、抗多径干扰能力强的优点。目前，OFDM技术在数字音频广播(DAB)、地面数字视频广播(DVB-T)、无

2、线局域网等领域得到广泛应用。它将是4G移动通信的核心技术之一。OFDM广泛用于各种数字传输和通信中，如移动无线FM信道，高比特率数字用户线系统(HDSL)，不对称数字用户线系统(ADSL)，甚高比特率数字用户线系统HDSL，数字音频广播(DAB)系统，数字视频广播(DVB)和HDTV地面传播系统。1999年，IEEE802.11a通过了一个SGHz的无线局域网标准，其中OFDM调制技术被采用为物理层标准，使得传输速率可以达54MbPs。这样，可提供25MbPs的无线ATM接口和10MbPs的以太网无线帧结构接口，并支持语音、数据、图像业务。这样的速率完全能满足室内、室外的各种应用场合。OFDM

3、由于技术的成熟性，被选用为下行标准很快就达成了共识。而在上行技术的选择上，由于OFDM的高峰均比(PAPR)使得一些设备商认为会增加终端的功放成本和功率消耗，限制终端的使用时间，一些则认为可以通过滤波，削峰等方法限制峰均比。不过，经过讨论后，最后上行还是采用了SC-FDMA方式。拥有我国自主知识产权的3G标准一一TD-SCDMA在LTE演进计划中也提出了TD-CDM-OFDM的方案B3G/4G是ITU提出的目标，并希望在2010年予以实现。B3G/4G的目标是在高速移动环境下支持高达100Mb/S的下行数据传输速率，在室内和静止环境下支持高达IGb/S的下行数据传输速率。而OFDM技术也将扮演

4、重要的角色。2. OFDM的原理研究与分析 2.1 OFDM的关键技术(1) 时域和频域同步OFDM系统对定时和频率偏移敏感，特别是实际应用中与FDMA、TDMA和CDMA等多址方式结合使用时，时域和频率同步显得尤为重要。(2) 信道估计在OFDM系统中，信道估计器的设计主要有两个问题：一是导频信息的选择，由于信道常常是衰落信道，需要不断对信道进行跟踪，因此导频信息也必须不断的发送；二是复杂度较低和导频跟踪能力良好的信道估计器的设计。(3) 信道编码和交织为了提高数字通信系统的性能，信道编码和交织是普遍采用的方法。对于衰落信道衰落中的随机错误，可以采用信道编码；对于衰落信道中的突发错误，可以采

5、用交织技术。(4) 降低峰值平均功率比由于OFDM信号在时域上表现为N个正交子载波信号的叠加，当这N个信号恰好均以峰值相加时，OFDM信号也将产生最大峰值，该峰值功率是平均功率的N倍。尽管峰值功率出现的概率较低，但为了不失真地传输这些高PAPR的OFDM信号，从而导致发送效率极低，接收端对前端放大器以及A/D变换器的线性度要求也很高。因此提出了基于信号畸变技术、信号扰码技术和基于信号空间扩展等降低OFDM系统PAPR的方法。(5) 自适应技术在OFDM系统中使用自适应技术，还应考虑频率分组、时间间隔、信道总延迟和信道估计误差等因素，其中信道估计误差对性能的影响较大。2.2 OFDM的优缺点OF

6、DM主要有下列一些优点:(1) OFDM在对抗干扰及衰落的优势：把高速数据流通过串并变换，使得每个子载波上的数据符号持续长度相对增加，可有效对抗信号波形间的干扰ISI，适用于多径信道存在频率选择性，而所有子载波都处于深衰落的概率极小，OFDM系统可通过动态比特分配和动态子信道分配的方法，充分利用信噪比较高的子信道，提高系统性能。(2) OFDM系统由于子载波之间存在正交性，允许子信道的频谱相互重叠，因此与常规的频分复用相比，OFDM可以最大限度地利用频谱资源。(3) 各子信道的正交调制和解调可通过离散傅利叶反变换(IDFT)和离散傅利叶变换(DFT)实现。对N很大(N>32)的系统，可以

7、通过快速傅立叶变换(FFT)来实现。基于FFT的系统在计算方面更有效，并且随着大规模集成电路技术和DSP的发展，IFFT和EFT都非常容易实现。(4) OFDM的开放灵活性：无论从无线数据业务的使用需求,还是从移动通信系统自身要求，都希望物理层支持非对称高速数据传输,而OFDM系统可以很容易使用不同数量的子载波来实现上下链路中不同的传输速率。OFDM技术的不足: 由于OFDM系统内存在多个正交子载波，而且其输出信号是多个子信道的叠加，因此对子信道的正交性有严格要求。而由于无线信道的时变性，还有发射机载波和本地振荡器的频率偏差，所以OFDM易受频率偏差的影响。如果多个子信号的相位一致时，所得到的

8、叠加信号的瞬时功率会远大于信号的功率，出现较大的峰值与均值功率比(PAR)，这个比值的增大会降低射频放大器的功率效率，使系统性能恶化。本文主要研究了OFDM调制解调系统，对噪声信道抗干扰能力的研究以及在实际中的基本应用。利用Matlab进行仿真和性能分析。2.3 原理及数学描述2.3.1 OFDM的载波调制正交频分复用OFDM（OrthogonalFrequencyDivisionMultiplex)是一种多载波调制方式，通过减小和消除码间串扰的影响来克服信道的频率选择性衰落。它的基本原理是将信号分割为N个子信号，然后用N个子信号分别调制N个相互正交的子载波。由于子载波的频谱相互重叠，因而可以

9、得到较高的频谱效率。下图是OFDM基带信号处理原理图。其中，（a）是发射机工作原理，（b）是接收机工作原理。图2.1 OFDM基带信号处理原理图2.3.2 OFDM的调制解调原理OFDM通过把需要发射的数据流分解为若干个并行的数据子流，这样每个数据子流在速率上就会降低很多，然后再进行相关调制，将它们调制到一组总数为N，频率之间的间隔相等，且又两两正交的子载波上。OFDM的调制可以用离散傅里叶反变化(IDFT)来实现，相应的有在解调端可以用离散傅里叶变换(DFT)来实现。而这两种傅里叶变换都有相应的快速算法，在系统的效率和相应时间上会有所提高。系统的实际应用中一般采用IFFT和FFT技术。2.3

10、.3串并转换数据传输的典型形式是串行数据流，符号被连续传输，每一个数据符号的频谱可占据整个可利用的带宽。但在并行守护据传输系统中，许多符号被同时传输，减少了那些在串行系统中出现的问题。在OFDM系统中，每个传输符号速率的大小大约在几十bps到几十Kbps之间，所以必须进行串并变换，将输入串行比特流转换为可以传输的OFDM符号。由于调制模式可以自适应调节，所以每个子载波的调制模式是可变化的，因而每个子载波可传输的比特数也是可以变化的，所以串并转换需要分配给每个子载波数据段的长度是不一样的。在接收端执行相反的过程，从各个子载波处来的数据被转换回原始的串行数据。2.3.4保护间隔和循环前缀应用OFD

11、M的一个重要原因在于它可以有效的对抗多径时延扩展。通过把输入数据流串并变换到N个并行的子信道中，使得每一个调制子载波的数据周期可以扩大为原始数据符号周期的N倍。为了最大限度的消除符号间干扰，还可以在每个OFDM符号间插入保护间隔(GI)，而且该保护间隔长度一般要大于无线信道中的最大时延扩展，这样一个符号的多径分量就不会对下一个符号造成干扰。在这段保护间隔内，可以不插入任何信号，即是一段空闲的传输时段。2.3.5 OFDM的基本参数各种OFDM参数的选择就是需要在多项要求冲突中进行折中考虑。通常来讲，如前所述，首先要确定3个参数：带宽、比特率以及保护间隔，按照惯例，保护间隔的时间长度应该为应用移

12、动环境信道的时延均方根值的24倍。一旦确定了保护间隔，则OFDM符号周期长度就可以确定。为了最大限度的减少由于插入保护间隔所带来的信噪比的损失，希望OFDM符号周期长度要远远大于保护间隔长度。但是符号周期长度又不可能任意大，否则OFDM系统中包括更多的子载波数，从而导致子载波间隔相应减少，系统的实现复杂度增加，而且还加大了系统的峰值平均功率比，同时使系统对频率偏差更加敏感。因此在实际应用中，一般选择符号周期是保护间隔长度的5倍，这样由于插入保护比特所造成的信噪比损耗只有1dB左右。2.3.6调制方式可以通过改变发射的射频信号的幅度、相位和频率来调制信号。对于OFDM系统来说，只能采用前两种调制

13、方法，而不能采用频率调制的方法，这是因为子载波是频率正交，而且携带独立的信息，调制子载波频率会破坏这些子载波的正交特性，这是频率调制不能在OFDM系统中采用的原因。2.3.7 16QAM调制16QAM 是用两路独立的正交 4ASK 信号叠加而成，4ASK 是用多电平信号去键控载波而得到的信号。它是 2ASK 体制的推广，和 2ASK 相比，这种体制的优点在于信息传输速率高。16QAM 信号采取正交相干解调的方法解调，解调器首先对收到的 16QAM 信号进行正交相干解调，一路与 cos c t 相乘，一路与 sin c t 相乘。然后经过低通滤波器，低通滤波器 LPF 滤除乘法器产生的高频分量，

14、获得有用信号，低通滤波器LPF 输出经抽样判决可恢复出电平信号。16 进制的正交振幅调制是一种振幅相位联合键控信号。16QAM 的产生有 2 种方法：（1）正交调幅法，它是有 2 路正交的四电平振幅键控信号叠加而成；（2）复合相移法：它是用 2 路独立的四相位移相键控信号叠加而成。2.3.8 QPSK调制QPSK是在2PSK（二相调相）的基础上发展而来的一种多进制相位调制。QPSK是利用载波的四种不同相位来表征传送的数字信息。在QPSK调制中，首先对输入的二进制数据按二位数字编成一组，以此构成双比特码元。其组合共有4种，即有4种不同状态。故可以用M=4种相位或相位差来表示。这里M=4，故称为四

15、相调相。同样，若采用八相调制方式，在一个码元时间内可传送3位码，其信息传送速率是二相调制方式的3倍。由此可见，采用多相调制的级数愈多，系统的传输速率愈高，但相邻载波之间的相位差愈小，接收时要区分它们的困难程度就愈大，将使误码率增加。3. 基于MATLAB的OFDM系统仿真实验3.1 MATLAB的特点与功能MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使

16、用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。3.2 MATLAB仿真系统3.2.1 仿真结构图图3.1 系统仿真原理图3.2.2 主要模块（1）串并转换 OFDM系统是多载波传输的一种特殊方式，而发送端产生的数据为串行的，因此需要进行串并转换。OFDM将高速输入的串行数据比特流转换成并行传输的低速数据流，而且串并转换之后会提高系统的抗干扰能力。串并转换的示意图如图3.2。图3.2 串并转换示意图（2）QPSK调制QPSK调

17、制的方法有两种，一种是数字方法，另一种是相位选择法。本程序中选择的是利用数字的方法来调制信号。这里通过调用一个子函数qpskmod，程序如下：function iout,qout=qpskmod(paradata,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

18、+ 2. ( m2 - ii ) .* paradata2(1:para),m2+ii+count2); end iout(1:para),jj)=isi; qout(1:para),jj)=isq; count2=count2+ml;end（3）IFFT/FFT运算通过N点的IDFT运算，把频域数据符号变换为时域数据信号，经调制后发送到信道中去。在实际应用中，一般采用FFT/IFFT运算，因为它可以显著的降低运算的复杂度。在Matlab软件里可以使用函数fft()和ifft()来对数据进行FFT/IFFT运算，可以省去很多复杂的运算。（4）保护间隔和循环前缀 在OFDM系统中，保护间隔是一种

19、循环复制增加了符号的波形长度，在符号的数据部分，每一个子载波内有一个整数倍的循环，此种符号的复制产生了一个循环的信号，即将每个OFDM符号的后一段时间的样点复制到OFDM符号的前面，形成前缀，在交接点没有任何的间断。当信道的最大多径时延扩展小于保护间隔的时间时，OFDM系统可以完全克服ISI的影响。具体的m语言实现为： ich3=ich2(fl-gl+1:fl,:);ich2; qch3=qch2(fl-gl+1:fl,:);qch2; （5）并串转换 这一过程是串并转换的逆过程，将N个子载波的数据传送到一个载波信道中去，将并行数据转换为串行数据序列进行传输。示意图如图3.3。图3.3 并串转

20、换示意图（6）加入高斯噪声 白噪声是根据噪声的功率谱密度是否均匀来定义的，而高斯噪声则是根据它的概率密度函数呈正态分布来定义的，高斯型白噪声称高斯白噪声，在通信系统的理论分析中，特别是在分析、计算系统抗噪声性能时，经常假定系统中信道噪声（即前述的起伏噪声）为高斯型白噪声。在matlab软件里产生高斯噪声的函数一般有两个WGN和AWGN。Awgn的用法为： y = awgn(x,SNR) 在信号x中加入高斯白噪声。信噪比SNR以dB为单位。x的强度假定为0dBW。如果x是复数，就加入复噪声。y = awgn(x,SNR,SIGPOWER) 如果SIGPOWER是数值，则其代表以dBW为单位的信号

21、强度；如果SIGPOWER为''measured''，则函数将在加入噪声之前测定信号强度。本程序正是采用的ReData=awgn（TrData，SNR，measured）来给复数发射数据TrData加入高斯噪声。（7）QPSK解调在进行解调之前，信号要完成去掉保护间隔和进行FFT运算的任务。去掉保护间隔也就可以去掉符号间的干扰。这里通过调用一个子函数qpskdemod，具体的m语言实现为： function demodata=qpskdemod(idata,qdata,para,nd,ml)demodata=zeros( para,ml*nd );demodat

22、a(1:para),(1:ml:ml*nd-1)=idata(1:para),(1:nd)>=0;demodata(1:para),(2:ml:ml*nd)=qdata(1:para),(1:nd)>=0;（8）接收信号 完成解调之后下面就是进行解调信号的判决得到接收信号，程序中将发送段的原始信号和接收端的信号的波形输了出来。（9）MATLAB仿真源程序：clear;clc;SNR=10;%信噪比为10fl=128;%FFT长度为128Ns=6;%OFDM信号的个数为6para=128;%子载波个数为128sr=250000;%符号速率为250000br=sr.*2;%子载波的比特

23、率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,para,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);qc

24、h2=imag(fy);%保护间隔设定ich3=ich2(fl-gl+1:fl,:);ich2; qch3=qch2(fl-gl+1:fl,:);qch2;%衰减的计算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);%通过形成复数发射数据TrData=ich4+qch4.*sqrt(-1);%在此加入一个高斯白噪声ReData=awgn(TrD

25、ata,SNR,'measured');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:gl+fl,:);%FFTRex=idata2+qdata2.*sqrt(-1);ry=fft(Rex);ReIChan=real(ry);ReQChan=imag(ry);ReIchan=ReIChan/kmod;ReQchan=ReQChan/k

26、mod;%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;endeop;nop=nop+1;fprintf('%dt%et%dn',iii,noe2/nod2,eop)end%输出结果per=eop/nop;ber=no

27、e/nod;figure(1)subplot(2,1,1),stem(Signal),grid minor;title('signal');xlabel('x'),ylabel('y');subplot(2,1,2),stem(ReSig),grid minor;title('Resignal')figure(2)subplot(2,1,1),stem(ich2),grid minor;xlabel('x'),ylabel('y');subplot(2,1,2),stem(qch2),grid mi

28、nor;xlabel('x'),ylabel('y');figure(3)subplot(2,1,1),stem(idata),grid minor;xlabel('时间'),ylabel('幅度');title('加噪声后I路波形')subplot(2,1,2),stem(qdata),grid minor;xlabel('时间'),ylabel('幅度');title('加噪声后Q路波形')figure(4)subplot(2,1,1),stem(ReIChan)

29、,grid minor;subplot(2,1,2),stem(ReQChan),grid minor;xlabel('x'),ylabel('y');%QPSK星座图figure(5)for alfa=0:0.001*pi:2*pi plot(cos(alfa),sin(alfa),'b') hold onendgrid;xlabel('I路');ylabel('Q路');title('调制后的星座图');hold off4. 仿真实验结果分析及总结4.1 实验分析通过图3.4可以看出，通过将图中发送数据和接收到的数据进行对比发现，经过O