2DPSK信号数字化解调技术研究

1、2DPSK信号数字化解调技术研究2009-08-11 11:49:43 | 天际的荒草 发表于分类: Documents 标签：2DPSK, 2DPSK， MATLAB, DFT, digital demodulate, Discrete Fourier Transform, MATLAB, SDR, software radio, 数字化解调, 软件无线电 0目  录摘要.IABSTRACT II第一章  绪论11.1 软件无线电概述.11.1.1软件无线电的基本概念及由来11.1.2 软件无线电的主要特点及优点11.1.3 软件无线电的基本结构21.2 软件无线电中的调

2、制解调技术.51.2.1 软件无线电中的调制解调问题51.2.2 正交调制原理.61.2.3 正交解调原理.81.3 课题设计内容.10第二章  几种基于DFT的数字化解调技术.112.1 DFT的基本原理.112.2 AM信号解调.112.2.1 解调方法.112.2.2 解调失真度122.3 QDPSK信号数字化解调.122.3.1 解调方法.132.3.2 位同步方法132.3.3 QDPSK信号的数字化解调.14第三章  2DPSK信号数字化解调算法.163.1 2PSK及2DPSK信号原理.163.2 基于DFT的2DPSK信号数字化解调方法.173.3 解调中的

3、位同步方法.183.4 相关说明与讨论19第四章  2DPSK信号数字化解调的仿真研究.214.1 仿真软件工具的选择.214.2 程序流程图.224.3 仿真主要代码介绍.264.4 仿真输出结果.27第五章  毕业设计总结305.1 基于DFT的2DPSK信号数字化解调算法的优点.305.2 设计中遇到的问题和不足之处305.3 心得体会31结束语.32参考文献33摘  要软件无线电（SDR）自从1992年被正式提出以来，受到了越来越多的关注，应用也越来越广泛。软件无线电技术一个很重要的内容就是信号的数字化调制与解调。本文首先阐述了SDR的概况，然后介绍了几种

4、信号的数字化解调算法。并且提出了一种在软件无线电中基于离散傅立叶变换（DFT）算法的二进制相对相移键控（2DPSK）信号的数字化解调算法，通过对载波周期内的采样值进行DFT来提取相位信息并恢复原始调制信号。用MATLAB对这个算法进行了仿真研究的结果表明：采用该方法可以正确的实现2DPSK信号的解调，与传统的解调方法相比,不仅解调过程简单，易于实现，计算量小，而且抗干扰性能得到了明显的改善。将该方法用于2DPSK信号方式的数字化接收机设计中具有实际的意义。关键词：  软件无线电， DFT， 2DPSK，数字化解调， MATLABABSTRACTSince the software r

5、adio (SDR) has been formally proposed in 1992, it has been attentioned more and more, and its application increased extensivly. An important element of the SDR is the digital signal modulation and demodulation. This paper first describes the SDR profiles, and then introduces several digital signals

6、demodulation algorithm. And presents a DFT- based digital demodulation scheme for 2DPSK signal in software radio.The idea is to do DFT on samples in every one (or several) carrier circle to get the phase of the 2DPSK signal and resume the modulated signal. MATLAB algorithm for the simulation results

7、 show that: 2DPSK signal can be demodulated accurately with DFT algorithm. Comparing with traditional methods, not only the process is simple and easy to realize with less computation, but also noise- resisting property is improved obviously. This method has practical significance in design of 2DPSK

8、 signal digital receiver.KEY WORDS： software radio（SDR），  2DPSK，  MATLAB，Discrete Fourier Transform（DFT），  digital  demodulate第一章  绪论1.1软件无线电概述1.1.1软件无线电的基本概念及由来软件无线电（SDRsoftware definded radio）是指在无线电系统中用软件处理的方式来实现智能天线、调制解调等功能，其突出特点是使产品的生产和维护简化，可靠性提高，能兼容多种信号方式，可通过更新软件来升级系统功能

9、等。通信装备从传统电台发展到软件无线电，中间有一个从量变到质变的过程。最开始人们在设计中尝试将所需要支持的多种不同的设备集成在一个结构里，由硬件模块完成在不同设备之间的切换。如果使用这种方法，系统的体积和重量会随着要支持的设备种类增加而增加，显然，这种做法难以解决多少问题。为了避免系统体积无限增大，随着数字信号处理技术的发展，出现了可编程数字无线电，开始尝试硬件模块复用以降低成本。可编程数字无线电可以进行部分重构，简单地修改电台参数，重构功能主要集中在基带处理部分，比如选择前向纠错编码的方式等。电台功能，尤其是发射和接收操作，依然主要由硬件控制，限制了电台向多频段多功能发展的能力。随着越来越多

