# Block Least Mean Squares Algorithm over Distributed Wireless Sensor Network

Authors: T. Panigrahi; P. M. Pradhan; G. Panda; B. Mulgrew
Journal: Journal of Computer Networks and Communications
 (2012)

---

## Abstract

In a distributed parameter estimation problem, during each sampling instant, a typical sensor node communicates its estimate either by the diffusion algorithm or by the incremental algorithm. Both these conventional distributed algorithms involve significant communication overheads and, consequently, defeat the basic purpose of wireless sensor networks. In the present paper, we therefore propose two new distributed algorithms, namely, block diffusion least mean square (BDLMS) and block incremental least mean square (BILMS) by extending the concept of block adaptive filtering techniques to the distributed adaptation scenario. The performance analysis of the proposed BDLMS and BILMS algorithms has been carried out and found to have similar performances to those offered by conventional diffusion LMS and incremental LMS algorithms, respectively. The convergence analyses of the proposed algorithms obtained from the simulation study are also found to be in agreement with the theoretical analysis. The remarkable and interesting aspect of the proposed block-based algorithms is that their communication overheads per node and latencies are less than those of the conventional algorithms by a factor as high as the block size used in the algorithms.


---

## Body

## 1. Introduction

A wireless sensor network (WSN) consists of a group of sensors nodes which perform distributed sensing by coordinating themselves through wireless links. Since the nodes operate in a WSN function with limited battery power, it is important to design the networks with a minimum of communication among the nodes to estimate the required parameter vector [1, 2]. In the literature, a number of research papers have appeared which address the energy issues of sensor networks. According to the energy estimation scheme based on the 4th power loss model with Rayleigh fading [3], the transmission of 1 kb of data over a distance of 100 m, operating at 1 GHz using BPSK modulation with 10-6 bit-error rate, requires 3 J of energy. The same energy can be used for executing 300 M instructions in a 100 MIPS/watt general purpose processor. Therefore, it is of great importance to minimize the communication among nodes by maximizing local estimation in each sensor node.Each node in a WSN collects noisy observations related to certain desired parameters. In the centralized solution, every node in the network transmits its data to a central fusion center (FC) for processing. This approach has the disadvantage of being nonrobust to the failure of the FC and also needs a powerful central processor. Again the problem with centralized processing is the lack of scalability and the requirement for a large communication resource [1]. If the intended application and the sensor architecture allow more local processing, then it would be more energy efficient compared to communication extensive centralized processing. Alternatively, each node in the network can function as an individual adaptive filter to estimate the parameter from the local observations and by cooperating with the neighbors. So there is a need to search for new distributed adaptive algorithms to reduce communication overhead for low-power consumption and low-latency systems for real-time operation.The performance of distributed algorithms depends on the mode of cooperation among the nodes, for example, incremental [4, 5], diffusion [6], probabilistic diffusion [7], and diffusion with adaptive combiner [8]. To improve the robustness against the spatial variation of signal-to-noise ratio (SNR) over the network, recently an efficient adaptive combination strategy has been proposed [8]. Also a fully distributed and adaptive implementation to make individual decisions by each node in the network is dealt with in [9].Since in block filtering technique [10], the filter coefficients are adjusted once for each new block of data in contrast to once for each new input sample in the least mean square (LMS) algorithm, the block adaptive filter permits faster implementation while maintaining equivalent performance as that of widely used LMS adaptive filter. Therefore, the block LMS algorithms could be used at each node in order to reduce the amount of communications.With this in mind, we present a block formulation of the existing cooperative algorithm [4, 11] based on the distributed protocols. Distinctively, in this paper, the adaptive mechanism is proposed in which the nodes of the same neighborhood communicate with each other after processing a block of data, instead of communicating the estimates to the neighbors after every sample of input data. As a result, the average bandwidth for communication among the neighboring nodes decreases by a factor equal to the block size of the algorithm. In real-time scenarios, the nodes in the sensor network follow a particular protocol for communication [12–14], where the communication time is much more than the processing time. The proposed block distributed algorithm provides an excellent balance between the message transmission delay and processing delay, by increasing the interval between two messages and by increasing the computational load on each node in the interval between two successive transmissions. The main motivation here is to propose communication-efficient block distributed LMS algorithms (both incremental and diffusion type). We analyze the performance of the proposed algorithms and compare them with existing distributed LMS algorithms.The reminder of the paper is organized as follows. In Section2, we present the BDLMS algorithm and its network global model. The performance analysis of BDLMS and its learning characteristics obtained from a simulation study are presented in Section 3. Performance analysis of the BILMS and its simulation results are presented in Section 4. The performance of the proposed algorithms in terms of communication cost and latency is compared with the conventional distributed adaptive algorithms in Section 5. Finally, Section 6 discusses the conclusions of the paper.


## 2. Block Adaptive Distributed Solution

