Hi,
I'm modeling rate equations describing an optical feedback effect back into laser cavity with a target in harmonic motion, to get its steady-state laser output. (lags are in harmonic motion) But the problem shows that dynamic models cannot provide the same tilt fringes same as steady-state equations derived from dynamic equations.
I use MATLAB solvers (dde23, ddesd) to numerically solve these dynamic equations to get its steady-state solutions. However, it fails to obtain the same simulated waveform as one I obtain from solving its steady-state equations directly.
To attempt to obtain its steady-state solutions from the dynamic models of time-dependent lags, I have tried MATLAB and Julia but fails:
dde23 of MATLAB: to update its initial conditions in a for loop for each constant lags;
ddesd of MATLAB: directly set time-dependent lags;
Julia's dde package: fails to get laser lasing;
For Mathematica, I find it doesn't support time-dependent lags till now;
For Python, Pydelay package use the same Bogacki-Shampine method in dde23;
Could anyone give any suggestion on this issue? I have spent a bunch of time to deal with this issue but no progress. TAT
Any suggestion is appreciated.
Thank you very much!
Xiao