基于DSP软件锁相环文档格式.docx

基于DSP软件锁相环文档格式.docx

《基于DSP软件锁相环文档格式.docx》由会员分享，可在线阅读，更多相关《基于DSP软件锁相环文档格式.docx（22页珍藏版）》请在冰点文库上搜索。

数字鉴相器的Z域模型如下：

实现数字鉴相器的方法之一就是借助信号的正交分解，图4是该方法的原理框图。

其中：

LPF表示低通滤波，A是低通滤波器带来的常数增益。

因此鉴相器的输出：

经过反Z变换得到数字环路滤波器的时域表达式为：

通过反Z变换得到数控振荡器的时域表达式：

变量uc（nT）数值较小且变化不会太快,因此式（14）成立：

综合式（11）、式（12）和式（13），得到NCO输出信号的表达式：

如果以数字频率描述数控振荡器，则称其数字中心频率为ω0T，数字偏置频率为ωn2·

uc（nT）·

T。

因此，该数控振荡器的灵敏度与数字灵敏度分别为

·

3多速率条件下的软件锁相环

　　在数字化接收机中，经常碰到多速率条件下的抽样率转换问题。

所谓多速率系统是指在一个数字系统中存在2个或2个以上的抽样率［4］。

构成软件锁相环鉴相器的混频器通常工作在系统采样频率上。

在满足奈奎斯特采样定律的前提下，数字化接收机的系统采样率一般高达数10M。

而数字鉴相器组成部分的反正切表，由于混频之后的数据经过多倍抽取，工作频率已经下降到与信号波特率相近的水平。

数据抽取同时也降低DSP的运算量，由DSP完成的环路滤波的处理速度近似等于信号波特率。

此外由于软件锁相环中的数控振荡器需要给混频器提供同样速率的正交载波，其工作速率与混频器相等，需要进行插来调整速率。

　　为了合理利用DSP有限的计算资源，总是在满足同步需要的前提下尽可能地降低环路滤波的工作速率，也就是通常所说的环路频率。

环路频率是软件锁相环的一个重要参数，他同时决定着锁相环算法的计算量与捕获速度。

环路频率过高将带来额外的计算负担，环路频率太低又不能满足捕获速度的需要，在应用常取系统波特率作为环路频率的大小。

该扩展模型对应的线性相位Z域模型如图5所示。

D表示数据抽取，I表示数据插。

数据在抽取之前先要进行抗混叠滤波，可用于抗混叠滤波的FIR滤波器有CIC滤波器、半带滤波器等。

　　实际情形中，由于零阶保持插几乎不需要额外的运算量，因此经常被采用。

实际上环路频率fL总是能够跟上环路滤波器输出信号的变化速率。

换言之，对环路滤波器输出信号按照fL的速率进行采样保留了他的全部信息，因此零阶保持插对系统性能不会有太大影响［4］。

　　结合上述软件锁相环的基本原理，下面借助Matlab仿真观察软件锁相环的系统响应。

本文以频率阶跃信号作为输入，观察软件锁相环的系统响应从而进一步验证本文建立的一系列软件锁相环模型。

　　设定系统采样频率为1MHz，仿真时间0.1s，信号中心频率125kHz，起始相位-π／4，输入频率阶跃100Hz，起始点为0.02s，抽取因子为8。

一般情况下，都希望环路工作在欠阻尼状态，取阻尼因子ξ=0.707，ωn由2πΔF（快捕带宽）决定［2］，分别取2π*40，2π*50，2π*100。

仿真出相位误差响应曲线、NCO偏置频率曲线和频率阶跃信号的相位曲线，如图6所示。

　　从图6可以看出，软件锁相环在［0，0.02］区间相位误差为0，处于锁定状态。

在t=0.02s时刻，输入信号频率产生了大小为100Hz的阶跃，导致软件锁相环进入捕获过程。

由于软件锁相环的校正作用，当ωn=2π*50时，系统在t=0.05s时刻重又进入同步状态，相位误差依旧为0。

由相位误差响应曲线可以看到，锁相环可以无相差的跟踪频率阶跃信号，同时表明虽然锁相环鉴相误差为0，但是由于环路滤波器的理想积分作用其输出的控制信号并不为0，由该控制信号产生的100Hz偏置频率保证了NCO输出与输入信号的同步。