10、的基本功能改由软件完成，系统的灵活性日益增强。另外，器件技术的进步驱使模拟/数字接口逐渐靠近天线，可以由软件实现更多的功能，模拟电路对射频特性的影响不断降低。由此逐渐产生一种全新的体系结构，诞生出软件无线电。1992年5月，MILTRE公司的Jeo Mitola首次明确提出了软件无线电（soft radio）的概念。其中心思想是：构造一个具有开放性、标准化、模块化的通用硬件平台，将各种功能，如工作频段、调制解调类型、数据格式、通信协议等用软件来完成，并使宽带A/D（模拟/数字）和D/A（数字/模拟）转换器尽可能靠近天线，以研制出具有高度灵活性、开放性的新一代无线通信系统。可以说这种电台是可用软

11、件控制和再定义的电台。选用不同的软件模块就可以实现不同的功能，而且软件可以升级更新，其硬件也可以像计算机一样不断地升级换代。由于软件无线电的各种功能是用软件实现的，如果要实现新的业务或调制方式只要增加一个新的软件模块即可。同时，由于它能形成各种调制波形和通信协议，故还可以与旧体制的各种电台通信，大大延长了电台的使用周期，也节约了开支。有些人也把软件无线电称为“超级计算机”。1.1.2 软件无线电的主要特点及优点软件无线电的主要特点可以归纳如下:（1）具有很强的灵活性  软件无线电可以通过增加软件模块，很容易增加新的功能。可以与其他任何电台进行通信，并可以作为其他电台的射频中继。可以通

12、过无线加载来改变软件模块或更新模块。为了减少开支，可以根据所需功能的强弱，取舍选用的软件模块。（2）具有较强的开放性  软件无线电由于采用了标准化、模块化的结构，其硬件可以随着器件和技术的发展而更新或扩展，软件也可以随需要而不断升级。软件无线电不仅能和新体制电台通信，还能和旧体制电台兼容。这样，既延长了旧体制电台的使用寿命，也保证了软件无线电本身有很长的生命周期。软件无线电的主要优点有:(1)简少了通信设备的硬件电路, 使系统的可靠性大大提高, 生产和维护因此变得十分简单;(2)可以通过更新软件来实现系统功能和性能指标的升级;(3)可以以很低的代价来扩展系统的功能, 达成对多种信号体

13、制兼容的工作能力;(4)通过改进数字信号处理算法, 可以获得比基于电路的传统通信设备更好的性能指标。软件无线电这一新概念一经提出，就得到了全世界无线电领域的广泛关注。由于软件无线电所具有的灵活性、开放性等特点，使其不仅在军、民无线通信中获得应用，而且将在其他领域例如电子战、雷达、信息化家电等领域得到推广，这将极大促进软件无线电技术及其相关产业（集成电路）的迅速发展。1.1.3 软件无线电的基本结构软件无线电的基本思想1是以一个通用、标准、模块化的硬件平台为依托，通过软件编程来实现无线电台的各种功能，从基于硬件、面向用途的电台设计方法中解放出来。功能的软件化实现势必要求减少功能单一、灵活性差的硬

14、件电路，尤其是减少模拟环节，把数字化处理（A/D和D/A变换）尽量靠近天线。软件无线电强调体系结构的开放性和全面可编程性，通过软件的更新改变硬件的配置结构，实现新的功能。软件无线电采用标准的、高性能的开放式总线结构，以利于硬件模块的不断升级和扩展。理想软件无线电的组成结构如图1.1所示。图1.1  软件无线电结构框图软件无线电主要有天线、射频前端、宽带A/D-D/A转换器、通用和专用数字信号处理器以及各种软件组成。软件无线电的天线一般要覆盖比较宽的频段，例如1MHz2GHz，要求每个频段的特性均匀，以满足各种业务的需求。例如，在军事通信中，可能需要VHF/UHF的视距通信、UHF卫星

15、通信，HF通信作为备用通信方式。为便于实现，可在全频段甚至每个频段使用几付天线，并采用智能化天线技术。射频前端在发射时主要完成上变频、滤波、功率放大等任务，接受时实现滤波、放大、下变频等功能。在射频变换部分，宽带、线性、高效射频放大器的设计和电磁兼容问题的处理是比较困难的。当然，如果采用射频直接数字化方式，射频前端的功能可以进一步简化，但对数字处理的要求提高。要实现射频直接带通采样，要求A/D转换器有足够的工作带宽（2GHz以上），较高的采样速率（一般在60MHz以上），而且要有较高的A/D转换位数，以提高动态范围。目前8位A/D转换器的工作带宽已做到1.5GHz以上。模拟信号进行数字化后的处

16、理任务全由DSP软件承担。为了减轻通用DSP的处理压力，通常把A/D转换器传来的数字信号，经过专用数字信号处理器件（如数字下变频器DDC）处理，降低数据流速率，并把信号变至基带后，再把数据送给通用DSP进行处理。通用DSP主要完成各种数据率相对较低的基带信号的处理，例如信号的调制解调，各种抗干扰、抗衰落、自适应均衡算法的实现等，还要完成经信源编码后的前向纠错（FEC）、帧调整、比特填充和链路加密等算法。由于DSP技术和器件的发展，高速、超高速的数字信号处理器不断涌现，如TMS320C6X,ADSP21160等，DSP已能基本满足软件无线电的技术需求。如果采用多芯片并行处理的方法，其处理能力还将

