# Patent

## Title
Predictive optimization of wavelength division multiplexed systems

## Abstract

The performance of a WDM system is improved by adjusting launch powers of each channel responsive to a prediction of system performance degradation. The prediction of system performance degradation could relate to the effect of Stimulated Raman Scattering (SRS) in which case, given launch power, a value representative of the effect of SRS is estimated for each channel in a WDM system using a method described herein. Using this SRS effect estimate, a value is computed for the estimated output power of each channel. These estimated output power values are then used to refine the estimate of the SRS effect on each channel. Subsequently, the refined SRS effect estimates are used to refine the estimated output power values. These steps are repeated and eventually an estimate of SRS effect on each wavelength channel is settled upon that is as accurate as that which results from other, more computationally complex, methods.

## Background
FIELD OF THE INVENTION
The present invention relates to communication over optical fibers and, more particularly, to predictive optimization of wavelength division multiplexed systems.
BACKGROUND OF THE INVENTION
Optical networks, which use optical fiber as a transmission medium, are well known for use in transmission of digital data signals. To make efficient use of a single optical fiber, many unique data signals may be transmitted over the same fiber so long as each data signal modulates an optical signal with a wavelength different from the other optical signals on the same fiber. This transmission scheme is called Wavelength Division Multiplexing (WDM). When the wavelengths of the different optical signals are only marginally different from one another, the transmission scheme may be called Dense Wavelength Division Multiplexing (DWDM). In a network using DWDM, two elements connected by a single physical link (optical fiber) may communicate using a number of signal channels, where each signal channel is an optical signal with a distinct wavelength. However, each signal channel may not be independent of the other signal channels transmitted on the same fiber, even though the other signal channels are transmitted at different wavelengths. An effect called Stimulated Raman Scattering (SRS) may allow crosstalk from one or more signal channels to cause changes (amplification or attenuation) to occur in a signal channel of interest.
SRS is an important non-linear effect in optical fibers named for the scientist who discovered a phenomenon in the scattering of light, called the Raman Effect, in 1928. Raman scattering occurs in a silica fiber when light propagates through the fiber. Stimulated Raman Scattering is an inelastic scattering process in which an incident photon loses its energy to create another photon of reduced energy at a lower frequency. That is, energy of a first signal at one wavelength may effect a second signal at a longer wavelength. SRS is often portrayed as a plot of Raman gain, gR, as a function of frequency (or wavelength). In such a plot, the magnitude of the Raman gain may be shown to vary over a Raman bandwidth (i.e., the bandwidth over which gR is non-zero).
As WDM systems move towards densely packed wavelengths in DWDM and related high channel count, SRS becomes a critical factor in system performance. Accordingly, modeling SRS accurately has been a research topic recently.
It has been shown that a WDM system is optimized when all channels have equal power or, alternatively, Optical Signal to Noise Ratio (OSNR) at the output of the fiber. Several schemes have been proposed for equalizing channel output power. In one scheme, disclosed by Chraplyvy in U.S. Pat. No. 5,225,922, issued Jul. 6, 1993, the output power of each channel is measured at the receiver. Based on these measurements, adjustments are made to the optical input signal power for each channel. In another scheme, disclosed by Khaleghi, et al, in U.S. Pat. No. 6,040,933, issued Mar. 21, 2000, optical power measurements of each channel are taken at inputs of optical amplifiers along the transmission path of the system. The channels are equalized by adjusting the optical input signal powers of each channel based on the measurements. In each of these schemes, there is a requirement for accurate power measurement and a feedback path over which instructions may be received for making adjustments based on the measurements. However, neither scheme may be used in the absence of a feedback path. Further, accurate optical power measurement generally involves expensive hardware. It would be desirable, then, to optimize a WDM system by equalizing output power for each channel in the system even in the absence of a feedback path.
SUMMARY OF THE INVENTION
A prediction of system performance degradation is determined. Based on this prediction, output power for each channel in a WDM system is estimated and, based on the estimated output power, the launch power of each channel may be adjusted. The prediction of system performance degradation could relate to the effect of Stimulated Raman Scattering (SRS) in which case, given launch power, a value representative of the effect of SRS is estimated for each channel in a WDM system using a method described herein. Using this SRS effect estimate, an intermediate value is computed for the output power of each channel. These intermediate output power values are then used to refine the estimate of the SRS effect on each channel. Subsequently, the refined SRS effect estimates are used to refine the intermediate power values. These steps are repeated and eventually an estimate of SRS effect on each channel is settled upon that is within engineering tolerance of that which results from other, more computationally complex, methods. This estimate of SRS effect on each channel may then be used to optimize the WDM system.
Advantageously, the method quantifies SRS accurately yet involves low computational complexity and is especially well suited to practical implementation.
In accordance with an aspect of the present invention there is provided a method for improving the performance of a wavelength division multiplexed system, the wavelength division multiplexed system including a plurality of wavelength channels. The method includes determining a prediction of system performance degradation, based on the prediction, determining an estimated output power for each wavelength channel of the plurality of wavelength channels and, based on the estimated output power, generating a control signal, the control signal for adjusting a launch power for each wavelength channel. In another aspect of the invention a system control unit is provided for performing this method. In a further aspect of the present invention, there is provided a software medium that permits a general purpose computer to carry out this method.
In accordance with another aspect of the present invention there is provided an optical transmission system including a plurality of channel transmitters, each for generating a corresponding wavelength channel, a multiplexer for combining the wavelength channels into a wavelength multiplexed signal sent on an optical fiber and a system control unit for: determining a prediction of system performance degradation; based on the prediction, determining an estimated output power for each wavelength channel at an output of the optical fiber; and, based on the estimated output power, generating a control signal, the control signal for transmitting to the channel transmitter for adjusting a launch power for each wavelength channel.
In accordance with a further aspect of the present invention there is provided a method for estimating an effect of Stimulated Raman Scattering (SRS) in a wavelength division multiplexed system. The method includes determining an estimate of the effect of SRS on each wavelength channel in the wavelength division multiplexed system based on an indication of launch power, determining a plurality of estimated output power values, one of the plurality of estimated output power values corresponding to each wavelength channel, based on the estimate of the effect of SRS on each wavelength channel and improving the estimate of the effect of SRS on each wavelength channel based on the plurality of estimated output power values to result in a plurality of improved SRS estimates. In another aspect of the invention a system control unit is provided for performing this method. In a further aspect of the present invention, there is provided a software medium that permits a general purpose computer to carry out this method.
Other aspects and features of the present invention will become apparent to those of ordinary skill in the art upon review of the following description of specific embodiments of the invention in conjunction with the accompanying figures.



