Robust Distributed Collaborative Beamforming for Wireless Sensor Networks with Channel Estimation Impairments.

We propose a new collaborative beamforming (CB) solution robust (i.e., RCB) against major channel estimation impairments over dual-hop transmissions through a wireless sensor network (WSN) of K nodes. The source first sends its signal to the WSN. Then, each node forwards its received signal after multiplying it by a properly selected beamforming weight. The latter aims to minimize the received noise power while maintaining the desired power equal to unity. These weights depend on some channel state information (CSI) parameters. Hence, they have to be estimated locally at each node, thereby, resulting in channel estimation errors that could severely hinder CB performance. Exploiting an efficient asymptotic approximation at large K, we develop alternative RCB solutions that adapt to different implementation scenarios and wireless propagation environments ranging from monochromatic (i.e., scattering-free) to polychromatic (i.e., scattered) ones. Besides, in contrast to existing techniques, our new RCB solutions are distributed (i.e., DCB) in that they do not require any information exchange among nodes, thereby dramatically improving both WSN spectral and power efficiencies. Simulation results confirm that the proposed robust DCB (RDCB) techniques are much more robust in terms of achieved signal-to-noise ratio (SNR) against channel estimation errors than best representative CB benchmarks.

1. Introduction

Collaborative beamforming (CB) stands out today as a key technique that offers tremendous capacity, coverage, and power gains [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]. Using CB, K autonomous and independent sensor nodes relay the information from a desired source to a target destination through a two-hop communication link by estimating then transmitting weighted replicas of the desired signal in the first and second time slots, respectively. The beamforming weights are designed so as to optimize an objective function while satisfying some practical constraints. Due to its numerous merits, CB has gained the attention of the research community. Ref. [2] introduced the CB concept and analyzed its performance in wireless sensor networks (WSNs). Ref. [12] evaluated the CB’s beampattern characteristics while [13] designed a technique that narrows down its mainbeam and minimizes its sidelobe effect. Refs. [14,15] proposed CB solutions that improve WSN energy efficiency and reduce its nodes collaboration time while [17,18,19,20,21,22,23,24] extended the CB applicability range to scattered environments.

Nevertheless, CB’s implementation require all other beamforming nodes’s information such as their locations. This can entail an excessively large amount of information exchange among the nodes. The optimal beamforming weights may then be efficiently approximated by another set of weights that can be computed based solely on locally available information at each individual node, thereby creating a distributed CB (DCB) solution capable of relieving the nodes from the cumbersome information exchange procedure.

Despite its advantages, CB (or DCB) inevitably suffers in practice from channel estimation errors. Indeed, the beamforming weights often depend on and, hence, require the estimation of channel state information (CSI) parameters locally at each node. Unfortunately, such a process could result in several estimation errors that may cause severe channel mismatch and, hence, dramatically hinder the CB (or DCB) performance. To overcome this shortcoming, Refs. [25,26,27,28,29,30,31,32,33,34,35,36,37,38] developed new robust CB (RCB or RDCB) techniques against such estimation errors. These techniques could be roughly divided into two categories: worst-case and stochastic. The former are designed to handle the worst-case scenario when errors reach their maximum and, hence, can be extremely inefficient when actually these errors are in most real-world conditions subject to random perturbations and/or unbounded [38,39]. The latter are more robust since their design accounts for random errors. Nevertheless, they have some drawbacks of their own. Indeed, they rely very often on iterative greedy suboptimal search approaches that explore a daunting number of potential solutions. Unfortunately, WSN nodes find their extremely limited computing and power capabilities severely burdened and quickly exhausted or depleted. Besides, their robustness very often deteriorate drastically in the presence of large channel estimation errors and, hence, become unsuitable for hostile wireless environments. More importantly, almost all existing stochastic CB techniques suffer from another major drawback: the key fact that they do not offer distributed solutions. Indeed, the weights depend on other node information which are locally unavailable. Although robust to small errors, their implementation requires in real-world operating conditions huge information exchange among all nodes. The required overwhelming data overhead could starve to "death" the very limited computing and power capabilities of WSN nodes, very often found already exhausted and depleted (cf. above), and, if not enough, could dramatically degrade their spectrum efficiency. [32] proposed a stochastic CB solution that is both robust and distributed (i.e., RDCB) where the total transmission power of the relays is minimized while being subject to an outage probabilistic quality of service (QoS) constraint. Nevertheless, the RDCB in [32] assumes that estimation errors relatively are much smaller than the channel estimates and, hence, suffers from severe performance degradation in harsh environments that characterize real-world operating conditions.

