Optimization of SignaltoNoiseandDistortion Ratio for Dynamic Range Limited Nonlinearities
Abstract
Many components used in signal processing and communication applications, such as power amplifiers and analogtodigital converters, are nonlinear and have a finite dynamic range. The nonlinearity associated with these devices distorts the input, which can degrade the overall system performance. Signaltonoiseanddistortion ratio (SNDR) is a common metric to quantify the performance degradation. One way to mitigate nonlinear distortions is by maximizing the SNDR. In this paper, we analyze how to maximize the SNDR of the nonlinearities in optical wireless communication (OWC) systems. Specifically, we answer the question of how to optimally predistort a doublesided memoryless nonlinearity that has both a “turnon” value and a maximum “saturation” value. We show that the SNDRmaximizing response given the constraints is a doublesided limiter with a certain linear gain and a certain bias value. Both the gain and the bias are functions of the probability density function (PDF) of the input signal and the noise power. We also find a lower bound of the nonlinear system capacity, which is given by the SDNR and an upper bound determined by dynamic signaltonoise ratio (DSNR). An application of the results herein is to design predistortion linearization of nonlinear devices like light emitting diodes (LEDs).
Nonlinear distortion, dynamic range, clipping, predistortion, optical wireless communication.
1 Introduction
In addition to being nonlinear, many components in a signal processing or communication system have a dynamic range constraint. For example, light emitting diodes (LEDs) are dynamic range constrained devices that appear in intensity modulation (IM) and direct detection (DD) based optical wireless communication (OWC) systems [1][2]. To drive an LED, the input electric signal must be positive and exceed the turnon voltage of the device. On the other hand, the signal is also limited by the saturation point or maximum permissible value of the LED. Thus, the dynamic range constraint can be modeled as twosided clipping. The same situation may happen in other applications such as digital audio processing [3].
Both nonlinearity and clipping result in distortions which may cause system performance degradation. SNDR is a commonly used metric to quantify the distortion that is uncorrelated to the signal [4][7]. Previous work in this area mainly concentrated on a family of amplitudelimited nonlinearities that is common in radio frequency (RF) system design involving nonlinear components such as power amplifiers (PAs) and mixers.
Different from the previous work, our study discusses the class of nonlinearities with a twosided dynamic range constraint that is more commonly found in optical and acoustic systems. Authors in [8][12] illustrated the impact of LED nonlinearity and clipping noise in OWC systems. Some predistortion strategies were proposed in [13][15]. However, to the best of our knowledge, the optimal nonlinear mapping under the twosided dynamic range constraint has not been studied.
There are two major differences from the amplitudelimited nonlinearity. First, the signal will be subject to turnon clipping and saturation clipping to meet the dynamic range constraint. Second, DC biasing must be used to shift the signal to an appropriate level to minimize distortion. In this paper, we will show that the ideal linearizer that maximizes the SNDR is a doublesided limiter that has an affine response. The parameters of the response can be calculated from the distribution of the input signal and the noise power.
In additional to deriving the SNDRoptimal predistorter, we also relate a lower bound on channel capacity to the SNDR, further motivating the SNDR considerations. Finally, we employ another common distortion metric, dynamic signaltonoise ratio (DSNR) to provide an upper bound on the doublesided clipping channel.
The remainder of this paper is organized as follows: Section II introduces the system model for dynamic range limited nonlinearity and the corresponding SNDR definition. In Section III, we derive the optimal nonlinear mapping that maximizes the SNDR and illustrate some examples. In section IV, we related the SNDR to the capacity of the nonlinear channel. Finally, Section VII concludes the paper. The detailed proofs of this paper are deferred to the Appendices.
2 System Model and SNDR Definition
2.1 System Model
Let us consider a system modeled by
(1) 
where is a realvalued signal with mean and variance ; is a zeromean additive noise process with variance ; is a memoryless nonlinear mapping with dynamic range constraint .
For notational simplicity, we omit the dependence in the memoryless system and replace and by and . Then we have an equivalent system modeled by
(2) 
where is a memoryless nonlinear mapping with dynamic range constraint and is a zeromean signal with variance .
2.2 SNDR Definition
According to Bussgang’s Theorem [16], the nonlinear mapping in (2) can be decomposed as
(3) 
where is the distortion caused by and is a constant, selected so that is uncorrelated with , i.e., . Thus
(4) 
The distortion power is given by
(5) 
The signaltonoiseanddistortion ratio (SNDR) is defined as
(6) 
The definition of SNDR here is a little bit different from that in [7], because all the signals are real and the distortion contains DC biasing. Thus, the distortion power is modeled as variance rather than the secondary moment.
We see from (6) that the SNDR is related to the distribution of , the noise power and the nonlinear mapping . Our aim in the next section is to determine the function that maximizes the SNDR given a signal distribution and the twosided clipping constraint.
3 SNDR Optimization and Examples
3.1 Optimization of SNDR
Similar to [7], let us use a function to normalize the nonlinear mapping :
(7) 
where . Let and substitute (7) into (6), we obtain
(8) 
where is the variance of and .
The SNDR optimization problem can be stated as follows:
(9)  
(10) 
for a given distribution of , dynamic range and noise power .
Fig. 1 illustrates an example of the . The region of is divided into three sets , and .
(11)  
(12)  
(13) 
Thus, to determine a nonlinear mapping , we need to find the sets , , and the shape of the function in .
We will solve this problem with the following steps:

find the optimal given , , ;

show that should be as large as possible;

determine and for the optimal solution.
Lemma 1
Assume that the sets , and are known, and . The function that maximizes the SNDR expression in (8) is of the form
(14) 
where
(15) 
(16) 
with
(17) 
and is the indicator function:
(18) 
This lemma holds if and only if satisfies for all .

See Appendix missing6.
This result rules out the functions whose shape over is nonlinear. Fig 2 demonstrates examples of functions that may satisfy Lemma 1. Here, the slope of the linear curve in can be either positive or negative.
Lemma 1 answered the question pertaining to the best shape of the function with given , and . The remaining question is how to determine the optimal sets , and so that the SNDR is maximum. This turns out to be a very challenging problem since we are seeking joint optimization over multiple sets. Let us consider first.
Lemma 2
Given sets , and , if can be enlarged to such that or , then a higher SNDR can be achieved.

See Appendix missing7.
Fig. 3 shows how Lemma 2 works. can be enlarged by occupying the subsets of and . The larger the set , the better the SNDR that can be achieved. Just as Lemma 1, Lemma 2 holds if and only if satisfies for all , that is, or .
Even with the set determined, we still need to determine and .
Lemma 3
If , the that maximizes the SNDR satisfies and ; if , the that maximizes the SNDR satisfies and .

Let us compare the SNDR between Fig. missinga and Fig. b. For , if there is a subset of in or a subset of in , which is illustrated in Fig. b, then we see that is decreased while the variance of is increased. Thus, the of Fig. b is less than the SDNR of Fig. a. Similarly, we can draw the same conclusion for the case with .
In the final analysis, Lemma 1, Lemma 2 and Lemma 3 imply that the optimal , and , in the sense of maximizing the SNDR, are , and if ; or , and if .
Theorem 1
Within the class of satisfying , the following maximizes the SNDR expression in (8):
(19) 
for , or
(20) 
for , where the and are found by solving the following transcendental equations:
(21) 
(22) 
with
(23) 
(24) 
(25) 
(26) 
(27) 
and is the probability density function (PDF) of . The optimal SNDR is found as
(28) 
where
(29) 
and
(30) 