BRIEF DESCRIPTION OF THE DRAWINGS
In the figures which illustrate example embodiments of this invention:
FIG. 1 schematically illustrates a DWDM system;
FIG. 2 illustrates steps of a WDM system performance improving method in accordance with an embodiment of the present invention;
FIG. 3 illustrates steps of an SRS effect estimating method in accordance with an embodiment of the present invention;
FIG. 4 is a plot comparing the iterative method of the present invention to an analytical method for 60 channels; and
FIG. 5 is a plot comparing the iterative method of the present invention to an analytical method for 100 channels.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS
In FIG. 1, a DWDM system 100 includes a transmitting node 102, an optical fiber 104 and a receiving node 106. The transmitting node 102 further includes N channel transmitters 110A, 110B, 110C, . . . , 110N (generally referred to as a channel transmitter 110) that each provide a wavelength channel to a multiplexer 108. The multiplexer 108 multiplexes the channels received from the channel transmitters 110 and sends the multiplexed signal to the receiving node 106 over the fiber 104. A launch power given to each channel by each channel transmitter 110 is controlled by a system control unit 112. As such, each channel transmitter 110 may include a laser source and an attenuator/amplifier (not shown) to adjust the launch power for the respective channel. The system control unit 112 may be loaded with transmitter control software for executing methods exemplary of this invention from software medium 114 which could be a disk, a tape, a chip or a random access memory containing a file downloaded from a remote source.
In overview, operation of the system control unit 112 may be considered in view of FIG. 2. Initially, a launch power indication is received from each channel transmitter 110 (step 202). Based on the launch power, a prediction of performance degradation of each channel in the DWDM system 100 is determined (step 204). Based on this prediction, output power for each channel is estimated (step 206). It is then determined whether the estimated output power for each channel is substantially the same (step 208). If not, the launch power of each channel is adjusted (step 210) and the adjusted launch powers are provided to the system control unit 112 (step 202). The method is complete when the estimated output power for each channel, as determined in step 208, is substantially the same. This adjustment serves to optimize the DWDM system 100.
In particular, the prediction of system performance degradation could relate to the effect of Stimulated Raman Scattering (SRS). There is, then, a requirement for a method to be used in determining an estimated value representative of the effect of SRS for each channel in the WDM system. The SRS effect is determined using an equation called the coupled-equation. The coupled-equation may be expressed as 



  
    
      
        
          
            
              d
              
                
                  P
                  i
                
                ⁡
                
                  (
                  z
                  )
                
              
            
            
              d
              z
            
          
          =
          
            
              
                -
                
                  α
                  i
                
              
              ⁢
              
                
                  P
                  i
                
                ⁡
                
                  (
                  z
                  )
                
              
            
            +
            
              
                
                  P
                  i
                
                ⁡
                
                  (
                  z
                  )
                
              
              ⁢
              
                
                  ∑
                  
                    j
                    =
                    1
                  
                  
                    i
                    -
                    1
                  
                
                ⁢
                
                   
                
                ⁢
                
                  
                    g
                    ij
                  
                  ⁢
                  
                    
                      
                        P
                        j
                      
                      ⁡
                      
                        (
                        z
                        )
                      
                    
                    
                      
                        K
                        ij
                      
                      ⁢
                      
                        A
                        ij
                      
                    
                  
                
              
            
            +
            
              
                
                  P
                  i
                
                ⁡
                
                  (
                  z
                  )
                
              
              ⁢
              
                
                  ∑
                  
                    j
                    =
                    
                      i
                      +
                      1
                    
                  
                  N
                
                ⁢
                
                   
                
                ⁢
                
                  
                    g
                    ij
                  
                  ⁢
                  
                    
                      
                        λ
                        j
                      
                      ⁢
                      
                        
                          P
                          j
                        
                        ⁡
                        
                          (
                          z
                          )
                        
                      
                    
                    
                      
                        λ
                        i
                      
                      ⁢
                      
                        K
                        ij
                      
                      ⁢
                      
                        A
                        ij
                      
                    
                  
                
              
            
          
        
        ,
      
    
    
      
        (
        1
        )
      
    
  