17、大大提高。软件无线电的结构基本可以分为3种：射频低通采样数字化结构、射频带通采样数字化结构和宽带中频带通采样数字化结构，如图1.2图1.4所示。图1.2  射频低通采样数字化的理想软件无线电结构射频低通采样数字化的软件无线电，其结构简单，把模拟电路的数量减少到最低程度，如图1.2所示。从天线进来的信号经过滤波放大后就由A/D进行采样数字化，这种结构不仅对A/D转换器的性能如转换速率、工作带宽、动态范围等提出了非常高的要求，同时对后续DSP或ASIC(专用集成电路)的处理速度要求也特别的高，因为射频低通采样所需的采样速率至少是射频工作带宽的2倍。例如，工作在1MHz1GHz的软件无线电

18、接收机，其采样速率至少需要2GHz，这样高的采样率，A/D能否达到暂且不说，后续的数字信号处理器也是难以满足要求的。图1.3所示的射频带通采样软件无线电结构可以较好的解决上述射频低通采样软件无线电结构对A/D转换器、高速DSP等要求过高，以至无法实现的问题。图1.3  射频带通采样软件无线电结构这种射频带通采样软件无线电结构与低通采样软件无线电结构的主要不同点是，A/D前采用了带宽相对较窄的电调滤波器，然后根据所需的处理带宽进行带通采样。这样对A/D采样速率的要求就不高了，对后续DSP的处理速度也可以随之大大降低。但是需要指出的是，这种射频带通采样软件无线电结构对A/D工作带宽的要求

19、（实际上是对A/D中采样保持器的速度要求）仍然是比较高的。宽带中频带通采样软件无线电结构与目前的中频数字化接收机的结构是类似的，都采用了多次混频体制或叫超外差体制，如图1.4所示。这种宽带中频带通采样软件无线电结构的主要特点是中频带宽更宽（例如20MHz），所有调制解调等功能全部由软件加以实现。中频带宽是这种软件无线电与普通超外差中频数字化接收机的本质区别。显而易见，这种宽带中频带通采样软件无线电结构是上述三种结构中最容易实现的，对器件的性能要求最低，但它离理想软件无线电的要求最远，可扩展性、灵活性也是最差的。图1.4  软件无线电的中频数字化结构1.2软件无线电中的调制解调技术1.

20、2.1 软件无线电中的调制解调问题为了能够进行无线传输以及达到多路复用、提高抗噪声的目的, 在发射端需要用代表信息的基带信号(模拟或数字的低通信号)去控制载波的参数变化, 这就是调制; 在接收端则需要从调制信号中恢复原来的基带信号, 这就是解调。已经有许多不同类型的调制方式可以适应不同的需要。在软件无线电系统中, 调制和解调都是用程序来实现的(也称为全数字化调制解调)。要编写出各种类型调制信号的调制解调软件, 关键是确定信号处理算法。可以利用 FPGA(现场可编程逻辑器件)来实现需要的调制解调算法, 其计算速度比 DSP 更快, 但是灵活性及控制功能较差, 需要与 DSP 或单片机配合使用。建

21、立调制解调算法及程序的一条途径是把模拟电路的工作原理软件化2。比如要对 AM信号进行相干解调, 或建立载波同步乘法器、低通滤波等软件模块的做法虽然可行, 但是计算量很大。实际上, 根据软件无线电的特点, 可以建立与调制解调电路工作原理有所不同的调制解调算法。图 1.5为 SDR 接收机中广泛使用的数字正交解调方案。这是一个具有通用性的解调模型, 对不同方式的调制信号只需要设计相应的基带解调算法。对于 AM 信号, 基带解调算法为 A(n) =(I(n) +Q(n) ) 。对 LPF 的输出进行数据抽取是因为基带信号 I、Q 需要的采样率远低于对调制信号的采样率。这种解调方案利用软件中可以实现的

22、平方和开方运算而免去了复杂的载波同步过程,不仅减少了计算量, 也避免了因载波同步误差而引起的解调误差 (相位同步误差和比较小的频率同步误差都不影响解调效果)。因为仍然是相干解调, 所以这种解调方案具有良好的抗干扰性能。图1.5  正交数字化解调器但是这种解调方法的计算量还是比较大的,因为对每一次采样值都要分两路进行乘法和阶数较高的低通滤波。由于计算量大(或者说对计算速度要求高), 正交数字化解调目前也难以用价位合理的 DSP 来实现, 于是就出现了数字下变频器(DDC)这样的专用芯片。DDC 通过载波跟踪用数控振荡器产生 Cos 和 Sin 数字序列, 分别与调制信号的采样序列相乘,