See the proofs of Lemma 1, Lemma 2 and Lemma 3.
Theorem 1 establishes that the nonlinearity in the shape of Fig. 5 is optimal.
Figure 5: Illustration of optimal functions to maximize the SNDR Predistortion is a wellknown linearization strategy in many applications such as RF amplifier linearization. For the dynamic range constrained nonlinearities like LED electricaltooptical conversion, predistortion has been proposed to mitigate the nonlinear effects. Specifically, given a system nonlinearity , it is possible to apply a predistortion mapping so the overall response is linear. According to Theorem 1, it is best to make equal to the function given in (19) or (20) if is normalized with dynamic range constraint . Using the analytical tools presented above, we can answer the questions regarding the selection of the gain factor , DC biasing and the clipping regions on both sides, or equivalently, the sets and . Theorem 1 shows that these optimal parameters (in terms of SNDR) depend on the PDF of and the dynamic signaltonoise ratio . Thus, our work can serve as a guideline for the system design. In the next subsection, examples are given to illustrate the calculations of the optimal factors and .
3.2 Examples for selections of optimal parameters
In the last subsection, we learned that the optimal factors and can be calculated by solving two transcendental equations (21) and (22). However, there may not be closedform expressions for the solutions. Additionally, solving (21) and (22) may result in multiple solutions, but we only keep the realvalued ones since all the signals here are realvalued.
Here, let us take into account a specific class of input signals whose distributions exhibit axial symmetry, such as uniform distribution and Gaussian distribution. When the distribution of the input signal is axial symmetric, the optimal clipping regions and are also symmetric. Thus, , and . Then the factors and can be calculated:
(31) (32) We see that the DC biasing will be the midpoint of the dynamic range. When the gain factor , it can be further expressed as:
(33) When the gain factor , it can expressed as:
(34) There is still no closedform expression for gain factor . Next, as examples, let us consider the calculations for uniform distribution and Gaussian distribution specifically.
Example 1
When the original signal is uniformly distributed in the interval , we infer that the normalized signal is uniformly distributed in the interval with the PDF
(35) For the case with , it is straightforward to calculate
(36) (37) (38) Equation (38) can be rewritten as a quadratic equation
(39) Thus, we can obtain a closedform solution for the optimal :
(40) We know that there should be two solutions for equation (39). In fact, the other solution is , which means that both and are 0. Thus, the solution given by (40) is the unique optimal selection for the gain factor . If is desired, the optimal solution is
(41) Example 2
When the original signal is Gaussian distributed, then the normalized signal has a standard Gaussian distribution with the PDF
(42) For the case with , we have
(43) (44) where is the error function with the definition
(45) (46) Here the optimal does not have a closedform expression but can be easily calculated numerically. We can draw the similar conclusion for the case with .
3.3 Numerical results
Fig. 6 shows the optimal as a function of DSNR for the above examples.
Next, we illustrate the SNDR of two different nonlinear mappings. is the optimal solution chosen by Theorem 1. is a fixed mapping given below:
(47) The corresponding SNDR curves are shown in Fig. 7. This example illustrates that the nonlinearity yields a higher SNDR as compared to the other nonlinearity, as expected according to Theorem 1.
4 Relationship Between SNDR and Capacity
4.1 Lower Bound on Capacity
The capacity is given by
(48) where is the mutual information between and [18]. To obtain the capacity of the dynamic range constrained channel, we need to solve the following optimization problem:
(49) for a specific zeromean noise with variance . Moreover, it can be simplified as:
(50) which means that we need to find an input distribution in the interval to maximize the mutual information. Specially, when the noise is Gaussian, the issue is similar to Smith’s work in [17]. In this case, if DSNR is low, the capacity is achieved by an equal pair of mass points at and ; if DSNR is high, the asymptotic capacity is the same as the information rate due to a uniformly distributed input in [17].
However, in most cases, we are most interested in the achievable data rate given a nonlinear channel mapping with any input and any noise. Similar to the work in [7], we obtain a lower bound on the information rate:
(51) (52) by referring to (8). Since for any input distribution , by setting to be the PDF of a zeromean Gaussian r.v., we obtain
(53) with the SNDR evalutated for a Gaussian .
4.2 Upper Bound on Capacity
In this subsection, we find an upper bound for the capacity. Similar to [7], supposing is the PDF of that maximizes the capacity, i.e.,
(54) We can write the capacity as
(55) Next, we bound the entropy with the entropy of a Gaussian , yielding
(56) where with . Specifically, if the noise is Gaussian, we have the upper bound:
(57) 4.3 Example of Bounds
Since SNDR is determined by DSNR and the distribution of signal, we plot the bounds as functions of DSNR for Gaussian distributed signal, which is shown in Fig. 8. We also compare the lower bounds given by two different nonlinear mappings and , which are introduced in the last section. This example illustrates that the nonlinearity chosen according to Theorem 1 yields a tighter lower bound as compared to the other nonlinearity. In addition, we can see that the capacity of Gaussian channel as determined by Smith [17] is between the lower bounds and upper bound that we have.
5 Conclusion
The main contribution of this paper is the SNDR optimization within the family of dynamic range constrained memoryless nonlinearities. We showed that, under the dynamic range constraint, the optimal nonlinear mapping that maximizes the SNDR is a doublesided limiter with a particular gain and a particular bias level, which are determined based on the distribution of the input signal and the DSNR. In addition, we found that provides a lower bound on the nonlinear channel capacity, and serves as the upper bound. The results of this paper can be applied for optimal linearization of nonlinear components and efficient transmission of signals with doublesided clipping.
6 Proof of Lemma 1
Since we are solving the optimization problem w.r.t. a function, the functional derivative is introduced here [7][19]. By using the Dirac delta function as a test function, the notion of functional derivative is defined as:
(59) Just as the variable derivative operation, the linear property, product rule and chain rule hold for functional derivative. In addition, from (59), we infer that
(60) (61) To maximize the SNDR w.r.t , we need
(62) We infer that
(63) Similarly,
(64) (65)