发明名称 基于三维星座图簇调制的二维Chirped OFDM通信方法
摘要 本发明涉及基于三维星座图簇调制的二维Chirped OFDM通信方法,基于三维星座图簇调制的二维多载波通信技术,属于宽带无线数字通信领域。可以用于提高适用于双弥散信道传输的Chirped OFDM符号传输的准确性和安全性。本发明的方法根据正四面体的特殊性质,给出了一种从标准的MQAM星座图直接构造的三维星座图的设计方案,继承了MQAM的诸多优良性质,并且在此基础上增加了符号点间的距离,在新型的二维Chirped OFDM通信系统中减小误符号率,增加了通信的可靠性和安全性。
申请公布号 CN103685132B 申请公布日期 2016.07.06
申请号 CN201310641930.0 申请日期 2013.12.03
申请人 北京理工大学 发明人 郇浩;时鹏飞;陶然;王腾;常诚
分类号 H04L27/34(2006.01)I 主分类号 H04L27/34(2006.01)I
代理机构 代理人
主权项 基于三维星座图簇调制的二维Chirped OFDM通信方法,其特征在于步骤为:步骤1,依据通信原理中的数字带通调制技术关于正交振幅调制的设计方式,将标准的MQAM星座图中的M个星座点以QPSK星座图的点为单位,由外圈沿逆时针方向标记,划分为M/4个QPSK的星座图组合形式;每个QPSK星座图点和三维正四面体星座图点对应,得到M/4个正四面体星座图组合;步骤2,根据MQAM星座图点能量分布特点,即是距离原点的距离,将M/4个划分的QPSK星座图组合中相同能量的划归为一个组,共划分为K组,将它们与原点的K个不同距离,由小到大记为R<sup>k</sup>,k=1,2,…,K;步骤3,三维空间中构建K个同心球面,其半径为R<sup>k</sup>,k=1,2,…,K;在每个球面上构造相同方位的正六面体,取其4个不相邻的顶点构成正四面体星座图;根据步骤2中,划分的K组QPSK星座图组合,将正四面体星座图分组;步骤4,根据能量层由内而外,R<sup>k</sup>由小到大的顺序构建相应层上正四面体组合的星座图;选取内层上R<sup>1</sup>的正六面体的4个不相邻顶点构成一个正四面体,相邻层上的正六面体在同样的方位的4个不相邻顶点构造正四面体星座图,并沿着z轴逆时针旋转π/2;如果位于同一层R<sup>k</sup>上有两组正四面体,首先根据R<sup>k‑1</sup>层的正四面体,旋转π/2构成其中一个正四面体,然后在此基础上沿着z轴逆时针旋转π/2,构成一个正六面体星座图;所有由内而外的M/4个正四面体组旋转角度表示为:(π/2)<sup>m</sup>,m=0,1,…,(M/4)‑1;步骤5,根据正六面体坐标参数化M点三维星座图的坐标为S<sub>0</sub>,S<sub>1</sub>,…,S<sub>M‑1</sub>,在能量归一化准则<img file="FDA0000961177900000011.GIF" wi="799" he="110" />和星座图点间最小欧氏距离最大化准则以及外能量层之间星座点间距不小于最内层R<sup>1</sup>上正四面体星座点间的距离,求出最优的星座图点坐标;步骤6,将初始的3N长度的二进制数据通过比特映射到所求出的M点三维星座图上,由于每个点3bit,得到N个映射后的星座点为Q<sub>n</sub>=(x<sub>n</sub>,y<sub>n</sub>,z<sub>n</sub>)<sup>T</sup>、n=0,1,…,N‑1,其中T表示转置;通过和传统Chirped OFDM系统的融合,将Q<sub>0</sub>,Q<sub>1</sub>,…,Q<sub>N‑1</sub>作为新的二维分数阶域的Chirped OFDM符号:<maths num="0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>S</mi><mo>=</mo><mrow><mo>(</mo><msub><mi>Q</mi><mn>0</mn></msub><mo>,</mo><msub><mi>Q</mi><mn>1</mn></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>Q</mi><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msub><mfenced open = "(" close = ")"><mtable><mtr><mtd><msub><mi>x</mi><mn>0</mn></msub></mtd><mtd><msub><mi>x</mi><mn>1</mn></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>x</mi><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mn>0</mn></msub></mtd><mtd><msub><mi>y</mi><mn>1</mn></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>y</mi><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>z</mi><mn>0</mn></msub></mtd><mtd><msub><mi>z</mi><mn>1</mn></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>z</mi><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced><mrow><mn>3</mn><mo>&times;</mo><mi>N</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000961177900000021.GIF" wi="621" he="357" /></maths>步骤7,对Chirped OFDM系统的符号S做二维的逆分数阶傅里叶变换,行和列对应的分数阶旋转角度分别为α<sub>1</sub>和α<sub>2</sub>,得到的发送端的基带调制信号s=(s<sub>0</sub>s<sub>1</sub>…s<sub>N‑1</sub>),其中s<sub>n</sub>=s(n<sub>1</sub>,n<sub>2</sub>),n=0,1,…,N‑1表示为:<maths num="0002"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>s</mi><mrow><mo>(</mo><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>,</mo><msub><mi>n</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow><mo>=</mo><msqrt><mfrac><mrow><mn>1</mn><mo>-</mo><mi>j</mi><mi> </mi><msub><mi>cot&alpha;</mi><mn>1</mn></msub></mrow><mn>3</mn></mfrac></msqrt><mo>&CenterDot;</mo><msqrt><mfrac><mrow><mn>1</mn><mo>-</mo><mi>j</mi><mi> </mi><msub><mi>cot&alpha;</mi><mn>2</mn></msub></mrow><mi>N</mi></mfrac></msqrt><mo>&times;</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mfrac><mrow><mi>j</mi><mi> </mi><msub><mi>cot&alpha;</mi><mn>1</mn></msub></mrow><mn>2</mn></mfrac><msubsup><mi>n</mi><mn>1</mn><mn>2</mn></msubsup><mo>&CenterDot;</mo><msup><mi>&Delta;t</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow><mo>&times;</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mfrac><mrow><mi>j</mi><mi> </mi><msub><mi>cot&alpha;</mi><mn>2</mn></msub></mrow><mn>2</mn></mfrac><msubsup><mi>n</mi><mn>2</mn><mn>2</mn></msubsup><mo>&CenterDot;</mo><msup><mi>&Delta;t</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>&times;</mo><munderover><mi>&Sigma;</mi><mrow><msub><mi>k</mi><mn>1</mn></msub><mo>=</mo><mn>0</mn></mrow><mn>2</mn></munderover><munderover><mi>&Sigma;</mi><mrow><msub><mi>k</mi><mn>2</mn></msub><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>S</mi><mrow><mo>(</mo><mrow><msub><mi>k</mi><mn>1</mn></msub><mo>,</mo><msub><mi>k</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mfrac><mrow><mi>j</mi><mi> </mi><msub><mi>cot&alpha;</mi><mn>1</mn></msub></mrow><mn>2</mn></mfrac><msubsup><mi>k</mi><mn>1</mn><mn>2</mn></msubsup><mo>&CenterDot;</mo><msubsup><mi>&Delta;u</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><mi>j</mi><mfrac><mrow><mn>2</mn><msub><mi>&pi;n</mi><mn>1</mn></msub><msub><mi>k</mi><mn>1</mn></msub></mrow><mn>3</mn></mfrac></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mfrac><mrow><mi>j</mi><mi> </mi><msub><mi>cot&alpha;</mi><mn>2</mn></msub></mrow><mn>2</mn></mfrac><msubsup><mi>k</mi><mn>2</mn><mn>2</mn></msubsup><mo>&CenterDot;</mo><msubsup><mi>&Delta;u</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mi>j</mi><mfrac><mrow><mn>2</mn><msub><mi>&pi;n</mi><mn>2</mn></msub><msub><mi>k</mi><mn>2</mn></msub></mrow><mi>N</mi></mfrac></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000961177900000022.GIF" wi="1974" he="333" /></maths>其中,0≤n<sub>1</sub>,k<sub>1</sub>≤2,0≤n<sub>2</sub>,k<sub>2</sub>≤N‑1,j是复数单位,Δt是时域的分辨率,Δu<sub>1</sub>=2π|sinα<sub>1</sub>|/(3Δt)和Δu<sub>2</sub>=2π|sinα<sub>2</sub>|/(NΔt)分别是行和列的分数域分辨率;步骤8,接收端下变频后,接收到的包含高斯噪声基带信号为r=(r<sub>0</sub>,r<sub>1</sub>,…,r<sub>N‑1</sub>),其中r的第n<sub>1</sub>行n<sub>2</sub>列元素表示为r(n<sub>1</sub>,n<sub>2</sub>),r<sub>n</sub>=s<sub>n</sub>+n<sub>n</sub>,n=0,1,…,N‑1,n<sub>n</sub>=(n<sub>n</sub>,n<sub>n</sub>,n<sub>n</sub>)<sup>T</sup>是三维加性高斯白噪声;r经过二维分数阶傅里叶变换进行解调后,得到的信号为R=(R<sub>0</sub>,R<sub>1</sub>,…,R<sub>N‑1</sub>),其中R<sub>k</sub>=R(k<sub>1</sub>,k<sub>2</sub>)的表达式如下:<maths num="0003"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>R</mi><mrow><mo>(</mo><mrow><msub><mi>k</mi><mn>1</mn></msub><mo>,</mo><msub><mi>k</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow><mo>=</mo><msqrt><mfrac><mrow><mn>1</mn><mo>-</mo><mi>j</mi><mi> </mi><mi>cot</mi><mrow><mo>(</mo><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow></mrow><mn>3</mn></mfrac></msqrt><mo>&CenterDot;</mo><msqrt><mfrac><mrow><mn>1</mn><mo>-</mo><mi>j</mi><mi> </mi><mi>cot</mi><mrow><mo>(</mo><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow></mrow><mi>N</mi></mfrac></msqrt><mo>&times;</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mfrac><mrow><mi>j</mi><mi> </mi><mi>cot</mi><mrow><mo>(</mo><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><msubsup><mi>k</mi><mn>1</mn><mn>2</mn></msubsup><mo>&CenterDot;</mo><msup><mi>&Delta;t</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow><mo>&times;</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mfrac><mrow><mi>j</mi><mrow><mo>(</mo><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><msubsup><mi>k</mi><mn>2</mn><mn>2</mn></msubsup><mo>&CenterDot;</mo><msup><mi>&Delta;t</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>&times;</mo><munderover><mi>&Sigma;</mi><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>=</mo><mn>0</mn></mrow><mn>2</mn></munderover><munderover><mi>&Sigma;</mi><mrow><msub><mi>n</mi><mn>2</mn></msub><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>r</mi><mrow><mo>(</mo><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>,</mo><msub><mi>n</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mfrac><mrow><mi>j</mi><mi> </mi><mi>cot</mi><mrow><mo>(</mo><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>1</mn></msub></mrow><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><msubsup><mi>n</mi><mn>1</mn><mn>2</mn></msubsup><mo>&CenterDot;</mo><msup><mi>&Delta;t</mi><mn>2</mn></msup><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><msub><mi>&pi;n</mi><mn>1</mn></msub><msub><mi>k</mi><mn>1</mn></msub></mrow><mn>3</mn></mfrac></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>-</mo><mfrac><mrow><mi>j</mi><mi> </mi><mi>cot</mi><mrow><mo>(</mo><mrow><mo>-</mo><msub><mi>&alpha;</mi><mn>2</mn></msub></mrow><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><msubsup><mi>n</mi><mn>2</mn><mn>2</mn></msubsup><mo>&CenterDot;</mo><msup><mi>&Delta;t</mi><mn>2</mn></msup><mo>-</mo><mi>j</mi><mfrac><mrow><mn>2</mn><msub><mi>&pi;n</mi><mn>2</mn></msub><msub><mi>k</mi><mn>2</mn></msub></mrow><mi>N</mi></mfrac></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000961177900000023.GIF" wi="1974" he="318" /></maths>步骤9,求取解调后得到的信号R<sub>k</sub>与M点三维星座图S<sub>0</sub>,S<sub>1</sub>,…,S<sub>M‑1</sub>中各点的最小距离,选择距离最近的点作为要恢复的信号,最后将恢复的信号逆映射回原始3N长的二进制序列,则完成了整个二维通信过程。
地址 100081 北京市海淀区中关村南大街5号