where Pi(z) is the power in the ith channel as a function of propagation distance, z, αi is a fiber attenuation coefficient for the ith channel and is usually considered to be the same (α) for all channels. gij is a Raman gain coefficient between channel i and channel j that is calculated as gij=g′(fj−fi), where 



  
    g
    ′
  
  =
  
    
      d
      
        g
        R
      
    
    
      d
      f
    
  




and represents the slope of the Raman gain at frequency fi. The Raman gain coefficient, gij, is usually considered to be the same for all channels. Aij is the effective cross sectional area of the fiber and is usually considered to be the same (i.e., Ae) for all channels. 1/Kij is a polarization factor and is usually considered to be 0.5 for all channels. λi/λj is a wavelength ratio between channel i and channel j and is usually considered to be unity for all channels.

Estimating the SRS effect using equation 1 is a mathematical challenge. A common approach to simplify the problem is to take the so-called “undepleted-pump” assumption (see Govind P. Agrawal, “Nonlinear Fiber Optics,” Second Edition, Academic Press, 1995, hereby incorporated herein by reference) and solve the equation for two channels first and then extend the result to all channels. However, the result becomes unacceptably inaccurate as the channel count becomes high, as it does in DWDM systems. Recently, Christodoulides published a letter (Christodoulides, D. N. and Jander; R. B., “Evolution of Stimulated Raman Crosstalk in Wave-Length Division Multiplexed Systems,” IEEE Photonics Technology Letters, 1996, hereby incorporated herein by reference) detailing a closed form solution to the coupled-equation (equation 1).
The Christodoulides method solves the coupled-equation (equation 1) directly to obtain an accurate yet complex closed-form solution for the maximum crosstalk impairment, XT, of a channel as a function of propagation distance, z, given equal channel loading, 



  
    
      
        
          XT
          ⁡
          
            (
            z
            )
          
        
        =
        
          1
          -
          
            
              
                Ne
                
                  
                    -
                    
                      cNP
                      0
                    
                  
                  ⁢
                  
                    Z
                    e
                  
                  ⁢
                  
                    
                      (
                      
                        N
                        -
                        1
                      
                      )
                    
                    2
                  
                
              
              ⁡
              
                [
                
                  
                    sinh
                    ⁡
                    
                      (
                      
                        
                          
                            cNP
                            0
                          
                          ⁢
                          
                            Z
                            e
                          
                        
                        2
                      
                      )
                    
                  
                  
                    sinh
                    ⁡
                    
                      (
                      
                        
                          
                            cN
                            2
                          
                          ⁢
                          
                            P
                            0
                          
                          ⁢
                          
                            Z
                            e
                          
                        
                        2
                      
                      )
                    
                  
                
                ]
              
            
            .
          
        
      
    
    
      
        (
        2
        )
      
    
  