Consider a sensor network withN number of sensor nodes randomly distributed over the region of interest. The topology of a sensor network is modeled by an undirected graph. Let 𝒢 be an undirected graph defined by a set of nodes 𝒱 and a set of edges ℰ. Nodes i and j are called neighbors if the are connected by an edgey that is, (i,j)∈ℰ. We also considered a loop which consists of a set of nodes i1,i2,…,iN such that the node ik is ik+1’s neighbor, k=1,2,…,N, and i1 is iN’s neighbor. Every node in the network i∈𝒱 is associated with noisy output di to the input data vector ui. We have assumed that the noise is independent of both input and output data; therefore, the observations are spatially and temporally independent. The neighborhood of node i is defined as the set of nodes connected to node i which is defined as 𝒩i=j∣(i,j)∈ℰ [15].Now, the objective is to estimate anM×1 unknown vector w∘ from the measurements of N nodes. In order estimate this, every node is modeled as a block adaptive linear filter where each node updates its weights using the set of errors observed in the estimated output vector, and broadcasts that to its neighbors. The estimated weight vector of the kth node at time n is denoted as ŵk(n). Let uk(n) be the input data of kth node at time instant n, then the input vector to the filter at time instant n is(1)uk(n)=[uk(n),uk(n-1),…,uk(n-M+1)]T.
The corresponding desired output of the node for the input vector uk(n) is modeled as [16, 17]
(2)dk(n)=ukT(n)w°+υk(n),
where υk(n) denotes a temporally and spatially uncorrelated white noise with variance συ,k2.Theblock indexj is related to the time index n as
(3)n=jL+i,i=0,1,2,…,L-1,j=1,2,…,
where L is the block length. The jth block contains time indices n=[jL,jL+1,…,jL+L-1]. Combining input vectors of kth node for block j to form a matrix given by
(4)Xkj=[uk(jL),uk(jL+1),…,uk(jL+L-1)]T,
the corresponding desired response at jth block index of kth node is represented as
(5)dkj=[dk(jL),dk(jL+1),…,dk(jL+L-1)]T.
Let ekj represent the L×1 error signal vector for jth block of kth node and is defined as
(6)ekj=dk(j)-Xkjŵkj,
where ŵkj estimated weight vector of the filter when jth block of the data is input at the kth node and of the order of M×1.The regression input data and corresponding desired responses are distributed across all the nodes and are represented in two global matrices:(7)Xbgj=col{X1j,X2j,…,XNj},dbgj=col{d1j,d2j,…,dNj}.
The objective is to estimate the M×1 vector w from the above quantities, those collected the data across N nodes. By using this global data, the block error vector for the whole network is
(8)ebgj=dbgj-Xbgjŵ.
Now, the vector ŵ can be estimated by minimizing MSE function as
(9)min⁡ŵE‖dbg-Xbgŵ‖2.
The time index is dropped here for simple mathematical representation. Since the quantities are collected data across the network in block format; therefore, the block mean square error (BMSE) is to be minimized. The BMSE is given by [17, 18]
(10)BMSE=1L[E[dbgTdbg]-E[dbgTXbg]ŵ-ŵTE[XbgTdbg]-ŵTE[XbgTXbg]ŵ].
Let the input regression data u be Gaussian and defined by the correlation function r(l)=σ2α|l| in the covariance matrix, where α is the correlation index and σ2 is the variance of the input regression data, then the relation between correlation and cross-correlation quantities among blocked and unblocked data can be denoted as [10]
(11)RXbg=LRUg,RdXbg=LRdug,Rdbg=LRdg,
where RXbg=E[XbgTXbg], RdXbg=E[XbgTdbg], and Rdbg=E[dbgTdbg], which are the autocorrelation and cross-correlation matrices for global data in blocked form. Similarly, the correlation matrices for unblocked data are defined as RUg=E[UgTUg], Rdug=E[UgTdg], and Rdg=E[dgTdg] where the global distribution of data across the network is represented as Ug=[u1,u2,…,uN]T and dg=[d1,d2,…,dN]T. These relations are also valid for node data in individual nodes.Now, the block mean square error (BMSE) in (10) is reduced to
(12)BMSE=1L[LRdg-LRdugŵ-LŵTRdug-LwTRugŵ]=Rdg-Rdugŵ-ŵTRdug-wTRugŵ=MSE.
Comparing (12) with the MSE of conventional LMS for global data [17, 19], it can be concluded that the MSE in both the cases is same. Hence, block LMS algorithm has similar properties as that of the conventional LMS algorithm. Now, (9) for blocked data can be reduced to a form similar to that of unblocked data as
(13)min⁡ŵE‖dg-Ugŵ‖2.
The basic difference between blocked and unblocked LMS lies in the estimation of the gradient vector used in their respective implementation. The block LMS algorithm uses a more accurately estimated gradient because of the time averaging. The accuracy increases with the increase in block size. Taking into account the advantages of block LMS over conventional LMS, the distributed block LMS is proposed here.


### 2.1. Adaptive Block Distributed Algorithms

In adaptive block LMS algorithm, each nodek in the network receives the estimates from its neighboring nodes after each block of input data to adapt the local changes in the environment. Two different types of distributed LMS in WSN have been reported in literature, namely, incremental and diffusion LMS [6, 19]. These algorithms are based on conventional LMS for local learning process which in terms needs large communication resources. In order to achieve the same performance with less communication resource, the block distributed LMS is proposed here.


#### 2.1.1. The Block Incremental LMS (BILMS) Algorithm

In an incremental mode of cooperation, information flows in a sequential manner from one node to the adjacent one in the network after processing one sample of data [4]. The communications in the incremental way of cooperation can be reduced if each node need to communicate only after processing a block of data. For any block of data j, it is assumed that node k has access to the ŵk-1j estimates from its predecessor node, as defined by the network topology and constitution. Based on these assumptions, the proposed block incremental LMS algorithm can be stated by reducing the conventional incremental LMS algorithm ((16) in [19]) to a blocked data form as follows,
(14)ŵ0j=ŵj-1,ŵkj=ŵk-1j+μkL∑q=0L-1uk(jL+q)×(dk(jL+q)-ukT(jL+q)ŵk-1j)=ŵk-1j+μkLXjkT(dkj-Xkjŵk-1j),fork=1,2,…N,ŵj=ŵNj,
where μk is the local step size, and L is the block size.


#### 2.1.2. The Block Diffusion LMS (BDLMS) Algorithm

Here, each nodek updated its estimate by using a simple local rule based on the average of its own estimates plus the information received from its neighbor 𝒩k. In this case, for every jth block of data at the kth node, the node has access to a set of estimates from its neighbors 𝒩k. Similar to block incremental LMS, the proposed block diffusion strategy for a set of local combiners ckl and for local step size μk can be described as a reduced form of conventional diffusion LMS [6, 20] as
(15)θkj-1=∑l∈Nk,j-1cklŵkj-1,θk(-1)=0,ŵkj=θkj-1+μkL∑q=0L-1uk(jL+q)×(dk(jL+q)-ukT(jL+q)θkj-1).
The weight update equation can be rewritten in more compact form by using the data in block format given in (4) and (5) as
(16)ŵkj=θkj-1+μkLXkjT(dkj-Xkjθkj-1).
Comparing (15) with (19) in [21], it is concluded that the weight update equation is modified into block format.


## 2.1. Adaptive Block Distributed Algorithms

In adaptive block LMS algorithm, each nodek in the network receives the estimates from its neighboring nodes after each block of input data to adapt the local changes in the environment. Two different types of distributed LMS in WSN have been reported in literature, namely, incremental and diffusion LMS [6, 19]. These algorithms are based on conventional LMS for local learning process which in terms needs large communication resources. In order to achieve the same performance with less communication resource, the block distributed LMS is proposed here.


### 2.1.1. The Block Incremental LMS (BILMS) Algorithm

In an incremental mode of cooperation, information flows in a sequential manner from one node to the adjacent one in the network after processing one sample of data [4]. The communications in the incremental way of cooperation can be reduced if each node need to communicate only after processing a block of data. For any block of data j, it is assumed that node k has access to the ŵk-1j estimates from its predecessor node, as defined by the network topology and constitution. Based on these assumptions, the proposed block incremental LMS algorithm can be stated by reducing the conventional incremental LMS algorithm ((16) in [19]) to a blocked data form as follows,
(14)ŵ0j=ŵj-1,ŵkj=ŵk-1j+μkL∑q=0L-1uk(jL+q)×(dk(jL+q)-ukT(jL+q)ŵk-1j)=ŵk-1j+μkLXjkT(dkj-Xkjŵk-1j),fork=1,2,…N,ŵj=ŵNj,
where μk is the local step size, and L is the block size.


### 2.1.2. The Block Diffusion LMS (BDLMS) Algorithm

Here, each nodek updated its estimate by using a simple local rule based on the average of its own estimates plus the information received from its neighbor 𝒩k. In this case, for every jth block of data at the kth node, the node has access to a set of estimates from its neighbors 𝒩k. Similar to block incremental LMS, the proposed block diffusion strategy for a set of local combiners ckl and for local step size μk can be described as a reduced form of conventional diffusion LMS [6, 20] as
(15)θkj-1=∑l∈Nk,j-1cklŵkj-1,θk(-1)=0,ŵkj=θkj-1+μkL∑q=0L-1uk(jL+q)×(dk(jL+q)-ukT(jL+q)θkj-1).
The weight update equation can be rewritten in more compact form by using the data in block format given in (4) and (5) as
(16)ŵkj=θkj-1+μkLXkjT(dkj-Xkjθkj-1).
Comparing (15) with (19) in [21], it is concluded that the weight update equation is modified into block format.