23、 并完成低通滤波和数据抽取, 输出 I、Q 信号供 DSP 进行位同步和码元解调。由于 DDC 是电路实现的, 可以采用并行处理结构 (同时进行滤波所需要的多次乘法运算), 所以能够解调计算量大的问题。目前已经产品化的软件无线电接收机, 在结构上大多是采用的是“中频采样数字下变频器DSP 基带解调”方案, DDC也就成了软件无线电的关键技术之一。另一种解决方案就是零中频解调, 也称为直接变换法, 即用与发射载波同频率的模拟 Cos 和 Sin 信号分别与接收到的调制信号相乘, 经过 LPF 输出 I、Q 信号, 采样后供 DSP 进行基带解调。比如在基于雷达工作原理的第二代射频身份识别系统(R

24、FID,需要接收和解调受到射频 IC 卡控制的反射电磁波)中, 由于收、发是在同一个设备中, 可以把发射载波直接用于接收解调, 就可以采用“直接变换法”解调方案。DDC 器件价格较高, 功耗大, 功能和结构缺乏灵活性, 直接变换法性能不及中频解调方案, 所以, 如果研发出计算量小、解调质量好,能够用“中频采样DSP 软件解调”方案实现的解调算法, 是很有意义的。1.2.2 正交调制原理为了采用统一的硬件结构来实现多种调制方式，我们必须寻求一种通用的解调方法，正交调制的理论很好地解决了这一问题。尽管调制样式有多种多样，但实质上不外乎用调制信号去控制载波的一个或几个参数，使这个参数按照调制信号的规

25、律而变化的过程。调制信号的数学表达式为：S（n）A（n）cos(n)n+(n)调制信号可以分别“寄生”在已调信号的振幅、频率和相位中，相应的调制就是调幅、调频及调相这三人调制方式。将上式改写为：S(n)=A(n)cos n+(n)式中， 表示载波的角频率。所以，   S(n)=A(n)cos（n）cos( n)- A(n)sin（n） sin( n)=X (n) cos( n)-X (n) sin( n)式中，  X (n)= A(n)cos（n）X (n)= A(n)sin（n）这就是我们所希望获得的同相和正交两个分量，由上式可以看出它们包含了信号的幅度和相位信

26、息，根据X (n)、X (n)，就可以对各种调制信号进行解调。通用调制算法如图1.6所示。图 1.6  数字正交调制通用模型现以MSK信号正交调制为例来说明，二进制频移键控（MSK）信号：Smsk=cosn+2kTbn+k应用三角恒等式展开该式，并考虑到k±1，k0或，则可得：Smskn=coskcosn2Tbcosn-kcosksinn2Tbsinn  Ikcosn2Tbcosn+Qksinn2Tbsinn式中，  Ik=cosk ,  Qk= -kcosk 。DMSK（数字MSK）正交调制器实现框图如图1.7所示。图1.7  数字

27、MSK正交调制框图1.2.3  正交解调原理多种调制信号在理论上都可以通过正交解调算法实现解调3。正交解调必须首先实现正交分解，图1.8示出了数字正交下变频法实现正交分解的框图。图1.8  数字正交解调的通用模型经过数字正交下变频和低通滤波后形成I 、Q 两路正交基带信号，解调算法主要是用两路正交信号计算出信号的幅度和相位。解调算法如下：调幅解调：A(n)I2n+Q2(n)                 

28、              (1.1)调相信号解调：n=tan-1Q(n)I(n)                                 (

29、1.2)调频信号解调：f(n)=n=n-n-1                           (1.3)由式（1.2）和式（1.3）可见，对于调频调相信号的解调算法需要求出瞬时相位n，即需做除法与反正切运算，这在硬件实现上比较困难，若用FPGA实现则要消耗很多门资源，故需寻找简化算法。假设In和Qn为n时刻I路和Q路采样值，In-1

30、和Qn-1为n1时刻I路和Q路采样值，则：Dot(n)=InIn-1+QnQn-1cosncosn-1+sinnsinn-1cos(n-n-1)cos()和    cross（n）=InQn-1-QnIn-1=cosnsinn-1-sinncosn-1=sin(n-n-1)=sin()分别为前后两个数据采样点的相位差的正弦和余弦值，我们分别称Dot（n）为点积，cross（n）为叉积。通过点积值或符号可判断相差大小，可用于DPSK和DMSK的解调；通过叉积的值和符号可判断频差大小，可用于GMSK解调；通过点积与叉积符号可组合判决相差大小，从而对DQPSK、/4 D

