[发明专利]一种盲源分离的无线传感器网络节点选择方法

摘要

本发明提出一种盲源分离的无线传感器网络节点选择方法，基于传感器节点的观测噪声和剩余能量两个因素，首先通过数值拟合的方法确定节点平均观测噪声和盲源分离性能之间的关系；在此基础上，将节点选择问题建模为基于节点剩余能量和盲源分离性能折中的组合优化问题，在汇聚节点处采用基于凸优化的启发算法求解该优化问题以确定被选中的节点；最后，被选中的节点将采集到的混合信号传输到汇聚节点，在汇聚节点出进行盲源分离恢复出各个独立的源信号。采用该节点优化选择算法能有效提高节点剩余能量，延长无线传感器网络的生存周期，同时保证良好的信号分离性能。

主权项

一种盲源分离的无线传感器网络节点选择方法，其特征在于，该方法包含以下步骤：步骤1、确定参与盲源分离的传感器节点个数最优值k，具体步骤如下：(1)、设信号源个数为q，令k'＝q+1，定义效用函数：G=μ(Σi=1NEtoli-Σj∈SkETxj)+C;]]>其中，第一项为选中的k'个传感器节点传输数据后，WSN即Wireless sensor networks为无线传感器网络的全部节点的总剩余能量；μ为加权系数；N为WSN节点总数；为节点i传输数据之前的能量；Sk为选中的k'个节点构成的子集；为被选中的节点j(j∈Sk)传输混合信号的能耗；第二项C为k'个节点参与的信号盲分离的性能参数；(2)、给定参数μ的一个取值，采用FastICA即Fast Independent component analysis即快速独立成分分析的盲源分离方法，进行3000次独立的语音信号盲分离实验：每次实验随机选取参与盲源分离的k'个传感器节点，其对应的混叠系数用服从高斯分布的随机变量表示，每个被选中传感器节点的观测信噪比SNR是随机的，因此平均信噪比MSNR取值也是随机的；计算每次实验n对应的Gn，最后取均值作为k'为某一特定值时所对应的参数G；(3)、k'＝k'+1，判断k'是否大于N，如果是，则转(4)，否则转(2)；(4)、μ分别取100、200、400和600时，重复(2)(3)的方法，可以得到一组k'‑G关系曲线，当曲线达到峰值时对应的横坐标的值就是最优传感器节点个数k；步骤2、定义被选中传感器节点的平均观测信噪比MSNR和基于相关系数的盲源分离性能参数C之间的函数关系：C＝f(MSNR)＝a2MSNR2+a1MSNR+a0；其中，二阶多项式系数a2，a1，a0通过最小二乘拟合的方法确定，a2＝‑0.06627,a1＝1.783,a0＝85.23；步骤3、对参与盲源分离的传感器节点进行优化选择，将实现盲源分离的无线传感器网络节点优化选择建模为如下组合优化问题：ming(b1,...,bN)=μΣi=1Nbi[(BEelec+Bϵmpdi4)-Etoli]-[a2(1kΣi=1NbiSNRi)2+a1(1kΣi=1NbiSNRi)]]]>上式中，bi∈{0,1},i＝1,...N其中，bi为0‑1指示变量，当bi值为1时，表示传感器节点i被选中，将其采集到的语音混合信号发送到融合中心sink，当bi值为0时，表示传感器节点i未被选中，其采集到的混合信号将不被发送；为传感器节点i传输数据之前的能量；N是无线传感器网络节点总数；BEelec+Bεmpd4为传感器节点传输能耗模型，B为传输一帧数据包含的总比特数，Eelec＝50nJ/bit为常数，d为传输距离，参数εmp＝0.0013pJ/bit/m4；SNRi为传感器节点i的观测信噪比；k为步骤1中确定的最优节点个数；a2，a1为步骤2中确定的拟合系数；μ为加权系数；步骤4、通过基于凸优化的启发式方法求解步骤3中的组合优化问题，求解步骤如下：(1)将步骤3中的非凸约束条件bi∈{0,1},i＝1,...N松弛为凸约束条件0≤bi≤1,i＝1,...N；(2)初始化松弛后的问题的可行点0≤bi≤1，i＝1,...,N；λ2N×1＞0；v；其中，k为步骤1中确定的最优节点个数；N为WSN节点总数；v、λ2N×1为对偶变量；数μ＝20，εfeas＝10‑8，ε＝10‑6；(3)计算代理对偶间隙计算惩罚因子其中m＝2N，N为WSN中节点总数，f(b)＝[b1‑1,...,bN‑1,‑b1,...,‑bN]T，0≤bi≤1，i＝1,...N；(4)通过求解下述方程组求解牛顿步长的方程组确定当前搜索方向(Δb,Δλ,Δv)HDf(b)T1N×1-diag(λ)Df(b)-diag(f(b))011×N00ΔbΔλΔv=-rdualrcentrpri]]>其中，b＝[b1，...，bN]T，SNR＝[SNR1，...，SNRN]T，1N×1和11×N分别为全1矩阵，Df(b)＝[I；‑I]，IN×N为单位阵，rdual＝Hb+Df(b)Tλ+v,rpri＝1Tb‑k，rcent＝‑diag(λ)f(b)‑(1/t)12N；(5)smax＝min{1,min{‑λi/Δλi|Δλi＜0}}，初始化前进步长s＝0.99smax；用回溯法迭代求前进步长s，直到满足：||rt(b+sΔb,λ+sΔλ,v+sΔv||2≤(1‑αs)||rt(b,λ,v)||2，回溯法参数取值为α＝0.01,β＝0.5,rt(b,λ,v)计算公式如下：rt(b,λ,v)=Hb+Df(b)Tλ+v-diag(λ)f(b)-(1/t)12N1Tb-k]]>(6)判断是否满足迭代终止条件||rpri||2≤εfeas,||rdual||2≤εfeas,and如果不满足，返回(3)；如果满足，停止迭代，设得到的最优解为(7)对进行降序排序，将前k个置为1，后面N‑k个置为0，k为步骤1中确定的最优节点个数，对应bj＝1的节点j即为被选中的节点；步骤5、被选中的节点将观测数据发送到融合中心sink进行盲源分离。