用MATLAB实现OFDM仿真分析

1、3.1计算机仿真仿真实验是掌握系统性能的一种手段。它通过对仿真模型的实验结果来确定实际系统 的性能。从而为新系统的建立或系统的改进提供可靠的参考。通过仿真，可以降低新系统 失败的可能性，消除系统中潜在的瓶颈。优化系统的整体性能，衡量方案的可行性。从中 选择最后合理的系统配置和参数配置。然后再应用于实际系统中。因此，仿真是科学研究 和工程建设中不可缺少的方法。3.1.1仿真平台硬件CPU Pentium III 600MHz内存：128M SDRAM软件操作系统：Microsoft Win dows2000 版本 5.0仿真软件：The Math Works Inc. Matlab 版本6.5包

2、括MATLAB 6.5的 M文件仿真系统。Matlab是一种强大的工程计算软件。目前最新的 6.x版本（windows环境）是一种功能 强、效率高、便于进行科学和工程计算的交互式软件包。其工具箱中包括：数值分析、矩阵运算、通信、数字信号处理、建模和系统控制等应用工具程序，并集应用程序和图形于 一便于使用的集成环境中。在此环境下所解问题的Matlab语言表述形式和其数学表达形式相同，不需要按传统的方法编程。Matlab的特点是编程效率高，用户使用方便，扩充能力强，语句简单，内涵丰富，高效方便的矩阵和数组运算，方便的绘图功能。3.1.2基于MATLAB勺OFDMS统仿真链路根据OFDM基本原理，本

3、文给出利用MATLA编写OFD系统的仿真链路流程。串行数据 经串并变换后进行QDPS数字调制，调制后的复信号通过 N点IFFT变换，完成多载波调制， 使信号能够在N个子载波上并行传输，中间插入10训练序列符号用于信道估计，加入循环 前缀后经并串转换、D /A后进入信道,接收端经过N点FFT变换后进行信道估计，将QDPS解 调后的数据并串变换后得到原始信息比特。本文采用MATLA语言编写M文件来实现上述系统。M文件包括脚本M文件和函数M文件， M文件的强大功能为MATLA的可扩展性提供了基础和保障，使MATLA能不断完善和壮大，成 为一个开放的、功能强大的实用工具。M文件通过in put命令可以

4、轻松实现用户和程序的交 互，通过循环向量化、数组维数预定义等提高 M文件执行速度，优化内存管理，此外，还 可以通过类似C+H语言的面向对象编程方法等等。图3-1仿真链路模型3.2循环前缀及信道估计对系统误码率的改善分析本文主要研究循环前缀对系统误码率的改善分析3.2.1循环前缀OFDM系统中，每个并行数据支路都是窄带信号，可近似认为每个支路都经历平坦衰落，这样就减小了频率选择性衰落对信号的影响。同时，每路子数据流速率的降低，减小了符 号间干扰(ISI)。此外，还可以通过加保护间隔的办法完全消除符号间干扰。假设每个OFDM符号由丫个样值组成，由于时延扩展，接收端将会有和信道冲激响应持续时间相对应

