- This topic has 20 replies, 4 voices, and was last updated 9 years, 8 months ago by
mohammed hosni.
-
AuthorPosts
-
-
April 13, 2015 at 8:19 am #19631
mohammed hosni
Participantplease i need to design WDM communication system by using PMMA polymer fiber with the following parameter attenuation 150 db/km , length 20 m,temperature 70 Celsius,core size 980 micro meter , cladding size 1mm.
as there isn’t PMMA polymer fiber component in optisystem so i need to model PMMA polymer fiber with the mentioned parameter using matlab component. if there is another way to use optical fiber component in optisystem please tell me it if not please tell me how i can model first glass optical fiber using matlab component to compare the results with the optical fiber component which already exist in optisystem software .
then i will model PMMA polymer fiber.
i need your help for my master degree study.
thanks. -
April 13, 2015 at 10:42 am #19643
Damian Marek
ParticipantIf you have more information about the fiber you can do it directly in OptiSystem, like core refractive index, cladding index etc.
-
April 13, 2015 at 3:18 pm #19665Participant
what about the other mentioned parameters?
what about modeling optical fiber with matlab component? -
April 14, 2015 at 1:48 am #19680
Abdallah Ismail
ParticipantWe are modeling fiber using split step Fourier method that solve Schrodinger equation.
In case that your fiber is a special case of SMF or MMF, you can use Measured components of theses fibers with your parameters that Damian mentioned above to model your fiber. -
April 14, 2015 at 12:13 pm #19709
Dr Rk Sethi
ParticipantAbdallah Hi,
Can you elaborate, how are you modelling and simulating fiber.
thanks and regards,
Sethi Rk -
April 14, 2015 at 2:55 pm #19713Participant
Refractive index=1.4906
Chromatic dispersion=dn/dλ = -0.043329 µm-1
Group velocity dispersion
GVD = 81.808 fs2/mm
D = -446.31 ps/(nm km)
Dispersion formula
n2−1=0.99654λ2λ2−0.00787+0.18964λ2λ2−0.02191+0.00411λ2λ2−3.85727 -
April 20, 2015 at 4:57 am #19938Participant
Is this dispersion value is real?
-
April 21, 2015 at 3:14 pm #20017Participant
yes
-
April 24, 2015 at 1:18 pm #20104Participant
please find the attached files matlam m.file for modeling the attenuation and the dispersion of fibe that is just example but there is an error about the dispersion please help me to model the dispersion by matlab
-
April 24, 2015 at 1:19 pm #20105Participant
that is optisystem file
-
April 24, 2015 at 2:17 pm #20108Participant
I hope the attached files, helps you Mohamed.
-
April 24, 2015 at 2:18 pm #20109Participant
There is a problem with the attachments
-
April 24, 2015 at 2:19 pm #20111Participant
also this one
-
April 24, 2015 at 2:21 pm #20112Participant
that is it
-
April 24, 2015 at 2:22 pm #20113Participant
hope optiwave solve it
-
April 24, 2015 at 2:50 pm #20114Participant
thanks Abdallah Ismail
but how can i use this file with optisystem matlab component -
April 24, 2015 at 3:05 pm #20115Participant
it’ll not work without the function
-
April 24, 2015 at 3:06 pm #20116Participant
hope it attached for the last time
-
April 24, 2015 at 3:07 pm #20117Participant
copy this code to a function file with name sym_ssf
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55
function [to,fo] = sym_ssf(M,h,gamma,Dh,uf0)
% Symmetrized Split-Step Fourier Algorithm
%
% ==Inputs==
% M = Simulation step number ( M*h = simulation distance )
% h = Simulation step
% gamma = Nonlinearity coefficient
% Dh = Dispersion operator in frequency domain
% uf0 = Input field in the frequency domain
%
% ==Outputs==
% to = Output field in the time domain
% fo = Output field in the frequency domain
%
% written by Jong-Hyung Lee
for k = 1:M
%=============================================================
% Propagation in the first half dispersion region, z to z+h/2
%=============================================================
Hf = Dh.*uf0;
%==========================================================
% Initial estimate of the nonlinear phase shift at z+(h/2)
%==========================================================
% Initial estimate value
ht = ifft(Hf); % time signal after h/2 dispersion region
pq = ht.*conj(ht); % intensity in time
u2e = ht.*exp(h*1i*gamma*pq); %Time signal
%=============================================================
% Propagation in the second Dispersion Region, z+(h/2) to z+h
%=============================================================
u2ef = fft(u2e);
u3ef = u2ef.*Dh;
u3e = ifft(u3ef);
u3ei = u3e.*conj(u3e);
%========================================================
% Iteration for the nonlinear phase shift(two iterations)
%========================================================
u2 = ht.*exp((h/2)*1i*gamma*(pq+u3ei));
u2f = fft(u2) ;
u3f = u2f.* Dh;
u4 = ifft(u3f);
u4i = u4.*conj(u4);
u5 = ht.*exp((h/2)*1i*gamma*(pq+u4i));
u5f = fft(u5);
uf0 = u5f.*Dh;
u6 = ifft(uf0); u6i = u6.*conj(u6);
%=============================================================
% Maximum allowable tolerance after the two iterations
etol = 1e-5;
if abs(max(abs(u6i))-max(abs(u4i)))/max(abs(u6i)) > etol
disp(‘Peak value is not converging Reduce Step Size’),break
end
%=============================================================
end
to = u6; fo = uf0;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
April 24, 2015 at 3:15 pm #20118Participant
I think that code will not help you, where your fiber can not be modeled with SMF Schrodinger’s equation.
Schrodinger’s equation is not used for modeling multimode fiber as well as PMMA polymer fiber.
There are conditions for Schrodinger’s equation to be valid.
Try to look for modeling PMMA polymer fiber, and i think you can use another tools of optiwave for this purpose.
Regards, -
April 25, 2015 at 1:03 am #20121Participant
if i want this code for smf fiber how can i use this code with optisytem matlab component
-
-
AuthorPosts
- You must be logged in to reply to this topic.