31、QPSK、D8PSK等调制方式解调。点积与叉积只有和、差与乘法运算，用FPGA很容易实现。而且叉积所提取的频差可直接用于载波跟踪。1.3课题设计内容数字化调制是指用软件产生出调制信号的采样序列，再通过D/A转换得到模拟的调制信号，数字化解调则是指对已调波信号进行A/D转换，再通过数据处理来实现对信号的解调。数字化调制、解调是软件无线电技术（SDR）中的一个重要内容。SDR主要依靠软件来完成接收系统的各项功能，如智能天线、信号识别、调制解调等，其优点在于可以使产品的硬件大大简化，可靠性大大提高，便于生产和维护，可以通过更新软件来实现产品的功能升级等。/4QDPSK信号相对一般的QPSK信号具有频

32、谱更加集中，更有利于实现位同步的优点。数字化调制的基本要求是产生性能好的调制信号波形，计算量小。本课题要求：（1）对SDR的基本概念进行研究，重点是数字化调制、解调技术。（2）设计一个基于离散傅立叶变换的2DPSK信号的数字化解调算法。（3）用MATLAB语言编程产生出具有典型性的2DPSK信号。（4）对2DPSK信号进行数字化解调。（5）对解调的误码率进行研究。第二章 几种基于DFT的解调技术2.1 DFT的基本原理离散傅立叶变换（DFT）对于有限长序列是一种非常重要的数学变换。因为其实质上是有限长序列傅立叶变换的有限点离散采样，从而开辟了频域离散化的道路，使数字信号处理论可以再频域采用数字

33、运算的方法进行，这样就大大增加了数字信号处理的灵活性。更重要的是DFT有多种快速算法，统称为快速傅立叶变换（FFT），从而使信号的实时处理和设备的简化得以实现。因此，时域离散系统的研究于应用在许多方面取代了传统的连续时间系统。所以说，DFT不仅在理论上有重要意义，而且在各种信号的处理中亦起着核心作用。一般DFT被定义为：Xk=DFTxn=n=0N-1x(n)e-j2nk/N其中， k=0,1,2，N-1它的逆变换IDFT定义为：x(n)=IDFTX(k)=1Nk=0N-1X(k)ej2nk/N其中， n=0,1,2，N-1其中，e-jcos+jsin 。可以证明，离散傅立叶变换的逆变换是唯一的

34、。DFT的快速算法FFT的出现，使DFT在数字通信、语音信号处理、图象处理、功率谱估值、仿真、系统分析、雷达理论、光学、医学、地震以及数值分析等各个领域都得到广泛应用。2.2 AM信号解调2.2.1 解调方法文献4中提出了一种基于DFT（离散傅立叶变换）的AM信号解调算法，要点是对低中频AM信号进行整周期采样（比如取采样频率为载波频率的8倍），对每一个载波周期内的采样点（记为x1x8）进行DFT，计算出载波的幅值A(n)：I(n)=k=18xkcos(2k8)=0.707(x1-x3-x5+x7)-x8+x4Q(n)=k=18xksin(2k/8)=0.707(x1+x3-x5-x7)+x2-

35、x6A(n)=(I(n)2+Q(n)2)1/2显然，去除直流成分后，A(n)序列便是需要的解调输出。与一般的正交解调算法相比较, 由于省去了低通滤波和数据抽取过程, 对采样数据基本上只做加减运算, 每 8 个采样点才做一次平方、开方运算, 计算量大大降低, 为采用“中频采样DSP 解调”方案创造了条件。采用较低的采样频率( 比如每个载波周期采样 4 个点) 也可以正常解调, 当然较高的采样频率对抑制噪声是有利的。2.2.2 解调失真度要正确的完成对AM信号的解调，必须考虑输出的失真度，只有失真度限制在一定的范围内，才可能不失真的恢复原始调制信号（通常应小于0.1）。失真度的表达式如下式：DPx

36、Ps其中，Px表示谐波功率，Ps表示输出信号功率。在不存在定时误差的情况下，失真度很小，通过计算得出D3.4*10-7，可以完全忽略不计，改变调幅系数对解调输出失真没有明显的影响。在存在定时误差时，会对失真度有一定的影响。所谓定时误差是指由于AM信号载波频率偏离设计值（或接收端对载波频率测量存在误差）和采样定时精度有限等原因，使得对AM信号在一个载波周期内的平均采样点数不是整数。若对AM信号每20个载波周期采159.5点（定时误差大约为0.0031），得到失真度为D0.0081。每20个载波周期采159点，得到D=0.02.可见定时误差会使失真度增大，但在通常的定时误差条件下，失真度仍然是可以

37、忽略的。2.3 QDPSK信号数字化解调按照定义产生的PSK信号在码元切换时会发生相位跳变，需要经过限带滤波后才能够发射（或者先对基带信号进行限带滤波，再进行正交调制），在接收机的解调器之前也会对信号进行带通滤波，所以实际接收到的PSK信号的码元波形会分成两种区域，码元的中间部分是稳定区，前、后部分位过渡区。稳定区内的波形接近正弦波，过渡区内的波形则不是正弦波，并且幅度明显降低。因为信道噪声在整个码元内是相同的，所以在过渡区内信噪比很低。提出数字化解调方案应该考虑上述条件。为了表述方便，下面先以一个QDPSK信号为例来对解调方法进行说明，然后再就一般情况进行讨论。设进入解调器的中频QDPSK信