## 2.1.1. The Block Incremental LMS (BILMS) Algorithm

In an incremental mode of cooperation, information flows in a sequential manner from one node to the adjacent one in the network after processing one sample of data [4]. The communications in the incremental way of cooperation can be reduced if each node need to communicate only after processing a block of data. For any block of data j, it is assumed that node k has access to the ŵk-1j estimates from its predecessor node, as defined by the network topology and constitution. Based on these assumptions, the proposed block incremental LMS algorithm can be stated by reducing the conventional incremental LMS algorithm ((16) in [19]) to a blocked data form as follows,
(14)ŵ0j=ŵj-1,ŵkj=ŵk-1j+μkL∑q=0L-1uk(jL+q)×(dk(jL+q)-ukT(jL+q)ŵk-1j)=ŵk-1j+μkLXjkT(dkj-Xkjŵk-1j),fork=1,2,…N,ŵj=ŵNj,
where μk is the local step size, and L is the block size.


## 2.1.2. The Block Diffusion LMS (BDLMS) Algorithm

Here, each nodek updated its estimate by using a simple local rule based on the average of its own estimates plus the information received from its neighbor 𝒩k. In this case, for every jth block of data at the kth node, the node has access to a set of estimates from its neighbors 𝒩k. Similar to block incremental LMS, the proposed block diffusion strategy for a set of local combiners ckl and for local step size μk can be described as a reduced form of conventional diffusion LMS [6, 20] as
(15)θkj-1=∑l∈Nk,j-1cklŵkj-1,θk(-1)=0,ŵkj=θkj-1+μkL∑q=0L-1uk(jL+q)×(dk(jL+q)-ukT(jL+q)θkj-1).
The weight update equation can be rewritten in more compact form by using the data in block format given in (4) and (5) as
(16)ŵkj=θkj-1+μkLXkjT(dkj-Xkjθkj-1).
Comparing (15) with (19) in [21], it is concluded that the weight update equation is modified into block format.


## 3. Performance Analysis of BDLMS Algorithm

The performance of an adaptive filter is evaluated in terms of its transient and steady-state behaviors, which, respectively provide the information about how fast and how well a filter is capable to learn. Such performance analysis is usually challenging in interconnected network because each nodek is influenced by local data with local statistics {Rdx,k,RX,k}, by its neighborhood nodes through local diffusion, and by local noise with variance συ,k2. In case of block distributed system, the analysis becomes more challenging as it has to handle data in block form. The key performance metrics used in the analysis are MSD (mean square deviation), EMSE (excess mean square error), and MSE for local and also for global networks and are defined as
(17)ηkj=E‖w̃k-1j‖2(MSD),ζkj=E‖ea,kj‖2(EMSE),ξkj=E‖ek(j)‖2=ζkj+συ,k2(MSE),
and the local error signals such as weight error vector and a priori error at kth node for jth block are given as
(18)w̃k-1j=w°-ŵk-1j,ea,kj=ukjw̃k-1j.
The algorithm described in (15) is looking like the interconnection of block adaptive filters instead of conventional LMS adaptive algorithm among all the nodes across the network. As shown in (12) that the block LMS algorithm has similar properties to those of the conventional LMS algorithm, the convergence analysis of the proposed block diffusion LMS algorithm can be carried out similar to the diffusion LMS algorithm described in [18, 21].The estimated weight vector forjth block across the network is defined as
(19)ŵj=[ŵ1j;…;ŵNj].
Let C be the N×N metropolis with entries [ckl], then the global transaction combiner matrix G is defined as G=C⊗IM. The diffusion global vector for jth block is defined as
(20)θj=Gŵj.
Now, the input data vector at jth block is defined as
(21)Xj=diag⁡{X1j,…,XNj}.
The desired block responses at each node k are assumed which have to obey the traditional data model used in literature [16–18], that is,
(22)dkj=Xkjw°+vkj,
where vkj is the background noise vector of length L. The noise is assumed to be spatially and temporarily independent with variance συ,k2. Using blocked desired response for single node (17), the global response for kth block can be modeled as
(23)dbgj=Xbgjwg°+vj,
where wg° is the optimum global weight vector defined for every node and is written as wg∘=[w∘;,…,;w∘] and
(24)vj=[v1j;,…,;vNj](LN×1)
is the additive Gaussian noise for jth block index.Using the relations defined above, the block diffusion strategy in (15) can be written in global form as
(25)ŵj=θj-1+1LSXj(dbgj-Xjθj-1),
where the step sizes for all the nodes are embedded in a matrixS,(26)S=diag⁡{μ1IM,μ2IM,…,μNIM}(NM×NM).
Using (20), it can be written as
(27)ŵj=Gŵj-1+1LSXj(dbgj-XjGŵj-1).


### 3.1. Mean Transient Analysis

The mean behavior of the proposed BDLMS is similar to diffusion LMS given in [18, 21]. The mean error vector signal is given as
(28)E[w̃j]=(INM-1LSRX)GE[w̃j-1],
where RX=diag{RX,1,RX,2,…,RX,N} is a block diagonal matrix and
(29)RX,k=E[XkTXk]=LE[UkTUk]=LRU.
Hence, (28) can be written as
(30)E[w̃j]=(INM-SRU)GE[w̃j-1].
Comparing (30) with that of diffusion LMS ((35) in [21]), we can find that both block diffusion LMS and diffusion LMS yield the same characteristic equation for the convergence of mean; and it can be concluded that block diffusion protocol defined in (15) has the same stabilizing effect on the network as diffusion LMS,


### 3.2. Mean-Square Transient Analysis

The variance estimate is a key performance indicator in mean-square transient analysis of any adaptive system. The variance relation for block data is similar to that of conventional diffusion LMS(31)E‖w̃j‖Σ2=E‖w̃j-1‖Σ′2+1L2E[vjTXjSΣSXjTvj],(32)Σ′=GTΣG-1LGTΣSE[XjTXj]G-1LGTE[XjTXj]SΣG+1L2GTE[XjTXj]SΣSE[XjTXj]G.
Using E[XjTXj]=LE[UjTUj] from the definition in (32), we obtain(33)Σ′=GTΣG-GTΣSE[UjTUj]G-GTE[U(j)TUj]SΣG+GTE[UjTUj]SΣSE[UjTUj]G,
which is similar to (45) in [21]. Using the properties of expectation and trace [18], the second term of (31) is solved as(34)1L2EvjTXjSΣSXjTvj=1L2E[tr⁡[XjSΣSXjT]]E[vjTvj]=E[njTUjSΣSUjTnj],
where the noise variance vector nj is not in block form, and it is assumed that the noise is stationary Gaussian. Equations (31) and (32) may therefore be written as
(35)E‖w̃j‖Σ2=E‖w̃j-1‖Σ′2+1L2EnjTUjSΣSUjTn(j),(36)Σ′=GTΣG-GTΣSE(UjTUj)G-GTE(UjTUj)SΣG+GTE(UjTUj)SΣSE(UjTUj)G.
It may be noted that variance estimate (36) for BDLMS algorithm is exactly the same as that of DLMS [21]. In the block LMS algorithm, the local step size is chosen to be L times that of the local step size of diffusion LMS in order to have the same level of performance. As the proposed algorithm and the diffusion LMS algorithm have similar properties, the evolution of their variances is also similar. Therefore, the recursion equation of the global variances for BDLMS will be similar to (73) and (74) in [21]. Similarly, the local node performances will be similar to (89) and (91) of [21].


