科斯塔斯环总结报告

1、Costas环研究总结报告一、研究目的:高速移动通信系统中，收、发送机的相对运动使得接收信号不可避免地形成快衰落 信道，从而引入较大的多普勒偏移和多普勒变化率，故此必须纠正或补偿多普勒频偏和多普勒变化率，以保证收发系统间的信号同步。由于较大的多普勒频偏和多普勒变化率的存在，造成接收载波与给定载波在频率和相位上均存在较大的差别，为正确解调需使用载波跟踪机制来提取相干载波，costas环路是高效、可靠的跟踪环路，在成熟系统中常用来来提取接收信号相干载波，以实现对接收信号进行精确相干解调，确保收发系统信号间高可靠和高精度的同步。二、性能要求:依系统指标需求，要求该环路能够狗在 100 US （约合2

2、00个复值信号点）之内达到求解稳定，输出相干载波的实时频偏估测值，并纠正频偏，纠偏范围依系统指标需求应达到1khz-30khz （载波2Mhz）;频率抖动范围在频偏的 10%以内（噪声w 1Odb）;相位抖动范围 在土 50 （噪声w 10db）。基于costas环路的高效性和可靠性以及较低的实现代价，本课题经过仿真验证，选择该环路作为该功能模块的实现环路。三、costas环路基本原理概述：Costas环路在系统架构中的接入示意图如图1所示：图1costas环路系统架构中的接入图Costas环路功能结构图如图2所示:图2costas环路功能结构图m(t)是基带调制信号 cos(Wct)是调制载