当快捕带宽发生变化导致改变时，锁相环的捕获速度也发生了变化，快捕带宽越宽，捕获速度越快。

4软件锁相环的DSP实现

　　在宽带数字化接收机的实现中，数字下变频采用通用可编程下变频器HSP50214B。

在实现载波同步、码元同步软件锁相环的整个反馈环路中，数控振荡器、鉴相器由HSP50214B完成，环路滤波在TMS320C6X中完成。

DSP实现框图如图7所示。

环路延时是一个应该重视的因素。

带来软件锁相环环路延时主要有以下2种原因：

（1）环路FIR滤波器带来的延时；

（2）数据等待处理带来的额外延时。

　　在数字化接收机中，采用粗同步与细同步两级。

粗同步环路时延大，反应速度慢；

细同步环路时延小，反应速度快，粗同步保证有效信号落在滤波器的通带之，细同步可以在粗同步基础上获得较大捕获带和同步带。

此外还采用抛弃若干采样点，消除不必要的环路延时。

　　可以看出，软件锁相环具有处理灵活的优点，他摆脱了复杂的硬件电路设计，解决了许多模拟环遇到的难题。

目前，由于DSP功能越来越强大，工作速度越来越高，也为软件锁相技术的发展创造了必要的条件。

基于定点DSP的软件锁相环的设计和实现

论文报告

摘要：

软件锁相环是软件接收机中执行载波恢复功能的关键部分。

提出了一种48位定点扩展精度的算法，可以有效地实现软件锁相环。

与浮点算法比较，能极降低DSP的运算量，降低功耗，同时保证动态围运算精度。

低轨道卫星软件接收机软件锁相环定点扩展精度算法

关键词：

低轨小卫星通信是近年来卫星通信应用中一个方兴未艾的重要领域，"

创新一号"

小卫星是我国研制的具有完全自主知识产权的存储与转发通信小卫星，cascom手持终端是专门为这颗小卫星研制的低功耗地面手持通信终端，支持调制数据速率达76．8kbps的BPSK窄带信道。

基于TI公司的低功耗16位定点数字信号处理器TMS320VC5510（最高运算能力为200MIPS），完全用软件实现低中频数字接收机，其中包括执行载波恢复功能的软件锁相环SPLL（Soft-warePhase-LockedLoops）。

在用浮点算法实现软件锁相环时，由于TMS320VC5510是一个定点DSP处理器；

没有浮点处理单元，只能用编译器产生模拟浮点运算的指令，运算量需要67．2MIPS，效率很低，因此需要一种能在TMS320VC5510上执行的定点算法，有效地降低运算量。

本文提出了一种48位定点扩展精度算法实现SPLL，提高了效率，减少了运算量，同时保证了环路计算的精度和动态围。

1软件锁相环

1．1软件锁相环的结构

图1表示软件数字接收机中的解调器。

它包括由改进的costas环路构成的载波跟踪环路。

采样后的中频信号经过数字混频，滤掉高频分量，通过改进的costas环路产生控制信号，控制数控振荡器（NCO）得到新的本振参考信号。

其中的相位检测器和环路滤波器结构如图2所示。

1．2软件锁相环的设计

由图2可见，SPLL的计算由计算相位误差和更新环路中间变量、输出控制信号两部分组成。

算法描述可用伪码白表示：

//中断发生

{

//读取基带数据

I_baseband=I（n）

Q_baseband=Q（n）

//计算相相位误差d（n）

d（n）=sign（I_baseband）×

K_norm**×

Q_baseband//更新环路中间变量s_pll（n）

s_pll（n）=C2×

d（n）+s_pll（n-1）

//输出控制信号Δf（n）

Δf（n）=C1×

d（n）+s_pll（n）

//设置NCO载波频率f（n）

f（n）=Δf（n）+f0***

}//end

注：

**K_norm是归一化因子，由I_baseband和Q_baseband和初始值决定：

***f0是固定的NCO中心频率

算法描述中的中断周期就是环路采样时间间隔。

中断发生后，第一步读取基带同相项数据和正交项数据

I_baseband=I（n）=Acosθe 

（1）

Q_baseband=Q（n_=Asinθe 

（2）

A是基带信号幅度，θe是相位误差。

第二步计算硬判决的同相数据乘以相位误差。