### 3.3. Learning Behavior of BDLMS Algorithm

The learning behavior of BDLMS algorithm is examined using simulations. The characteristic or variance curves are plotted for block LMS and are compared with that of DLMS. The row regressors with shift invariance input [18] are used with each regressor having data as
(37)uk(i)=[uk(i),uk(i-1),…,uk(i-M+1)]T.
In block LMS, the regressors for L=3 and M=3 are given as
(38)Xk(1)=(uk(1)00uk(2)uk(1)0uk(3)uk(2)uk(1)),Xk(2)=(uk(4)uk(3)uk(2)uk(5)uk(4)uk(3)uk(6)uk(5)uk(4)).
The desired data are generated according to the model given in literature [18]. The unknown vector w° is set to [1,1,…,1]T/M.The input sequence{uk(i)} is assumed to be spatially correlated and is generated as
(39)uk(i)=ak⋅uk(i-1)+bk⋅nk(i),i>-∞.
Here, ak∈[0,1) is the correlation index, and nk(i) is a spatially independent white Gaussian process with unit variance and bk=σu,k2·(1-ak2). The regressors power profile is given by {σu,k2}∈(0,1]. The resulting regressors have Toeplitz covariance with corelation sequence rk(i)=σu,k2·(ak)|i|,i=0,1,2,…,M-1.Figure1 shows an eight-node network topology used in the simulation study. The network settings are given in Figures 2(a) and 2(b).Figure 1
Network topology used for block diffusion LMS.Network statistics used for the simulation of BDLMS. (a) Network corelation index per node. (b) Regressor power profile.
(a)(b)


### 3.4. The Simulation Conditions

The algorithm is valid for any block of length greater than one [10], while L=M is the most preferable and optimum choice.The background noise is assumed to be Gaussian white noise of varianceσυ,k2=10-3, and the data used in the study is generated using dk(n)=uk(n)w∘+vk(n). In order to generate the performance curves, 50 independent experiments are performed and averaged. The results are obtained by averaging the last 50 samples of the corresponding learning curves. The global MSD curve is shown in Figure 3. This is obtained by averaging E∥w̃kj-1∥2 across all the nodes over 100 experiments. Similarly, the global EMSE curve obtained by averaging E∥eja,k∥2, where eja,k=xkjw̃kj-1, across all the nodes over 100 experiments is displayed in Figure 4. The global MSE is depicted in Figure 5. It shows that in both the cases the MSE is exactly matching.Figure 3
Global mean-square deviation (MSD) curve for diffusion and block diffusion LMS.Figure 4
Global excess mean-square deviation (EMSE) curve for diffusion and block diffusion LMS.Figure 5
Global mean-square deviation (MSE) curve for diffusion and block diffusion LMS.Since the weights are updated and then communicated for local diffusion after everyL data samples, the number of communications between neighbors is reduced by L times compared to that of the diffusion LMS case where the weights are updated and communicated after each sample of data.The global performances are the contributions of all individual nodes, and it is obtained by taking the mean performance of all the nodes. The simulation results are provided to compare with that obtained by diffusion LMS for individual node. The local MSD evolution at node 1 is given in Figure6(a) and at node 5 is given in Figure 6(b). Similarly, the local EMSE evolution at nodes 1 and 7 is depicted in Figure 7. The convergence speed is nearly the same in both MSD and EMSE evolution, but the performance is slightly degraded in case of BDLMS. The loss of performance in case of BDLMS could be traded for the huge reduction in of communication bandwidth.Local mean-square deviation (MSD) comparison between block diffusion LMS and diffusion LMS. (a) MSD curve at node 1. (b) MSD curve at node 7.
(a)(b)Figure 7
Local EMSE at node 7 for the same network.


## 3.1. Mean Transient Analysis

The mean behavior of the proposed BDLMS is similar to diffusion LMS given in [18, 21]. The mean error vector signal is given as
(28)E[w̃j]=(INM-1LSRX)GE[w̃j-1],
where RX=diag{RX,1,RX,2,…,RX,N} is a block diagonal matrix and
(29)RX,k=E[XkTXk]=LE[UkTUk]=LRU.
Hence, (28) can be written as
(30)E[w̃j]=(INM-SRU)GE[w̃j-1].
Comparing (30) with that of diffusion LMS ((35) in [21]), we can find that both block diffusion LMS and diffusion LMS yield the same characteristic equation for the convergence of mean; and it can be concluded that block diffusion protocol defined in (15) has the same stabilizing effect on the network as diffusion LMS,


## 3.2. Mean-Square Transient Analysis

The variance estimate is a key performance indicator in mean-square transient analysis of any adaptive system. The variance relation for block data is similar to that of conventional diffusion LMS(31)E‖w̃j‖Σ2=E‖w̃j-1‖Σ′2+1L2E[vjTXjSΣSXjTvj],(32)Σ′=GTΣG-1LGTΣSE[XjTXj]G-1LGTE[XjTXj]SΣG+1L2GTE[XjTXj]SΣSE[XjTXj]G.
Using E[XjTXj]=LE[UjTUj] from the definition in (32), we obtain(33)Σ′=GTΣG-GTΣSE[UjTUj]G-GTE[U(j)TUj]SΣG+GTE[UjTUj]SΣSE[UjTUj]G,
which is similar to (45) in [21]. Using the properties of expectation and trace [18], the second term of (31) is solved as(34)1L2EvjTXjSΣSXjTvj=1L2E[tr⁡[XjSΣSXjT]]E[vjTvj]=E[njTUjSΣSUjTnj],
where the noise variance vector nj is not in block form, and it is assumed that the noise is stationary Gaussian. Equations (31) and (32) may therefore be written as
(35)E‖w̃j‖Σ2=E‖w̃j-1‖Σ′2+1L2EnjTUjSΣSUjTn(j),(36)Σ′=GTΣG-GTΣSE(UjTUj)G-GTE(UjTUj)SΣG+GTE(UjTUj)SΣSE(UjTUj)G.
It may be noted that variance estimate (36) for BDLMS algorithm is exactly the same as that of DLMS [21]. In the block LMS algorithm, the local step size is chosen to be L times that of the local step size of diffusion LMS in order to have the same level of performance. As the proposed algorithm and the diffusion LMS algorithm have similar properties, the evolution of their variances is also similar. Therefore, the recursion equation of the global variances for BDLMS will be similar to (73) and (74) in [21]. Similarly, the local node performances will be similar to (89) and (91) of [21].


