On-line version ISSN 2179-1074
J. Microw. Optoelectron. Electromagn. Appl. vol.11 no.1 São Caetano do Sul June 2012
Increasing energy efficiency in OCDMA network via distributed power control
Fábio R. DurandI; Bruno A. AngélicoI; Taufik AbrãoII
In this work, we investigate the utilization of transmission power control as mechanism to increase the energy efficiency in optical code division multiplexing access (OCDMA) networks. We have modeled the energy efficiency considering the optical fiber transmission and network infrastructure as encoders, decoders and star coupler. In the analyzed scenario, this model confirms that the energy consumption of the network infrastructure is larger than the energy consumption of the transmission infrastructure; On the other hand, we have proposed a scheme to define the energy efficiency according to the BER level accomplished with the quality of service (QoS) requirements. The main results showed that it is possible to save 70% of the transmitted energy per bit with the penalty of one order of magnitude of BER.
Index Terms: optical code division multiplexing access, distributed power control algorithm.
The informatics and telecommunications industries have become a component of the economies from the nations. The growing of Internet, information and communications technology (ICT) infrastructures have been showing an exponential behavior and the electrical power consume has increasing in the same way. However, the energy resources available in the world are limited and the energy efficiency in all sort of industry has been encouraged. The energy efficiency in the ICT is represented by the energy consumption per bit of data transported and/or processed . Nowadays, the ICT sector is responsible for approximately 5 % of the total electrical power consumption in developed national economies and especially the Internet consumes 1 % of this total electrical power consumption . Recent studies have showed the importance of the consideration of energy consumption in optical communications design, considering the transmission infrastructure (transmitters, receivers, fibers and amplifiers)  and network infrastructure (switchers and routers)  aspects. In a global scale network, the energy consumption of the switching infrastructure is larger than the energy consumption of the transport infrastructure . In this context, it is necessary to improve the energy efficiency of switching and to optimize the network design to reduce the quantity of switching and overheads. Low energy efficiency is observed mainly at access network, because the access equipment consumes almost 50% of energy consumed in the core or metro equipment . However, the amount of bits transported in core or metro network is larger than in the access network.
Optical-code-division multiple access (OCDMA) based networks  have attracted a lot of interests due to various advantages including asynchronous operation, high network flexibility, protocol transparency, simplified network control and potentially enhanced security . In OCDMA each different code defines a user and different codeusers can share a common channel. In a common channel, the interference that may arise between different codeusers is known as multiple access interference (MAI) and can limit the number of codeusers utilizing the channel simultaneously. OCDMA can be divided into noncoherent (unipolar) and coherent (bipolar) systems. The noncoherent systems are based only on intensity modulation of optical power , and coherent systems are based on modulation of amplitude and phase . The coherent code is trueorthogonal and noncoherent code is pseudoorthogonal. As a consequence, the performance of coherent codes is higher than noncoherent codes when analyzing the signaltointerference plus noise ratio (SINR) . However, the main drawback of coherent OCDMA lies in the technical implementation difficulties, concomitant with the utilization of phaseshifted optical signals . The noncoherent codes can be classified in one dimensional (1D) and two dimensional (2D) codes. In 1D codes, the bits are subdivided temporally into many short chips with a designated chip pattern representing a user's code. On the other hand, in 2D codes the bits are subdivided into individual time chips, and each chip is assigned to an independent wavelength from a discrete set of wavelengths . The 2D codes have better performance than 1D codes and can significantly enhance the number of active and potential users . In the OCDMA, the nearfar ratio increases the multiple access interference effects, because each active node in transmission mode contributes with MAI and the power penalty varies dynamically. Therefore, the static power budget design does not solve this problem . Furthermore, if the distances between the nodes are quite different, as would be the case in practice, the power received from different nodes will be significantly different. Thus, considering as reference one OCDMA node, the performance of closer nodes is many orders of magnitude better than that of far nodes.
Hence, an efficient power control is needed to overcome this problem and to enhance the performance and throughput of the optical network . This could be achieved via signal-to-interference (SINR) optimization . In this case, analogous to a wireless CDMA cellular system, the power control, centralized or distributed, is one of the most important issues since it has a significant impact on both performance and capacity; it is the most effective way to avoid the nea-rfar problem and to increase capacity. Previously, this problem has been investigated in OCDMA networks to solve the nearfar problem  and to obtain quality of service (QoS) at physical layer ; however this issue has not been investigated with focus on the energy efficiency.
In this work, we study the utilization of distributed power control as mechanism to obtain energy efficiency improvement in OCDMA access network with non-coherent 2-D codes. The performance of this architecture networks is analyzed considering the loss characteristics of encoders, decoders and star coupler. Our objective consists in investigating the viability of power control in order to increase the energy efficiency of the OCDMA network and to determine the best tradeoff between the energy efficiency and near-far mitigation.
This paper is organized as following: Section II describes the architecture of the network for 2D codes. Section III discusses the proposed methodology for increasing OCDMA energy efficiency. The performance of proposed methodology systems is developed in Section IV. Section V shows representative numerical results for the proposed approach; finally, in Section VI the main conclusions are presented.
II. NETWORK ARCHITECTURE
The OCDMA architecture considered in this work is formed by K nodes interconnected by passive star coupler, in a broadcastandselect pattern as shows Fig. 1. For viability characteristics, we consider network equipment, such as codeprocessing devices (encoders and decoders at the transmitter and receiver) and star coupler. Such devices could be made using robust, lightweight, and lowcost technology platforms with commercial-off-the-shelf technologies , .
The transmitting and receiving nodes create virtual path based on the code and the total link length is given by dij = + , where is the link length from the transmitting node to the star coupler and is the link length from the receiving node to the star coupler. The received power at a jth node is given by Prj = astar pi exp (−α f dij ) , where pi is the transmitted power by ith transmitter node, αf is the fiber attenuation (km1) and astar is the star coupler attenuation (linear units). Considering deciBell units, astar = 10 log (K )−[10log 2 (K )log 10 δ], where, δ is the excess loss ratio.
The 2D wavelengthhopping/ timespreading code sequence is illustrated in Fig. 2. This 2D code is transmitted and its destiny in the network is determined by a particular code sequence. The 2D codes can be represented by Nλ× NT matrices, where Nλ is the number of rows, that is equal to the number of available wavelengths, and NT is the number of columns, that is equal to the code length. The code length is determined by the bit period TB which is subdivided in small units called chips, each of duration Tc = TB/ NT. In each code, there are w short pulses of different wavelength, where w is called the weight of the code. An (Nλ× NT, w, Ia, Ic) code is the collection of binary Nλ× NT matrices each of code weight w; Ia and Ic are nonnegative integers and represent the constraints on the autocorrelation and crosscorrelation , respectively.
The 2D OCDMA utilizes multiwavelength sources such lightemitting diodes, amplified spontaneous emission noise from erbiumdoped fiber amplifier (EDFA), gain switched Fabry Pérot lasers, and supercontinuum generation . These sources avoid the need for rapidly wavelength hop according the wavelengthhopping pattern. The encoder essentially creates a combination of two patterns: a wavelengthhopping pattern and a timespreading pattern. The common technology applied in encoders/decoders with delay lines are array waveguide gratings (AWGs), thin-film filters (TFFs), simultaneously as fiber Bragg gratings (FBGs), holographic Bragg reflectors (HBRs), chirped Moire gratings (CMGs) . The schemes of encoders commonly utilized are based on AWG, TFF and FBGs, as are showed in Fig. 3.
The losses associated with the encoders/ decoders are given by 
where aAWG is the AWG loss, aDelay is the delay line loss, aTFF is the TFF loss, aBragg is the FBG loss and aCirulator is the circulator loss. The loss usual value for these equipments areaAWG= 2.5 dB, a Delay = 1dB, a TFF = 0.5 dB, a Bragg = 0.5 dB and a Cirulator = 3dB. The AWG encoder has approximately uniform loss (6 dB) independently of the number of wavelengths (Nλ) . For this characteristic, we will consider the AWG as encoder/decoder in the rest of this work.
III. THEORETICAL ANALYSIS
a. Energy efficiency and power allocation
The necessary energy in each i-th node for the transmission of 1 bit can be expressed as ,
where K is the number of nodes, Tbit is the time to transmit one bit over the network given by Tbit=1/R, where R is the bit rate (bits/seconds) and pi is the transmitted power. In order to determine the energy per bit values it is necessary to define the individual node transmitting power (pi). The pi is obtained by power budget or power control and it is associated to a QoS, a SINR and a maximum bit error rate (BER) tolerated by the optical nodes. Under a power control situation, each optical node adjusts its transmitting power in an attempt to maximize the number of transmitted bits with minimum energy consumption. This concept can be formulated by the energy efficiency definition ,
where M is the number of bits in each transmitted packet, L is the number of information bits contained in each data packet and f (γi) is the efficiency function, which approximates the probability of error-free packet reception. It can be approximated by f (γi )=(1− BERi )M , with γi being the SINR for the i-th node, given by ,
where gij are the total loss in the path that connects i-th-Tx node to j-th-Rx node, Ii is the total interference power level comes from the others transmitters nodes and Ni is the receiver power noise level. Both terms Ii + Ni will be described in details in the next subsection, eq. (8).
In the same way this concept is used to define a metric named utility that is the number of bits received per energy expended or the relation of the throughput and power dissipation . Hence, for each ith node, the maximum number of transmitted bits occurs at power level for which the partial derivative of εi with respect to pi equals to zero, ∂εi/∂pi = 0 . So, the derivative of energy efficiency can be obtained referring to efficiency function f (γi ) and eq. (4), setting:
From (5) assuming pi > 0, the necessary condition to maximize the energy efficiency is immediately obtained:
So, adopting the widely accepted approximation for BER performance with uncoding system (M=L) and binary modulation, the bit error rate can be approximated as BERi = e −γi. Hence, we have: As a result, we can conjecture that the optimum power allocation criterion in terms of energy efficiency is given by
where κ is a constant of proportionality. Besides, further simplification can be obtained if the condition M=1 could be assumed:
b. Proposed scheme
Note that in order to satisfies (6) it is necessary adjust the SINR at each received node equals to the target SINR, γi = . In this context, we propose the utilization of power control iterative algorithm in order to establish the lower energy per bit according to the OCDMA network QoS requirements. Fig. 4 shows the flowchart with the adopted scheme to assign the energy per bit transmitted in each node. This scheme aims to save energy per bit at each node and, as a consequence, decreases the total energy consumption overall the network.
This scheme is based on the power control with the restriction of QoS requirement based on BER level. The entries are the BER level and the network topology, comprising the following parameters: node distances, links lengths, fiber parameters, EDFA preamplifiers value, OCDMA code parameters, and so forth.
The iterative power control algorithm defines the transmitted power according to the number of active transmitting nodes aiming to establish individual target value of BER, and taking into account as low level energy per bit as possible. The outputs of the scheme are the necessary number of iterations of the iterative power control algorithm and the optimized value of energy per bit necessary for each node. In the next subsections we will illustrate BER computation, the main characteristics of the power control procedure, as well as the way to apply the power control in order to obtain a more energy efficient OCDMA network.
IV. PERFORMANCE EVALUATION
a. Signaltointerference plus noise ratio (SINR) and Bit error rate (BER)
In , Inaty et al. derive an expression for optical SINR in a multirate chipsynchronous 2DOCDMA system. Assuming singlerate case, the SINR for the i-th node of a 2-D code OCDMA system can be precisely written as:
where gij and Gamp are, respectively, the total loss and gain of the preamplifier in the path that connects ithTx node to j-th-Rx node. σ 2 is the average variance of the Hamming aperiodic cross-correlation amplitude , PN is the spontaneous amplified emission (ASE) noise power in the preamplifier considering the two polarization mode presented in a single mode fiber and pj is the other nodes transmitted power. Note that the usual receiver noise power includes thermal noise, shot noise and optical preamplifier noise. However, ASE in the optical preamplifier will be the main limiting factor (in addition to the MAI), compared to thermal and shot noise at the receiver . In our work, the receiver noise power is represented by
where nsp is the spontaneous emission factor, typically around 2 to 5, h is Planck's constant, f is the carrier frequency, G is the amplifier gain and Bo is the optical bandwidth. Ideally, in order to reduce the ASE noise power, the optical bandwidth can be set to a minimum of Bo = 2R. Assuming Gaussian noise approximation, the BER associated with the ith node is related to its SINR by where erfc(.) is the complementary error function.
b. Power control
The power budget analysis establishes the transmitted power necessary to reach the photodiode sensibility. The power budget in dB must satisfy the follow inequality,
where ps is the receiver power sensibility, NMAI is the MAI interference variance, CEnc is the encoder losses and CDec is the decoder losses. The losses associated with the encoders/ decoders are given by (1).
We observe in (10) the static characteristics of the power budget design, mainly because the MAI varies with the number of active nodes transmitting simultaneously. The design based on static power budget results in a transmitting power higher or lower than the necessary as a result, an increase in the nearfar problem occurs. To solve this problem, it is adequate to apply dynamic power control that defines the transmitting power according the number of active transmitting nodes.
The power control in optical networks is an optimization problem aiming to establish individual target value of SINR, denoted by ; hence, each node unit transmitting with power pi, has to be controlled in order to satisfy the relation γi ≥ . The maximum achievable SINR at the receiving nodes is equal to the maximum achievable carrier to interference ratio (CIR) at the receiver output times NT2/σ2. Denoting Γi as the CIR at the required decoder input, in order to get a certain QoS which is associated to a maximum tolerable BER at the i-th optical node, and defining the K-dimensional column vector of the transmitted optical power p = [p1, p2,…… pK]T , then the optical power control optimization problem consists in finding the optical power vector p* that minimizes the following cost function :
subject to the constraints:
where 1T = [1, ..., 1] and Γ* is the minimum CIR to achieve a desired QoS. The elements gij that represent the connections of transmitterreceiver pairs constitute the network attenuation matrix:
Using matrix notations, (12) can be written as,
where I is a K-dimension identity matrix, H is the normalized interference matrix, whose elements can be evaluated by
and the ith element of vector u is given by:
Note that in (15) we have a scaled version of the noise power. Solving (14), substituting inequality by equality, we get the optimized power vector solution through matrix inversion:
This optical power vector represents the case of power equilibrium at the receiver node, and is the optimal power required achieving the target CIR. Increasing the value of the CIR would result in higher optical power values that could result in a maximum power higher than the allowed. In this case, a solution would be fixing or decreasing the target CIR value or removing (switching off) some users from the network. The centralized power control is obtained by matrix inversion, as illustrated in (17), and corresponds to an existence of a central node. The central node storages information about all physical network architecture like fiber length between nodes, amplifier position, regular updating about transmission establishment and dynamic of the optical traffic. These characteristics are the drawback of centralized control . On the other hand, the distributed power control algorithm (DPCA) synthesis consists of the development of a systematic procedure for the vector p evolution in order to reach the optimum value, based on the , γi and pi values. The optimum solution for the power allocation problem satisfies the following iterative process  (For more details about the derivation of (18), see Appendix I):
where n is the number of iterations and α is the convergence. Note that the convergence factor α in (18) is the numerical integration step to solve an ordinary differential equation, which with some minor alterations is also considered in many other studies of power control . The algorithm represented by eq. (18) was evaluated for positive and no greater than 1 value of α, since it was shown in  that the algorithm divergence outside this interval. This parameter is responsible for the convergence speed: values close to 1 indicate fast convergence; however, the quality of the solution after convergence is not excellent when compared with values close to 0 . There are more details about several aspects that we do not discuss here, such as, convergence, proximity to the optimum value, and sensibility to estimation errors as well, that are discussed in details in .
Observing the scalar equations in (18) we can infer that the transmitted power actualization of each node depends only on its own parameters. This apparent uncoupling among nodes follows from the fact that we are estimating the SINR directly and not calculating it from the equation (8). It is obvious that for obtaining the SINR from (8) it would demand an expensive procedure in order to estimate all channel gains (interferers and desired user), imposing a centralized implementation of the DPCA, and finally increasing the algorithm implementation complexity. So, we must avoid this approach. However, the recursion in eq. (18) can be
effectively implemented in each (i-th) optical node because all necessary parameters (i.e., αi, the QoS level given by , and the transmitted power pi[n] ), except γi [n], can be considered known at the i-th node. Under the distributed power control approach, the SINR γi [n] is obtained at the correspondent destination optical node that demodulates the signal from the optical user i-th. In this way, the destination optical node estimates γi [n] , quantizes it in a convenient number of bits, and transmits this information to the optical user ith through a return channel. Indeed, it is possible to measure γi [n] without the effective knowledge of the information from all the optical interferer nodes, only sensing the ith optical signal for a period of time at the destination optical node, following by a sensing of all (sum) interferer signals (in the absence of ith transmitted signal) in a second period of time. In summary, (18) depends only on local parameters allowing the power control to work in a distributed manner, i.e., each one of the K links would be able to carry out separately the respective power control procedure. This brief explanation justifies the name distributed power control algorithm (DPCA) for the set of recursions implicit in (18).
To cope with the SINR optimization, the minimum power constraint (which is also called sensitivity level) assures detection of the optical signal by all optical devices. On the other hand, the maximum power constraint guarantees the minimization of nonlinear physical impairments, because it makes the aggregated power on a link to be limited to a maximum value.
The convergence characteristics of the transmitted power given by DPCA present an asymptotic behavior for any strictly positive initial condition (pi > 0) . Therefore, the energy efficiency presents the same asymptotic behavior of the transmitting power, and then an increase in the number of iterations results in the convergence for the transmitting energy .
V. NUMERICAL RESULTS
The scenario considered in our study was presented in Fig. 1. The adopted distances between Tx nodes and star coupler are shown in Fig. 5.a, while distances between Rx nodes and star coupler considering 31 nodes are presented in Fig. 5.b. The nodes are uniformly distributed over an area with a radius between 2 km and 50 km; hence, the range of the total link length is within [4; 100] km. We consider the transmission of an OCDMA signal over a nonzerodispersionshifted (NZD) fiber (ITUT G.655) with attenuation αf = 0.2 dB/km and nonlinear parameterγ= 2 (W.km)1 utilizing AWG encoders/decoders. Applying typical parameters values: nSP = 2, h = 6.63× 10−34 J/Hz, f = 193.1THz, G = 20 dB, Bo = 30 GHz, immediately we obtain ASE power noise PN = 15 × 10−7 A2. The transmission rate considered is 2.5 Gbps, SINR target = 27 dB for 2D codes (4 × 101, 4, 1, 0).
Fig. 6.a and 6.b present the energy per bit assignment per node for OCDMA system considering only the transmission infrastructure (transmitters, receivers, fibers and amplifiers) and transmission plus network infrastructure (encoders, decoders and switchers), respectively.
As one can observe from Fig. 6.a and 6.b, there is a considerable difference between the values of energy per bit per node considering only the transmission infrastructure and considering the overall system, i.e. aggregating losses and multiple access interference in transmission plus network. Indeed, in accordance with , , our numerical results show the lower energy per bit necessary for the transmission system even we do not use optically amplified links and our system adopt only EDFA preamplifier. The losses due to the OCDMA encoder/decoder, as well as the loss introduced by the star coupler, have been resulted in a substantial difference of three orders of magnitude in the energy per bit regards to other conventional OCDMA systems with same parameters considered herein. This result is confirmed by static power budget considering losses of 6 dB for encoder, 6 dB for decoder and 14 dB for star coupler. These numerical results for the energy per bit were obtained applying the centralized control strategy, i.e. solving (17) by matrix inversion. However, the drawback of centralized control is the necessity of central node with parameters knowledge of the physical network architecture.
Fig. 7 shows the energy per bit per node as a function of the number of iterations for the classical iterative DPCA  with α= 0.8 , considering the same scenario (transmission and network infrastructure) presented in Fig. 5. This adopted value for the convergence velocity factor allows to combine a relative fast convergence velocity with acceptable quality of the solution after 4050 iterations . A decrease of 75% in the number of iterations required to convergence was observed in  when the value of α increase from 0.2 to 0.8. Indeed, from Fig. 7, one can see when iteration evolving, the energy per bit per node of the DPCA approaches the values calculated by centralized control strategy, eq. (7), represented by horizontal dot lines. Furthermore, we observe an increase in the energy per bit per node with the number of iterations and this characteristic could be utilized to control the energy efficiency. In our proposed scheme, the number of iterations is the parameter in which the transmitted energy per bit per node can be controlled according to the SINR penalty. However, the SINR penalty affects the BER and could result in an unrealistic OCDMA system. This issue is directly or indirectly considered in Figs. 8 and 9.
Fig. 8 shows the BER performance as a function of the number of active optical nodes for a code weight of 4 and a code length of 101, applying the proposed distributed DPCA scheme. The scenario is the same adopted in Fig. 5 with α = 0.8. The BER system with DPCA considering different percentage of energy per bit of 5%, 27% and 100% is obtained; these percentages are obtained with, 3, 12 and 60 iterations, respectively. Note that in this scenario, 100% energy per bit is achieved after 60 iterations (total DPCA convergence), in a similar way showed in Figure 7.
It can be seen from Fig. 8 the impact of power control procedure (in terms of transmitting energy per bit) on the bit error rate performance. The deployment of DPCA with transmission of 100% of energy per bit approaches the BER considering the ideal system , i.e. equal distance for all nodes. Indeed, with lower energy per bit (5%), the SINR declines because the MAI and nearfar effects are strongly influenced by an increase in the number of active nodes. Under this situation, an error occurs when crosscorrelational pulses from the (K 1) interfering reaches a level higher than the autocorrelation peak, changing a bit zero to a bit one. As a consequence, under the nonoptimized SINR, the bits are more vulnerable to distortion, which reduces the likelihood of their successful reception. It is worth mentioning that with the utilization of only 27 % of the maximal energy per bit, the BER performance approaches the result with 100% of energy per bit. Besides, the errorfree limit (BER=10-9) is exceeded when more than 4 or 8 nodes are actives with 5% or 100% of the energy per bit, respectively.
On the other hand, a forward error correction (FEC) limit of 10-5 is exceeded when more than 8 or 14 optical nodes become actives with 5% or 100% of the energy per bit, respectively. However, the FEC scheme presents high energy consumption of 250 pJ/b with low energy efficiency, so that it is difficult to perform at high bit rates . In this context, further studies including FEC schemes are necessary in order to evaluate the global energy efficiency impact on the OCDMA system.
Fig. 9 presents the BER performance against the percentage of transmitting energy per bit for the same scenario adopted in Fig. 8. We show the BER with DPCA power control for the number of active nodes of 5, 10 and 15. As expected, the low percentage of energy per bit results in a high BER, mainly for higher number of active nodes. Furthermore, one can control the optical BER level performance, according the system requirement, by adjusting the percentage of energy per bit transmitted. Thus, there is a trade-off between transmitting energy and BER. For example, with 5 active nodes it is possible to obtain a BER level of 109 with 9% of the energy per, which results in 91% of reduction in the energy consumption when compared with a system without any energy control mechanism. Furthermore, the general behavior of this figure suggests that it is possible to save up to 70% of energy per bit transmitted with the penalty of one order of magnitude of BER.
In this work we have investigated the utilization of distributed power control algorithm (DPCA) as mechanism to increase the energy efficiency in OCDMA access network with non-coherent 2-D codes. In our proposed scheme, based on the classical DPCA, we have modeled the energy efficiency considering optical fiber transmission and network infrastructure as encoders, decoders and star coupler. This model confirms that the energy consumption of the network infrastructure is larger than the energy consumption of the transmission infrastructure for our numerical example. In the proposed scheme it is possible to obtain the energy efficiency according to the level of BER and the quality of service (QoS) requirements. The obtained results have shown economy of 70% in the transmitted energy per bit, at the penalty of one order of magnitude of BER.
In the proposed scheme it was possible to define the energy efficiency according to the BER level accomplished the quality of service (QoS) requirements.
 K. Hinton, J. Baliga, M. Z. Feng, R.W.A. Ayre and R. S. Tucker, "Power consumption and energy efficiency in the internet". IEEE Network, vol. 25, no. 2, pp. 612, March/April 2011. [ Links ]
 S. Yoo, "Energy Efficiency in the Future Internet: the Role of Optical Packet Switching and Optical Label Switching", IEEE J Selected Topics in Quantum Electronics, vol. 17, no. 2, pp. 406 418, March-April 2011. [ Links ]
 Rodney S. Tucker, "Green Optical Communications Part I: Energy Limitations in Transport", IEEE J Selected Topics in Quantum Electronics, vol. 17, no. 2, pp. 245 260, March-April 2011. [ Links ]
 Rodney S. Tucker, "Green Optical Communications Part II: Energy Limitations in Networks", IEEE J Selected Topics in Quantum Electronics, vol. 17, no. 2, pp. 261 274, March-April 2011. [ Links ]
 A. Stok and E. H. Sargent, "Lighting the local area: Optical codedivision multipleaccess and quality of service provisioning," IEEE Network, vol. 14, no. 6, pp. 4246, Nov./Dec. 2000. [ Links ]
 G.C. Yang and W.C. Kwong, Prime codes with applications to CDMA optical and wireless networks, Artech House, Boston, MA, 2002. [ Links ]
 K. Fouli and M. Maier, "Ocdma and optical coding: Principles, applications, and challenges," IEEE Communications Magazine, vol. 45, no. 8, pp. 2734, August 2007. [ Links ]
 H. Yin and D. J. Richardson, Optical code division multiple access communication networks: theory and applications. Berlin: SpringerVerlag and Tsinghua University Press, 2009. [ Links ]
 F R. Durand, L. Galdino, L. H. Bonani, F. R. Barbosa, M. L. F. Abbade and Edson Moschim, "The Effects of Polarization Mode Dispersion on 2D WavelengthHopping Time Spreading Code Routed Networks", Photonics Network Communications, vol. 20, no. 1, pp. 27 32, Aug. 2010. DOI 10.1007/s1110701002426. [ Links ]
 N. Tarhuni, T. Korhonen, M. Elmusrati and E. Mutafungwa, "Power Control of Optical CDMA Star Networks", Optics Communications, vol. 259, pp. 655 664, Mar. 2006. [ Links ]
 E. Inaty, H. Shalaby, P. Fortie, and L. Rusch, Multirate "Optical Fast Frequency Hopping CDMA System Using Power Control", Journal of Lightwave Tech., vol. 20, n. 2, pp.166 177, March 2003. [ Links ]
 F. R. Durand and T. Abrão, "Distributed SNIR Optimization Based on the Verhulst Model in Optical Code Path Routed Networks With Physical Constraints", Journal of Optical Communications and Networking, vol. 3, no. 9, pp. 683691, Sep. 2011. doi:10.1364/JOCN.3.000683 [ Links ]
 F. R. Durand, M. S. Filho and T. Abrão, "The effects of power control on the optical CDMA random access protocol", Optical Switching and Networking, (In press) doi:10.1016/j.osn.2011.06.002 [ Links ]
 T. J. Gross, T. Abrão and P. J. E. Jeszensky. "Distributed power control algorithm for multiple access systems based on Verhulst model". AEU -- International Journal of Electronics and Communications. vol. 65, n.4, April 2011, pp. 361372. [ Links ]
 Y. Huang, V. Baby, I. Glesk, C. Bres, C. Greiner, D. Iazikov, T. Mossberg, and P. Prucnal, "Novel multicodeprocessing platform for wavelength-hopping timespreading optical cdma: A path to device miniaturization and enhanced network functionality," IEEE Journal of Selected Topics in Quantum Electronics, vol. 13, no. 5, pp. 14711479, Sept./Oct. 2007. [ Links ]
 C. Brès and P. Prucnal, "Code-empowered lightwave networks," Journal of Lightwave Technology, vol. 25, no. 10, pp. 29112921, Oct. 2007. [ Links ]
 K. O. Hill and G. Meltz, "Fiber bragg grating technology fundamentals and overview", Journal of Lightwave Technology, vol. 15, no. 8, pp. 12631276, Aug. 1997. [ Links ]
 D. Goodman and Narayan Mandayan "Power control for wireless data", IEEE Personal Communications, vol. 7, no. 2, pp. 48 54, April 2000. [ Links ]
 G. Foschini and Z. Miljanic, "A simple distributed autonomous power control algorithm and its convergence", IEEE Trans. Veh. Technol., vol. 42, n. 4, pp. 641646, Nov. 1993. [ Links ]
 Z. Uykan and H. Koivo, "Sigmoidbasis nonlinear power-control algorithm for mobile radio systems", IEEE Trans. Veh. Technol. vol. 53, n.1, pp. 265271, Jan 2004. [ Links ]
Received 29 Aug. 2011; for review 2 Sept. 2011; accepted 24 Feb. 2012
The derivation of (18) is obtained considering that the optimum solution for the power allocation problem satisfies the following set of equations :
where . Now, if z (x) is an arbitrary odd function that has null value only at origin, equation (A.1) is the only equilibrium point of the following continuous time dynamical system :
This control scheme is general and constitutes a systematic form to describe the dynamical behavior of various DPCAs. The sigmoidal , the Foschini , and recently the Verhulst  algorithm can be seen as particular structures of the general form. In the interval t to t + αi, the secant line to the curve pi(t) is indistinguishable from the tangent line in t for small values of αi if pi(t) is a continuous curve. In this case, we can write the approximation:
for αi > 0, i = 1, ...,K, which is the wellknown Euler integration method. From the above exposition, the sequence pi[n] = pi(nαi), can be recursively obtained by:
for i = 1, ...,K. From definitions in (8) and (14) results:
The DPCA developed by Foschini and Miljanic  and illustrated by (18), can be stated defining: zi (x)= x , αi= α, i = 1, ...,K.