智能優(yōu)化算法應用：基于鵜鶘算法3D無線傳感器網(wǎng)絡(WSN)覆蓋優(yōu)化 - 附代碼

摘要：本文主要介紹如何用鵜鶘算法進行3D無線傳感器網(wǎng)(WSN)覆蓋優(yōu)化。

1.無線傳感網(wǎng)絡節(jié)點模型

本文主要基于0/1模型，進行尋優(yōu)。在二維平面上傳感器節(jié)點的感知范圍是一個以節(jié)點為圓心，半徑為$R_n$的圓形區(qū)域，該圓形區(qū)域通常被稱為該節(jié)點的“感知圓盤”，$R_n$稱為傳感器節(jié)點的感知半徑，感知半徑與節(jié)點內(nèi)置傳感器件的物理特性有關(guān)，假設節(jié)點$n$的位置坐標為$(x_n,y_n,z_n)$在0-1感知模型中，對于平面上任意一點$p(x_p,y_p,z_p)$,則節(jié)點$n$監(jiān)測到區(qū)域內(nèi)點$p$的事件發(fā)生概率為：
$$P_r(n,p)=\{\begin{array}{l}1,d(n,p)\le R_n\\ 0,esle\end{array}\text{\hspace{1em}(1)}$$
其中$d(n,p)=\sqrt{(x_n?x_p{)}^2+(y_n?y_p{)}^2+(z_n?z_p{)}^2}$為點和之間的歐式距離。

2.覆蓋數(shù)學模型及分析

現(xiàn)假定目標監(jiān)測區(qū)域為二維平面，在區(qū)域$Area$上投放同型結(jié)構(gòu)傳感器節(jié)點的數(shù)目為N,每個節(jié)點的位置坐標值假設已被初始化賦值，且節(jié)點的感知半徑r。傳感器節(jié)點集則表示為：
N o d e { x 1 , . . . , x N } (2) Node\{x_1,...,x_N\} \tag{2} Node{x1?,...,xN?}(2)
其中 n o d e i = { x i , y i , z i , r } node_i=\{x_i,y_i,z_i,r\} nodei?={xi?,yi?,zi?,r},表示以節(jié)點$(x_i,y_i,z_i)$為圓心,r為監(jiān)測半徑的球，假定監(jiān)測區(qū)域$Area$被數(shù)字化離散為$m?n?l$個空間點，空間點的坐標為$(x,y,z)$，目標點與傳感器節(jié)點間的距離為：
$$d(nod{e}_i,p)=\sqrt{(x_i?x{)}^2+(y_i?y{)}^2+(z_i?z{)}^2}\text{\hspace{1em}(3)}$$
目標區(qū)域內(nèi)點被傳感器節(jié)點所覆蓋的事件定義為$c_i$。則該事件發(fā)生的概率$P{c}_i$即為點$(x,y,z)$被傳感器節(jié)點$nod{e}_i$所覆蓋的概率:
$$P_{cov}(x,y,z,nod{e}_i)=\{\begin{array}{l}1,ifd(nod{e}_i,p)\le r\\ 0,esle\end{array}\text{\hspace{1em}(4)}$$
我們將所有的傳感器節(jié)點在目標監(jiān)測環(huán)境中的區(qū)域覆蓋率$CoverRatio$定義為傳感器節(jié)點集的覆蓋面積與監(jiān)測區(qū)域的面積之比，如公式所示：
$$CoverRatio=\frac{\sum P_{cov}}{m?n?l}\text{\hspace{1em}(5)}$$
那我們的最終目標就是找到一組節(jié)點使得覆蓋率最大。

3.鵜鶘算法

鵜鶘算法原理請參考：https://blog.csdn.net/u011835903/article/details/124809854
鵜鶘算法是尋找最小值。于是適應度函數(shù)定義為未覆蓋率最小，即覆蓋率最大。如下：
$$fun=argmin(1?CoverRatio)=argmin(1?\frac{\sum P_{cov}}{m?n?l})\text{\hspace{1em}(6)}$$

4.實驗參數(shù)設定

無線傳感器覆蓋參數(shù)設定如下：

%% 設定WNS覆蓋參數(shù),
%% 默認輸入?yún)?shù)都是整數(shù)，如果想定義小數(shù)，請自行乘以系數(shù)變?yōu)檎麛?shù)再做轉(zhuǎn)換。
%% 比如范圍1*1，R=0.03可以轉(zhuǎn)換為100*100，R=3；
%區(qū)域范圍為AreaX*AreaY*AreaZ
AreaX = 100;
AreaY = 100;
AreaZ = 100;
N = 20 ;%覆蓋節(jié)點數(shù)
R = 15;%通信半徑

鵜鶘算法參數(shù)如下：

%% 設定鵜鶘優(yōu)化參數(shù)
pop=30; % 種群數(shù)量
Max_iteration=30; %設定最大迭代次數(shù)
lb = ones(1,3*N);
ub = [AreaX.*ones(1,N),AreaY.*ones(1,N),AreaZ.*ones(1,N)];
dim = 3*N;%維度為3N，N個坐標點

5.算法結(jié)果

從結(jié)果來看，覆蓋率在優(yōu)化過程中不斷上升。表明鵜鶘算法對覆蓋優(yōu)化起到了優(yōu)化的作用。

6.參考文獻

[1] 史朝亞. 基于PSO算法無線傳感器網(wǎng)絡覆蓋優(yōu)化的研究[D]. 南京理工大學.

7.MATLAB代碼