## 3.3. Learning Behavior of BDLMS Algorithm

The learning behavior of BDLMS algorithm is examined using simulations. The characteristic or variance curves are plotted for block LMS and are compared with that of DLMS. The row regressors with shift invariance input [18] are used with each regressor having data as
(37)uk(i)=[uk(i),uk(i-1),…,uk(i-M+1)]T.
In block LMS, the regressors for L=3 and M=3 are given as
(38)Xk(1)=(uk(1)00uk(2)uk(1)0uk(3)uk(2)uk(1)),Xk(2)=(uk(4)uk(3)uk(2)uk(5)uk(4)uk(3)uk(6)uk(5)uk(4)).
The desired data are generated according to the model given in literature [18]. The unknown vector w° is set to [1,1,…,1]T/M.The input sequence{uk(i)} is assumed to be spatially correlated and is generated as
(39)uk(i)=ak⋅uk(i-1)+bk⋅nk(i),i>-∞.
Here, ak∈[0,1) is the correlation index, and nk(i) is a spatially independent white Gaussian process with unit variance and bk=σu,k2·(1-ak2). The regressors power profile is given by {σu,k2}∈(0,1]. The resulting regressors have Toeplitz covariance with corelation sequence rk(i)=σu,k2·(ak)|i|,i=0,1,2,…,M-1.Figure1 shows an eight-node network topology used in the simulation study. The network settings are given in Figures 2(a) and 2(b).Figure 1
Network topology used for block diffusion LMS.Network statistics used for the simulation of BDLMS. (a) Network corelation index per node. (b) Regressor power profile.
(a)(b)


## 3.4. The Simulation Conditions

The algorithm is valid for any block of length greater than one [10], while L=M is the most preferable and optimum choice.The background noise is assumed to be Gaussian white noise of varianceσυ,k2=10-3, and the data used in the study is generated using dk(n)=uk(n)w∘+vk(n). In order to generate the performance curves, 50 independent experiments are performed and averaged. The results are obtained by averaging the last 50 samples of the corresponding learning curves. The global MSD curve is shown in Figure 3. This is obtained by averaging E∥w̃kj-1∥2 across all the nodes over 100 experiments. Similarly, the global EMSE curve obtained by averaging E∥eja,k∥2, where eja,k=xkjw̃kj-1, across all the nodes over 100 experiments is displayed in Figure 4. The global MSE is depicted in Figure 5. It shows that in both the cases the MSE is exactly matching.Figure 3
Global mean-square deviation (MSD) curve for diffusion and block diffusion LMS.Figure 4
Global excess mean-square deviation (EMSE) curve for diffusion and block diffusion LMS.Figure 5
Global mean-square deviation (MSE) curve for diffusion and block diffusion LMS.Since the weights are updated and then communicated for local diffusion after everyL data samples, the number of communications between neighbors is reduced by L times compared to that of the diffusion LMS case where the weights are updated and communicated after each sample of data.The global performances are the contributions of all individual nodes, and it is obtained by taking the mean performance of all the nodes. The simulation results are provided to compare with that obtained by diffusion LMS for individual node. The local MSD evolution at node 1 is given in Figure6(a) and at node 5 is given in Figure 6(b). Similarly, the local EMSE evolution at nodes 1 and 7 is depicted in Figure 7. The convergence speed is nearly the same in both MSD and EMSE evolution, but the performance is slightly degraded in case of BDLMS. The loss of performance in case of BDLMS could be traded for the huge reduction in of communication bandwidth.Local mean-square deviation (MSD) comparison between block diffusion LMS and diffusion LMS. (a) MSD curve at node 1. (b) MSD curve at node 7.
(a)(b)Figure 7
Local EMSE at node 7 for the same network.


## 4. Performance Analysis of BILMS Algorithm

To show that the BILMS algorithm has guaranteed convergence, we may follow the steady-state performance analysis of the algorithm using the same data model as the one which is commonly used in the conventional sequential adaptive algorithms [5, 22, 23].The weight-energy relation is derived by using the definition of weighteda priori and a posteriori error [18]
(40)‖w̃kj‖Σ2+|eja,kΣ|2‖Xkj‖Σ2=‖w̃k-1j‖Σ2+|ejp,kΣ|2‖Xkj‖Σ2.
Since (40) is similar to that of (35) in [19]. Thus, the performance of BILMS is similar to that of ILMS. The variance expression is obtained from the energy relation (40) by replacing a posteriori error by its equivalent expression and then averaging both the sides
(41)E[‖w̃kj‖Σ2]=E[‖w̃k-1j‖Σ′2]+|μkL|2E[VjkTXkjΣXjkTVkj]Σ′=Σ-μkL(ΣXjkTXkj+XjkTXkjΣ)+|μkL|2XjkTXkjΣXjkTXkj.
The variance relation in (41) is similar to the variance relation of ILMS in [19]. The performance of ILMS is studied in detail in literature. It is observed that the theoretical performance of block incremental LMS and conventional incremental LMS algorithms are similar because both have the same variance expressions. Simulation results provide the validation of this analysis.


### 4.1. Simulation Results of BILMS Algorithm

For the simulation study of IBLMS, we have used the regressors with shift-invariance as with the same desired data used in the case of BDLMS algorithm. The time-correlated sequences are generated at every node according to the network statistics. The same network has been chosen here for simulation study as defined for block diffusion network in Section3.3. In incremental way of cooperation, each node receives information from its previous node, updates it by using own data, and sends the updated estimate to the next node. The ring topology used here is shown in Figure 8. We assume that the background noise to be temporarily and spatially uncorrelated additive white Gaussian noise with variance 10-3. The learning curves are obtained by averaging the performance of 100 independent experiments, generated by 5,000 samples in the network. It can be observed from figures that the steady-state performances at different nodes of the network achieved by BILMS matche very closely with that of ILMS algorithm. The EMSE plots which are more sensitive to local statistics are depicted in Figures 9(a) and 9(b). A good match between BILMS and ILMS is observed from these plots. In [19], the authors have already proved the theoretical matching of steady-state nodal performance with simulation results. As the MSE roughly reflects the noise power and the plot indicates the good performance of the adaptive network, it may be inferred that the adaptive node performs well in the steady state.Figure 8
Network topology.Network nodal performance. (a) Mean-square error (MSE) versus node. (b) excess mean-square error (EMSE) versus node. (c) Mean-square deviation (MSD) versus node.
(a)(b)(c)The global MSD curve shown in Figure10 is obtained by averaging E∥ψ̃k(j-1)∥2 across all the nodes and over 50 experiments. Similarly, the global EMSE and MSE plots are displayed in Figures 11 and 12, respectively. These are obtained by averaging E∥ea,k(j)∥2, where ea,k(j)=xk,jψ̃k(j-1) across all the nodes over 50 experiments.Figure 10
Global mean-square deviation (MSD) curve for incremental LMS and block incremental LMS.Figure 11
Global excess mean-square deviation (EMSE) curve for incremental LMS and block incremental LMS.Figure 12
Global mean-square deviation (MSE) curve for incremental LMS and block incremental LMS.If the weights are updated afterL data points and then communicated for local diffusion, the number of communications between neighbors is reduced by L times that of ILMS where the weights are updated after processing each sample of data. Therefore, similar to BDLMS, the communication overhead in BILMS also gets reduced by L times that of ILMS algorithm.The performance comparison between two proposed algorithms BDLMS and BILMS for the same network is shown in Figures13–15. One can observe from Figure 13 that the MSE for BILMS algorithm is converging faster than BDLMS. Since the same noise model is used for both the algorithms, therefore after convergence, the steady-state performances are the same for both of them. But in case of MSD and EMSE performances in Figures 14 and 15, little difference is observed. It is due the different cooperation scheme used for different algorithms. However, the diffusion cooperation scheme is more adaptive to the environmental change compared to the incremental cooperation. But a higher number of communication overhead are required for BDLMS than BILMS algorithm.Figure 13
Global mean-square error curve for BILMS and BDLMS.Figure 14
Global mean-square deviation curve for BILMS and BDLMS.Figure 15
Global excess mean-square error curve for BILMS and BDLMS.


