发明名称 一种多入多出-正交频分复用系统同步方法
摘要 本发明公开了一种MIMO-OFDM同步方法,它通过利用IFFT变换的特性,将发射端各天线插入的训练序列,在各发射天线到达接收天线信号延迟不同时,叠加后仍然可以得到两个相同的半段序列。在接收端的时域进行时间粗同步,频域进行时间精同步;收端时间同步取得后,可容易的取得频率同步。由于本发明通过在发端插入偶位均为零,奇位均为训练序列的同步符号,可以区分不同时延,适用于分布式MIMO系统。因此它具有易于实现,实用性强,应用面广等多种优点;同时它又具有带宽占用少,频谱利用率高的优点。
申请公布号 CN100561999C 申请公布日期 2009.11.18
申请号 CN200610020797.7 申请日期 2006.04.26
申请人 电子科技大学 发明人 肖悦;彭小勇;雷霞;李少谦
分类号 H04L27/26(2006.01)I;H04B7/26(2006.01)I 主分类号 H04L27/26(2006.01)I
代理机构 代理人
主权项 1、一种多入多出-正交频分复用同步方法,它由发端和收端两部分组成,具体步骤如下:发端对发射信号的处理步骤如下:步骤1:定义发射天线的个数为M,M取正整数,且M的数目大于1;选择M个长均为N<sub>1</sub>的PN序列b<sub>m</sub>[k],N<sub>1</sub>取正整数;将这个PN序列b<sub>m</sub>[k]尾部加一个“1”,构成长度为Q的序列c<sub>m</sub>(k),k∈[0,Q-1],m∈[1,M];此时的c<sub>m</sub>(k)取值为复数形式,即c<sub>m</sub>(k)∈{1+j,-1-j};步骤2:将步骤1得到的序列c<sub>m</sub>(k)进行插零处理,生成每个天线发射的训练序列T<sub>m</sub>(i),i∈[0,N-1],N为OFDM系统的FFT点数,且取正整数;插零处理的具体做法是:序列c<sub>1</sub>(k)的第一个信息c<sub>1</sub>(1),插入到第一个发射天线的训练序列T<sub>1</sub>(i)的第一位;序列c<sub>1</sub>(k)的第二个信息c<sub>1</sub>(2),插入到训练序列T<sub>1</sub>(i)的第2M+1位;第三个信息c<sub>1</sub>(3)插入到T<sub>1</sub>(i)的第4M+1位,以此类推,直到第k个信息c<sub>1</sub>(k)插入到T<sub>1</sub>(i)的第2(Q-1)M+1位,训练序列T<sub>1</sub>(i)的其余位均插零;对第二个发射天线而言,序列c<sub>2</sub>(k)的第一个信息c<sub>2</sub>(1),插入到对应训练序列T<sub>2</sub>(i)的第三位;序列c<sub>2</sub>(k)的第二个信息c<sub>2</sub>(2),插入到训练序列T<sub>2</sub>(i)的第2M+3位;第三个信息c<sub>2</sub>(3)插入到训练序列T<sub>2</sub>(i)的第4M+3位,以此类推,直到第k个信息c<sub>2</sub>(k)插入到T<sub>2</sub>(i)的第2(Q-1)M+3位,训练序列T<sub>2</sub>(i)的其余位均插零;以此类推,序列c<sub>M</sub>(k)的第一个信息c<sub>M</sub>(1),插入到第一个发射天线的训练序列T<sub>M</sub>(i)的第2M-1位;序列c<sub>M</sub>(k)的第二个信息c<sub>M</sub>(2),插入到训练序列T<sub>M</sub>(i)的第2M+2M-1位;第三个信息c<sub>M</sub>(3)插入到T<sub>M</sub>(i)的第4M+2M-1位,以此类推,直到第k个信息c<sub>M</sub>(k)插入到T<sub>M</sub>(i)的第2(Q-1)M+2M-1位,训练序列T<sub>M</sub>(i)的其余位均插零;这里的N和发射天线数M之间必须满足:N=2MQ;<img file="C2006100207970002C1.GIF" wi="1322" he="137" />步骤3:将步骤2得到的训练序列T<sub>m</sub>(i)对应插入数据序列1到数据序列m中;具体插入方法是:训练序列T<sub>1</sub>(i)插入到第一个天线上的数据序列1中,训练序列T<sub>2</sub>(i)插入到第二个天线上的数据序列2中,以此类推,直到训练序列T<sub>M</sub>(i)插入到第M个天线上的数据序列M中;T<sub>m</sub>(i)插入各天线数据序列的位置均是相同的,然后将插入后的序列,进行IFFT运算,得到的结果再进行加循环前缀CP处理,这样得到的结果再进行组帧处理后,由各射频天线发射出去;收端对接收信号的处理步骤如下:步骤4:在第p路接收天线上,将接收到的发端发射的序列,通过一个大小为N的滑动窗,p∈[1,P],P是接收天线数,取正整数,且大于1;将窗口中的数据r<sub>i</sub>分为前后两个半段,前半段是r<sub>i</sub>,0≤i≤N/2-1,后半段是r<sub>i</sub>,N/2≤i≤N-1;将这两个半段序列前面τ长度的信息都去掉,用后面的信息对应求相关,得到相关值φ(d),d表示时间偏移,取正整数:<maths num="0001"><![CDATA[<math><mrow><mi>&phi;</mi><mrow><mo>(</mo><mi>d</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>&tau;</mi></mrow><mrow><mi>N</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn></mrow></munderover><mrow><mo>(</mo><msubsup><mi>r</mi><mrow><mi>d</mi><mo>+</mo><mi>i</mi></mrow><mo>*</mo></msubsup><mo>&CenterDot;</mo><msub><mi>r</mi><mrow><mi>d</mi><mo>+</mo><mi>N</mi><mo>/</mo><mn>2</mn><mo>+</mo><mi>i</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤5:将步骤4得到的相关值φ(d)进行归一化处理,得到一个相关峰M(d):<maths num="0002"><![CDATA[<math><mrow><mi>M</mi><mrow><mo>(</mo><mi>d</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>2</mn><mo>&CenterDot;</mo><mo>|</mo><mi>&phi;</mi><mrow><mo>(</mo><mi>d</mi><mo>)</mo></mrow><mo>|</mo></mrow><msup><mrow><mo>(</mo><mi>R</mi><mrow><mo>(</mo><mi>d</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,R(d)为进行公式(2)运算的信息功率之和<maths num="0003"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>(</mo><mi>d</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>&tau;</mi></mrow><mrow><mi>N</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>|</mo><msub><mi>r</mi><mrow><mi>d</mi><mo>+</mo><mi>i</mi></mrow></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>&tau;</mi><mo>+</mo><mi>N</mi><mo>/</mo><mn>2</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>|</mo><msub><mi>r</mi><mrow><mi>d</mi><mo>+</mo><mi>i</mi></mrow></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤6:然后设置一个硬判门限和长度为N的滑动窗,将从相关峰M(d)中超过门限的时间偏移点所对应的滑动窗中的序列开始,连续将L个长度为N的序列选出,设这段序列为g<sub>i</sub>(t),i=0,1,...,L-1,t=0,1,...,N-1;记录超过门限的时间偏移点为<img file="C2006100207970003C4.GIF" wi="79" he="88" />,L取正整数,且其取值大于1;步骤7:将步骤6得到的L个长度为N的序列g<sub>i</sub>(t),分别进行FFT运算,得到G<sub>i</sub>(k),i=0,1,...,L-1,k=0,1,...,N-1:<maths num="0004"><![CDATA[<math><mrow><msub><mi>G</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>g</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><mo>-</mo><mn>2</mn><mi>&pi;jnk</mi><mo>/</mo><mi>N</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤8:然后将步骤7得到的序列G<sub>i</sub>(k),按照发端步骤2中插入训练序列的方式,将对应位置的信息抽取出来,具体的抽取方式为:检测第一路发射天线发射的信息到达第p路接收天线的时间精同步点,将G<sub>i</sub>(k)中第一位,第2M+1位,第4M+1位,以此类推到第2(Q-1)M+1位信息抽取出来;检测第二路发射天线发射的信息到达第p路接收天线的时间精同步点,将G<sub>i</sub>(k)中第三位,第2M+3位,第4M+3位,以此类推到第2(Q-1)M+3位信息抽取出来;以此类推,检测第M路发射天线发射的信息到达第p路接收天线的时间精同步点,将G<sub>i</sub>(k)中第2M-1位,第2M+2M-1位,第4M+2M-1位,以此类推到第2(Q-1)M+2M-1位信息抽取出来;和本地序列c<sub>m</sub>(k)进行相关相乘,就可以得到第m路发射天线信号的时间精同步点<img file="C2006100207970004C2.GIF" wi="118" he="70" /><maths num="0005"><![CDATA[<math><mrow><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mrow><mn>2</mn><mo>,</mo><mi>m</mi><mo>,</mo><mi>p</mi></mrow></msub><mo>=</mo><munder><mrow><mi>arg</mi><mi>max</mi></mrow><mi>i</mi></munder><mo>{</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>Q</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>G</mi><mi>i</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mn>2</mn><mi>m</mi><mo>-</mo><mn>1</mn><mo>+</mo><mn>2</mn><mi>jM</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>c</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,G<sub>i</sub><sup>*</sup>(2m-1+2jM)表示对G<sub>i</sub>(2m-1+2jM)求共轭的结果;步骤9:将步骤6得到的时间偏移点<img file="C2006100207970004C4.GIF" wi="61" he="69" />与步骤8得到的时间精同步点<img file="C2006100207970004C5.GIF" wi="98" he="69" />求和,得到第m路发射天线发射的信号到达第p路接收天线的时间同步点:<maths num="0006"><![CDATA[<math><mrow><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mrow><mi>m</mi><mo>,</mo><mi>p</mi></mrow></msub><mo>=</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mrow><mn>1</mn><mo>,</mo><mi>p</mi></mrow></msub><mo>+</mo><msub><mover><mi>&theta;</mi><mo>^</mo></mover><mrow><mn>2</mn><mo>,</mo><mi>m</mi><mo>,</mo><mi>p</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤10:由步骤9得到时间同步点后,得到时间调整后的信息;再将时间调整后的信息进行频率同步,得到频率同步处理后的信息;然后将频率同步处理后的信息进行去循环前缀CP处理;将去循环前缀CP处理后的信息再进行FFT处理。
地址 610054四川省成都市建设北路二段四号