5、的前 L(L L就可完全消除ISI。最初的保护间隔是用空数据填充的，这虽然消除了 ISI，但却破坏了信道间的正交性。 后来,Peled和Ruiz提出了用循环前缀填充保护间隔的方法，即把丫个样值的最后M个复制到个OFDM符号的前端作为保护间隔，利用循环卷积的概念，只要循环前缀的长度大 于信道的冲激响应，信道间仍是正交的。符号周期由T增加至T = T + T， T是保护时隙，增加保护时隙会降低频谱利用率,所以T 一般小于等于T/4。3.2.2系统仿真结果分析用上述OFD系统的实现方法进行OFD传输系统的仿真,系统参数见表3-1并行子信道个数64FFT的长度64载波个数64符号速率25000循环次数

6、100调制方式QPSK循环一次OFDM符号个数10首先基本参数设置成如上表所示，再分别设置保护间隔为18和32，然后改变输入信号的信噪比为：3、4、5、6、7、8、9、10、11、12、13。para=iE4图3-2相同信道（64）插入不同保护间隔下误码率比较图3-3相同信道（128）插入不同保护间隔下误码率比较图3-4相同信道（256）插入不同保护间隔下误码率比较如图3-2可知，误码率随着信噪比的增加而逐渐减小，直至趋于理想化的误码率为0。paiA=2,5Edatal （下方曲线）是插入保护间隔为18时的误码率曲线，data2 （上方曲线）是插入保护 间隔为32时的曲线，可见在插入保护间隔过

7、大时，会造成在信噪比低时的高误码率，但随 着信噪比的增加，误码率会趋于相同。改变并行子信道的个数为200,如图3-3，3-4，做相同的比较，可得同样的结论。如图3-5与图3-6可知，在插入相同保护间隔的情况下，并行子信道的个数越多，信道传输信息的误码率越低。并且大大提高了频谱的利用率。这就充分说明OFDM系统的优点，可以把串行高速数据流成若干路并行低速数据流，每路低速数据流被调制在彼此正交 的子载波上，然后所有子载波叠加在一起构成发送信号。这种技术具有在杂波干扰下传送 信号的能力，可以利用在容易受外界干扰或者抵抗外界干扰能力较差的传输介质中。因此,OFDM的应用前景十分光明。图3-5相同保护间

8、隔下插入不同信道（64、128、256）误码率比较图3-6 相同保护间隔（32）下插入不同信道（128、200）误码率比较OFDM技术在实现的过程中，需要根据相应的信道条件和系统要求进行合理设计，才能 发挥其优势。系统的参数选择，导频和同步方案的设计，均衡和编码技术的结合使用，都 需要在实现之前进行优化设计，所以仿真工作是必不可少的。OFD啲系统仿真是一个非常 复杂的问题。根据研究问题的侧重点，可以采用不同的仿真工具和方法分开进行。目前许 多系统级的仿真工具都支持 OFDM系统的仿真，对于具体的算法则可以使用更加灵活的程 序语言工具进行仿真。OFDM系统的仿真是可以优化整个系统的参数和指标，缩

9、短开发的周 期。结合我们的实践经验，系统的分析了OFDM实现中的关键技术，通过一个实例给出了OFDM系统仿真的基本框架，但在具体的系统设计中，还有很多更复杂的问题需要解决，尤 其是同步技术，本文给出的只是系统设计时需要宏观考虑的问题。第4章结论随着数字时代的到来，OFD技术给通信领域注入了新的活力。正如绪论中所述，OFDM技术良好的性能在计算机网络、数字音频/视频广播、移动通信、卫星通信等多个领域中 都有体现。本文首先介绍了 OFD系统的关键技术，然后又分析了 OFD系统的原理，在OFD系统原 理的基础上，详细设计了一个基于 MATLA的系统仿真链路模型，并进行了计算机仿真，同 时详细分析了仿

10、真结果。通过对该系统的仿真及分析，可以得到如下结论：a) IFFT/FFT变换可以有效的实现OFD信号的调制解调。b) OFD可以有效的提供系统频带利用率。c)信息流的误码率随着信噪比的增加而逐渐减小。d)信息传输时，并行信道的个数越多信息流的误码率越小。本文所做的只是系统传输误码率的仿真，在实际应用中，大多需要将OFD同相和反相两路信号调制成一路信号后，进行射频调制。这时，系统的同步问题将是仿真的重点。需 要充分考虑系统的频率同步问题，频率的偏移将导致系统子载波得正交性变差，严重影响 系统的正常工作。OFDM同步也是当前通信领域里的一个研究热点。事物是在不断发展中，OFDM技术也需要发展。通

11、过与新技术的结合，OFD技术将不 断完善。其良好的性能预示着它将拥有一个广阔的前景。致 谢在论文的结束，特附上此致谢，感谢在写本文期间给予了我热心帮助的老师、同学和 朋友：感谢在本科学习期间给予我热心帮助的老师、同学和朋友。首先要感谢指导老师：常怡萍，感谢这一个多学期来常老师在学习、工作中给予的关 心和帮助。本文的研究工作也是在常老师的悉心指导下完成的。常老师严谨求实的工作作 风都使我受益匪浅。尤其是她对当今科技前沿的把握，使我所做的工作具有非常广泛的应 用前景。在她耐心细致的指导下，我的专业知识得以扩充，学习能力得以提高，论文的工 作也得以顺利地进行。所以，再次向常老师表达我最诚挚的谢意。在

