(TigerSHARC处理器技术及其应用)第9章TS系列DSP系统设计实例.ppt

9.1.1 基于高阶循环累积量的载频估计　　在移动无线通信信道中，由于存在传播延时及多普勒效应，信号的载波会出现偏移。如果检测器是相干的，那么接收机就必须估计这种载波偏移。盲前向载波偏移估计方法对于解决突发模式传输和高效率调制中的同步问题是非常有潜力的，因为它不需要很长的捕获时间，而且不需要那些浪费带宽的训练序列。因此，盲前向载波偏移估计方法经常应用于广播网络的同步，以及那些匹配滤波器前面的粗载频偏差估计环节。 　　这里利用接收序列的高阶循环统计量来完成数字通信信号的盲前向载波偏移估计。接收机接收的信号可以表示为 9.1.2 通信信号参数分析的硬件及软件实现　　通信信号参数分析的硬件组成如图9.1-1所示。这是一个零中频接收机，接收的射频信号通过正交混频器下变频到直流，本振是一个固定的自由振荡器。由于没有使用锁相本振，因此接收信号和本振的频率与相位必然存在一定的偏差。IQ两路基带信号送到AD进行量化，经过FPGA缓存后，送给DSP进行后续处理，其中包括对残留的载波频率偏差估计等。 图9.1-1 通信信号参数分析硬件组成 　　下面给出在VisualDSP++5.0环境中实现基于高阶循环累积量的载频估计算法。　　在实际硬件实现中IQ数据由AD量化得到，在这里为了程序的完整，首先产生了一段模拟数据。程序如下所示：　　#include stdlib.h 　　#include stdio.h 　　#include math.h　　#include complex.h　　#include FFTDef.h 　　#define M 1024 //数据长度　　#define RAND_RANGE 2147483647.0 //rand函数产生随机数的范围　　#define PI 3.1415927　　#define SNR 5 //信噪比dB　　#define N_FFT 1024 //FFT的点数　　　　#pragma align 4096　　complex_float c40[M]; //四阶共轭时变相关　　#pragma align 4096　　complex_float s[M]; //已调复信号 　　#pragma align 16　　float twiddles[MAX_FFT_SIZE] = {　　 #include twid1k_m.dat //FFT旋转因子　　 }; 　　float C40_Amp[M]; //循环四阶累积量幅度，为了方便绘图，将以下三个变量设为全局变量　　float alfa[N_FFT]; //循环频率　　float S_Amp[M]; //s的FFT幅度　　 　　extern void FFT32(float inbuf[], float outbuf[], int size, int realOrComplex); //FFT函数声明，参见VisualDSP++?　　//平台自带的例子　　 　　int main( int argc, char *argv[] )　　{ 　　 int N;　　 int *sc; //二进制符号序列 　　 float fs; //采样率Hz　　 float fc; //载频Hz　　 float sr; //符号速率Hz　　 int i, j;　　 float theta, rayleigh; //相位，Rayleigh噪声　　 complex_float c20;　　 float sigma; //高斯白噪声方差　　 int N_FFT_2;　　 complex_float S[M]; //s的FFT结果　　 complex_float C40[M]; //循环四阶累积量　　 float fMax;　　 float fe; //载频偏差的估计结果 　　 fs = 320; //采样率Hz　　 sr = 45; //符号速率Hz　　 fc = 20; //载频 Hz　　 　　 N = M*sr/fs+1; //计算二进制符号的个数　　 sc = malloc (N*size