We propose a new RDCB solution robust against major channel estimation impairments, namely phase synchronization, localization, direction-of-arrival (DoA), and/or channel scatterers/coefficients estimation errors over dual-hop transmissions from a source to a destination communicating through a WSN of K nodes. In the first time slot, the source sends its signal to the WSN while, in the second, each node forwards its received signal after multiplying it by a properly selected beamforming weights. The latter aims to minimize the received noise power while maintaining the desired power equal to unity. These weights depend a priori on some CSI parameters. Hence, they have to be estimated locally at each node, thereby resulting in channel estimation errors that could severely hinder the CB performance. Exploiting an efficient asymptotic approximation at large K, we develop alternative solutions that not only account for estimation errors, but also adapt to different implementation scenarios and wireless propagation environments ranging from monochromatic (i.e., scattering-free) to polychromatic (i.e., scattered) ones. Besides, in contrast to existing techniques, our new RCB solutions are also distributed (i.e., RDCB) in that they do not require any information exchange among nodes, thereby dramatically improving the WSN spectral and power efficiencies. Simulations results confirm that the proposed RDCB techniques are much more robust in terms of achieved signal-to-noise ratio (SNR) against channel estimation errors than the nominal optimal CB solution (i.e., optimized without being aware of and, hence, accounting for impairments) and RDCB in [32] benchmarks, yet at much lower complexity, power cost, and overhead, making them suitable for WSN deployment in the harsh environments that characterize operation in real-world conditions.

The paper is organized as follows. Section 2 describes the dual-hop communication system model. The proposed monochromatic RDCB (M-RDCB) and polychromatic RDCB (P-RDCB) techniques are derived in Section 3. Section 4 analyzes theoretically the performance gains that the new M-RDCB and P-RDCB solutions could potentially achieve by integrating robustness in their designs in monochromatic and polychromatic environments. Simulation results are discussed in Section 5 and Section 6 draws out concluding remarks.

Uppercase and lowercase bold letters denote matrices and column vectors, respectively. and are the entry of a matrix and entry of a vector, respectively. The operators , , and denote the complex conjugate, the transpose, and the complex conjugate transpose or Hermitian, respectively. is the 2-norm of a vector and is the absolute value. The operator ⊙ is the element-wise product. stands for the statistical expectation and is the first-order Bessel function of the first kind.

2. System Model

As illustrated in Figure 1, the system of our interest consists of a wireless sensor network (WSN) comprised of K nodes each of which is equipped with a single isotropic antenna and uniformly and independently distributed on , the disc with center at O and radius R, a receiver , and a source S both located in the same plane containing . We assume that there is no direct link from the source to the receiver due to high pathloss attenuation [2,15,17,18,19,20,21]. Moreover, let denote the source’s polar coordinates and s its narrow-band. In this paper, we assume that the signal bandwith’s reciprocal is large with respect to the time delays of all rays. For this reason, the time notion is ignored when denoting the source signal. unit power signal. Without any loss of generality, S is assumed to be at . Let , , and denote the k-th node’s polar coordinates, backward, and forward channel, respectively. is assumed to be a zero-mean unit-variance circular Gaussian random variable. Since the WSN nodes are independent and completely autonomous, we consider here that the k-th WSN is only aware of its coordinates and both its backward and forward channels while being obvious to those of all other nodes in the network. A dual-hop communication is established from the source S to the receiver . In the first time slot, the source sends its signal s to the WSN. Let denote the received signal vector at the sensor nodes given by where and is the nodes’ noise vector. In the second time slot, the k-th node multiplies its received signal with the complex conjugate of the beamforming weight and forwards the resulting signal to the receiver. It follows from (1) that the received signal at O is where is the beamforming vector, , , and n is the receiver noise. Let and denote the received power from the source, and the aggregate noise power due to the thermal noise at the receiver and the forwarded noises from the terminals, respectively. It holds from (2) that where . Although several approaches can be adopted to properly design the beamforming weights [35], we are only concerned in this paper with minimizing the aggregate noise power while maintaining the beamforming response , and hence, the power received from the source equal to unity. The latter is nothing else but the well-known minimum variance distortionless response (MVDR) beamformer [40,41] with a relaxed distortionless response constraint. Mathematically speaking, we have to solve the following optimization problem:where denotes the ideal optimal beamforming vector. The optimization problem in (5) can be rewritten as or, equivalently as

