Advanced
Fuzzy Adaptive Modified PSO-Algorithm Assisted to Design of Photonic Crystal Fiber Raman Amplifier
Fuzzy Adaptive Modified PSO-Algorithm Assisted to Design of Photonic Crystal Fiber Raman Amplifier
Journal of the Optical Society of Korea. 2013. Jun, 17(3): 237-241
Copyright ©2013, Optical Society of Korea
  • Received : January 01, 2013
  • Accepted : April 04, 2013
  • Published : June 25, 2013
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Majid Akhlaghi
Young Researchers and Elite Club, Omidieh Branch, Islamic Azad University, Omidieh, Iran
m.akhlaghi@sutech.ac.ir
Farzin Emami
Department of Opto Electronic, Shiraz University of Technology, Shiraz, Iran
Abstract
This paper presents an efficient evolutionary method to optimize the gain ripple of multi-pumps photonic crystal fiber Raman amplifier using the Fuzzy Adaptive Modified PSO (FAMPSO) algorithm. The original PSO has difficulties in premature convergence, performance and the diversity loss in optimization as well as appropriate tuning of its parameters. The feasibility and effectiveness of the proposed hybrid algorithm is demonstrated and results are compared with the PSO algorithm. It is shown that FAMPSO has a high quality solution, superior convergence characteristics and shorter computation time.
Keywords
I. INTRODUCTION
Nowadays a Raman amplifier is considered as a key component in realizing the wide-band amplification due to the low-noise operation and optical amplification in any wavelengths [1] . Over the past few years, photonic crystal waveguides have been widely studied, modeled and fabricated due to their peculiar properties such as endlessly single-moded [2] , high-nonlinearity [3 , 4] , overall controllable dispersion properties [5 , 6] , bandwidth improvement [7 , 8] and enhanced coupling properties [9] . A new fiber called the holey fiber is a category of photonic crystal fibers (PCFs) and has been used in Raman amplifiers as the gain medium [10 , 11] . Such fibers consist of a pure silica core surrounded by a regular array of longitudinal air holes and can offer tight modal confinement. Thus they can provide an effective nonlinearity per unit length and can have the same order or higher than a conventional fiber. Their Raman efficiency defined as the Raman gain coefficient can be divided by the pump effective area and is expected to be much higher than that of a standard single-mode fiber (SMF). In [12] , authors report the experimental study of the effect of a nonlinear PCF on the noise characteristics of a distributed Raman amplifier. There are two types of reported optimization work about the PCF RAs. In the first type, Raman gain is inherently flat or may possibly have optical amplification just by using one laser pump [13] ; whereas in the second type, because of the gain non-uniformity, gain flatness can achieve by application of many pumps [14] . In triangular PCFs, gain can be increased by varying the geometrical parameters in the fiber cross section such as hole distances, pitch wavelength and their ratio and higher GeO2 percentages [15] . The present study aimed to propose a more powerful algorithm called fuzzy adaptive modified particle swarm optimization (FAMPSO) to choose pump power and pump wavelength for a 10 Km long PC fiber with backward eight-pump configuration in order to create a flattened gain in the C and L bands. This method not only enhances the accuracy of design Raman amplifier parameters but also increases the reliability and speed of design parameters.
II. METHODS
Taking into account the loss, noise and Raman interaction, the optical power evaluation along the waveguide (z-direction) is determined by [15] :
PPT Slide
Lager Image
Raman gain coefficient of PCF with different Pitch and Germanium mole fraction.
PPT Slide
Lager Image
PPT Slide
Lager Image
Where p ± ( z , λi ) are the forward or backward signal/pump powers at the wavelength, λi , α is the attenuated factor and n ± ( z , λi ) defines the forward or backward noise power which is due to the amplified Rayleigh backscattering and spontaneous Raman scattering effects. Raman gain γR of this amplifier is shown in Fig. 1 [15] .
- 2.1. Classical PSO
PSO is an optimization algorithm inspired by a habitation such as birds. This method is based on two points: artificial live and evolution. In the extension of PSO algorithm, the possible solutions of the desired optimization problem are regarded as non-bulky birds and any element of this colony called a particle. They can fly in an n-dimensional space and correct their trajectory based on the previous experiments as well as the neighboring particles. For particle number i, the displacement vector Xi is defined as:
PPT Slide
Lager Image
Where S is the searching space. This particle consists of the velocity vector Vi in the form of:
PPT Slide
Lager Image
The best former position Pi in the search space of S , is defined as:
PPT Slide
Lager Image
Using these definitions, we can find the new positions by the aid of some weighting and learning factors which affect the algorithm searching procedure severely [16] .
PPT Slide
Lager Image
PPT Slide
Lager Image
Where Vik is velocity of particle i at iteration k, w is inertia weight factor, c 1 , c 2 are acceleration constants, xik is the position of particle i at iteration k, Pik is the best position of particle i until iteration k, Pgk is the best position of group until iteration k and finally r1, r2 are random numbers between 0 and 1. There are three tuning parameters; c 1 , c 2 and w which can control the algorithm behavior strongly. We can use w to affect the former velocity on the next velocity. Increasing w causes an increment in searching ability and decreasing w can continue an accurate searching in the previously scanned areas. So, a proper selection provides a compromise between the local and the global searches. In many works, the selected w was big at first and after an initial search, this value would be lowered; a linearized model [17] . c 1 and c 2 are the best private and global positions, respectively. Since c 1 is the amount of particle trust in the previous experiments, it is called cognitive parameter, whereas c 2 describes the trust amount on the group and hence it is entitled the social parameter. For c 1 > c 2 the particle is attracted to the best private position and the reverse is true for the other case. Many optimization procedures utilize c 1 = c 2 = 2. To improve the PSO convergency, we reform and modify this algorithm by combining three algorithms; modified PSO, adaptive PSO and fuzzy PSO.
- 2.2. Modified PSO (MPSO)
One of the disadvantages of the original PSO is falling in a local optimum point, so in order to improve the convergence property and accuracy of the PSO algorithm, this paper presented a modified particle swarm optimization method. Premature convergence can take place under different situations such as the population has converged to local optima, the population has lost its diversity, and the search algorithm has proceeded slowly or hasn’t proceeded at all. Mutation is a powerful strategy to increase the population diversity and improve the PSOs performance [18] . This paper employed a new mutation operator described as follows: In each iteration, three vectors are selected randomly ( Z1 , Z2 , Z3 ) from the initial population after these vectors are mutated. It is worthwhile noting that for seeking the search space uniformly, the selected vectors should be different ( Z1 Z2 Z3 ). A mutant vector ( Xmut ) is generated as: [19] :
PPT Slide
Lager Image
Where β * is the mutation constant with the value of: 0.93< β * <1. The target vector is mixed with the mutated vector, using the following scheme to yield the trial vector:
PPT Slide
Lager Image
Usual crossover is selected in the range of: 0.1 < crossover < 0.9. In this paper β * = 0.93 and crossover = 0.1 [16] .
- 2.3. Adjusting the Learning Coefficients
Learning coefficients, C 1 and C 2 , are used for determination of the local and overall best response in relation to particle rates, respectively. C 1 and C 2 in conventional PSO are considered as the constant value or change linearly during each iteration. The results show that these characteristics cannot obtain optimal results in some cases. So in this paper, these coefficients are computed as follows [17] .
PPT Slide
Lager Image
We chose n = 2 and G 0 as the G_best_value in the first iteration.
- 2.4. Fuzzy Formulation for Tuning of w
Weight coefficient determines the effect of the past velocity experience on the present one. The right choice of this parameter results in the balance between global and local search of PSO algorithm. Great weight coefficient enhances the global search performance and small weight coefficient improves the local search. The parameter w is often held constant or linearly changed for the entire run of a PSO but this method cannot obtain the global optima in many cases. The best choice for adjusting the computed weight coefficient is based on changes in G_best_value and in this regard the fuzzy logic method is applied. The inputs of the fuzzy method are NFV and w while the output of this method is Δw which are computed as follows [20] :
PPT Slide
Lager Image
In Eq. (11) FV is a current best performance evaluation and the calculated value of FV from Eq. 11 at the first iteration may be used as FVmin for the next iterations. FVmax is a very large value which is greater than any feasible solution. Usual inertia weight value is about 0/4 < w < 0/9 . In the fuzzy set, both positive and negative corrections are essential for the inertia weight. Therefore, for the inertia weight correction, a range of _0.1 to 0.1 has been selected [20] .
PPT Slide
Lager Image
For achieving the optimal point, objective functions should be described by membership functions. In this paper, for simplicity, all the membership functions are modeled with triangular shape and also they can take one of the three following status: Small (S), Medium (M) and Large (L). In Table 1 , output variables of fuzzy sets are presented in three forms of linguistic values; NE (Negative), ZE (Zero), and PE (Positive) with related membership functions, as shown in Fig. 2 . Rules of fuzzy sets are shown in Table 1 . There are nine possible rules for two input variables and three linguistic values for all input variables. Each fuzzy rule as from an “IF-THEN” statement like the following rule, for example: If NFV is ( L ) and w is ( M ) then Δw is ( ZE ). The fuzzy rules are shown in Table 1 used to modify the inertia weight correction ( Δw ). Each rule represents a relation between the input and output spaces. The membership functions are depicted in Fig. 2 are utilized for fuzzifying each input and output variable.
Based on the above theoretical description, we simulated our optimization method by using FAMPSO as shown in the flowchart of Fig. 3 .
Fuzzy rules of the input and output variables
PPT Slide
Lager Image
Fuzzy rules of the input and output variables
PPT Slide
Lager Image
The membership functions.
PPT Slide
Lager Image
Flowchart of FAMPSO algorithm.
III. NUMERICAL RESULTS
Consider a WDM system with 16 -channels. There is a grid of 3.75 nm with signal powers of about 20 mW. In this simulation we assumed that the lossy fiber has a length of about 10 Km. Runge-Kutta and a shooting method is used to solve the pump and signal interaction equations [21] . Our goal is to minimize the gain ripple in the C-L band by utilizing of the FAMPSO method to introduce optimizing sixteen parameters; the wavelengths and power levels for the 8 pumps range are from 1420 nm - 1520 nm and 0 mW - 120 mW respectively. For Optimizing gain ripple, the FAMPSO method should be minimizing the following function:
PPT Slide
Lager Image
In this equation is a predefined average gain for which minimization must be done around it. In this work, the particle numbers are three times of the optimization parameters. Initially any increases in the number of particles can improve the convergence speed, but accuracy and convergence speed will be reduced for more iterations. Application of the FAMPSO algorithm to optimize the Raman gain of a PCFRA is plotted in Fig. 4 for a system with eight backward pumps. As shown, variations of the gains restricted to a range of about 0.1 dB, which is very much smaller than other reported methods such as the genetic algorithm [22] and PSO algorithm [23]. The optimization results for eight-pump PCFRA are shown in Table 2 . As shown in this Table we can design a
PPT Slide
Lager Image
Optimized Raman gain for eight backward pump using FAMPSO.
Optimization results for eight-pumps of PCFRA
PPT Slide
Lager Image
Optimization results for eight-pumps of PCFRA
uniform On-Off Raman gain with eight pumps for which the maximum powers is less than 120 mW. Fig. 5 shows the gain ripples versus the number of iterations for the standard PSO and the proposed FAMPSO method. Indeed, the time consumption of the proposed method comparing with the classical PSO is shown in this figure. As shown in this figure, FAMPSO algorithm is performing better than PSO in terms of convergence speed and accuracy. Finally, Fig. 6
PPT Slide
Lager Image
Gain ripples versus the number of iteration for the standard PSO and the proposed FAMPSO.
PPT Slide
Lager Image
Gain ripples versus the number of the independent algorithm running in (a) FAMPSO method, (b) PSO method.
shows the development of the minimum gain ripples for 20 independent algorithms running which are obtained by PSO and FAMPSO algorithms. As it is seen, the FAMPSO method is more reliable than the classical PSO.
IV. CONCLUSION
We simulated the optimum solution for a system of RA with eight backward pumps, for photonic crystal fiber we used a new optimization algorithm named FAMPSO. Designing a PCF-RA with a flattened gain done in the C and L band and the pump powers and the pump wavelengths were calculated so that a gain ripple of about 0.1dB reached. This method not only enhances the accuracy of design Raman amplifier parameters but also increases the reliability and speed of design parameters.
References
Wang D. , Zou N. , Kang Z. , Liu J. , Namihira Y. , Zhang Y. 2011 “Design and analysis of gain-flattened raman amplifierswith novel highly nonlinear photonic crystal fibers” in Proc. 7th Int. Conf. on WiCOM shanghai, China 1 - 4
Birks T. A. , Knight J. C. , Russell P. S. J. (1997) “Endlesslysingle-mode photonic crystal fiber” Opt. Lett. 22 961 - 963
Yusoff Z. , Lee J. H. , Belardi W. , Monro T. M. , Teh P. C. , Richardson D. J. (2002) “Raman effects in a highly nonlinear holey fiber: amplification and modulation” Opt. Lett. 27 424 - 426
Saitoh K. , Koshiba M. (2004) “Highly nonlinear dispersionflattened photonic crystal fibers for supercontinuum generation in a telecommunication window” Opt. Express 12 2027 - 2032
Saitoh K. , Koshiba M. , Hasegawa T. , Sasaoka E. (2003) “Chromatic dispersion control in photonic crystal fibers: application to ultra-flattened dispersion” Opt. Express 11 843 - 852
Sinha R. K. , Varshney S. K. (2003) “Dispersion properties of photonic crystal fibers” Microw. Opt. Technol. Lett. 37 129 - 132
Danaie M. , Kaatuzian H. (2011) “Bandwidth improvement for a photonic crystal optical Y-splitter” J. Opt. Soc. Korea 15 283 - 288
Aghababaeian H. , Vadjed-Samiei M. H. , Granpayeh N. (2011) “Temperature stabilization of group index in silicon slotted photonic crystal waveguides” J. Opt. Soc. Korea 15 398 - 402
Bozorgi M. , Granpayeh N. (2011) “Directional emission from photonic crystal waveguide output by terminating with CROW and employing the PSO algorithm” J. Opt. Soc. Korea 15 187 - 195
de Matos C. J. S. , Hansen K. P. , Taylor J. R. (2003) “Experimental characterization of Raman gain efficiency of holey fiber” Electron. Lett. 5 424 - 425
Zhao C. L. , Li Z. , Yang X. , Lu C. , Jin W. , Demokan M. S. (2005) “Effect of a nonlinear photonic crystal fiber on the noise characterization of a distributed Raman amplifier” IEEE Photon. Technol. Lett. 3 561 - 563
Varshney S. K. , Saitoh K. , Fujisawa T. , Koshiba M. 2006 “Design of gain-flattened highly nonlinear photonic crystal fiber Raman amplifier using a single pump: a leakage loss approach” in Proc. OFC Conference Anaheim, California, USA paper OWD4.
Chang J. , Wang T. , Tao Z. 2009 “Optimal design of flattenedgain spectrum of Raman amplifiers based on photonic crystal fibers” in Proc. Photonics and Optoelectronics, SOPO Symposium Shanghai, China 1 - 4
Bottacini M. , Poli F. , Cucinotta A. , Selleri S. (2004) “Modeling of photonic crystal fiber Raman amplifiers” J. Lightwave Technol. 7 1707 - 1713
Emami F. , Akhlaghi M. (2012) “Gain ripple decrement of S-band Raman amplifier” IEEE Photon. Technol. Lett. 24 1349 - 1352
Chaoli S. , Jianchao Z. , Jeng P. (2011) “A modified particle swarm optimization with feasibility-based rules for mixedvariable optimization problems” Int. J. Innovative Computing, Information and Control 7 3081 - 3096
Niknam T. , Mojarrad H. , Nayeripour M. (2011) “A new fuzzyadaptive particle swarm optimization for non-convex economicdispatch” Int. J. Innovative Computing, Information and Control 7 1764 - 1778
Niknam T. , Mojarrad H. D. , Meymand H. Z. (2011) “A novel hybrid particle swarm optimization for economic dispatch with valve-point loading effects” Applied Energy Conversion and Management 52 1800 - 1809
Bajpai P. , Singh S. N. (2007) “Fuzzy adaptive particle swarm optimization for biddinrg strategy in uniform price spot market” IEEE Trans. Power Sys. 22 2152 - 2160
Dashtestani Z. L. , Korushavi F. , Rahmani M. H. (2008) “Anefficient shooting method for fiber amplifiers and lasers” J. Optics & Laser Technology 40 1041 - 1046
Ferreira G. C. M. , Cani S. P. N. , Pontes M. J. , Segatto M. E. V. (2011) “Optimization of distributed Raman amplifiers using a hybrid genetic algorithm with geometric compensation technique” IEEE J. Photonics 3 390 - 399
Jiang H. M , Xie K. , Wang Y. F (2010) “Shooting algorithm and particle swarm optimization based Raman fiber amplifiers” Opt. Commun. 283 3348 - 3352