12、这两年的学习生活中，还得到了学校多位老师的热情帮助。他们经常给我讲解当 今世界前沿技术，给了我许多启发。在此一并感谢。最后还要感谢我的同学以及本科学习 期间的舍友们，和他们一起度过两年和谐融洽的宝贵时间，和他们建立的友谊和我所学到 的知识同样弥笃可贵。这两年来的每一点成绩和进步都得益于我身边的很多人，在此谨向所有关心过我， 帮助过我的老师和朋友表示衷心的感谢。电讯技术1998.06期移动通信1997.06 期1998.04 期1998.04 期1998.05 期1998.02 期参考文献佟学俭、罗涛.OFDM移动通信技术原理与应用.第一版.人民邮电出版社.2003樊昌信，张甫栩，徐炳祥，吴成柯

13、.通信原理.第五版.国防工业出版社.2001郑君里，应启衍，杨为理信号与系统第二版北京高等教育出版社 2000陈亚勇等编著MATLAB信号处理详解人民邮电出版.2001.李引新，赵姚同正交频分多路（OFDM的实现方式及性能分析张永生，李道本正交频分复用码分多址技术的特点及其应用梁兵，董在望OFDM仿真软件系统研究广播与电视技术郭坚正交频分复用技术极其在 CDMA中的应用数字通信高焕英 OFDM调制技术研究及硬件实现方案无线电技术张俊峰 正交频分复用（OFDM）调制技术天津通信技术王文博，郑侃宽带无线通信OFDM技术人民邮电出版社尹长川，罗涛，乐光新多载波宽带无线通信技术北京邮电大学出版社123

14、4789101112131415张云峰，彭军用MATLA取现OFDM技术的仿真电脑与信息技术14卷1期 2006.2阮沈勇王永利桑群芳编MATLAB程序设计电子工业出版社2004.1（ MATLAB应用与提高系列）李建东、杨家玮，个人通信，人民邮电出版社，1998.Q%*准备部分*para=128;fftle n=128;noc=128;n d=6;ml=2;sr=2500;br=sr.*ml;gile n=18;ebn 0=3;o%*主循环部分*nloop=100;noe = 0;nod = 0; eop=0;no p=0;for iii=1: nloopQ%k*发射部分Q%*产生数据*附M

15、ATLA仿真程序seldata=ra nd(1,para* nd*ml)0.5;Q%k*串并转*paradata=reshape(seldata,para ,n d*ml);*Q%*调制*Q%*离散傅里叶反变换）*Q%*插入保护间隔*Q%*接收部分ich,qch=qpskmod(paradata,para, nd,ml); kmod=1/sqrt(2);ich1=ich.*kmod;qch1=qch.*kmod;x=ich1+qch1.*i; y=ifft(x); ich2=real(y); qch2=imag(y);ich3,qch3= gii ns(ich2,qch2,fftle n, g

16、ile n,n d); fftle n2=fftle n+gile n;Q%k*衰减计算 *spow=sum(ich3.A2+qch3.A2)/nd./para;att n=0.5*spow*sr/br*10.A(-eb nO/1O);att n=sqrt(att n);%*高斯白信道*ich4,qch4=comb(ich3,qch3,att n);*Q%k*去除保护间隔 *中原工学院毕业设计论文ich5,qch5= girem(ich4,qch4,fftle n2,gile n,n d);(%* FFT(离散傅里叶变换）*rx=ich5+qch5.*i; ry=fft(rx); ich6=r

17、eal(ry); qch6=imag(ry);(%* QPSK角年调 *ich7=ich6./kmod;qch7=qch6./kmod;demodata=qpskdemod(ich7,qch7,para, nd,ml);%*并串变换 *demodata仁reshape(demodata,1,para* nd*ml);%*比特误码率 (BER) *% in sta ntan eous nu mber of error and data no e2=sum(abs(demodata1-seldata);no d2=le ngth(seldata);% cumulative the nu mber of error and data in noe and nod