OFDM系统之利用L-STF和L-LTF序列进行同步

OFDM系统之利用L-STF和L-LTF序列进行同步

什么是同步？为什么要进行同步?

同步，简而言之就是在接收端找到OFDM符号的起始点，以对它进行FFT，从而恢复出原始信号。发送端不断发送OFDM信号，而接收机接收到这些信号会有一定的延时，如果不进行同步而直接进行FFT运算，是无法恢复出原始信号的。

什么是L-STF和L-LTF序列？

L-STF（短训练序列）和L-LTF（长训练序列）是IEEE 802.11ac协议中用于进行同步的序列。在此协议中，同步分为两个步骤——粗同步（也称数据包检测）和细同步（符号同步）。其中，粗同步用于判断数据包的到来，细同步用于将数据包的位置精确到bit位。L-STF和L-LTF序列分别用于粗同步和细同步。

如图所示为802.11ac中数据包的结构，L-STF由10个重复的短码元组成，L-LTF由一段符号间隔和2个重复码元构成，如图所示：

如何生成L-STF序列和L-LTF序列？

以下代码利用matlab生成L-STF序列和L-LTF序列
首先定义这样一个长64的序列L_STF

L_STF = [0,0,0,0,0,0,...     
         0,0,1+1i,0,0,0,-1-1i,0,0,0,1+1i,0,0,0,-1-1i,...    
         0,0,0,-1-1i,0,0,0,1+1i,0,0,0,0,0,0,0,-1-1i,...
         0,0,0,-1-1i,0,0,0,1+1i,0,0,0,1+1i,0,0,0,1+1i,...
         0,0,0,1+1i,0,0,...
         0,0,0,0,0];

对其进行IFFT运算，得到L_STF_time，L_STF_time为16个元素一重复的长64的序列，即包含4个重复序列

L_STF_time = ifft(L_STF,64);

对这个序列进行拼接，得到含10个重复码元的序列，即为最终的L-STF序列。

L_STF_time=[L_STF_time,L_STF_time,L_STF_time(1:32)];

对L-LTF序列，也有类似操作

L_LTF = [0,0,0,0,0,0,...   
             1,1,-1,-1,1,1,-1,1,-1,1,1,1,1,1,1,-1,-1,1,1,-1 ... 
             1,-1,1,1,1,1,0,1,-1,-1,1,1,-1,1,-1,1,-1,-1,-1,-1 ...
             -1,1,1,-1,-1,1,-1,1,-1,1,1,1,1,...
             0,0,0,0,0];     
L_LTF_time = ifft(L_LTF,64);
L_LTF_time=[L_LTF_time(33:64),L_LTF_time,L_LTF_time];

其中，在最后拼接时，将原来的L_LTF_time序列中的尾端序列复制到首端，形成符号间隔GI，然后重复拼接，得到最终的L-LTF序列。

如何实现同步？

实现同步有多种方法——自相关算法，互相关算法，最大似然机制等等。这里介绍利用自相关算法实现粗同步和利用互相关算法实现细同步。

粗同步

粗同步是对接收到的序列进行自相关。由于包头内含有多个重复的S序列，可以想象，即使是失真后，序列之间的相似度仍然比一般序列要高，根据这一特点，我们取两个长16的窗口，延序列滑动，对窗口内的数据进行自相关，如果得到的结果能在一段时间内维持在一个临界值以上（该临界值通过经验得出），则可以认为找到了包头，即确认了数据包的到来。但此时得到的位置可能是在序列中的任意位置，我们接下来要通过细同步来确定OFDM符号的起始点。

function [index,peak,keep,p,energy]=Coarse_Syn(rcv_data_noise,L)
    for i=1:L-48
        for m=1:16
            c(m)=rcv_data_noise(i+m-1+16)*conj(rcv_data_noise(i+m-1+32));
    %       c(m)=noise_premble_data(i+m-1)*conj(noise_premble_data(i+m-1+16));
            d(m)=rcv_data_noise(i+m-1)*conj(rcv_data_noise(i+m-1+16));
        end
        p(i)=sum(c);
%         co(i)=sum(con);
%         absp(i)=abs(real(p(i)))+abs(imag(p(i)));
        absp(i)=real(p(i)).*real(p(i))+imag(p(i)).*imag(p(i));
        r(i)=sum(d);
%         absr(i)=abs(real(r(i)))+abs(imag(r(i)));
        absr(i)=real(r(i)).*real(r(i))+imag(r(i)).*imag(r(i));
        peak(i)=absp(i)-absr(i); 
        keep(i)=absr(i)/absp(i);
    end

细同步

细同步的思想跟粗同步类似，只不过细同步是将接收到的序列与本地L-LTF序列进行互相关，找到互相关的结果最大的位置，显然，只有当窗口起点在L-LTF序列的起点时，得到的互相关的结果最大。由于L-LTF序列长度已知，就可以得到OFDM符号开始的位置。

Search_Window = 300;
Data_For_FineSyc = Tx_Continue_AddStfLtf(Coarse_Index_Est:Coarse_Index_Est+Search_Window+160);
for i= 1:Search_Window
    for j=1:160
        Xcor(j)=Data_For_FineSyc(j+i)*conj(L_LTF_time(j));
    end
    Sum_Xocr(i)=abs(sum(Xcor));
end
fine_index=find(Sum_Xocr==max(Sum_Xocr));
figure(1)
plot(Sum_Xocr);
title('Find Syc Peak');
Data_Start=Coarse_Index_Est+fine_index+160

由此，我们就实现了同步。
以上是我自己的理解，如有不对，敬请指正。

参考资料：IEEE802.11ac定时同步与信道估计算法研究链接.