38、号一个码元包含10个载波周期，中间5个周期为稳定区，信号上叠加了高斯型噪声，载波频率已知（由收发双方约定，或在接收机中利用某种方法测定），采样频率取为载波频率的8倍，码元数据与相位跳变的关系是：000,  10/2，  11，  013/22.3.1解调方法先假设已经实现位同步，取码元稳定区内的40个采样值（计为x x ）进行DFT，求出5次谐波（即载波）谱值的实部I和虚部Q，再求出相位作为本码元的相位，用本码元相位减去前一码元相位，便得到本码元的相位跳变值T。有关的计算公式如下：I=k=140xkcos5*2k/40=140k=140xkcosk/4Q=140k

39、=140xksink/4                       (2.1)=tan-1|Q/I|当 I0,Q0 时    ；   当I0,Q<0时      -  ；当 I<0,Q<0时   +；  当 I<0

40、,Q0时    - ；           （2.2）码元的相位跳变为：Tb-a+d      （mod 2）                         （2.3）其中b是本

41、码元相位，a是前一码元相位，d是调整位同步点时的附加相位（见下段内容）。从T到码元解调数据的判决条件为:/4T-/4  判为 00；3/4T/4  判为015/4T3/4  判为 11；7/4T5/4  判为10         （2.4）2.3.2位同步方法：正交数字化解调是对基带信号I、Q执行位同步算法，而本方案是直接用QDPSK信号的采样值进行位同步，其原理是：按载波周期（对应连续8个采样值）进行DFT求出载波幅值A1，在码元切换处A1将出现极小值，所以通过寻找A1的极

42、小值点，就可以实现位同步。初始位同步： 为了提高位同步的定位精度，每接受到4个采样值就用最新的8个采样值（记为x1x8）计算一次A1：I1=18k=18xkcosk/4;Q1=18k=18xksink/4;A1=I12+Q12                               

43、(2.5)对前后的A1值进行比较，就可以找出极小值点，如果该点的A1值也明显小于A1的正常值，则该点就是一个码元切换点（而不是由于干扰引起的）。继续寻找下一个码元切换点，如果前后两个码元切换点的距离（用采样点数表示）近似等于80的整数倍，便认为位同步成功，可以开始进行码元解调。否则就继续寻找下一个码元切换点。位同步维持：  对每一个码元，在计算过码元的相位后，便从该码元的倒数第2个载波周期开始按（2.5）式计算A1值，并寻找码元切换点，如果在4个载波周期内没有出现或者找到的码元切换点与原来的位同步点相同，则位同步点不调整；否则，就要进行位同步点调整：如果找到的码元切换点超前于现有的位

44、同步点，就把位同步点前移一个采样点，如果找出的码元切换点落后于现有的位同步点，就把位同步点后移一个采样点。如果调整了位同步点，就需要在后一个码元的相位跳变值上加上一个修正值d，见（2.3）式。位同步点后移时d/4，位同步点前移时d-/4 。2.3.3 QDPSK信号的数字化解调文献5中提到了一种基于DFT的QDPSK信号解调算法。接收到的限带 QDPSK 信号的码元波形分为两个区域, 中间部分是稳定区, 前、后部分为过渡区。信号波形在码元稳定区内基本上是正弦波, 在过渡区内幅度明显降低。设进入解调器的中频QDPSK信号的一个码元包含10个载波周期，中间7个周期为稳定区，采样频率取为载波频率的8

45、倍。解调方法：用码元稳定区内的56个采样值（记为x1x56）进行DFT，求出7次谐波（即载波）谱值的实部I和虚部Q，在求出本码元的相位，用本码元相位减去前一码元相位，即得到本码元的相位跳变值，进而可判决得到本码元的数据。DFT部分的算式如下：I=156k=156xkcosk/4Q=156k=156xksink/4=tan-1 (I+j*Q)QDPSK信号解调的位同步可以利用DFT算法实现。按载波周期（连续8个采样值）进行DFT求出载波幅值A1,因为在码元切换出A1将出现极小值，所以通过在码元过渡区内寻找A1的极小值点就可以实现初始位同步及位同步维持。第三章2DPSK信号数字化解调算法3.1 2

46、PSK及2DPSK信号原理6二进制移相键控（2PSK）方式是受键控的载波相位按基带脉冲而改变的一种数字调制方式。设二进制符号及其基带波形与以前假设的一样，那么，2PSK信号的形式一般表示为：e0t=nang(t-nTs)cosct这就是说，在一码元持续时间内，当发送二进制符号0时，e0t取0相位；发送二进制符号1时，e0t取相位。这种以载波的不同相位直接去表示相应数字信息的相位键控，通常被称为绝对移相方式。但我们看到，如果采用绝对移相方式，由于发送端是以某一个相位做基准的，因而在接收系统中也必须有这样一个固定基准相位做参考。如果这个参考相位发生变化（0相位变为相位或者相位变成0相位），则恢复的