Figure 1

System model.

The solution of the above convex optimization problem can be expressed as and, hence, the k-th node’s weight is given by

It follows from (9) that in order to implement , the k-th node must estimate both its backward and forward channels . Unfortunately, in practice, such a process results in channel estimation errors which may hinder the beamforming performance. As such, is only valid in ideal conditions where implementation impairments do not exist. In real-world conditions, is substituted by where and are the k-th estimates of the backward and forward channels, respectively. Another drawback of the nominal optimal CB solution (as referred to in the remainder of the paper) (i.e., optimized without being aware of and, hence, accounting for channel parameter estimation impairments), which must be emphasized herein, is that the k-th node must be aware of the channel estimates of all other nodes in the WSN. To this end, each node must broadcast its channel information through the network and, hence, ’s implementation requires a huge overhead. The latter might not only cause the depletion of the WSN nodes scarce energy resources and the deterioration of their spectral efficiency, but also from severely amplified performance losses due to the cumulatively increasing impact of CSI feedback errors. In what follows, we propose new RDCB techniques robust against channel parameter estimation errors both in monochromatic and polychromatic environments. We focus in this work on phase synchronization, localization, DoA, and/or channel scatterers/coefficients estimation errors.

3. Proposed Robust Distributed Collaborative Beamforming (RDCB) Techniques

In order to overcome channel estimation impairments, one should start by taking an in-depth look into the backward channel structure. The latter falls under two main categories: (i) single-ray (i.e., monochromatic) channels that ignore the scattering phenomenon to assume only a unique line-of-sight ray; and (ii) multi-ray (i.e., polychromatic) ones that account for the scattering present in most real-world environments. As far as the forward channel is concerned, we make no particular a priori assumption on its structure (i.e., whether it is monochromatic or polychromatic). And so we assume that its estimate includes an additive channel estimation error term.

3.1. Monochromatic (i.e., Scattering-Free) Environments

In such environments, can be expressed as where is the k-th node’s initial phase. Hence, the corresponding beamforming weight can be estimated from a pilot signal received at each node from the source S, assumed in this case to be far from the nodes (i.e., ) along two options: (1) either each node estimates the initial phase ; or (2) it estimates both the direction-of-arrival (DoA) and the coordinates . The first option requires the implementation of phase synchronization techniques [42,43] while the second relies both on DoA estimation and localization algorithms [44,45]. Nevertheless, both options incur estimation errors of different nature that hinder the accuracy of and, hence, the performance of CB.

3.1.1. Implementation Option 1 (Phase Synchronization)

This implementation option results in a phase jitter due to synchronization and phase offset estimation errors among nodes. Therefore, the k-th node’s backward channel estimate is given as where and is the k-th node’s phase jitter that depends on its local oscillator characteristics. We will show later that ’s performance deteriorates as increases due to channel mismatch (i.e., ) it causes. To overcome this challenging issue, we propose in this paper to anticipate the inevitable phase jitter by accounting for its impact in the CB design. Actually, one could modify the optimization problem in (7) as where , , and is the error incurred when estimating the k-th node’s backward channel without a priori knowledge of its structure (i.e., whether monochromatic or polychromatic) using a training sequence sent from the receiver. The robust proposed beamforming vector is then given by

As can be observed from (14), depends on both the actual and estimated channel values. Because nodes are of course unaware of the actual channel, we need to substitute by an equivalent quantity that depends only on known parameters. To this end, we propose to investigate the asymptotic expression for this term at large K. It follows from the definitions of , , and that By resorting to the law of large numbers, we have for large K that , , , , and . Assuming that and are independent and uniformly distributed over and , respectively, one could obtain for large K where and are the variances of and , respectively. As the number of nodes in WSNs is typically large, we can substitute (16) in (14) to obtain

A straightforward inspection of (17) reveals that is exclusively dependent on , , , and . The first and second are locally estimated by the k-th node while and depend on its local oscillator characteristics and the adopted phase synchronization technique and, hence, could be stored in its local memory before WSN deployment. Furthermore, is independent of the forward and backward channels of all other nodes. This is an important DCB feature since it avoids any information exchange among WSN nodes thereby saving their scarce energy resources and improving the WSN spectral efficiency. It is worth noting in the absence of backward channel phase estimation errors (i.e., ) and forward channel estimation errors (i.e., ) that reduces to the nominal M-DCB solution given in [18] by thereby underlining unambiguously both the challenges and the merits of implementing DCB robustness to channel parameter estimation errors.