In equation 2, N is the number of channels and 



  
    
      
        c
        =
        
          
            
              g
              ′
            
            ⁢
            Δ
            ⁢
            
               
            
            ⁢
            f
          
          
            2
            ⁢
            
              A
              e
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




where Δf is the interchannel frequency spacing, Ae, is the effective cross sectional area of the fiber and Ze, is an effective propagation distance, given by 



  
    
      
        
          Z
          e
        
        =
        
          
            
              1
              -
              
                e
                
                  
                    -
                    α
                  
                  ⁢
                  
                     
                  
                  ⁢
                  z
                
              
            
            α
          
          .
        
      
    
    
      
        (
        4
        )
      
    
  




The result achieved using equation 2 is accurate, but very computationally intensive for practical implementation in embedded systems since it involves calculation of sinh, exponents, etc.
According to the following iterative method, considered in view of FIG. 3, the undepleted pump assumption is used to solve the coupled-equation (equation 1) and the solution is iteratively modified. After a few iterations, the solution from the iterative method quickly approaches the solution from the Christodoulides method, yet with a significantly reduced computational complexity.
The DWDM system 100 (FIG. 1) is assumed to be an equally-spaced N-channel DWDM system. Furthermore, channel 1 has the lowest wavelength, channel N has the highest wavelength and the wavelengths of the channels increase monotonically with reference number. For a generic ith channel, there will be an increase in power as each channel with a shorter wavelength transfers power to the ith channel due to SRS. Similarly, the ith channel transfers power to those channels with longer wavelengths due to SRS.
The amount of power transfer (crosstalk), x, over the length of the fiber 104 due to SRS between any two channels may be computed (step 302) using the undepleted pump assumption as: 



  
    
      
        
          x
          
            [
            1
            ]
          
        
        =
        
          
            
              g
              R
            
            2
          
          ⁢
          
            (
            
              Δλ
              
                Δλ
                R
              
            
            )
          
          ⁢
          
             
          
          ⁢
          
            (
            
              
                L
                e
              
              
                A
                e
              
            
            )
          
          ⁢
          
            P
            0
          
        
      
    
    
      
        (
        5
        )
      
    
  