## 4.1. Simulation Results of BILMS Algorithm

For the simulation study of IBLMS, we have used the regressors with shift-invariance as with the same desired data used in the case of BDLMS algorithm. The time-correlated sequences are generated at every node according to the network statistics. The same network has been chosen here for simulation study as defined for block diffusion network in Section3.3. In incremental way of cooperation, each node receives information from its previous node, updates it by using own data, and sends the updated estimate to the next node. The ring topology used here is shown in Figure 8. We assume that the background noise to be temporarily and spatially uncorrelated additive white Gaussian noise with variance 10-3. The learning curves are obtained by averaging the performance of 100 independent experiments, generated by 5,000 samples in the network. It can be observed from figures that the steady-state performances at different nodes of the network achieved by BILMS matche very closely with that of ILMS algorithm. The EMSE plots which are more sensitive to local statistics are depicted in Figures 9(a) and 9(b). A good match between BILMS and ILMS is observed from these plots. In [19], the authors have already proved the theoretical matching of steady-state nodal performance with simulation results. As the MSE roughly reflects the noise power and the plot indicates the good performance of the adaptive network, it may be inferred that the adaptive node performs well in the steady state.Figure 8
Network topology.Network nodal performance. (a) Mean-square error (MSE) versus node. (b) excess mean-square error (EMSE) versus node. (c) Mean-square deviation (MSD) versus node.
(a)(b)(c)The global MSD curve shown in Figure10 is obtained by averaging E∥ψ̃k(j-1)∥2 across all the nodes and over 50 experiments. Similarly, the global EMSE and MSE plots are displayed in Figures 11 and 12, respectively. These are obtained by averaging E∥ea,k(j)∥2, where ea,k(j)=xk,jψ̃k(j-1) across all the nodes over 50 experiments.Figure 10
Global mean-square deviation (MSD) curve for incremental LMS and block incremental LMS.Figure 11
Global excess mean-square deviation (EMSE) curve for incremental LMS and block incremental LMS.Figure 12
Global mean-square deviation (MSE) curve for incremental LMS and block incremental LMS.If the weights are updated afterL data points and then communicated for local diffusion, the number of communications between neighbors is reduced by L times that of ILMS where the weights are updated after processing each sample of data. Therefore, similar to BDLMS, the communication overhead in BILMS also gets reduced by L times that of ILMS algorithm.The performance comparison between two proposed algorithms BDLMS and BILMS for the same network is shown in Figures13–15. One can observe from Figure 13 that the MSE for BILMS algorithm is converging faster than BDLMS. Since the same noise model is used for both the algorithms, therefore after convergence, the steady-state performances are the same for both of them. But in case of MSD and EMSE performances in Figures 14 and 15, little difference is observed. It is due the different cooperation scheme used for different algorithms. However, the diffusion cooperation scheme is more adaptive to the environmental change compared to the incremental cooperation. But a higher number of communication overhead are required for BDLMS than BILMS algorithm.Figure 13
Global mean-square error curve for BILMS and BDLMS.Figure 14
Global mean-square deviation curve for BILMS and BDLMS.Figure 15
Global excess mean-square error curve for BILMS and BDLMS.


## 5. Performance Comparison

In this section, we present an analysis of communication cost and latency to have a theoretical comparison of the performances of distributed LMS with block distributed LMS.


### 5.1. Analysis of Communication Cost

Assuming that the messages are of fixed bit width, the communication cost is modeled as the number of messages transmitted to achieve the steady-state value in the network. LetN be the number of nodes in the network, and let M be the filter length. The block length L is chosen to be the same as the filter length. Let h be the average time required for the transmission of one message, that is, for one communication between the nodes [24–26].


#### 5.1.1. ILMS and BILMS Algorithms

In the incremental mode of cooperation, every node sends its own estimated weight vector to its adjacent node in a unidirectional cyclic manner. Since at any instant of time, only one node is active/allowed to transmit to only one designated node, the number of messages transmitted in one complete cycle isN-1. Let K be the number of cycles required to attain the steady-state value in the network. Therefore, the total number of communications required to converge the system to steady state is given by
(42)CILMS=(N-1)K.
In case of BILMS also, at any instant of time, only one node in the network is active/allowed to transmit to one designated follower node, as in the case of ILMS. But, in case of BILMS, each node sends its estimated weight vector to its follower node in the network after an interval of L sample periods after processing a block of L data samples. Therefore, the number of messages sent by a node in this case is reduced to K/L, and accordingly, the total communication cost is given by
(43)CBILMS=(N-1)KL.


#### 5.1.2. DLMS and BDLMS Algorithms

The diffusion-based algorithms are communication intensive. In DLMS mode of cooperation, in each cycle, each node in the network sends its estimated information to all its connected nodes in the network. So the total number of messages transmitted by all the nodes in a cycle is(44)c=∑i=1Nni,
where ni is the number of nodes connected to the ith node, and the total communication cost to attain convergence is given by
(45)CDLMS=cK.
In this proposed block diffusion strategy, the number of connected nodes ni and the total size of the messages remain the same as that of DLMS. But, in case of BDLMS algorithm, each node distributes the message after L data samples. Therefore the communication is reduced by a factor equal to the block length, and the total communication cost in this case is given by
(46)CBDLMS=cKL.


### 5.2. Analysis of Duration for Convergence

The time interval between the arrival of input to a node and the time of reception of corresponding updates by the designated node(s) may be assumed to be comprised of two major components. Those are processing delay to perform the necessary computations in a node to obtain the estimates to be updated and the communication delay involved in transferring the message to the receiver node(s). The processing delay will very much depend on the hardware architecture of the nodes to perform the computation which could be widely varying. But, without losing much of the generality of analysis, we can assume that each node hasM parallel multipliers and one full adder to implement the LMS algorithm. Let TM and TA be the time required for executing a multiplication and an addition, respectively. Therefore, the processing delay needed for single update in LMS is
(47)D=2TM+(M+1)TA.
The communication delay is mostly due to the implementation of protocols for transmission and reception, which remains almost the same for different nodes. The location of nodes will not have any major contribution to the delay unless the destination node is far apart, and a relay node is required to make the message reach the destination. In this backdrop, we can assume that the same average delay h is required to transfer each message for all receiver-transmitter pairs in the network.