3、波，假定环路锁定，且不考虑噪声影响，则振荡器输 出的两路互为正交的本地载波分别为=cos(wct d )v2 =sin (wct)(t)= (t)为本地载波与输入相干载波之间的相位差。由图1所示输入信号 m(t)cos(wct)分别与锁定后的本地载波即式(2)、(3)相乘得到1v3 = m(t)cos(wct) cos(wct r(t) m(t)cos&(t) cos(2wct (t)】 2(4)1v4 二 m(t)cos(w ct) sin( Wet 二)m(t) I si n(珂t) si n(2w ct 二)】2(5)经过低通滤波器滤除 2倍分量得到1v 5 m (t ) c o

4、s (t ()1v 6 m (t ) s pn (t ()V5和V6相乘得到vd(t)二 KdSin(2珥t)(8)可以把上面部分看做是一个鉴相的过程，vd (t)就是鉴相结果。环路滤波器的功能结构图如图3所示：图3环路滤波器(10)则vt)经过环路滤波器得到 Vd (t)二v门(t) K1 Vd (t - 1)K2环路滤波器输出相角 vd(t) , vd (t)控制vco的相位和频率，使 Vd (t) 一直变小直到接近 于0,此时锁相环进入稳定状态，此时 V1就是所需的相干载波，而 V5就是相干解调输出。四、costas环路具体算法实现im = m(t)sin (wc _ wd)t:(13)

5、QPS信号S(t)Q路+I路SinQ路I路imQ路输出90 相孑移环路滤L 右求相角振荡器-波器""cosI路输出re图4 QPSK信号相干载波提取与相干解调实现F变频后得到两路信号(假设在理想状态下点的偏转只和频偏有关忽略噪声的影响)(9):为初相，w c等于频偏和载波的比值。.j (Wd t)vco输出信号可以表示为 e令(9)与(10)相乘m(t)ej(wc)訂心=m(t)ej(WcWd)t"(10)(11)re = m(t)cos(wc - wd )t 1(12)鉴相：1）反正切函数法二 arctan(m/re)二 arctansi n( Wc-Wd)t

6、)(14)cos(wc - Wd )二(-Wd)t + ® - k (k= 0,1,2)(21 + 1)基于qpsk调制(n =1,2,3” 4Qpsk调制所有的4点都位于45,135, 225,315，若想正确解调信号，则解调后的所有点都应位于旋转45的坐标轴上。相角为(15)= (Wc-Wd)t + d4(k = 0,1,2)(n =0,123) 42)符号函数法,_ sign(im) re _ sign(re) im；re2 im2(16)3）乘法器+低通滤波器法令已调信号 m(t)cos(wct) vco输出 sin(wct J1v7 二 m(t)cos(Wct) sin (

7、Wet ) m(t)si n( :) cos(2w c t )】 2(17)经过低通滤波器滤除 2倍分量得到= 1m(t )sin()(18)综上所述三种鉴相方法，（2）和都用到v - sin（可（这公式只在二很小时成立，所以（1）精确度最高，本研究采用反正切函数法。环路滤波器：环路滤波器对输入信号的噪声起抑制作用， 并调节环路的矫正速度。 常用的环路滤波器 有无源比例积分滤波器，有源比例积分滤波器，理想比例积分滤波器。选择理想比例积分滤 波器进行讨论：图3理想比例积分滤波器kY(k)二 K1X(k) K2' X(n)n=0(19)K1和K2是此环路滤波器的两个系数,1KoKd1KoK

8、d8匚 WnTs二2_44 WnTs (WnTs)4( WnTs)244 WnTs (WnTs)2(20)K 二 KoKd 为环路增益，u是锁相环阻尼系数,8 bWn,为环路固有角频率，Bl是环路噪声带宽，Ts是环路滤波器的采样周期（即NCO相位调整时间间隔）根据（20）式计算环路滤波器系数：1）确定阻尼系数（工程取值一般为0.707）。2）根据锁相环的跟踪精度要求及跟踪范围要求确定其等效噪声带宽（通常取BL：0.1Rb Rb为信息数据速率）。3）由量化的阶数及 NCO输出的量化幅度值计算环路增益。VCO用于根据给定的相角产生正弦和余弦三角函数。根据接收信号得到相角，经过环路滤波器后去控制VC

9、O的输出，使得相角逐渐变小，最终使相角减小到很小的数值，得到稳定的wd可以估计出频率偏移，re和im就是信号的相干解调后的实部和虚部。五、VHDL仿真实现已下变频调制后的 QPSK信号I和Q两路，载波2Mhz，噪声10dB，设定环路滤波器系 数 6=12，c2=1/2A7。VHDL各模块端口描述：图4 VHDL结构图顶层模块COMPONENT costas_PLL3PORT(clk200m : IN std_logic;200M时钟用于高速计算clk2m : IN std_logic;2M时钟 信号输入时钟rst : IN std_logic;复位信号高有效IMn : IN stdogic_v

10、ector(15 downto 0);输入信号虚部RE_in : IN std_logic_vector(15 dow nto 0);输入信号实部chu_enable : IN std_logic; chu 序列有效标志位 data_e nable : IN std_logic; 有效数据有效标志位IM_out : OUT stdo gic_vector(15 dow nto 0); 虚部输出 RE_out : OUT stdogic_vector(15 downto 0);_-实部输出 chu_v : OUT std_logic; -chu 序列标志位 data_v : OUT stdo g

11、ic; -有效数据标志位 clk_out : OUT std_logic; -输出时钟);END COMPONENT;复乘模块Comp onent complexport( clk:in std_logic;输入 200M 高速时钟ar:in stdogic_vector(15 downto 0); 输入数据实部 ai:in stdogic_vector(15 downto 0); 输入数据虚部br:in stdogic_vector(15 downto 0);-相角余弦输入 bi:in stdogic_vector(15 downto 0);-相角正弦输入pr:out std_logic_v

12、ector(31 downto 0); 输出实部pi:out std_logic_vector(31 dow nto 0);e nd component; 输出虚咅 B鉴相模块comp onent phase_out port(clk:in std_logic; 输入 200M 高速时钟reset:i n std_logic; 复位信号enable:in std_logic;鉴相使能标志位x_in:in stdogic_vector(15 downto 0); 输入复乘器的实部结果y_in:in stdogic_vector(15 downto 0); 输入复乘器的虚部结果phase_out:

13、out std_logic_vector(15 downto 0); 相位角输出(相位响应输出)end comp onent;环路滤波器模块comp onent loop_filter port(clk:in std_logic; 200M 高速工作时钟reset:i n std_logic; 复位信号coun t2:in stdo gic_vector(15 dow nto 0);计数器(循环计数周期 100)phase_in:in stdogic_vector(15 downto 0); 输入鉴相器输出的相角角 freq_part1:out std_logic_vector(15 down

14、to 0); 频率响应输出phase_out1:out stdogic_vector(15 downto 0);-输出一个经环路滤波器平滑后的相位角 end comp onent;VCO模块 component cordic-用 cordic 模块作为 VCOport(clk:in std_logic;-200M 高速时钟phase_i n:in stdo gic_vector(15 dow nto 0); 环路滤波器输出的相位角si n:out std_logic_vector(15 dow nto 0);输入的相角的正弦函数输出cosi ne:out stdo gic_vector(15

15、dow nto 0); 输入的相位角的余弦函数输出 end comp onent;仿真结果输出入下图：进入锁相环的信号星座图如下图:进入的锁村!坏数星座图据图&信号经过锁相环后输出的星座图（此星座图是锁相环锁定后估计出平均频偏后逐点纠频所得，由于初相角未知固此星座图有个固定偏移）得出结论：在130点后锁相环进入锁定状态， 估计出频偏大概为6Khz,频偏抖动范围是500hz，相位在-5范围内抖动，星座聚拢在四个点上，锁相环成功锁定。锁相环的锁定频偏的范围曲n-raPE!n 口1000-1500-200025002000 -1000 0 1000 2000In-Phase图5 :进入锁相环的信号星座图信号经过锁相环后其频率响应曲线如下图箍相环换率响应曲绒图6:锁相环频率响应曲线图信号经过锁相环的相位响应曲线如下图:图7锁相环相位响应曲线图信号经过锁相环后的星座图如下图经过锁相环后输出| Sav侶 FigiHE|amroJPEno可以通过改变程序中的环路滤波器系数来改变，本程序的锁定范围大概在800hz-30Khz。锁相环锁定时间和其精度是负相关的，若要缩短锁定时间必然要牺牲精度，反之