3.1.2. Implementation Option 2 (Localization and direction of arrival (DoA) Estimation)

With option 2, each node must perform both self-localization and DoA estimation with, once again, some inevitable estimation errors that hinder the channel information accuracy. In such a case, the estimated backward channel can be written as where is the error on the radial coordinate and is the combined error on the angle coordinate and ( = 0). Adopting similar steps as in Section 3.1.1, one can prove that the proposed RDCB can be expressed in this scenario as where expectation is taken over , , , , and . As could be observed from (17), each node is able to compute its own weight using only its local information, thereby avoiding any information exchange that may dramatically deteriorate the WSN power and spectral efficiencies. However, every node needs to compute the expectation in the right hand sight (RHS) of (17), thereby burdening the proposed beamformer’s implementation complexity. In what follows, we prove owing to the adoption of a mild assumption that its possible to derive this expectation term in (20) in closed form. Assuming that and are statistically independent, we have . The pdf of can be determined as

Therefore, its average is determined by where . Please note in the second line that we resort to the variable change . We also remove the imaginary part of the equation as it is a sinus function which is odd and, hence, its integral over a zero-centered interval is null. Besides, we have

Injecting (22) and (23) in (20) yields (24) in which depends only on the coefficients of own estimated channels, , , and . Since each terminal can locally estimate its own channel, the proposed M-RDCB solution does not incur any noticeable overhead, computation, or power costs. Whereas , , and can be easily broadcast over the WSN at very negligible increase in such three cost items. It is also worth noting in absence of localization and DoA estimation errors (i.e., ) and forward channel estimation errors (i.e., ) that the proposed M-RDCB, , reduces once again to the nominal M-DCB solution given in Equation (18) above.

3.2. Polychromatic Environments

We assume here that the source is scattered by a given number of scatterers located in the same plane containing . These scatterers generate out of the transmit signal L rays or “spatial chromatics” (with reference to their angular distribution) that form a polychromatic propagation channel. The l-th ray or chromatic is characterized by its angle deviation from the source direction and its complex amplitude . In such a case, the backward channel of the k-th node is given by

It is noteworthy that (25) reduces to (11) when there is no scattering (i.e., and ). It follows from (25) that each node must estimate in polychromatic environments its polar coordinates and the l-th ray’s DoA and its amplitude . This would often result in errors which may cause a channel mismatch, thereby hindering the proposed beamforming performance. The backward channel estimate of the k-th node is then given by where and are the errors on and the combined phase , respectively. It follows then from (26) that where and . Exploiting the law of large numbers and assuming that and are statistically independent, we obtain at large K where is the variance of and is determined as follows where , , and

Injecting (28) in (14) yields the following new P-RDCB:

Obviously, we observe from (31) that the k-th node requires knowledge of both and to be able to derive its corresponding weight . However, to avoid the costly calculations of the integrals in (29) and (30), we derive in the sequel their expressions in closed form.

Let us first focus on . Assuming that is sufficiently small to satisfy , could be rewritten as where is a random variable whose pdf is where , , and . Please note that we assume both and in (33) to be uniformly distributed random variables over . Furthermore, if is picked large enough so that holds, we have

Finally, the closed form expression of the above expectation could be easily found using the primitive G of given by where and are the Fresnel S and C functions, respectively, and

As far as is concerned, its closed form expression could also be obtained following similar steps. Indeed, for small and large , we have and, hence, it can be solved using the primitive in (35). Please note that the results obtained so far hold for any given angular spread (AS) which is nothing but the variance of the random variables s characterizing the scattering’ strength. However, when AS is relatively small (i.e., low scattering effects), one could derive much more compact and simple expressions both for and . In such a case, s are small and so are s. Consequently, (29) boils down to after exploiting in the second line the Taylor series expansion around 0 of the Bessel function. Besides, for small AS, could be rewritten as

After performing the variable change and exploiting the equivalence between the Bessel and hypergeometric functions, we obtain

It follows from (31), (34), (35), (38), and (40) that depends solely on locally available information at each node, thereby lending itself to a power- and spectrum-efficient RDCB implementation over WSNs even in polychromatic (i.e., scattered) environments. It is worth noting once more, in the absence of localization, DoA, and scatterers estimation errors (i.e., ), that the proposed P-RDCB, , reduces to the nominal P-DCB solution in given [18] by: where and thereby highlighting yet again both the faced hurdles and the gained advantages of implementing DCB robustness to channel parameter estimation errors.