#### 5.2.1. Estimation of Delays for the ILMS and BILMS Algorithms

In case of ILMS, the duration of each updating cycle by all the nodes is(48)ND+(N-1)h,
and the total duration for convergence of the network is given as
(49)LILMS=[ND+(N-1)h]K.
If the same hardware as that of ILMS is used for the implementation of BILMS, the delay for processing one block of data is 2MTM+M(M+1)TA=MD. Then the duration of one cycle of update by the block incremental LMS is N{2MTM+M(M+1)TA}+(N-1)h, and the duration of convergence of this algorithm is
(50)LBILMS=[NMD+(N-1)h]KL.
For L=M, the above expression could be reduced to
(51)LBILMS=[ND+(N-1)hL]K.
Comparing (51) with (49), we can find that in BILMS the processing delay remains the same as that in ILMS, but the communication overhead is reduced by L times.


#### 5.2.2. Estimation of Delays for the DLMS and BDLMS Algorithms

Similar to ILMS, it is also assumed here that the updates of a node reaches all the connected nodes after the same average delayh. Therefore, the communication delay remains the same as that of ILMS, but in this case, it needs more processing delay to process the unbiased estimates received from the connected neighboring nodes. The total communication delay in a cycle in this case can be given by cTA+NTM, where c is the total number of messages transferred in a cycle given by (44). Now, the total duration of a cycle in diffusion LMS with the same hardware constraints is given by
(52)LDLMS=[cTA+NTM+ND+(N-1)h]K.
In case of DBLMS, the total communication delay per cycle is reduced by a factor of L, which can be expressed as
(53)LBDLMS=[cTA+NTM+NMD+(N-1)h]KL.
The mathematical expressions of communication cost and latency for the distributed LMS and the block distributed LMS algorithms are summarized in Table 1. A numerical example is given in Table 2 to show the advantage of block-distributed algorithms over the sequential-distributed algorithms. The authors have simulated the hardware for 8-bit multiplication and addition in TSMC 90 nm. The multiplication and addition time are found to be TA=10-5ns,TM=10-3ns. We assume the transmission delay h=10-2s. Looking at the convergence curves obtained from the simulation studies, we can say that the network attains steady state after 250-input data in DLMS and 50-input data in ILMS case. The filter length M as well as the block size L are taken to be 10 in the numerical study.Table 1
Comparison of performances of distributed algorithms and block distributed algorithms.
ParameterILMSBILMSDLMSBDLMSCommunication Cost(N-1)K(N-1)K/LcKcK/LDuration of convergence[ND+(N-1)h]K[NMD+(N-1)h]K/L[cTA+NTM+ND+(N-1)h]K[cTA+NTM+ND+(N-1)h]K/LTable 2
Numerical comparison of performances of sequential and block distributed adaptive algorithms.
ParameterILMSBILMSDLMSBDLMSCommunication cost950954500450Duration of convergence9.5s0.95s9.5s0.95s


## 5.1. Analysis of Communication Cost

Assuming that the messages are of fixed bit width, the communication cost is modeled as the number of messages transmitted to achieve the steady-state value in the network. LetN be the number of nodes in the network, and let M be the filter length. The block length L is chosen to be the same as the filter length. Let h be the average time required for the transmission of one message, that is, for one communication between the nodes [24–26].


### 5.1.1. ILMS and BILMS Algorithms

In the incremental mode of cooperation, every node sends its own estimated weight vector to its adjacent node in a unidirectional cyclic manner. Since at any instant of time, only one node is active/allowed to transmit to only one designated node, the number of messages transmitted in one complete cycle isN-1. Let K be the number of cycles required to attain the steady-state value in the network. Therefore, the total number of communications required to converge the system to steady state is given by
(42)CILMS=(N-1)K.
In case of BILMS also, at any instant of time, only one node in the network is active/allowed to transmit to one designated follower node, as in the case of ILMS. But, in case of BILMS, each node sends its estimated weight vector to its follower node in the network after an interval of L sample periods after processing a block of L data samples. Therefore, the number of messages sent by a node in this case is reduced to K/L, and accordingly, the total communication cost is given by
(43)CBILMS=(N-1)KL.


### 5.1.2. DLMS and BDLMS Algorithms

The diffusion-based algorithms are communication intensive. In DLMS mode of cooperation, in each cycle, each node in the network sends its estimated information to all its connected nodes in the network. So the total number of messages transmitted by all the nodes in a cycle is(44)c=∑i=1Nni,
where ni is the number of nodes connected to the ith node, and the total communication cost to attain convergence is given by
(45)CDLMS=cK.
In this proposed block diffusion strategy, the number of connected nodes ni and the total size of the messages remain the same as that of DLMS. But, in case of BDLMS algorithm, each node distributes the message after L data samples. Therefore the communication is reduced by a factor equal to the block length, and the total communication cost in this case is given by
(46)CBDLMS=cKL.


## 5.1.1. ILMS and BILMS Algorithms

In the incremental mode of cooperation, every node sends its own estimated weight vector to its adjacent node in a unidirectional cyclic manner. Since at any instant of time, only one node is active/allowed to transmit to only one designated node, the number of messages transmitted in one complete cycle isN-1. Let K be the number of cycles required to attain the steady-state value in the network. Therefore, the total number of communications required to converge the system to steady state is given by
(42)CILMS=(N-1)K.
In case of BILMS also, at any instant of time, only one node in the network is active/allowed to transmit to one designated follower node, as in the case of ILMS. But, in case of BILMS, each node sends its estimated weight vector to its follower node in the network after an interval of L sample periods after processing a block of L data samples. Therefore, the number of messages sent by a node in this case is reduced to K/L, and accordingly, the total communication cost is given by
(43)CBILMS=(N-1)KL.


## 5.1.2. DLMS and BDLMS Algorithms

The diffusion-based algorithms are communication intensive. In DLMS mode of cooperation, in each cycle, each node in the network sends its estimated information to all its connected nodes in the network. So the total number of messages transmitted by all the nodes in a cycle is(44)c=∑i=1Nni,
where ni is the number of nodes connected to the ith node, and the total communication cost to attain convergence is given by
(45)CDLMS=cK.
In this proposed block diffusion strategy, the number of connected nodes ni and the total size of the messages remain the same as that of DLMS. But, in case of BDLMS algorithm, each node distributes the message after L data samples. Therefore the communication is reduced by a factor equal to the block length, and the total communication cost in this case is given by
(46)CBDLMS=cKL.


## 5.2. Analysis of Duration for Convergence

