发明名称 复数神经网络信道预测方法
摘要 本发明公开了一种复数神经网络信道预测方法,主要解决MIMO系统中由于信道时变而导致的信道衰落问题。其技术方案是,1.基站对信道进行测量得到含有估计误差的信道系数训练序列;2.根据得到的信道系数序列得到相应的训练样本和期望输出;3.输入训练样本进行复数小波神经网络训练,得到最终的网络权值;4.基站利用训练后的复数小波神经网络进行信道系数预测。本方法简便易行,效果良好,适用于减少由信道时变对MIMO系统信道的影响。
申请公布号 CN105142177A 申请公布日期 2015.12.09
申请号 CN201510473741.6 申请日期 2015.08.05
申请人 西安电子科技大学 发明人 刘祖军;李兴旺;孙德春;汪嘉曦
分类号 H04W24/06(2009.01)I;G06N3/02(2006.01)I 主分类号 H04W24/06(2009.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种复数神经网络信道预测方法,包括如下步骤:(1)在正式通信前基站BS对信道进行测量,得到含有估计误差的信道系数训练序列h(1),h(2),h(r),…,h(N),r为1到N的整数,N为训练序列总个数;(2)初始化复数神经网络的输入层节点个数p=4、隐含层节点个数m=10、输出层节点个数q=1,学习概率η=0.01,动量因子α=0.935,隐含层节点j与输出层节点k之间的网络权值w<sub>jk</sub>为0到1之间随机复数,输入层节点i与隐含层节点j之间的网络权值v<sub>ij</sub>为0到1之间随机复数,其中i为1到p的整数,j为1到m的整数,k为1到q的整数,神经元函数Ψ<sub>C</sub>(t)采用复数变化后的Morlet小波函数,初始化小波伸缩平移参数{a,b}为0到1之间随机数;(3)根据信道系数训练序列,分别得到训练样本集<img file="FDA0000775602180000013.GIF" wi="200" he="68" />和期望输出集<img file="FDA0000775602180000014.GIF" wi="228" he="68" />x(n)={h(n),h(n+1),…,h(n+p‑1)},d(n)=h(n+p),其中,n为1到w的整数,w为训练样本总个数;(4)输入训练样本x(n),利用神经元函数Ψ<sub>C</sub>(t)分别计算隐含层节点j的输出y<sub>j</sub>(n)和输出层节点k的输出O<sub>k</sub>(n):<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>y</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&psi;</mi><mi>C</mi></msub><mrow><mo>(</mo><mrow><msub><mi>s</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>&psi;</mi><mi>C</mi></msub><mrow><mo>(</mo><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><msub><mi>v</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msub><mi>&psi;</mi><mi>C</mi></msub><mrow><mo>(</mo><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover><mrow><mo>&lsqb;</mo><mrow><mrow><mo>(</mo><mrow><msub><msup><mi>v</mi><mi>R</mi></msup><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><msup><mi>x</mi><mi>R</mi></msup><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><msub><msup><mi>v</mi><mi>I</mi></msup><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><msup><mi>x</mi><mi>I</mi></msup><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>+</mo><mi>i</mi><mrow><mo>(</mo><mrow><msub><msup><mi>v</mi><mi>R</mi></msup><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><msup><mi>x</mi><mi>I</mi></msup><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><msub><msup><mi>v</mi><mi>I</mi></msup><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><msup><mi>x</mi><mi>R</mi></msup><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow><mo>&rsqb;</mo></mrow></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000775602180000011.GIF" wi="1652" he="320" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>O</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&psi;</mi><mi>C</mi></msub><mrow><mo>(</mo><mrow><msub><mi>L</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>&psi;</mi><mi>C</mi></msub><mrow><mo>(</mo><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>w</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><mi>y</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msub><mi>&psi;</mi><mi>C</mi></msub><mrow><mo>(</mo><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><mrow><mo>&lsqb;</mo><mrow><mrow><mo>(</mo><mrow><msub><msup><mi>w</mi><mi>R</mi></msup><mrow><mi>j</mi><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><msup><mi>y</mi><mi>R</mi></msup><mi>j</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><msub><msup><mi>w</mi><mi>I</mi></msup><mrow><mi>j</mi><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><msup><mi>y</mi><mi>I</mi></msup><mi>j</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mo>+</mo><mi>i</mi><mrow><mo>(</mo><mrow><msub><msup><mi>w</mi><mi>R</mi></msup><mrow><mi>j</mi><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><msup><mi>y</mi><mi>I</mi></msup><mi>j</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><msub><msup><mi>w</mi><mi>I</mi></msup><mrow><mi>j</mi><mi>k</mi></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><msub><msup><mi>y</mi><mi>R</mi></msup><mi>j</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow><mo>&rsqb;</mo></mrow></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000775602180000012.GIF" wi="1743" he="328" /></maths>其中,ψ<sub>C</sub>(s<sub>j</sub>(n))为s<sub>j</sub>(n)输入到神经元函数得到的输出,ψ<sub>C</sub>(L<sub>k</sub>(n))为L<sub>k</sub>(n)输入到神经元函数得到的输出,s<sub>j</sub>(n)为隐含层节点j的输入,L<sub>k</sub>(n)为输出层节点k的输入,v<sub>ij</sub>(n),w<sub>jk</sub>(n)为输入训练样本x(n)时的网络权值,参数上标I表示复数虚数部分,R表示复数实数部分,x<sub>i</sub>(n)为输入层第i节点的输入信道系数;(5)根据输出层节点k的输出O<sub>k</sub>(n)和期望输出d<sub>k</sub>(n)得到误差信号:e<sub>k</sub>(n)=d<sub>k</sub>(n)‑O<sub>k</sub>(n);(6)根据误差信号得到误差能量和ε(n),并将该误差能量和ε(n)作为网络训练的代价函数:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&epsiv;</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>q</mi></munderover><mo>&lsqb;</mo><msup><mrow><mo>(</mo><msup><msub><mi>d</mi><mi>k</mi></msub><mi>R</mi></msup><mo>(</mo><mi>n</mi><mo>)</mo><mo>-</mo><msup><msub><mi>O</mi><mi>k</mi></msub><mi>R</mi></msup><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msup><msub><mi>d</mi><mi>k</mi></msub><mi>I</mi></msup><mo>(</mo><mi>n</mi><mo>)</mo><mo>-</mo><msup><msub><mi>O</mi><mi>k</mi></msub><mi>I</mi></msup><mo>(</mo><mi>n</mi><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA0000775602180000021.GIF" wi="1131" he="148" /></maths>(7)根据得到的代价函数构建隐含层节点j与输出层节点k之间的网络权值调整公式Δw<sub>jk</sub>(n)和输入层节点i与隐含层节点j之间的网络权值调整公式Δv<sub>ij</sub>(n);(8)根据得到的权值调整公式分别对隐含层节点j与输出层节点k之间的网络权值w<sub>jk</sub>(n)和输入层节点i与隐含层节点j之间的网络权值v<sub>ij</sub>(n)进行修改,修改完成后输入训练样本x(n+1),返回步骤(4),直到所有w个训练样本输入完毕,得到最终修改后的隐含层节点j与输出层节点k之间的网络权值w<sub>jk</sub>′(n)和输入层节点i与隐含层节点j之间的网络权值v<sub>ij</sub>′(n);(9)在<img file="FDA0000775602180000022.GIF" wi="50" he="59" />时刻,输入信道系数序列<img file="FDA0000775602180000023.GIF" wi="854" he="86" />根据最终修改后的隐含层节点j与输出层节点k之间的网络权值w<sub>jk</sub>′(n)和输入层节点i与隐含层节点j之间的网络权值v<sub>ij</sub>′(n),得到<img file="FDA0000775602180000024.GIF" wi="51" he="58" />时刻后的q个时刻的信道系数<img file="FDA0000775602180000025.GIF" wi="215" he="83" /><img file="FDA0000775602180000026.GIF" wi="581" he="146" />其中,<img file="FDA0000775602180000027.GIF" wi="131" he="78" />为隐含层节点j的输出,<img file="FDA0000775602180000028.GIF" wi="605" he="162" /><img file="FDA0000775602180000029.GIF" wi="120" he="74" />为输入层第i节点的输入信道系数。
地址 710071 陕西省西安市太白南路2号