新·TD-SCDMA基带信号分析处理算法.docVIP

TD-SCDMA上行基带信号分析处理算法 1 算法流程及总体框架 我们的信号分析处理算法可以概括为：先进行下行同步，再搜寻SYNC_UL并进行上行Midamble同步，然后利用Midamble校正频偏、相偏，最后再解码。流程图如下： 具体各个模块的功能和算法原理在下文详细介绍。 2 算法详解 2.1 下行同步 下行同步模块的功能是查找下行导频时隙（DwPTS）的位置并确定小区使用的下行同步码号。具体分两步实现： 利用DwPTS附近特殊的功率结构查找DwPTS的大致位置。（初搜） 2、利用SYNC_DL良好的自相关性来确定SYNC_DL的码号及其精确位置。（精搜） 2.1.1 下行同步初搜 TS0(848chips) GP(16chips) GP(32chips) 下行同步码 GP(96chips) 上行同步码 由于DwPTS附近特殊的功率结构（DwPTS两边的GP在实际信道中功率近乎为0，而下行同步码则以全功率发射，如上图所示），则用下行同步码的64码片的功率和除以其两边各32码片的功率和，所得商值要远远大于在其他位置用相同方法所得的商值。 据此，我们可以构造一个长为128码片的滑动窗，对接收到的数据进行逐码片滑动，每次滑动时，求出其中间64码片的功率和与两边各32码片功率和的比值，最后求取这些比值中的最大值，其最大值所对应的位置即为DwPTS的大致位置。 具体实现过程如下所示： W=128*N;%滑动窗长，N为采样倍率 Wm=64*N;%滑动窗中间部分长度 Wk=32*N;%滑动窗的一侧长度 Ap=sum(abs(Mux_Data(1:6400*N)).^2)/(6400*N); %Ap为信道平均功率 KK=zeros(1,6400*N+W); %KK用来存放每次滑动所得的商值 for n = 1:(6400*N+W) B(n) = sum(abs(Mux_Data((n+Wk):(n+Wk+Wm-1))).^2); %B为滑动窗中间部分功率和 if B(n)3*Ap*Wm %设定一个门限 A(n)=sum(abs(Mux_Data(n:(n+Wk-1))).^2)+... sum(abs(Mux_Data((n+W-Wk):(n+W-1))).^2); %A为滑动窗两边功率和 KK(n) = B(n)/A(n); else KK(n)=0; end end end [row,Position]= max(KK); %求最大值及其位置 注意： 程序中间有一个门限判定的过程，其作用是对一些明显不可能的位置直接略去，不予计 算，这样可大大减少计算量。 N为采样倍率，在此处使用的是一个4倍采样的数据。 对所有商值进行绘图，其结果如下： 可求得其最大值位置为15893，此值即为DwPTS 的大致起始位置。 2.1.2 下行同步精搜 SYNC_DL具有良好的自相关性，相同的SYNC_DL序列的相关值要远大于不同SYNC_DL序列的相关值。据此，我们在上一步计算出DwPTS 大致起始位置的基础上，以SYNC_DL的起始位置为中心向前、向后各延伸一定的距离（此处为16*N），作为SYNC_DL的可能起始位置。用码表中的32个SYNC_DL码分别与之进行相关计算，并求其最大值，则此最大值所对应的SYNC_DL码号即为基站所使用的SYNC_DL码号，所对应的位置即为SYNC_DL的起始位置。 具体实现过程如下所示： Sync_W=Mux_Data(Position+16*N:Position+112*N-1); %构造精搜窗，在初搜SYNC_DL位置前后延伸16*N作为SYNC_DL的可能起始位置进行相关 for n=1:32*N+1 %32*N+1个位置的循环 for n1=1:32 %32个SYNC_DL码的循环 res=abs(corrcoef(Sync_W(n:n+64*N-1),tempbin3(n1,:))); %tempbin3为用码表生成的SYNC_DL序列，corrcoef为相关计