θe 

（3）

其中硬判决函数如式（4），

以及θe≈sinθe. 

（5）

由式

（1）和

（2）得：

由式（3）、（5）和（6）得：

即把同相数据硬判决后结果乘以正交项数据后再乘以归一化因子K_norm。

K_norm初始值由I_baseband和Q_baseband的初始值决定，

由于定时恢复环路和AGC（自动增益控制）环路的作用，K_norm在解调过程中近似保持恒定。

第三步更新环路中间变量s_pll（n），

d（n）+s_pll（n） 

（9）

第五步设置NCO载波频率f（n）,

f（n）=Δf（n）+f0 

（10）

fo是固定的NCO中心频率。

至此，一次完整的SPLL计算完成。

在第三步和第四步计算中，环路滤波器系数Cl，C2可以通过环路采样时间间隔T（或者环路更新时间间隔）、环路自由频率ωn及环路阻尼系数ξ确定，如下两式：

C1=（1/K0Kd）（8ξωnT）（4+4ωnT+（ωnT）2 

（11）

C2=1/K0Kd（4（ωnT）2/（4+4ξωnT+（ωnT）2 

（12）

Kd为相位检测器的增益，由于在实现软件锁相环时，基带信号的同相和正交分量都经过归一化处理，故Kd=1；

K0为数控振荡器的增益，K0=2πT。

T为调制数据速率的倒数（1／76800），ξ一般取0．707。

在启动载波恢复之前有一个频率捕获过程，通过1024点的FFT，可保证接收的中频信号与本振信号之间的频率差△f0最大.

248位定点扩展精度算法

图3描述的算法在TMS320VC5510上用C语言直接利用浮点运算实现时，只能通过C编译器产生模拟浮点运算的定点指令。

这种方法效率很低，每次环路计算需花费875个指令周期。

在调制数据速率为76．8kbps的数字接收机中，需要67．2MIPS的运算量。

为了降低环路计算的运算量，同时保持浮点运算具有动态围大、精度高的优点，笔者提出了一种48位定点扩展精度计算的方法。

参加运算的每个操作数由三个16位定点数W2、W1、W0级联表示，其中高16位为二进制补码的整数部分，低32位为二进制补码的小数部分，符号位在最高位，又可称为Q15．32格式，如图3所示。

一个Q15．32数的表示围是（-32768，32768），小数分辨率是1/232（2﹒3283e-10），远远超过16位定点表示的精度，即3e-5（1/2^15-1）。

以下用加（ADD_48）、减（SUB_48）、乘（MULT_48）三种基本运算来说明定点扩展精度算法。

操作数X由X2、X1、X0构成，操作数Y由Y2、Y1、Y0构成，结果W由W2、W1、W0构成。

执行48位加法运算时，W2W1W0=X2X1X0十Y2Y1Y0，首先把小数部分X1X0和Y1Y0相加，结果保存到W1W0中，

产生的进位位CARRY与X2、Y2相加，结果保存到W2。

执行48位减法运算时，W2W1W0=X2X1X0－Y2Y1Y0，首先X1X0减去Y1Y0，结果保存到W1W0，产生借位位BORROW，再由X2减去Y2和借位位BORROW，结果保存到W2。

两个Q15．32数相乘时，乘积是一个Q30．64数，出于前面实现SPLL时对动态围和计算精度的要求，该Q30．64数可以双向截位为Q15．32的48位定点数。

具体做法是保留符号位和整数部分的低15位以及小数部分的高32位。

48位定点数的乘法由图4所示。

除了以上加、减、乘三种基本运算外，48位窄点扩展精度算法还包括取负（NEC-48）、数据拷贝（MOVE_48）两种操作。

取负操作即将X1X0取负，结果保存到W1W0，产生借位位BORROW，再用0减去X2和借位位BORROW，结果保存到w2；

数据拷贝，即把X1X0拷贝到W1W0，X2拷贝到W2。

在TMS320VC5510可编程DSP的基础上，利用48位定点扩展精度算法实现SPLL。

在实现过程中，采取了模块化的思路。

首先，把SPLL整个环路计算封装成一个可调用的C语言函数。

函数参数包括C1、C2、K_norm、基带信号的I及Q分量、环路中间变量、调整频率。

