发明名称 一种基于SIM-OFDM的叠加编码调制方法
摘要 本发明涉及通信抗干扰技术领域,尤其涉及叠加编码调制技术通信系统,子载波索引调制及其相关时域检测技术。本发明借鉴SIM-OFDM的思想和空间调制SM的思想,将一串比特数据的一部分用于传输,另一部分用作索引比特,并且分别进行编码。在OFDM系统中传输的时候,根据索引比特编码后的比特数据,选择特定的子载波来传输数据。在接收端采用最大似然的方法,求出索引比特的软信息,进而对索引比特进行译码。本发明与之前的SCM相比,不但由于SCM的叠加层数变少,复杂度降低,还可以在一定程度上保持甚至提高误码率性能。同时,用SCM调制代替QAM调制,在相同的传输效率下,也能使得性能得到提高。
申请公布号 CN104113393A 申请公布日期 2014.10.22
申请号 CN201410351516.0 申请日期 2014.07.23
申请人 电子科技大学 发明人 谭佳;肖悦;王顺顺;柏慧荣;但黎琳;李少谦
分类号 H04L1/00(2006.01)I;H04L27/26(2006.01)I 主分类号 H04L1/00(2006.01)I
代理机构 成都宏顺专利代理事务所(普通合伙) 51227 代理人 李玉兴
主权项 一种基于SIM‑OFDM的叠加编码调制方法,其特征在于,包括以下步骤:S1、由数字化信源产生M<sub>1</sub>+M<sub>2</sub>×L位信息比特,经过一个串并转换器,前M<sub>1</sub>位作为索引比特,后M<sub>2</sub>×L位作为数据比特,将数据比特送入另一个串并转换器;S2、索引比特编码和交织:将S1所述M<sub>1</sub>位比特进行编码,对编码完成后的索引比特进行交织,得到索引数据{u<sub>j</sub>},其中,编码方式为Turbo码,编码码率为R,j=1,2,…,M<sub>1</sub>/R;S3、进行SCM调制:S1所述M<sub>2</sub>×L位比特经过一个串并转换器后,将数据分在L层上传输,第l层的数据为<img file="FDA0000541968090000011.GIF" wi="149" he="79" />经过扩频长度为S的扩频器,交织,BPSK调制,得到<img file="FDA0000541968090000012.GIF" wi="138" he="85" />各层分别经过加权值{ω<sup>(l)</sup>,ω<sup>(2)</sup>,…,ω<sup>(L)</sup>}进行加权,将加权后的数据进行叠加,得到<img file="FDA0000541968090000013.GIF" wi="301" he="141" />送入SIM调制器,等待发送,其中,m=1,2,…,M<sub>2</sub>,j=1,2,…,M<sub>2</sub>×S;S4、进行SIM调制:对S3得到的<img file="FDA0000541968090000014.GIF" wi="296" he="140" />采用子载波2选1的方式进行SIM调制,即S3所得M<sub>2</sub>×S位数据通过M<sub>2</sub>×S×2个子载波来传输;S5、将S4得到的M<sub>2</sub>×S×2位数据进行IFFT变换到时域<img file="FDA0000541968090000015.GIF" wi="439" he="135" />对所述<img file="FDA0000541968090000016.GIF" wi="408" he="135" />加循环前缀后,过信道加噪声,进行FFT变换到频域,去掉循环前缀得到<img file="FDA0000541968090000017.GIF" wi="462" he="141" />其中,N=2×M<sub>2</sub>×S,K代表实际有数据的子载波数,<img file="FDA0000541968090000018.GIF" wi="73" he="92" />是对应子载波的信道参数值,所述<img file="FDA0000541968090000019.GIF" wi="75" he="92" />是未知的,需要通过ML的方法来求得,<img file="FDA00005419680900000110.GIF" wi="74" he="84" />表示第j时刻第l层的比特,n<sub>j</sub>为0均值,方差为δ<sup>2</sup>=N<sub>0</sub>/2的高斯噪声;S6、ML准则检测索引比特得到信道参数,具体如下:S61、用<img file="FDA0000541968090000021.GIF" wi="77" he="84" />和<img file="FDA0000541968090000022.GIF" wi="80" he="84" />表示第j个比特为1和为0的时候,采用ML的方法,算出索引比特的似然比信息<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>L</mi><mrow><mo>(</mo><msub><mi>u</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>log</mi><mo>{</mo><munder><munder><mi>&Sigma;</mi><mrow><msup><mi>u</mi><mo>&prime;</mo></msup><mo>=</mo><msubsup><mi>u</mi><mi>j</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup></mrow></munder><mrow><msup><mi>x</mi><mo>&prime;</mo></msup><mo>&Element;</mo><mi>X</mi></mrow></munder><msup><mi>e</mi><mfrac><msup><mrow><mo>|</mo><mo>|</mo><mi>y</mi><mo>-</mo><msub><mi>h</mi><msup><mi>u</mi><mo>&prime;</mo></msup></msub><msup><mi>x</mi><mo>&prime;</mo></msup><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mrow><mo>-</mo><msup><mi>&delta;</mi><mn>2</mn></msup></mrow></mfrac></msup><mo>/</mo><munder><munder><mi>&Sigma;</mi><mrow><msup><mi>u</mi><mo>&prime;</mo></msup><mo>=</mo><msubsup><mi>u</mi><mi>j</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msubsup></mrow></munder><mrow><msup><mi>x</mi><mo>&prime;</mo></msup><mo>&Element;</mo><mi>X</mi></mrow></munder><msup><mi>e</mi><mfrac><msup><mrow><mo>|</mo><mo>|</mo><mi>y</mi><mo>-</mo><msub><mi>h</mi><msup><mi>u</mi><mo>&prime;</mo></msup></msub><msup><mi>x</mi><mo>&prime;</mo></msup><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mrow><mo>-</mo><msup><mi>&delta;</mi><mn>2</mn></msup></mrow></mfrac></msup><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA0000541968090000023.GIF" wi="1009" he="276" /></maths>其中,X表示所有发送符号的集合,h<sub>u'</sub>表示对应子载波的信道参数,所述h<sub>u'</sub>根据索引比特u'来确定;S62、根据S61所述L(u<sub>j</sub>)得到<img file="FDA0000541968090000024.GIF" wi="76" he="91" />的值,若L(u<sub>j</sub>)≥0,则所述<img file="FDA0000541968090000025.GIF" wi="67" he="84" />是第2j个子载波上的信道参数值,若L(u<sub>j</sub>)<0,则所述<img file="FDA0000541968090000026.GIF" wi="66" he="84" />是第2j‑1个子载波上的信道参数值;S63、根据S62所述,得到信道参数H={h<sub>j</sub>},其中,j=M<sub>2</sub>×S;S7、将S6得到的信道参数信息H送入SCM的检测译码,将S61所述索引比特的似然比信息L(u<sub>j</sub>)解交织后送入译码器译码,译码完成以后,得到检测完成的索引比特数据,其中,译码方式对应于发送端的编码器;S8、SCM检测,具体如下:S81、在GA模块中,采用高斯近似检测,加权系数为W={w<sup>(1)</sup>,w<sup>(2)</sup>,…,w<sup>(l)</sup>},第j时刻的接收信息为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>r</mi><mi>j</mi></msub><mo>=</mo><msub><mi>h</mi><mi>j</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msup><mi>w</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msup><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>+</mo><msub><mi>n</mi><mi>j</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0000541968090000027.GIF" wi="449" he="141" /></maths>S82、将接收信息送入GA模块,针对S3所述<img file="FDA0000541968090000028.GIF" wi="76" he="85" />可以得到<img file="FDA0000541968090000029.GIF" wi="408" he="85" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>&zeta;</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>=</mo><msub><mi>r</mi><mi>j</mi></msub><mo>-</mo><msub><mi>h</mi><mi>j</mi></msub><msup><mi>w</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msup><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>,</mo><mi>E</mi><mrow><mo>(</mo><msub><mi>r</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>h</mi><mi>j</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msup><mi>w</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msup><mi>E</mi><mrow><mo>(</mo><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>,</mo><mi>Var</mi><mrow><mo>(</mo><msub><mi>r</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>h</mi><mi>j</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msup><mi>w</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msup><mi>Var</mi><mrow><mo>(</mo><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>+</mo><msup><mi>&delta;</mi><mn>2</mn></msup><mo>,</mo></mrow>]]></math><img file="FDA00005419680900000210.GIF" wi="1637" he="142" /></maths>其中,<img file="FDA00005419680900000211.GIF" wi="81" he="84" />表示针对<img file="FDA00005419680900000212.GIF" wi="68" he="85" />的干扰;S83、初始化,<img file="FDA00005419680900000213.GIF" wi="517" he="84" />得到接收信号r<sub>j</sub>的均值和方差;S84、求得<img file="FDA00005419680900000214.GIF" wi="87" he="85" />的均值<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>E</mi><mrow><mo>(</mo><msubsup><mi>&zeta;</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>=</mo><mi>E</mi><mrow><mo>(</mo><msub><mi>r</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>h</mi><mi>j</mi></msub><msup><mi>w</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msup><mi>E</mi><mrow><mo>(</mo><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA00005419680900000215.GIF" wi="603" he="85" /></maths>和方差<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>Var</mi><mrow><mo>(</mo><msubsup><mi>&zeta;</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>=</mo><mi>Var</mi><mrow><mo>(</mo><msub><mi>r</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msub><mi>h</mi><mi>j</mi></msub><msup><mi>w</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msup><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA00005419680900000216.GIF" wi="665" he="85" /></maths>S85、根据S84所述求得第j时刻传送的第l个比特的对数似然比为<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>L</mi><mrow><mo>(</mo><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>=</mo><mi>log</mi><mrow><mo>(</mo><mfrac><mrow><mi>Pr</mi><mrow><mo>(</mo><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mo>+</mo><mn>1</mn><mo>|</mo><msub><mi>r</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><mrow><mi>Pr</mi><mrow><mo>(</mo><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mo>-</mo><mn>1</mn><mo>|</mo><msub><mi>r</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>2</mn><msub><mi>h</mi><mi>j</mi></msub><mrow><mo>(</mo><msub><mi>r</mi><mi>j</mi></msub><mo>-</mo><mi>E</mi><mrow><mo>(</mo><msubsup><mi>&zeta;</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow></mrow><mrow><mi>Var</mi><mrow><mo>(</mo><msubsup><mi>&zeta;</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000541968090000031.GIF" wi="1005" he="187" /></maths>S86、根据S85求得L层所有用户的传送数据对应的似然比信息L(c);S9、在SCM中的每一层,根据S86所述L(c)进行解交织得到L(a),用<img file="FDA0000541968090000032.GIF" wi="149" he="79" />表示第l层的第i个比特,将所述<img file="FDA0000541968090000033.GIF" wi="146" he="79" />送入解扩频模块进行解扩频,得到第l层的第1个比特的对数似然比<img file="FDA0000541968090000034.GIF" wi="477" he="147" />其中,S为发射机中的扩频长度,<img file="FDA0000541968090000035.GIF" wi="75" he="84" />为的扩频比特;S10、得到S9所述<img file="FDA0000541968090000036.GIF" wi="151" he="79" />的解交织后的比特外信息<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>Ext</mi><mrow><mo>(</mo><msubsup><mi>a</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>=</mo><msubsup><mi>s</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mi>L</mi><mrow><mo>(</mo><msubsup><mi>d</mi><mn>1</mn><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>-</mo><mi>L</mi><mrow><mo>(</mo><msubsup><mi>a</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000541968090000037.GIF" wi="627" he="85" /></maths>S11、将S10所述比特外信息经过交织,得到<img file="FDA0000541968090000038.GIF" wi="207" he="85" />将所述<img file="FDA0000541968090000039.GIF" wi="176" he="85" />再次送入GA模块,对<img file="FDA00005419680900000310.GIF" wi="72" he="85" />的均值和方差进行更新:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>E</mi><mrow><mo>(</mo><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>=</mo><mi>tanh</mi><mrow><mo>(</mo><mfrac><mrow><mi>Ext</mi><mrow><mo>(</mo><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>,</mo><mi>Var</mi><mrow><mo>(</mo><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mi>E</mi><mrow><mo>(</mo><msubsup><mi>c</mi><mi>j</mi><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA00005419680900000311.GIF" wi="1002" he="172" /></maths>S12、完成L层所有数据的均值和方差的更新,将更新后的发送数据的均值和方差代入S8‑S9,按照设置的迭代次数CH进行迭代,根据迭代进行到第CH次时得到的<img file="FDA00005419680900000312.GIF" wi="146" he="78" />可以求得<img file="FDA00005419680900000313.GIF" wi="111" he="83" />其中,<img file="FDA00005419680900000314.GIF" wi="87" he="83" />表示第k层第m个比特的估计;S13、将S12所述<img file="FDA00005419680900000315.GIF" wi="86" he="84" />和S7所述检测完成的索引比特数据送入串并转换器,得到最终输出比特。
地址 611731 四川省成都市高新区(西区)西源大道2006号