Search

PurposeThe purpose of this paper is to demonstrate an effective and robust numerical solution for Raman fiber amplifier (RFA) equations which have no explicit solution. MATLAB BVP solvers are addressed for the solution.Design/methodology/approachThe continuation method proposed for the solution of RFA equations using MATLAB BVP solvers is explained. Scripts for improving the power values at the boundaries with continuation, extending fiber length with continuation and calculation of the analytical partial derivatives using the MATLAB Symbolic toolbox are introduced. Comparisons among the different MATLAB BVP solvers have been made. Using the continuation method, signal evolutions for different kinds of RFA amplifier configurations are plotted.FindingsThe paper finds that MATLAB BVP solver with the continuation method can be used in the design of various kinds of RFAs for high powers/long gain fiber spans.Research limitations/implicationsThe paper will assist the fiber optic research community who suffer from two or more point boundary‐value problems. Moreover, the stiffness of the signal evolution which is faced with high pump powers and/or long fiber lengths can be solved with continuation. This superiority of the solver can be used to overcome any stiff changes of the signals for future studies.Practical implicationsThe increased research interests and practical demands for RFAs have been calling for reasonable and efficient means for the performance evaluation of RFAs before the real amplifiers are fabricated. The solution method presented in this paper will be an efficient means for the solution of this issue.Originality/valueMATLAB BVP solvers have been proven to be effective for the numerical solution of RFAs with multiple pumps and signal waves. Using the continuation method, in a distributed RFA with ten pump sources, 2,400 mW total input pump power is achieved. The improvement of the total power is about 1.4 times compared with those of the previously reported methods. Using the MATLAB BVP solvers, total power/fiber span can be improved further using the continuation process with the cost of computational time. This is a notable and promising improvement from a RFA designer's point of view.

PurposeThe purpose of this paper is to demonstrate an effective and faster numerical solution for nonlinear‐coupled differential equations describing fiber amplifiers which have no explicit solution. MATLAB boundary value problem (BVP) solver of bvp6c function is addressed for the solution.Design/methodology/approachCoding method with the bvp6c is introduced, signal evolution, threshold calculation method is introduced, gain and noise figure are plotted and superiority of the bvp6c solver is compared with the Newton‐Raphson method.Findingsbvp6c function appears to be an effective tool for the solution fiber amplifier equations and can be used for different pump configurations of BFAs and RFAs. The excellent agreement between the proposed and reported results shows the reliability of the proposed threshold power calculation method.Research limitations/implicationsThe paper eases the work of the fiber optic research community, who suffer from two point BVPs. Moreover, the stiffness of the signal evolution which is faced with high pump powers and/or long fiber lengths can be solved with continuation. This superiority of the solver can be used to overcome any stiff changes of the signals for the future studies.Practical implicationsThe main outcome of this paper is the numerically calculation of the threshold values of fiber amplifiers without the necessity of the experiment. The robustness improvement of the solution is that the solver is able to solve the equations even with the poor guess values and the solution can be obtained without the necessity of analytical Jacobian matrix.Originality/valueMATLAB bvp6c solver has proven to be effective for the numerical solution of nonlinear‐coupled intensity differential equations describing fiber amplifiers with two‐point boundary values. Beside the signal evolution, thresholds of Brillouin and Raman fiber amplifiers can also be calculated by using the proposed solver. This is a notable and promising improvement of the paper, at least from a fiber optic amplifier designer point of view.