DSP的中断例程（ISR）可以直接调用环路计算函数，而且通过输入不同的Cl、C2，适用于不同的载波恢复环路中。

另外，在函数部用汇编语言进行编程，以充分利用DSP的计算能力，把48位定点扩展精度算法的五个基本操作封装成用汇编指令写的宏（macro），对照计算流程，调用这些宏，完成SPLL的核心计算部分。

经统计，每次环路计算需132个指令周期，总的运算量10．1MIPS，是浮点算法运算量（67．2MIPS）的14％。

用48位扩展精度算法实现软件接收机中的SPLL，解决了浮点算法运算量大的问题，同时还具备浮点算法动态围大、精度高的优点，已经成功应用于"

小卫星地面手持低功耗通信终端中。

另外，本文提出的SPLL实现算法，通过修改环路滤波器系数，也可以应用在其他软件接收机中，具有很好的扩展性

数字下变频器HSP50214B在中频数字接收机中的应用

软件无线电技术在侦察接收机中的应用越来越广泛，新型雷达和通信侦察装备的数字接收机大多由软件无线电技术实现。

软件无线电的目标是将数字化技术应用于接近天线的器件上，从而在其输出端直接数字化射频信号，提高系统的灵活性。

超高速A／D采样器件为实现完全数字化的无线电技术提供了可能，但同时对DSP（数字信号处理器）的处理速度提出过高的要求。

因此，将射频信号变频到中频再进行采样，数字信号经专用数字下变频器DDC（DigitalDownConverter）完成信号下变频和抽样处理后再送给DSP处理，则可大大减轻DSP的信号处理负担。

这里介绍专用数字下变频器HSP50214B在中频数字接收机中的应用。

2数字下变频器HSP50214B的部结构

HSP50214B是Intersil公司生产的可编程数字下变频器件。

该器件可将数字信号的载波频率进一步降低，甚至降到基带（即载波频率为0Hz）；

还可抽取输入的数字信号，在允许限度可降低数据量，使信号能够被DSP实时有效的处理，功能强大。

其输出信号直接送至DSP进行后续处理，包括解调、解码和协议控制等。

图1为HSP50214B数字下变频器的结构框图。

HSP50214B数字下变频器的主要功能单元：

（1）输入单元其电平检测单元选取一定长度的数据，比较并累积该段数据的绝对值与预设门限电平的误差，该误差累积量可由外部接口读出，进行外部自动增益控制。

（2）载波NCO单元该单元输出具有一定频率和初始相位的正、余弦两路信号，实现对输入信号的正交混频。

产生I／O通道数据。

（3）滤波抽取单元级联积分梳状（CIC）滤波抽取组、半带（HB）滤波抽取组、255阶FIR滤波器组成抽取／低通滤波器，实现低通滤波和数据抽取；

重采样多相滤波器进行分数倍采样的转换，使整个抽取系统的输出速率能满足特殊场合的要求。

（4）增益控制单元AGC为部增益控制，检测DDC幅度输出和预设门限之间的误差，从而调整FIR滤波器的输出增益，以提高小信号的增益并减少多级抽取造成的幅度衰减。

（5）坐标变换单元直坐标到极坐标转换器和数字鉴频器共同完成对各类幅度、频率和相位调制信号的解调。

（6）输出单元有直接串行输出、直接并行输出和先入先出（FIFO）3种输出方式，可提供同相分量、正交分量、瞬时幅度、瞬时相位和瞬时频率等5种输出数据类型，可根据调制方式和DSP的接口方式选择合适的输出方式。

（7）控制接口单元通过接口数据总线C[7：

0]、地址总线A[2：

0]和读，写信号WR／RD实现控制字及相应参数的写入和部寄存器容及状态信息的读出。

3HSP50214B电路设计及接口配置

在侦察接收机中实现数字中频滤波，数字下变频器是关键部分，嵌入式处理器可配置中频数字接收机中的数字下变频器并接收数据。

接收机系统首先预处理中频信号，然后经A/D采样器进行中频带通采样，采样信号通过HSP50214B实现数字下变频，将中频信号搬移到基带，最后送给DSP处理器进行解调、解扩。

图2为数字下变频器HSP50214B的主要接口电路。

经前端预处理后的模拟中频信号（IF信号）送至A／D采样器进行转换采样，采用ANALOG公司的AD9245A／D采样器，其输出是14位，最高采样速率可达80MHz，输入围大，功耗低，性价比高。

