I run tutorial files for optical bistability. However, the calculation result is different from it of website.
I think the code for removing high frequency components doesn’t work, which is
Ey_through_w = 2*((1:length(w))<=(length(w)/2+0.1))*fft(pinch(Ey),1,zero_pad); Ey_through_t = invfft(pinch(Ey_through_w)*filter1); Ey_through_t = Ey_through_t(1:length(t)); Esy_through_w = 2*((1:length(w))<=(length(w)/2+0.1))*fft(pinch(Esy),1,zero_pad); Esy_through_t = invfft(pinch(Esy_through_w)*filter1); Esy_through_t = Esy_through_t(1:length(t));
and the result of calculation of tutorial file is
, which is different from website result.
I also has another question.
In tutorial, a bistability property is considered in the structure composed of a waveguide and a resonator.
In this case, the one time monitor located in center of waveguide is enough to get the right result. However, in free space case, it shows very different result depends on the monitor location.
So, I set many time monitors like below picture and sum of the all electric field from monitors. And then take fourier transform. Is it right method?
To sum it up,
- The result of the bistability example is different from the result of web site.
- How can I calculate the reflected or transmitted electric field in free space case?