where Δλ is wavelength spacing between the two channels, ΔλR is Raman bandwidth, Le is an effective length of the fiber 104 and P0 represents launch power. By default, launch power is equal for each channel. The superscript [1] denotes the result is from the 1st iteration. Notably, equation 5 is related to equation 3, in that 



  
    
      
        
          
            
              g
              ′
            
            ⁢
            Δ
            ⁢
            
               
            
            ⁢
            f
          
          
            2
            ⁢
            
              A
              e
            
          
        
        =
        
          
            
              g
              R
            
            2
          
          ⁢
          
            (
            
              Δλ
              
                Δλ
                R
              
            
            )
          
          ⁢
          
             
          
          ⁢
          
            
              (
              
                
                  L
                  e
                
                
                  A
                  e
                
              
              )
            
            .
          
        
      
    
    
      
        (
        6
        )
      
    
  




Equation 6 allows equation 5 to be simplified to

x[1]=cP0.  (7)

The power at the output of the ith channel in the N-channel system 100, given SRS, may then be estimated as follows (step 304): 



  
    
      
        
          P
          i
          
            [
            1
            ]
          
        
        =
        
          
            P
            0
          
          ⁢
          
            
              ∑
              
                k
                =
                1
              
              N
            
            ⁢
            
               
            
            ⁢
            
              [
              
                1
                -
                
                  
                    (
                    
                      k
                      -
                      i
                    
                    )
                  
                  ⁢
                  
                    x
                    
                      [
                      1
                      ]
                    
                  
                
              
              ]
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




For the second iteration, the power computed using equation 8, Pi[1], is substituted for P0 in the right-hand side of equation 5 and a second iteration crosstalk estimate, Xi [2], is computed for each channel i (step 306), leading to an improved crosstalk estimate for each channel, 



  
    
      
        
          x
          i
          
            [
            2
            ]
          
        
        =
        
          
            
              g
              R
            
            2
          
          ⁢
          
            (
            
              Δλ
              
                Δλ
                R
              
            
            )
          
          ⁢
          
             
          
          ⁢
          
            (
            
              
                L
                e
              
              
                A
                e
              
            
            )
          
          ⁢
          
            P
            i
            
              [
              1
              ]
            
          
        
      
    
    
      
        (
        9
        )
      
    
  
  
    
      
        
           
        
        ⁢
        
          =
          
            cP
            i
            
              [
              1
              ]
            
          
        
      
    
    
      
        (
        10
        )
      
    
  




Using the improved (and channel specific) crosstalk estimate, the output power for channel i may be refined (step 308) to 



  
    
      
        
          
            P
            i
            
              [
              2
              ]
            
          
          =
          
            
              P
              i
              
                [
                1
                ]
              
            
            ⁢
            
              
                ∑
                
                  k
                  =
                  1
                
                N
              
              ⁢
              
                 
              
              ⁢
              
                [
                
                  1
                  -
                  
                    
                      δ
                      ⁡
                      
                        (
                        
                          k
                          -
                          i
                        
                        )
                      
                    
                    ⁢
                    
                      x
                      k
                      
                        [
                        2
                        ]
                      
                    
                  
                
                ]
              
            
          
        
        ,
      
    
    
      
        (
        11
        )
      
    
  




where δ is a factor introduced to compensate for the fact that power transfer in each channel occurs over the length of the fiber. δ has an optimum range of between 0.575 and 0.625. A value of 0.6 has been used for δ with acceptable results.

Generally, then, for the mth iteration: 



  
    
      
        
          x
          i
          
            [
            m
            ]
          
        
        =
        
        ⁢
        
          
            
              g
              R
            
            2
          
          ⁢
          
            (
            
              Δλ
              
                Δλ
                R
              
            
            )
          
          ⁢
          
             
          
          ⁢
          
            (
            
              
                L
                e
              
              
                A
                e
              
            
            )
          
          ⁢
          
            P
            i
            
              [
              
                m
                -
                1
              
              ]
            
          
        
      
    
    
      
        
           
        
        ⁢
        
          (
          12
          )
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            cP
            i
            
              [
              
                m
                -
                1
              
              ]
            
          
          ⁢
          
             
          
          ⁢
          and
        
      
    
    
      
        (
        13
        )
      
    
  
  
    
      
        
          P
          i
          
            [
            m
            ]
          
        
        =
        
        ⁢
        
          
            P
            i
            
              [
              
                m
                -
                1
              
              ]
            
          
          ⁢
          
            
              ∑
              
                k
                =
                1
              
              N
            
            ⁢
            
               
            
            ⁢
            
              
                [
                
                  1
                  -
                  
                    
                      
                        δ
                        
                          m
                          -
                          1
                        
                      
                      ⁡
                      
                        (
                        
                          k
                          -
                          i
                        
                        )
                      
                    
                    ⁢
                    
                      x
                      k
                      
                        [
                        m
                        ]
                      
                    
                  
                
                ]
              
              .
            
          
        
      
    
    
      
        (
        14
        )
      
    
  