47、数字信息就会发生0变为1或者1变为0，从而造成错误的恢复。考虑到实际通信时参考基准相位的随机跳变是可能的，而且在通信中不易被觉察，比如，由于某种突然的骚动，系统中的分频器可能发生状态的转移、锁相环路的稳定状态也可能发生转移等等。这样，采用2PSK方式就会在接收端发生错误的恢复。这种现象，通常称为2PSK方式的“倒”现象或者“反向工作”现象。为此，实际中一般不采用2PSK方式，而是采用一种所谓的相对（差分）移相（2DPSK）键控方式。2DPSK方式即是利用前后相邻码元的相对载波相位值去表示数字信息的一种方式。例如，假设相位值用相位偏移表示（定义为本码元初相与前一码元初相之差），并设=0 

48、;                         数字信息“0”=                        数字信息“1”

49、则数字信息序列与2DPSK信号的码元相位关系可举例表示如下：数字信息：      0  0  1  1   1   0   0   1   0   12DPSK信号相位： 0  0  0    0         0   0  或者：    

50、;   0    0   0   0       02DPSK信号的波形与2PSK信号的波形不同，2DPSK信号的波形的同一相位并不对应相同的数字信息符号，而前后码元相对相位的差才是唯一决定信息符号的因素。这说明，解调2DPSK信号时并不依赖于某一固定的载波相位参考值，只要前后码元的相对相位关系不破坏，则鉴别这个相位关系就可正确恢复数字信息，这就避免了2PSK方式的倒现象发生。二进制移相键控系统在抗噪声性能以及信道频带利用率等方面比二进制FSK和OOK优越，因而被广泛应用于数字通信中。

51、考虑到2PSK方式有倒现象，故它的改进型2DPSK信号受到重视。目前，在话带内以中速传输数据时，2DPSK是CCITT建议选用的一种数字调制方式。3.2 基于DFT的2DPSK信号数字化解调方法实际中接收到的2DPSK信号在经过带通滤波后，由于码元跳变出的高频分量被过滤掉，滤波后的 2DPSK 信号波形分为稳定区和过渡区, 码元中间部分是稳定区, 前、后部分为过渡区7。稳定区内的信号基本无损失, 波形近似为正弦波, 而过渡区内的波形则不是正弦波, 并且幅度明显降低。调制信息基本上只存在于码元稳定区。从上述分析出发, 可以得到如图 4.1 所示的基于 DFT 的数字解调方案。图3.1 

52、  DFT解调算法框图对每个码元稳定区内的采样点按照下面式子做DFT：I=1Nk=1Nxkcos2k/nQ=1Nk=1Nxksin2k/n其中, n 代表每个载波周期的采样点个数, N 代表做 DFT 时使用的稳定区内的采样点个数（通常取多个载波整周期）。然后，提取出前后码元的相位跳变信息T来进行解调判决：计算tan-1（Q/I）,并根据Q和I的正负情况确定的取值范围。当 I0,Q0 时    ；   当I0,Q<0时      -  ；当 I<0,Q<0时

53、   +； 当 I<0,Q0时    - ；把本码元的相位记为b，前一码元的相位记为a，则Tb-a+d    （mod 2）                          （3.1）其中d是进行了位同步点调整时附加的相位。d±m*2/nm为位同步调整

54、时移动的采样点个数，n为每个载波周期的采样点个数。从T到码元解调数据的判决条件为：2T2   判为 0；2T32    判为 1；图3.2  2DPSK移相信号矢量图可见，在每个码元周期只需要计算一次相位值即本码元的相位，然后相减得到跳变相位，就可以依据判决条件恢复原始数据，而不需要像文献7中所提到的对每个码元要随着窗函数的移动多次计算谱值，因而大大减轻了计算量，非常适合于软件无线电的数字化实时解调。3.3 解调中的位同步方法位同步是数字化解调中不可缺少的一个重要环节。在本方案中，对采样点按照载波周期做DFT提取幅度信息则可以建立起初

55、始位同步。其原理是：由于稳定区内的信号基本无损失，而过渡区内信号的幅度则明显降低，如果按照载波周期对采样值进行DFT，在码元波形过渡区得到的载波幅值A1会比在码元稳定区内求得的A1值小，而且存在一个极小值，当该极小值小于给出的判决门限值时（防止噪声引起的误判），可确定为码元的切换点。因此通过对A1值进行检查就可以确定码元切换点。在实际的通信中，为方便初始位同步的建立，可以发射一连串的1序列，以保证前后码元存在相位的跳变，更有利于寻找码元切换点。在初始位同步建立后，由于码元定时误差和噪声干扰，位同步点可能有偏差，还必须进行调整。对每个码元计算过相位后，在过渡区仍旧按载波周期计算A1值，并检测幅值

