大学毕业论文---基于matlab的扩频通信系统仿真.doc

毕业论文基于MATLAB的扩频通信系统仿真基于MATLAB的扩频通信系统仿真摘 要本文阐述了扩展频谱通信技术的理论基础和实现方法，利用MATLAB提供的可视化工具SIMULINK建立了扩频通信系统仿真模型，详细讲述了各模块的设计，并指出了仿真建模中要注意的问题。在给定仿真条件下，运行了仿真程序，得到了预期的仿真结果。同时，利用建立的仿真系统，研究了扩频增益与输出端信噪比的关系，结果表明，在相同误码率下，增大扩频增益，可以提高系统输出端的信噪比，从而提高通信系统的抗干扰能力。关键词扩频通信,信噪比， 误码率 Simulation of the Spread Spectrum Communication System Based on MATLABAuthor: Guo JianTutor: Tang HailingAbstractThe theory base and realizing methods of the spread spectrum communication technology was presented in this study. The simulation model of the spread spectrum communication system was built by using SIMULINK, which is provided by MATLAB. In addition, each module of the simulation model was introduced in detail，and pointed out the problems that must be pay attention to in the system simulation. On the basis of the designed simulation conditions, the simulation program was run and the anticipant results were gained. Moreover, the relationship between the spread spectrum gain and the fan-out error rate was also studied by use of the simulation system. The results showed that on the base of the same error rate, if the spread spectrum gain was enlarged, the Signal-to-Noise of the system fan-out would be enhanced and the anti-jamming capability of the communication system would also be enhanced.Keywords: spread spectrum communication, Signal-to-Noise, error rate目录1.绪论11.1 扩展频谱简介11.2 扩展频谱技术特点11.3 研究扩频通信的目的和意义21.4 本文的主要内容22.扩展频谱技术32.1 理论基础32.2 频谱的扩展实现42.3 m序列52.5 本章小结83.SIMULINK模型建立93.1 随机整数发生器（Random Integer genarator）93.2 PN序列发生器（PN Sequence Generator）103.3 最小相移键控调制（M-PSK Modulator Passband）113.4 加性高斯白噪声信道（AWGN Channel）113.5 最小相移键控解调 （M-PSK Demodulator Passband）123.6 误码率分析仪（Error Rate Calculation）123.7 其他设置133.7.1 RELAY 的设置133.7.2 product的设置143.7.3 display 的设置143.7.4 频谱分析仪（spectrum scope）的设置144.仿真结果与分析164.1 Simulink仿真结果164.2 Simulink仿真结果分析184.3 m文件的仿真结果194.4 m文件的仿真结果分析205.结论21致谢22参考文献23附录241.绪论1.1 扩展频谱简介扩展频谱通信具有很强的抗干扰性能，其多址能力、保密、抗多径等功能也倍受人们的关注，被广泛地应用于军事通信和民用通信中。扩频通信系统利用了扩展频谱技术，将信号扩展到很宽的频带上，在接收端对扩频信号进行相关处理即带宽压缩，恢复成窄带信号。对干扰信号而言，由于与扩频信号不相关，则被扩展到一个很宽的频带上，使之进入信号通频带内的干扰功率大大降低，相应增加了相关器输出端的信号/干扰比，对大多数人为干扰而言，扩频通信系统都具有很强的对抗能力。本文利用MATLAB对扩频系统中的m序列的产生、频谱、相关函数，以及整个扩频系统工作原理进行了仿真，为今后扩频通信系统在各个领域的应用和研究提供了依据。1.2 扩展频谱技术特点由于扩频通信能大大扩展信号的频谱，发端用扩频码序列进行扩频调制，以及在收端用相关解调技术，使其具有许多窄带通信难于替代的优良性能，能在“军转民”后，迅速推广到各种公用和专用通信网络之中，主要有以下几项特点：1 易于重复使用频率，提高了无线频谱利用率 无线频谱十分宝贵，虽然从长波到微波都得到了开发利用，仍然满足不了社会的需求。在窄带通信中，主要依靠波道划分来防止信道之间发生干扰。 为此，世界各国都设立了频率管理机构，用户只能使用申请获准的频率。 2 抗干扰性强，误码率低 扩频通信在空间传输时所占有的带宽相对较宽，而收端又采用相关检测的办法来解扩，使有用宽带信息信号恢复成窄带信号，而把非所需信号扩展成宽带信号，然后通过窄带滤波技术提取有用的信号。这祥，对于各种干扰信号，因其在收端的非相关性，解扩后窄带信号中只有很微弱的成份，信噪比很高，因此抗干扰性强。在目前商用的通信系统中，扩频通信是唯一能够工作于负信噪比条件下的通信方式。1.3 研究扩频通信的目的和意义扩频通信是通信的一个重要分支和信道通信系统的发展方向。采用扩频信号进行通信的优越性在于用扩展频谱的方法可以换取信噪比的好处，即接收机输出的信噪比相对于输入的信噪比有很大改善，从而提高了系统的抗干扰能力。扩频技术还具有保密性好、易于实现多址通信等优点，因此该技术越来越受到人们的重视。近年来，随着超大规模集成电路技术、微处理器技术的飞速发展，以及一些新型元器件的应用，扩频通信在技术上已迈上了一个新的台阶，不仅在军事通信中占有重要地位，而且正迅速地渗透到了个人通信和计算机通信等民用领域，成为新世纪最有潜力的通信技术之一。因此研究扩频通信具有很深远的意义。1.4 本文的主要内容本文第一章对扩展频谱进行了简介，并对扩展频谱技术的特点，研究扩频通信的目的和意义进行了阐述。第二章对扩展频谱技术的理论基础，扩展频谱的理论来系统的介绍扩展频谱技术，并在此简单介绍了m序列移位寄存器。第三章主要内容为直序扩频通信系统，通过对直序扩频通信的简化框图以及扩展和解扩过程图来介绍直序扩频，并介绍了直序扩频的特点及其应用，最后利用程序得出直序扩频通信系统的仿真结果。第四章通过扩频通信的应用及未来发展方向来得出结论，并附上参考文献，方便查找相关资料。最后附录为直接序列扩展频谱仿真程序，以便更好的验证直接序列扩展频谱。2.扩展频谱技术2.1 理论基础扩频通信的基本特点是其传输信息所用信号的带宽远大于信息本身的带宽。除此以外，扩频通信还具有如下特征:1是一种数字传输方式；2带宽的展宽是利用与被传信息无关的函数(扩频函数)对被传信息进行调制实现的；3在接收端使用相同的扩频函数对扩频信号进行相关解调，还原出被传信息。Shannon定理指出：在高斯白噪声干扰条件下，通信系统的极限传输速率（或称信道容量）为 （2-1）式中：C为系统的信道容量（bit/s）； B为系统信道带宽（Hz）；S为信号的平均功率； N为噪声功率由式中可以看出。若白噪声的功率谱密度为n0，噪声功率N=n0B ，则信道容量C可表示为： （2-2）由上式可以看出，B、n0、S确定后，信道容量C就确定了。由Shannon第二定理知，若信源的信息速率 小于或等于信道容量C，通过编码，信源的信息能以任意小的差错概率通过信道传输。为使信源产生的信息以尽可能高的信息速率通过信道，提高信道容量是人们所期望的。由Shannon公式可以看出：（1）信道容量C为常数时，带宽B与信噪比S/N可以互换，即可以通过增加带宽B来降低系统对信噪比S/N 的要求。（2）要增加系统的信息传输速率，则要求增加信道容量。增加信道容量的方法可以通过增加传输信号带宽B，或增加信噪比S/N来实现。由式（1）可知， B与C成正比，而C与S/N成对数关系，因此，增加B比增加S/N更有效。扩频通信系统由于在发送端扩展了信号频谱，在接收端解扩还原了信息，这样的系统带来的好处是大大提高了抗干扰容限。理论分析表明，各种扩频系统的抗干扰性能与信息频谱扩展后的扩频信号带宽比例有关。一般把扩频信号带宽B与信息带宽F之比称为处理增益GP，即： （2-2）B为系统信道带宽（Hz）； Gp扩频处理增益； F信息带宽；由上式可以看出，B与F差别越大，Gp越大，也就是说，扩频的增益越大。它表明了扩频系统信噪比改善的程度。除此之外，扩频系统的其他一些性能也大都与Gp有关。因此，处理增益是扩频系统的一个重要性能指标。系统的抗干扰容限Mj定义如下： （2-3）式中：S/N 输出端的信噪比； Ls 系统损耗；Gp扩频处理增益； Mj系统的抗干扰容限；由此可见，抗干扰容限Mj与扩频处理增益Gp成正比，扩频处理增益提高后，抗干扰容限大大提高，甚至信号在一定的噪声湮没下也能正常通信。通常的扩频设备总是将用户信息(待传输信息)的带宽扩展到数十倍、上百倍甚至千倍，以尽可能地提高处理增益【1】。2.2 频谱的扩展实现频谱的扩展是用数字化方式实现的。在一个二进制码位的时段内用一组新的多位长的码型予以置换，新码型的码速率远远高出原码的码速率，由傅立叶分析可知新码型的带宽远远高出原码的带宽，从而将信号的带宽进行了扩展。这些新的码型也叫伪随机（PN）码，码位越长系统性能越高。通常，商用扩频系统PN码码长应不低于12位，一般取32位，军用系统可达千位。当选取上述任意一个序列后，如M序列，将其中可用的编码，即正交码，两两组合，并划分为若干组，各组分别代表不同用户，组内两个码型分别表示原始信息"1"和"0"。系统对原始信息进行编码、传送，接收端利用相关处理器对接收信号与本地码型相关进行相关运算，解出基带信号( 即原始信息)实现解扩，从而区分出不同用户的不同信息。微波无线扩频通信的原理见图1。 信息调制扩频解扩解调扩频序列发生器扩频序列发生器同步电路B1B2B2D信息接收图1扩频系统基本原理图根据扩展频谱的方式不同，扩频通信系统可分为：直接序列扩频（DS）、跳频（FH）、跳时（TH）、线性调频以及以上几种方法的组合（混频）3。2.3 m序列Shannon在证明编码定理的时候，提出了用具有白噪声统计特性的信号来编码。白噪声是一种随机过程，它的瞬时值服从正态分布，功率谱在很宽的频带内都是均匀的，它有及其优良的相关特性。但是至今无法实现对白噪声的放大、调制、检测、同步及控制等，而只能用具有类似带限白噪声统计特性的伪随机码来逼近它，并作为扩频系统的扩频码。m序列是最长线性移位寄存器，是最重要的伪随机序列之一，这种序列易于产生，有优良的自相关特性。m序列是由移位寄存器加反馈后形成的，其结构如图2所示。图中n-1 （ i=1，2，3， r）为移位寄存器中每位寄存器的状态； ci（ i=1，2，3，r ）为第 n位寄存器的反馈系数。当n =0时，表示无反馈，将反馈线断开；当ci=1时表示存在反馈，将反馈线连起来。在此结构中c0=cr=1，c0不能为0，c0为0就不能构成周期性序列，因为 c0=0意味着无反馈，为静态移位寄存器。cr也不能为0，即第r 位寄存器一定要参加反馈，否则，r级的反馈移位寄存器将减化为r-1级的或更低的反馈移位寄存器。不同的反馈逻辑，即 ci（ i=1，2，3， r-1）取不同的值，将产生不同的移位寄存序列2。2.4 直序扩频扩频通信与一般的通信系统相比，主要是在发射端增加了扩频调制，而在接收端增加了扩频解调的过程，扩频通信按其工作方式不同主要分为直接序列扩频系统、跳频扩频系统、跳时扩频系统、线性调频系统和混合调频系统。现以直接序列扩频系统为例说明扩频通信的实现方法。直接序列扩频工作方式，简称直扩方式（DS方式）。就是用高速率的扩频序列在发射端扩展信号的频谱，而在接收端用相同的扩频码序列进行解扩，把展开的扩频信号还原成原来的信号。 直接序列扩频方式是直接用伪噪声序列对载波进行调制，要传送的数据信息需要经过信道编码后，与伪噪声序列进行模2和生成复合码去调制载波。信息调制扩频解扩解调PN码发生器PN码发生器同步电路B1B2B2D 图2直接序列扩展频谱系统简化框图由上图看出，在发射端，信源输出的信号与伪随机码产生器产生的伪随机码进行模2加，产生一速率与伪随机码速率相同的扩频序列，然后再用扩频序列去调制载波，这样得到已扩频调制的射频信号。在接收端，接收到的扩频信号经高放和混频后，用与发射端同步的伪随机序列对扩频调制信号进行相关解扩，将信号的频带恢复为信息序列的频带，然后进行解调，恢复出所传输的信息。 调制后的待传信号PN码待传信息 图3 直序扩频的频谱扩展过程直接序列扩频（DS-Direct Scquency）的频谱扩展和解扩过程见图3和图4所示。直接序列扩频就是用高码率的扩频码序列在发端直接去扩展信号的频谱，在收端直接使用相同的扩频码序列对扩展的信号频谱进行解调，还原出原始的信息。解调信息PN码解扩后的信号 图4 直接序列的频谱解扩过程在图上我们可以看出：1，在发端，信息码经码率较高的PN码调制以后，频谱被扩展了。在收端，扩频信号经同样的PN码解调以后，信息码被恢复；2，信息码经调制、扩频传输、解调然后恢复的过程，类似与PN码进行了二次"模二相加"的过程。直序扩频系统的内容十分广泛。根据需要不同，实际直序扩频系统的扩频、调制、解扩、解调等部分可以采用不用的方案4。2.5 本章小结直接序列扩频(Direct Sequence Spread Spectrum）系统是将要发送的信息用伪随机码（PN码）扩展到一个很宽的频带上去，在接收端，用与发端扩展用的相同的伪随机码对接收到的扩频信号进行相关处理，恢复出发送的信息。直序扩频系统的特点有：1 抗干扰性强 2 隐蔽性好3 易于实现码分多址（CDMA）4 抗多径干扰5 直扩通信速率高6 抗衰落7 远近"效应"8 组网能力9 窄带系统的兼容3.SIMULINK模型建立基于MATLAB/Simulink所建立的扩频通信系统的仿真模型，能够反映扩频通信系统的动态工作过程，可进行波形观察、品剖分析和性能分析等，同时能根据研究和设计的需要扩展仿真模型，实现以扩频通信为基础的现代通信的模拟仿真，为系统的研究和设计提供强有力的平台。图2为基于MATLAB/Simulink的直接序列扩频通信系统的仿真模型。图5 基于Simulink 的直接序列扩频通信系统的仿真模型3.1 随机整数发生器（Random Integer genarator）随机整数发生器(Randon Integer generator)作为仿真系统的信源，随机整数发生器产生二进制随机信号，采用时间、初始状态可自由设置，从而满足扩频通信系统所需信源的要求。3.2 PN序列发生器（PN Sequence Generator）用PN序列发生器（PN Sequence Generator）产生伪随机码，对基带信号进行扩频。扩频过程通过信息码与PN码进行双极性变换后相乘加以实现。解扩过程与扩频过程相同，即将接收的信号用PN码进行第二次扩频处理。3.3 最小相移键控调制（M-PSK Modulator Passband）对基带信号进行数字调制，将频谱从基带搬移至频带，以实现远距离传输的效果。3.4 加性高斯白噪声信道（AWGN Channel）信号通过加性高斯白噪声信道，以达到对信号的传输后进行再进行解扩解调的目的。改信道信噪比参数设置为负数，可验证扩频通信系统具有良好的抗噪声能力。3.5 最小相移键控解调 （M-PSK Demodulator Passband）对解扩后的信号进行解调，恢复基带信号。3.6 误码率分析仪（Error Rate Calculation）在误码率计算中，接收到的信号，由于经过扩频解扩、调制解调、相关统计等处理，会存在一个延迟，在误码仪模块的对话框中要设置一个合适的延迟。3.7 其他设置3.7.1 RELAY 的设置3.7.2 product的设置3.7.3 display 的设置3.7.4 频谱分析仪（spectrum scope）的设置系统中有四个频谱分析仪器，分别命名为：baseband、mod、SS和DSS。其中，baseband能观察基带信号的频谱；mod能观察到基带信号调制到频带后的频谱；SS能观察扩频后的频谱；DSS能观察解扩解调后的接受端输出信号的频谱。这四个频谱分析仪的参数设置大致一致，其前方的零阶保持器的采样时间略有不同。频谱分析仪名称BasebandModSSDSS采样时间（s）0.0011/300000.001/310.0014.仿真结果与分析4.1 Simulink仿真结果l 基带信号的频谱如下l 经调制后频谱如下l 扩频后频谱如下l 接收端恢复后的输出信号频谱如下l 误码计算仪的结果4.2 Simulink仿真结果分析可以看出，基带信号的频率分量主要是在-100Hz到100Hz，这与信号源所设置的频率相符合。基带信号进过M-PSK调制后，正频域的频率分量集中在2.9KHz-3.1KHz处，可看出载波频率为3KHz，这个结果与M-PSK调制器所设置的载波频率一致。频带信号经过与PN码相乘扩频后，观察其信号的频谱可以看出，正频域的频率分量主要集中在0Hz到20KHz这主瓣中，带宽扩展了约20倍。此时信号将经过模拟信道的AWGN信道模块，来到接收端。在接收端处，先经过解扩（即与扩频PN码同步的PN码相乘），得到解扩后的频带信号，再经过解调模块，恢复成基带信号。观察恢复后的信号的频谱，与信源处的信号基本一致。再看系统的误码率计算分析。 输出的数据是一个n 行（与输入数据数目相等）3列的矩阵。 第1列是差错率， 第2列是差错码的数量， 第3列是码元总数即前述的n。 可看出，差错率为0.002，差错码的数量为2，码元总数为1000。由此可见，系统具有较好的抗噪声性能。企抗噪声原理如图3所示。图6 扩频系统抗噪声原理示意图4.3 m文件的仿真结果仿真结果如图7和图8所示。其中，图4是信号比和误码率的值；图5 是绘制的函数关系图像。 图7 仿真结果的值图8 信噪比和误码率的关系函数图象4.4 m文件的仿真结果分析通过观察可知，信噪比越大，系统的误码率越小。这符合通信原理的规律。从图中还可得知，在负的信噪比条件下，扩频通信系统还能有比较小的误码率。可知扩频系统的抗干扰性能很好，能在负信噪比条件下工作。这是因为扩频系统特有的频谱扩展特点和特有的解扩技术，能使有用信号淹没在伪噪声之中，通过相关解扩将之还原。5.结论扩频通信技术具有广泛的应用：在军事上，在现代科技的许多领域中，扩频技术得到了非常广泛的应用。如空间系统通信、航空及航海电子系统等等。在民用通信中：90年代以来，扩频通信开始向民用通信领域发展。一是随着国际形势的变化，原来从事军品开发的部门与企业都纷纷转入民用开发，以寻找出路；二是数字蜂窝移动通信、个人通信等新兴通信方式，要求采用能节约频带的技术，解决频带拥挤问题，使人们考虑到扩频通信可以与现有通信并存，是提高频带利用率的有效途径；三是市场需求的推动，采用新技术是为了占有国内市场，竞争国际市场。扩频通信的未来发展方向：随着扩频通信的迅速发展，特别是2.4GHz频段扩频系统的大量使用，扩频系统之间的干扰问题越来越突出总的来说，扩频通信是当代通信技术的重大突破，是当今高新科技的热点之一，发展变化非常快，有时真是：山穷水复疑无路，柳暗花明又一春。大概这也是扩频通信的魅力所在【6】。扩频通信以其较强的抗干扰、抗衰落、抗多径性能而成为第三代通信的核心技术，结果表明，扩频系统对正弦干扰有着良好的抗干扰性，增大信噪比可以有效抑制正弦信号的干扰。可以广泛应用于对抗干扰性和保密性要求较高的军用或民用通信。本人通过课程设计，进行深入地研究学习扩频通信技术及对它进行仿真应用，将所学的知识进行归纳与总结，从而巩固通信专业基础知识，为以后的个人学习和工作打下基础。致谢非常感谢唐海玲老师，唐海玲老师在我大学的最后学习阶段毕业论文阶段给自己的指导，从最初的定题，到资料收集，到写作、修改，到论文定稿，她给了我耐心的指导和无私的帮助。为了指导我们的毕业论文，她放弃了自己的休息时间，她的这种无私奉献的敬业精神令人钦佩，唐老师对学科前沿的敏锐把握能力、严谨的治学态度和丰富的实践经验以及一丝不苟的工作作风，使我受益匪浅。无论是在平时学习阶段，还是在论文写作过程中，唐老师对我的教诲和对我在做论文时各方面的支持，我会永远记在心里。她严谨的工作作风和治学态度，深刻地影响着我，是我今后学习和工作的楷模。在此我向唐老师表示我诚挚的谢意。正是由于他们，我才能在各方面取得显著的进步，在此向他们表示我由衷的谢意，并祝所有的老师培养出越来越多的优秀人才，桃李满天下！参考文献1曾兴雯，刘乃安，孙献璞.扩展频谱通信及其多址技术M.西安：西安电子科技大学出版社，20042李建新，刘乃安，刘继平.现代通信系统分析与仿真-MATALAB通信工具箱M.西安：西安电子科技大学出版社，20013徐明远，邵玉斌.MATLAB仿真在通信与电子工程中的应用M.西安：西安电子科技大学出版社，20054刘树棠.现代通信系统M.北京：电子工业出版社，20065 张蕾，郑实勤. 基于MATLAB的直接序列扩频通信系统性能仿真分析研究J. 电气传动自动化，2007年第3期.6 张蕾，郑实勤.基于MATLAB的直接序列扩频通信系统性能仿真分析研究M.电气传动自动化，20077 张大亮，黎勇.扩频通信系统及MATLAB仿真.自然科学报.2008.3.8 查光明，熊贤祚. 扩频通信M. 西安电子科技大学出版社，1990年.9 邵保华. 软件无线电在扩频通信发射中的应用研究D. 哈尔滨工程大学，2004年.附录选取m序列作为仿真系统的扩频码，产生函数如下：function mout = mseq(stg, taps, inidata, n) %变量含义说明 % stg m序列阶数% taps 线性移位寄存器的系数% inidata 序列的初始化% n 输出序列的数目 % mout 输出的m序列if nargin < 4 n = 1; end mout = zeros(n,2stg-1); fpos = zeros(stg,1); fpos(taps) = 1; for ii=1:2stg-1 mout(1,ii) = inidata(stg); % 输出数据的存储 num = mod(inidata*fpos,2); % 反馈数据的计算 inidata(2:stg) = inidata(1:stg-1); % 线形移位寄存器的一次移位 inidata(1) = num; % 返回反馈值 end if n > 1 for ii=2:n mout(ii,:) = shift(mout(ii-1,:),1,0); end end发射机部分发射机部分包括扩频和调制两部分。 QPSK调制部分此处采用QPSK调制，相比SIMULINK的M-PSK调制的实现要求要低一些。function iout,qout=qpskmod(paradata,para,nd,ml) %各变量含义说明 % paradata 输入数据% iout 输出的实部数据% qout 输出的虚部数据% 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; count2=count2+ml; end扩频部分function iout, qout = spread(idata, qdata, code1) %变量含义说明 % idata 输入序列实部 % qdata 输入序列虚部 % iout 输出序列实部 % qout 输出序列虚部% code1 扩频码序列switch nargin case 0 , 1 error('lack of input argument'); case 2 code1 = qdata; qdata = idata; end hn,vn = size(idata); hc,vc = size(code1); if hn > hc error('lack of spread code sequences'); end iout = zeros(hn,vn*vc); qout = zeros(hn,vn*vc); for ii=1:hn iout(ii,:) = reshape(rot90(code1(ii,:),3)*idata(ii,:),1,vn*vc); qout(ii,:) = reshape(rot90(code1(ii,:),3)*qdata(ii,:),1,vn*vc); end信道仿真部分实现高斯白噪声代码函数function iout, qout = comb2(idata, qdata, attn) %各变量含义说明% idata 输入序列实部% qdata 输入序列虚部 % iout 输出序列实部 % qout 输出序列虚部% attn 根据信噪比得到的信号衰减水平 v = length(idata); h = length(attn); iout = zeros(h,v); qout = zeros(h,v); for ii=1:h iout(ii,:) = idata + randn(1,v) * attn(ii); qout(ii,:) = qdata + randn(1,v) * attn(ii); end 接收机的解扩部分函数function iout, qout = despread(idata, qdata, code1) %各变量含义说明% idata 输入序列实部 % qdata 输入序列虚部 % iout 输出序列实部 % qout 输出序列虚部 % code1 扩频码序列 switch nargin case 0 , 1 error('lack of input argument'); case 2 code1 = qdata; qdata = idata; end hn,vn = size(idata); hc,vc = size(code1); vn = fix(vn/vc); iout = zeros(hc,vn); qout = zeros(hc,vn); for ii=1:hc iout(ii,:) = rot90(flipud(rot90(reshape(idata(ii,:),vc,vn)*rot90(code1(ii,:),3); qout(ii,:) = rot90(flipud(rot90(reshape(qdata(ii,:),vc,vn)*rot90(code1(ii,:),3); end接收机解调部分的函数function demodata=qpskdemod(idata,qdata,para,nd,ml) %各变量含义说明% idata 输入数据的实部 % qdata 输入数据的虚部% demodata 解调后的数据% para 并行的信道数% nd 输入数据个数% ml 调制阶数demodata=zeros(para,ml*nd); demodata(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; %误码率分析 noe2 = sum(sum(abs(data-demodata); nod2 = user * nd * ml; noe = noe + noe2; nod = nod + nod2; % fprintf('%dt%en',ii,noe2/nod2); end %数据文件 ber = noe / nod; fprintf('%dt%dt%dt%en',ebn0,noe,nod,noe/nod); fid = fopen('BER.dat','a'); fprintf(fid,'%dt%et%ft%ftn',ebn0,noe/nod,noe,nod); fclose(fid); err_rate_final(ebn0+6)=ber;end%性能仿真图 figuresemilogy(SNR,err_rate_final,'b-*'); xlabel('信噪比/dB')ylabel('误码率')axis(-5,10,0,1)grid on初始化设置以及仿真循环设置Main函数为clear all;clc sr = 256000.0; % 符号速率ml = 2; % 调制阶数 br = sr * ml; % 比特速率 nd = 100; % 符号数SNR=-5:1:10; % Eb/No %*扩频码初值设定* user = 1; % 用户数stage = 3; % 序列阶数 ptap1 = 1 3; % 第一个线性移位寄存器的系数disp('-start-');%*扩频码的产生*for ebn0=-5:1:10 % m序列code = mseq(stage,ptap1,regi1,user); code = code * 2 - 1; clen = length(code); %*