In what follows, we analyze the performance of the the proposed M-RDCB and P-RDCB solutions.

4. Theoretical Performance Analysis of Robustness Gains

In this section, we assess analytically the performance of the new RDCB solutions (i.e., M-RDCB and P-RDCB) in terms of achieved average SNR (ASNR) against the nominal optimal CB solution in (10) so as to assess the theoretical potential gains they could achieve by integrating robustness in their designs. Let denote the ASNR achieved by any CB where the expectation is taken over all node coordinates, forward and backward channels, and channel estimation errors. Unfortunately, the derivation of in closed form turns out to be a tedious task, if not impossible. In this work, we propose to study instead another practically appealing metric that is the average signal to average noise ratio (ASANR) where and . Please note that [17,18,19,20,21,22,23] have shown that and have approximatively the same behaviors. Let us first derive the average received power from any source located at using .

4.1. Implementation in Scattering-Free Environments—Option 1

Let us first derive the average beampattern achieved by . Exploiting the Taylor series expansion around 0 of the exponential function, we obtain where = . Besides, we know that where stands for the Bessel functions of first kind. Injecting (45) in (44) leads to

It follows from (46) that . Consequently, using the nominal optimal DCB, , the power received at decreases with due to the channel mismatch. This is not surprising since the design of does not account for such errors. On the other hand, the average beampattern achieved by the proposed M-RDCB can be calculated as

The above result verifies that for any given estimation errors. Consequently, the proposed beamformer is much more robust than the nominal optimal DCB.

4.2. Implementation in Scattering-Free Environments—Option 2

If Option 2 is adopted, the average beampattern achieved by can be expressed as

Let . It has already been shown in (22) that where is the the hypergeometric function which has a peak at 0 and decreases when x grows large. Besides, we have

Injecting (49) and (50) in (48) yields

It follows from (51) that decreases with and due to the channel mismatch. On the other hand, the average beampattern achieved by the proposed M-RDCB, which accounts for such errors, can be determined as

It follows from (52) that for any localization and DoA estimation errors, in contrast to the nominal optimal DCB, thereby validating again the robustness of the proposed M-RDCB against channel estimation errors. Furthermore, we observe from (51) and (52) that the proposed M-RDCB achieves an important gain over its counterparts in terms of the received desired power, a gain that substantially increases at higher channel estimation errors.

Now, let us turn our attention to the noise powers. Using either , , or , the average noise power can be calculated, respectively, as and with Option 1 or with Option 2. It could be readily shown from (53)–(55) that , making since .

4.3. Implementation in Scattered Environments

From (10) and (25), turns out to be a complex quotient of several random variables and, hence, deriving the closed form expression of its average is extremely difficult, if not impossible. To circumvent this daunting issue, we propose in this paper to derive instead its asymptotic expression for large K. When is implemented in polychromatic environments, the following theorem holds:

For large K and any given AS, we have

See Appendix A. □

Since has a maximum at 0 and decreases rapidly with x, we have and, hence, . This means that the nominal optimal DCB is unable to satisfy the constraint in (5) due to channel estimation errors. Now, the next theorem introduces the asymptotic expression for .

For large K and any given AS, we have

See Appendix A. □

It follows from (57) that when K is large enough. Therefore, the power received from the source S is substantially improved if is implemented in lieu of . Furthermore, we have

Since and both and are lower than 1, holds when K is large enough. Consequently, in such a condition, , thereby proving the superiority of the new P-RDCB against the nominal optimal DCB in polychromatic environments.

5. Numerical Evaluation Results

This section evaluates numerically the performance of the proposed M-RDCB and P-RDCB techniques and gauge them against the nominal optimal CB solution in (10) and to the RDCB in [32] to emphasize on one hand the need for and the benefits of implementing robustness in DCB, and to assess on the other hand the performance gains of the proposed solutions against the best representative RDCB benchmark. The empirical quantities are obtained by averaging over random realizations of of , , for and , for . In all simulations, we assume that the number of rays or chromatics is , , and the noises’ powers and are 10 dB below the source transmit power.

