[发明专利]一种具有信噪比约束的无线传感器网络滚动时域信噪比估计方法

权利要求书

1.一种具有信噪比约束的无线传感器网络滚动时域信噪比估计方法，其特征在于，具体步骤如下：

(1)、利用信噪比公式，功率控制算法和流速控制算法，得到如下无线传感器网络状态空间模型：

γ‾(k+1)=(1-α)γ‾(k)+αγ^(k)+n(k)]]>

式中，为dB尺度的无线传感器网络节的信噪比水平，为dB尺度的无线传感器网络节点期望信噪比水平，α为功率调节参数，n(k)为随机白噪声，为流速调节参数，c为网络阻塞率，d(k)为流速增量，，H为无线传感器网络信道带宽；

(2)、以为无线传感器网络模型的状态变量，设定无线传感器网络信道容量为F，则最大信噪比为设定信号发送阈值为得到无线传感器网络的信噪比约束为得到具有信噪比约束的无线传感器网络模型：

xk+1=Axk+wkyk=Cxk+vkxk∈[γ‾min,γ‾min]]]>

式中，

A和C为权重参数矩阵，A=1-αα01-n~c,]]>C=1001;]]>

w_k为无线传感器网络的系统噪声，wk=n(k)n‾d(k);]]>

y_k为无线传感器网络的测量输出信号；

v_k为测量噪声；

(3)、设定权重矩阵∏，Q和R，设定滚动时域窗口N＝1，状态变量x_k由初始状态x_k-1和扰动噪声w_k-1决定，而x_k-1与w_k-1不相关，将无线传感器网络滚动时域信噪比估计转化为等价的不等式约束最小化问题：

minθ1Ψk(θ1)]]>

s.t.Ψk(θ1)=||xk-1-x^k-1||Π-12+||wk-1||Q-12+||yk-1-Cxk-1||R-12]]>

x_k＝Ax_k-1+w_k-1

y_k-1＝Cx_k-1+v_k-1

γ‾min≤xk≤γ‾max]]>

式中，

θ₁＝(x_k-1，w_k-1)表示决策变量；

表示欧几里得范数；

为k-1时刻的先验估计值即参考值，为k-1时刻的估计误差，

表示对估计值的信任度；

表示对扰动信号的估计；

(4)、将步骤(3)设定的不等式约束最小化问题转化为近似等价的等式约束优化问题：

minθ1L]]>

s.t.L=Ψk(θ1)-μΣi=1nlogαi-μΣi=1nlogβi-λTΣi=1n(xk-1i-αi-γ-imin)]]>

-λTΣi=1n(γ-imax-αi-βi)]]>

x_k＝Ax_k-1+w_k-1

y_k-1＝Cx_k-1+v_k-1

xk-1-α=γ‾min]]>

α+β=γ‾max]]>

式中，μ为障碍参数，λ为Lagrangian乘子，αⁱ，βⁱ，λⁱ，和分别为向量α，β，λ，x_k-1，和的第i个元素；

(5)、通过LOQO内点法求解步骤(4)中设定的等式约束优化问题，具体步骤如下：

S1-1：初始化，设定测试时间长度K，在可行域的区间范围内，任意初始化k-1时刻的变量y_k-1和序列{x_k-1，w_k-1，α，β，λ}；

S1-2：根据一阶KKT最优化条件，以和{x_k-1，w_k-1，α，β，λ}为初始迭代点，计算出估计增量Δx_k-1和Δλ：

Δxk-1=T-1(-▿xk-1Ψ-(2A+B)-1Λrλ+2μ(2A+B)-1e)]]>

Δλ=HxT-1(2μ(2A+B)-1e-(2A+B)-1Λrλ)+(I-HxT-1)▿xk-1Ψ-λ]]>

式中，μ=0.05(2α+β)Tλk,]]>T＝H_x+(2A+B)^-1Λ，Hx=▿xk-12L=Π-1+CTR-1C]]>为Hessian矩阵，矩阵A，B和Λ分别表示以αⁱ，βⁱ和λⁱ元素的对角矩阵，e为元素都为1的列向量，▿xk-1Ψ=2Π-1(xk-1-x^k-1)-2CTR-1(yk-1-Cxk-1),]]>

rλ=▿λL=xk-1-2α-β-γ‾min+γ‾max;]]>

S1-3：计算估计增量Δw_k-1，Δα和Δβ：

Δwk-1=-Hw-1rw]]>

Δα＝-0.5Λ^-1A(r_α+2Δλ)

Δβ＝-Λ^-1B(r_β+Δλ)

式中，Hw=▿wk-12L=Q-1]]>为Hessian矩阵，rw=▿wk-1L=2Q-1wk-1,]]>

rβ=▿βL=μe-BΛe;]]>

S1-4：更新估计值{x_k-1，w_k-1，α，β，λ}：

x_k-1＝x_k-1+ρΔx_k-1

w_k-1＝w_k-1+ρΔw_k-1

α＝α+ρΔα

β＝β+ρΔβ

λ＝λ+ρΔλ

式中，ρ为LOQO内点法求解器设定的搜索步长且ρ＝0.35；

S1-5：采用Matlab中的LOQO内点求解器判断估计值是否满足精度要求，如果满足精度要求，则估计值{x_k-1，w_k-1，α，β，λ}为最优估计值，然后转到S1-6；如果不满足则转到S1-2；

S1-6：基于滚动优化原理，根据最优估计值x_k-1，w_k-1计算当前k时刻的最优估计

xk*=Axk-1+wk-1]]>

S1-7：更新x^k-1=xk*;]]>

S1-8：判断终止条件，如果k＝K，结束，得到信噪比估计最优值；否则，k＝k+1，转到S1-2。