56、极小点。没有出现幅值极小点（表明前后码元相位连续）或者找到的幅值极小点与原位同步点相同，则位同步点不调整；否则，根据本次计算出来的幅值极小点对已建立的位同步点进行调整，并根据调整结果对码元的跳变值附加修正值d，且d的大小由公式（3.2）决定。d±m*2/n                            

57、;                 (3.2)其中，m为位同步调整时移动的采样点个数，n为每个载波周期的采样点个数。3.4 相关说明与讨论（1）在进行码元解调时也可以把靠近稳定区的一些采样值用于DFT计算（采样值的数目仍需要与整数个载波周期相对应），以求进一步降低误码率，因为那里的信噪比也比较高。其实稳定区和过渡区并无一个明确的分界点。（2）在进行位同步计算时，如果一个码元包含的载波周期比较多（比如一个码元包含50个载波周期），也可以每

58、n个采样值（一般情况下，为了提高位同步定位精度，每接收到n/2个采样值就用最新的n个采样值计算一次A1）进行一次DFT，求出I1和Q1值。I1和Q1计算公式是：I1=1nk=1nxkcos2k/nQ11nk=1nxksin2k/n（3）在一般的举例中我们会假设一个码元周期刚好等于整数个载波周期，但是解调方法并不要求一定如此。比如一个码元包含10.25个载波周期，将有82个采样点（假设n=8），只要把码元过渡区算为5.25个载波周期，仍然使用稳定区内的N个采样值来进行DFT计算就可以了（计算码元相位跳变的公式需要做相应改变）。(4)在实际系统中不可避免会存在位同步误差（确定的码元起点与实际起点有

59、偏差），只使用码元稳定区内的采样值进行解调计算，位同步误差的影响也可以消除或减小。（5）一般情况下举例时，采样频率是取为载波频率8倍时候多，如果取为4、5或者9、10等整数，也可以实现解调（有关计算公式需做修改）。但是降低采样频率不利于抑制噪声，会使误码率有所增大。选择8或者4还可以减小计算量，比如为8时，则对码元采样点按照前面公式做DFT时，可求得：I=1Nk=1Nxkcos2k/n=140k=140xkcosk/40.025x1+x9+x17+x25+x33-x5-x13-x21-x29-x37+0.707(x2+x10+     

60、0;    x18+x26+x34-x4-x12-x20-x28-x36-x6-x14-x22-x30-x38+          x8+x16+x24+x32+x40)可以看出，几乎只需要做加、减运算就可以了。另外，考虑到需要进行载波解调和位同步，在倍数为8时每个码元采样80点已经是比较少的了。（6）如果存在采样定时误差（采样频率不等于载波频率的整数倍），也会使误码率增大，但是不一定很严重。（7）本解调算法对于各种MDPSK以及/4QDPSK信号都是适用的，只要修改从码元的

61、相位跳变值到调制数据的判决条件式就可以了，也可以推广到MQAM（幅相联合调制）信号的解调，因为通过DFT不仅可以求出码元的相位及相位跳变，也可以求出码元的载波幅值。图3.3  采用中频自动调整方案的SDR接收机结构（8）为了使采样频率为中频的整数倍，可以采用如下的电路方案：按照设定的采样频率进行采样，在解调程序中去检查在码元稳定区中的40个采样点是否刚好对应5个载波周期，如果有偏差，再去控制变频器的本振（VCO）频率，使中频向正确的方向改变。相应的软件接收机结构如图3.3所示：第四章  2DPSK信号数字化解调的仿真研究4.1 仿真软件工具的选择在目前工程界使用的各种仿真工

62、具中，MATLAB是应用最广泛，也是使用最简便的一个，相比于其他各种仿真工具，它具有很多独特的优点。因此，本次设计中的算法仿真我们选用MATLAB作为仿真工具。MATLAB是矩阵实验室（Matrix Laboratory）之意。除去本身卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。MATLAB的基本数据单位是矩阵，它的指令表达式与数学工程那个中常用的形式十分相似，故用MATLAB来解算问题要比用C语言、FORTRAN等语言完成相同简捷得多。当前流行的MATLAB5.3/Simulink3.0包括拥有数百个内部函数的主包和三十几种工具包（Toolbox）。工具包又可以分为功能性工具包和学科工具包。功能性工具包用来扩充MATLAB的符号计算，可视化建模仿真，文字处理及实时控制等功能。学科工具包是专业性比较强的工具包，控制工具包，信号处理工具包，通信工具包等都属于此类。开放性使MATLAB广受用户欢迎。除内部函数外，所有MATLAB主包文件和各种工具包都是可读可修改的文件，用户通过对源程序的修改或加入自己编写程序构造新的专用工具包。一种语言之所以能如此迅速的普及，显示出如此旺盛的生命力，是由于它有着不同于其他语言的特点，正如同FORTRAN和C等