After a few iterations of steps 306 and 308 (about 5 or 6), it has been found that the estimate of output power (Pi) for each channel is almost identical to the output power that may be obtained using the Christodoulides method. The number of iterations used in practice may be preset as a parameter of a system, based on a required degree of accuracy, or iterations may continue until the difference between successive iterations is suitably small, i.e., less than a threshold (step 310).
In operation, the control system 112 (FIG. 1) sends a control signal to adjust the launch power for each channel to the computed Pi.
A comparison of this iterative method with the Christodoulides (analytical) method is illustrated in FIG. 4 for 60 channels and FIG. 5 for 100 channels.
Advantageously, the effect of SRS is quantified iteratively. Each iteration involves minimum amount of computation as the method uses addition and multiplication only (which are straightforward to implement). Further, c is a constant and can be computed off-line and computation-intensive operations (such as exponents, etc.), present in the Christodoulides method, have been avoided while obtaining good accuracy. This iterative method is particularly well suited to embedded systems.
Other modifications will be apparent to those skilled in the art and, therefore, the invention is defined in the claims.

## Claims

1. A method for improving the performance of a wavelength division multiplexed system, said wavelength division multiplexed system including a plurality of wavelength channels, said method comprising:
iteratively determining a prediction of system performance degradation based, at least in part, on a launch power for each wavelength channel of said plurality of wavelength channels, an effective fiber length and an effective fiber cross sectional area; 
based on said prediction, determining an estimated output power for said each wavelength channel of said plurality of wavelength channels; and 
based on said estimated output power, generating a control signal, said control signal for adjusting said launch power for said each wavelength channel.

2. The method of claim 1 wherein said adjusting acts to equalize said estimated output power for said each wavelength channel.

3. The method of claim 1 further comprising receiving an indication of said launch power for said each wavelength channel and basing said prediction of system performance degradation on said indication.

4. The method of claim 1 wherein said determining said estimated output power for each wavelenghth uses: 



  
    P
    i
    
      [
      m
      ]
    
  
  =
  
    
      P
      i
      
        [
        
          m
          -
          1
        
        ]
      
    
    ⁢
    
      
        ∑
        
          k
          =
          1
        
        N
      
      ⁢
      
         
      
      ⁢
      
        [
        
          1
          -
          
            
              
                δ
                
                  m
                  -
                  1
                
              
              ⁡
              
                (
                
                  k
                  -
                  i
                
                )
              
            
            ⁢
            
              
                g
                R
              
              2
            
            ⁢
            
              (
              
                Δλ
                
                  Δλ
                  R
                
              
              )
            
            ⁢
            
               
            
            ⁢
            
              (
              
                
                  L
                  e
                
                
                  A
                  e
                
              
              )
            
            ⁢
            
              P
              k
              
                [
                
                  m
                  -
                  1
                
                ]
              
            
          
        
        ]
      
    
  




where:
Pi[m] is the power in ith channel after an mth iteration; 
δ is a compensating factor; 
gR is Raman gain; 
Δλ is wavelength spacing; 
ΔλR is Raman bandwidth; 
Le is said effective fiber length; and 
Ae is said effective fiber cross sectional area.