The time interval between the arrival of input to a node and the time of reception of corresponding updates by the designated node(s) may be assumed to be comprised of two major components. Those are processing delay to perform the necessary computations in a node to obtain the estimates to be updated and the communication delay involved in transferring the message to the receiver node(s). The processing delay will very much depend on the hardware architecture of the nodes to perform the computation which could be widely varying. But, without losing much of the generality of analysis, we can assume that each node hasM parallel multipliers and one full adder to implement the LMS algorithm. Let TM and TA be the time required for executing a multiplication and an addition, respectively. Therefore, the processing delay needed for single update in LMS is
(47)D=2TM+(M+1)TA.
The communication delay is mostly due to the implementation of protocols for transmission and reception, which remains almost the same for different nodes. The location of nodes will not have any major contribution to the delay unless the destination node is far apart, and a relay node is required to make the message reach the destination. In this backdrop, we can assume that the same average delay h is required to transfer each message for all receiver-transmitter pairs in the network.


### 5.2.1. Estimation of Delays for the ILMS and BILMS Algorithms

In case of ILMS, the duration of each updating cycle by all the nodes is(48)ND+(N-1)h,
and the total duration for convergence of the network is given as
(49)LILMS=[ND+(N-1)h]K.
If the same hardware as that of ILMS is used for the implementation of BILMS, the delay for processing one block of data is 2MTM+M(M+1)TA=MD. Then the duration of one cycle of update by the block incremental LMS is N{2MTM+M(M+1)TA}+(N-1)h, and the duration of convergence of this algorithm is
(50)LBILMS=[NMD+(N-1)h]KL.
For L=M, the above expression could be reduced to
(51)LBILMS=[ND+(N-1)hL]K.
Comparing (51) with (49), we can find that in BILMS the processing delay remains the same as that in ILMS, but the communication overhead is reduced by L times.


### 5.2.2. Estimation of Delays for the DLMS and BDLMS Algorithms

Similar to ILMS, it is also assumed here that the updates of a node reaches all the connected nodes after the same average delayh. Therefore, the communication delay remains the same as that of ILMS, but in this case, it needs more processing delay to process the unbiased estimates received from the connected neighboring nodes. The total communication delay in a cycle in this case can be given by cTA+NTM, where c is the total number of messages transferred in a cycle given by (44). Now, the total duration of a cycle in diffusion LMS with the same hardware constraints is given by
(52)LDLMS=[cTA+NTM+ND+(N-1)h]K.
In case of DBLMS, the total communication delay per cycle is reduced by a factor of L, which can be expressed as
(53)LBDLMS=[cTA+NTM+NMD+(N-1)h]KL.
The mathematical expressions of communication cost and latency for the distributed LMS and the block distributed LMS algorithms are summarized in Table 1. A numerical example is given in Table 2 to show the advantage of block-distributed algorithms over the sequential-distributed algorithms. The authors have simulated the hardware for 8-bit multiplication and addition in TSMC 90 nm. The multiplication and addition time are found to be TA=10-5ns,TM=10-3ns. We assume the transmission delay h=10-2s. Looking at the convergence curves obtained from the simulation studies, we can say that the network attains steady state after 250-input data in DLMS and 50-input data in ILMS case. The filter length M as well as the block size L are taken to be 10 in the numerical study.Table 1
Comparison of performances of distributed algorithms and block distributed algorithms.
ParameterILMSBILMSDLMSBDLMSCommunication Cost(N-1)K(N-1)K/LcKcK/LDuration of convergence[ND+(N-1)h]K[NMD+(N-1)h]K/L[cTA+NTM+ND+(N-1)h]K[cTA+NTM+ND+(N-1)h]K/LTable 2
Numerical comparison of performances of sequential and block distributed adaptive algorithms.
ParameterILMSBILMSDLMSBDLMSCommunication cost950954500450Duration of convergence9.5s0.95s9.5s0.95s


## 5.2.1. Estimation of Delays for the ILMS and BILMS Algorithms

In case of ILMS, the duration of each updating cycle by all the nodes is(48)ND+(N-1)h,
and the total duration for convergence of the network is given as
(49)LILMS=[ND+(N-1)h]K.
If the same hardware as that of ILMS is used for the implementation of BILMS, the delay for processing one block of data is 2MTM+M(M+1)TA=MD. Then the duration of one cycle of update by the block incremental LMS is N{2MTM+M(M+1)TA}+(N-1)h, and the duration of convergence of this algorithm is
(50)LBILMS=[NMD+(N-1)h]KL.
For L=M, the above expression could be reduced to
(51)LBILMS=[ND+(N-1)hL]K.
Comparing (51) with (49), we can find that in BILMS the processing delay remains the same as that in ILMS, but the communication overhead is reduced by L times.


## 5.2.2. Estimation of Delays for the DLMS and BDLMS Algorithms

Similar to ILMS, it is also assumed here that the updates of a node reaches all the connected nodes after the same average delayh. Therefore, the communication delay remains the same as that of ILMS, but in this case, it needs more processing delay to process the unbiased estimates received from the connected neighboring nodes. The total communication delay in a cycle in this case can be given by cTA+NTM, where c is the total number of messages transferred in a cycle given by (44). Now, the total duration of a cycle in diffusion LMS with the same hardware constraints is given by
(52)LDLMS=[cTA+NTM+ND+(N-1)h]K.
In case of DBLMS, the total communication delay per cycle is reduced by a factor of L, which can be expressed as
(53)LBDLMS=[cTA+NTM+NMD+(N-1)h]KL.
The mathematical expressions of communication cost and latency for the distributed LMS and the block distributed LMS algorithms are summarized in Table 1. A numerical example is given in Table 2 to show the advantage of block-distributed algorithms over the sequential-distributed algorithms. The authors have simulated the hardware for 8-bit multiplication and addition in TSMC 90 nm. The multiplication and addition time are found to be TA=10-5ns,TM=10-3ns. We assume the transmission delay h=10-2s. Looking at the convergence curves obtained from the simulation studies, we can say that the network attains steady state after 250-input data in DLMS and 50-input data in ILMS case. The filter length M as well as the block size L are taken to be 10 in the numerical study.Table 1
Comparison of performances of distributed algorithms and block distributed algorithms.
ParameterILMSBILMSDLMSBDLMSCommunication Cost(N-1)K(N-1)K/LcKcK/LDuration of convergence[ND+(N-1)h]K[NMD+(N-1)h]K/L[cTA+NTM+ND+(N-1)h]K[cTA+NTM+ND+(N-1)h]K/LTable 2
Numerical comparison of performances of sequential and block distributed adaptive algorithms.
ParameterILMSBILMSDLMSBDLMSCommunication cost950954500450Duration of convergence9.5s0.95s9.5s0.95s


## 6. Conclusion

We have proposed the block implementation of the distributed LMS algorithms for WSN. The theoretical analysis and the corresponding simulation results demonstrate that the performance of the block-distributed LMS algorithms is similar to that of the sequential-distributed LMS. The remarkable achievement of the proposed algorithms is that a node requiresL(block size) times of less communications compared to the conventional sequential-distributed LMS algorithms. This would be of great advantage in reducing the communication bandwidth and power consumption involved in the transmission and reception of messages across the resource-constrained nodes in a WSN. In the coming years, with continuing advances in microelectronics, we can accommodate enough computing resources in the nodes to reduce the processing delays in the nodes, but the communication bandwidth and communication delay could be the major operational bottlenecks in the WSNs. The proposed block formulation therefore would have further advantages over the sequential counterpart in the coming years.



---

*Source: 601287-2012-03-07.xml*