由于HSP50214B的输入为14位，所以将AD9245的14位输出接到HSP50214B的14位输入即可。

HSP50214B的CLKIN引脚与AD9245的时钟相连，每个CLKIN时钟到达就对其输入数据采样1次。

CLKIN的时钟信号可达65MHz，PROCLK时钟是混频后的处理时钟最高不能超过CLKIN时钟。

采用CPLD电路控制数据的读出、写入、寄存器选择等信号。

8位控制接口数据总线C[7：

0]接TMS320VC5402型DSP的D端口，3位寄存器地址线A[2：

0]接CPLD的P1．2～P1．0。

CPLD是Ahera公司MAX7000系列中的EPM7128器件，可将I／O设置在3．3V或5V电源下工作，从而实现电平转换，同时满足5V器件HSP50214B和3．3V器件TMS320VC5402的电压要求。

设置HSP50214B为并行直接输出模式，输出I、Q两路正交分量，16位并口AOUT[15：

0]输出同相分量，BOUT[15：

0]输出正交分量。

DATARDY丙引脚通过CPLD与TMS320VC5402的中断引脚INTO相连，当AOUT端口产生新的数据时，DATARDY产生的负脉冲触发DSP产生中断信号，由其中断服务程序接收并处理数据。

DSP响应数字下变频器产生的中断时，分时并行接收数据，依次在下变频器的AOUT和BOUT两个输出端接收同相与正交分量。

这两个端口分别被映射为I／O空间的两个地址，当DSP对这两个地址读操作，通过CPLD进行译码，分别使能AOUT和BOUT，数据即可顺利被DSP读取，从而进行各种后续的基带处理。

4HSP50214B控制字的设置

HSP50214B共有256个32位的控制字寄存器，其偏移地址是000H~0FFH，控制字决定各功能模块的参数设置，如载波中心频率、滤波器阶数、255阶FIR滤波器系数、抽取因子、输出格式与方式的选择等。

用户通过控制接口可写入控制字，也可读出某些部寄存器的容及状态信息，如电平检测是否完成等。

TMS320VC5402通过HSP50214B的微处理器接口设置部寄存器参数。

HSP50214B通过写入控制字实现解调，每个控制字是32位，每次只能通过数据总线接口C[7：

0]写入8位控制字，需通过地址线A[2：

0]选择写入控制字的高低位，连续写4次，图3为其加载时序。

对控制字写操作的步骤为：

（1）按字节由低到高的顺序，把32位的控制字分为4组，每组8位数据，在每次写选通（WR）信号的上升沿到来时，依次把4组数据装载到HSP50214B主控制器的相应位上，对应地址用二进制表示分别为000、001、010、011（000：

保持寄存器1．001：

保持寄存器2，010：

保持寄存器3，011：

保持寄存器4）；

（2）在写选通（WR）信号的第5个上升沿到来时，将主控制器中接收到的数据装入目标寄存器并锁存，此时，装载目标寄存器的地址到100，即A（2：

0）=100；

（3）写完1个控制字后，必须等待4个处理时钟才能写下一个控制字，以保证数据装载成功。

只有控制字被正确写入，HSP50214B参数选择合理，解调的信号才不会严重失真。

5测试分析

实际应用可根据侦察接收机接收信号类型需求，设置各寄存器的参数，包括CIC、AGC、串并口、输出格式、滤波器的系数等，灵活设置下变频。

模块加电时，DSP初始化过程中将各寄存器的参数设置到HSP50214B的目的寄存器。

可对不同调制方式信号实现灵活解调分析，具有适应性强、软件升级方便等优点。

通过DSP仿真器进行硬件仿真测试，其结果符合理论计算数值，在1．3MHz的带宽围，镜像抑制比均在110dB以上；

当镜像抑制比为110dB时，I、Q两路所对应幅度的相对误差小于0．25％，相位正交性误差小于0．1°

，可满足高端侦察接收机高性能信号处理的要求。

6结束语

HSPS0214B数字下变频器在软件无线电中具有结构开放、软件可编程及功能多样等特点，在军事及民用数字接收机中都有巨大的应用潜力，它可使接收机系统具有良好的灵活性及可扩展性。