5. The method of claim 4 wherein said performance degradation occurs due to Stimulated Roman Scattering (SRS) and said determining said prediction comprises determining an effect, x, of SRS on said each wavelength channel where 



  x
  =
  
    
      
        g
        R
      
      2
    
    ⁢
    
      (
      
        Δλ
        
          Δλ
          R
        
      
      )
    
    ⁢
    
       
    
    ⁢
    
      (
      
        
          L
          e
        
        
          A
          e
        
      
      )
    
    ⁢
    
      P
      0
    
  




and P0 is said launch power.

6. A method for improving the performance of a wavelength division multiplexed system, said wavelength division multiplexed system including a plurality of wavelength channels, said method comprising:
receiving an indication of an initial launch power for each wavelength channel of said plurality of wavelength channels; 
based on said indication of said initial launch power, determining a prediction of system performance degradation, wherein said performance degradation occurs due to Stimulated Raman Scattering (SRS) and said determining said prediction comprises includes determining an effect of SRS on said each wavelength channel; 
said determining said effect of SRS including: 
determining an estimate of said effect of SRS on said each wavelength channel based on said indication of said initial launch power; 
based on said estimate of said effect of SRS on said each wavelength channel, determining a plurality of estimated output power values, one of said plurality of estimated output power values corresponding to said each wavelength channel; and 
improving said estimate of said effect of SRS on said each wavelength channel based on said plurality of estimated output power values to result in a plurality of improved SRS estimates; 

based on said plurality of improved SRS estimates, determining a plurality of refined estimated output power values; and 
based on said plurality of refined estimated output power values, generating a control signal for said each wavelength channel, said control signal for adjusting said initial launch power for said each wavelength channel to result in an adjusted launch power.

7. The method of claim 6 further comprising further improving said plurality of improved SRS estimates based on said plurality of refined estimated output power values to result in a plurality of further improved SRS estimates.

8. The method of claim 7 further comprising further refining said plurality of refined estimated output power values based on said plurality of further improved SRS estimates to result in a plurality of further refined estimated output power values.

9. The method of claim 8 further comprising repeating said further improving and said further refining.

10. The method of claim 9 further comprising performing said repeating until a difference between successive refined estimated output power values corresponding to a given wavelength channel is less than a threshold.

11. The method of claim 4 where said determining said effect of SRS further comprises using the undepleted pump assumption.

12. An optical transmission system comprising:
a plurality of channel transmitters, each for generating a corresponding wavelength channel; 
a multiplexer for combining said wavelength channels into a wavelength multiplexed signal sent on an optical fiber; 
a system control unit for: 
iteratively determining a prediction of system performance degradation based, at least in part, on a launch power for each wavelength channel of said plurality of wavelength channels, an effective fiber length and an effective fiber cross sectional area; 
based on said prediction, determining an estimated output power for said each wavelength channel at an output of said optical fiber; and 
based on said estimated output power, generating a control signal, said control signal for transmitting to said channel transmitter for adjusting said launch power for said each wavelength channel.

13. A system control unit in a wavelength division multiplexed system operable to:
iteratively determine a prediction of system performance degradation based, at least in part, on a launch power for each wavelength channel of said plurality of wavelength channels, an effective fiber length and an effective fiber cross sectional area; 
based on said prediction, determine an estimated output power for said each wavelength channel of a plurality of wavelength channels in said wavelength division multiplexed system; and 
based on said estimated output power, generate a control signal, said control signal for adjusting said launch power for said each wavelength channel.

14. A computer readable medium containing computer-executable instructions which, when performed by a processor in a system control unit in a wavelength division multiplexed system, cause the processor to:
determine a prediction of system performance degradation; 
based on said prediction, determine an estimated output power for each wavelength channel of a plurality of wavelength channels in said wavelength division multiplexed system; 
improve said prediction of system performance degradation based on said estimated output power for said each wavelength channel to result in an improved prediction of system performance degradation; 
determine a revised estimated output power for said each wavelength channel based on said improved prediction of system performance degradation; and 
based on said revised estimated output power, generate a control signal, said control signal for adjusting a launch power for said each wavelength channel.