发明名称 一种提高解码转发协作系统保密性能的波束成形优化方法
摘要 一种提高解码转发协作系统保密性能的波束成形优化方法,属于通信技术领域,由以下通信系统来实现,该系统包括信源节点、中继节点、信宿节点和窃听节点。本发明方法通过将目标问题转化为一个凸优化问题和一个一维线性搜索问题的组合能够以较高的精度得出单个功率限制下基于DF方式工作的协作通信系统的保密容量及相对应的最优波束成形矢量。本发明的意义在于相比较于由信源直接发射信号到信宿的直接传输系统,在总发射功率相同的情况下,使用该方法进行波束成形优化后的DF协作通信系统的物理层保密能力有着显著的优势。
申请公布号 CN102611488B 申请公布日期 2014.04.30
申请号 CN201210090515.6 申请日期 2012.03.30
申请人 山东大学 发明人 马丕明;马晖
分类号 H04B7/06(2006.01)I;H04W12/00(2009.01)I 主分类号 H04B7/06(2006.01)I
代理机构 济南金迪知识产权代理有限公司 37219 代理人 许德山
主权项 1.一种提高解码转发协作系统保密性能的波束成形优化方法,由以下通信系统来实现,该系统包括信源节点、中继节点、信宿节点和窃听节点,所有的节点都是单天线的节点,预先设整个通信系统工作在DF即解码转发方式下,在通信的第一阶段,信源节点将信号发送至各个中继节点,在第二个阶段,所有的中继节点对接收到的有污损的信号进行解码,还原出信源节点发送的原始信号,而后相互协作,将原始信号发送至信宿节点,窃听节点在这个过程中对信号进行窃听,系统中信源节点、信宿节点与窃听节点都只有一个,参与协作通信的中继节点的个数是已知的,记作N;整个系统的信道是准静态的,平坦衰落的无记忆信道,记信宿信道h=[h<sub>0</sub>,....,h<sub>N-1</sub>]<sup>H</sup>其中h<sub>i</sub>∈C表示第i个中继节点与信宿节点间的信道衰落因数,窃听信道g=[g<sub>0</sub>,.....,g<sub>N-1</sub>]<sup>H</sup>其中g<sub>i</sub>∈C表示第i个中继节点与窃听节点间的信道衰落因数;所有的噪声都是平稳的高斯白噪声,信宿节点处的噪声功率为<img file="FDA0000460649570000013.GIF" wi="84" he="69" />窃听节点处的噪声功率为<img file="FDA0000460649570000014.GIF" wi="85" he="69" />将波束成形矢量记做w=[w<sub>0</sub>,....,w<sub>N-1</sub>]<sup>T</sup>,其中<img file="FDA0000460649570000015.GIF" wi="350" he="71" />为第i个中继节点处的波束成形因子,该方法步骤如下:1)利用信道估计,估计出<img file="FDA0000460649570000016.GIF" wi="170" he="73" />g,h的值,设定每个中继节点的功率限制pi的值和求解精度ε,该求解精度指的是保密容量的真实值与求解值之间差的绝对值,计算出信宿信道的自相关矩阵R<sub>h</sub>=hh<sup>H</sup>,窃听信道的自相关矩阵R<sub>g</sub>=gg<sup>H</sup>;2) 初始化t<sub>min</sub>,t<sub>max</sub>:<img file="FDA0000460649570000017.GIF" wi="179" he="73" />t<sub>max</sub>的初始化通过求解优化问题(1)进行,记(1)的最优值为o<sub>pt</sub>,<maths num="0001"><![CDATA[<math><mrow><msub><mi>t</mi><mi>max</mi></msub><mo>=</mo><msub><msup><mi>O</mi><mn>2</mn></msup><mi>pt</mi></msub><mo>+</mo><msubsup><mi>&sigma;</mi><mi>E</mi><mn>2</mn></msubsup></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><munder><mrow><mi>max</mi><mi>imize</mi></mrow><mi>w</mi></munder><mo>:</mo><mi>Re</mi><mrow><mo>(</mo><msup><mi>w</mi><mi>H</mi></msup><mi>g</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>subject to</mi><mo>:</mo><mi></mi><msup><mrow><mo>|</mo><msub><mi>w</mi><mi>i</mi></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>&le;</mo><msub><mi>p</mi><mi>i</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>t<sub>max</sub>、t<sub>min</sub>分别为<img file="FDA00004606495700000110.GIF" wi="302" he="77" />的上限和下限,其中w<sup>*</sup>为欲求解的最优波束成形矢量,(1)问题是标准的凸优化问题,可以通过内点法进行求解,其中maximize代表优化w使其后的目标式的值最大化,subject to代表变量w受到其后限制条件的约束;3)计算t<sub>left</sub>=t<sub>max</sub>-0.618(t<sub>max</sub>-t<sub>min</sub>),t<sub>right</sub>=t<sub>min</sub>+0.618(t<sub>max</sub>-t<sub>min</sub>);4)将t<sub>left</sub>代入<img file="FDA00004606495700000111.GIF" wi="440" he="76" />求解R<sub>1</sub>(t<sub>left</sub>),同理,计算R<sub>1</sub>(t<sub>right</sub>),其中f<sub>1</sub>(t)为问题(2)的最优值的平方;<maths num="0003"><![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><munder><mrow><mi>max</mi><mi>imize</mi><mo>:</mo></mrow><mi>w</mi></munder><mi>Re</mi><mrow><mo>(</mo><msup><mi>w</mi><mi>H</mi></msup><mi>h</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>subject to</mi><mo>:</mo><msubsup><mi>&sigma;</mi><mi>E</mi><mn>2</mn></msubsup><mo>+</mo><msup><mi>w</mi><mi>H</mi></msup><msub><mi>R</mi><mi>g</mi></msub><mi>w</mi><mo>&le;</mo><mi>t</mi></mtd></mtr><mtr><mtd><msup><mrow><mo>|</mo><msub><mi>w</mi><mi>i</mi></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>&le;</mo><msub><mi>p</mi><mi>i</mi></msub><mo>,</mo><mo>&ForAll;</mo><mi>i</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>随着t的变化(2)式的最优值的平方会不断的变化,因而(2)式的最优值的平方可以被记为一个关于t的函数,即f<sub>1</sub>(t),同时(2)问题同样是一个标准的凸优化问题;5)若R<sub>1</sub>(t<sub>left</sub>)&lt;R<sub>1</sub>(t<sub>right</sub>)则对t<sub>min</sub>进行赋值t<sub>min</sub>=t<sub>left</sub>,若R<sub>1</sub>(t<sub>left</sub>)&gt;R<sub>1</sub>(t<sub>right</sub>)则对t<sub>max</sub>进行赋值t<sub>max</sub>=t<sub>right</sub>,若R<sub>1</sub>(t<sub>left</sub>)=R<sub>1</sub>(t<sub>right</sub>)则对t<sub>min</sub>和t<sub>max</sub>进行赋值t<sub>min</sub>=t<sub>left</sub>,t<sub>max</sub>=t<sub>right</sub>;6)计算t<sub>max</sub>-t<sub>min</sub>,若t<sub>max</sub>-t<sub>min</sub>>ε则跳转至步骤3),否则继续步骤7);7)计算t<sub>1</sub>=(t<sub>max</sub>+t<sub>min</sub>)/2;8)计算通信系统的保密容量<img file="FDA0000460649570000021.GIF" wi="746" he="78" />并通过将t<sub>1</sub>带入(2)求解相对应的最优波束成形矢量w<sup>*</sup>;9)输出通信系统的保密容量C<sub>S</sub>和最优波束成形矢量w<sup>*</sup>。
地址 250100 山东省济南市历城区山大南路27号