Figure 2 plots the ASANR and ASNR achieved by the proposed beamformer in monochromatic (i.e., scattering-free) environments versus the variances of channel errors , , and for different values of K. Figure 2a considers implementation Option 1 which results in a phase jitter while Figure 2b considers implementation Option 2 which results in localization and DoA estimation errors. We observe from these figures that the analytical and empirical curves of and match perfectly thereby validating the correctness of the derivations in Section 4. Furthermore, we notice that the ASNR and ASANR curves remain very close for or coincide almost perfectly otherwise, thereby proving the insightfulness of the ASANR metric. As far as the proposed M-RDCB performance is concerned, it is able to achieve optimal performance, even for small K, when the channel estimation errors are relatively small to moderate (i.e., in Option 1 or and in Option 2). This confirms the robustness of the new M-RDCB. For extremely large errors, however, it looses only a fraction of a dB. Actually, with the advances made during the two last decades in the field of phase synchronization, localization, and DoA estimation, these channel estimation errors are often very small, making our beamformer’s performance optimal if advanced parameter estimation algorithms are adopted. Nevertheless, the latter naturally come with increased complexity and cost, which certainly burden those of WSN nodes. In this context, the proposed M-RDCB offers the possibility of using inaccurate but low-cost estimation algorithms at negligible performance losses, making it a more appealing and practical for cost-effective WSN deployment in real-world conditions. All these observations corroborate all the results of Section 3.

Figure 2

Average signal to average noise ratio (ASANR) and average signal to noise ratio (ASNR) of proposed monochromatic robust distributed collaborative beamforming (M-RDCB) in monochromatic environments under: (a) implementation Option 1, and (b) implementation Option 2.

Figure 3 displays the achieved ASANR gain of proposed M-RDCB against the nominal optimal DCB and the M-RDCB in [32] for different values of K. Figure 3a considers the first implementation option while Figure 3b considers the second. We observe from these figures that M-RDCB largely outperforms both benchmarks for any given K, , , and . For instance, when Option 1 is adopted, it achieves for ASANR gains of and dB against the nominal optimal DCB and RDCB in [32], respectively. If Option 2 is adopted, these ASANR gains increase, respectively, to as much as and dB when and . As could be observed from Figure 3, these gains increase rapidly with both K and channel estimation errors. These observations corroborate all the results of Section 4 and verify the net superiority of robust M-RDCB techniques in scattering-free environments.

Figure 3

ASANR gain of proposed M-RDCB in monochromatic environments under: (a) implementation Option 1, and (b) implementation Option 2.

Figure 4 shows the ASANR, ASNR, and the ASANR gain of the proposed P-RDCB in polychromatic (i.e., scattered) environments. Figure 4a,c,e plot its achieved ASANR and ASNR versus K, , , and while Figure 4b,d,f compare its achieved ASANR with those of the nominal optimal DCB and RDCB in [32]. As expected, approaches the optimal ASANR performance level even in polychromatic environments, and that is for all tested AS values. In such environments, achieves ASANR gains of until and dB against both benchmarks, respectively. As can be observed from Figure 4b,d,f these gains increase rapidly with both K and channel estimation errors. For instance, according to Figure 4f, the ASANR gain over the two benchmarks increases by approximatively against both when is twice as large or and , respectively, when is four times as large. These additional observations further verify the high robustness of the proposed P-RDCB against channel estimation errors, key practical feature that allows its cost-effective integration in real-world WSN applications. Again, according to Figure 4, we verify that the analytical and empirical values of match perfectly when or are very close otherwise. All these observations corroborate once more the discussion of Section 4.

Figure 4

ASANR, ASNR, and ASANR gain of proposed P-RDCB in polychromatic environments.

6. Conclusions

We have proposed a new DCB solution robust against major channel estimation impairments, namely, phase synchronization, localization, direction-of-arrival (DoA), and/or channel scatterers estimation errors over dual-hop transmissions from a source to a destination communicating through a WSN of K nodes. Exploiting an efficient asymptotic approximation at large K, we have developed alternative RCB solutions that not only account for estimation errors, but also adapt to different implementation scenarios and wireless propagation environments ranging from monochromatic to polychromatic ones. Besides, in contrast to existing techniques, our new RCB solutions are distributed in that they do not require any information exchange among nodes, thereby dramatically improving both WSN spectral and power efficiencies. Simulation results have confirmed that the proposed RDCB techniques are much more robust in terms of achieved SNR against channel estimation errors than best representative CB benchmarks, yet at much lower complexity, power cost, and overhead, making them suitable for WSN deployment in the harsh environments that characterize operation in real-world conditions.