- This topic has 20 replies, 4 voices, and was last updated 10 years 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 #19665
mohammed hosni
Participantwhat 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 #19713
mohammed hosni
ParticipantRefractive 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 #19938
Abdallah Ismail
ParticipantIs this dispersion value is real?
-
April 21, 2015 at 3:14 pm #20017
mohammed hosni
Participantyes
-
April 24, 2015 at 1:18 pm #20104
mohammed hosni
Participantplease 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 #20105
mohammed hosni
Participantthat is optisystem file
-
April 24, 2015 at 2:17 pm #20108
Abdallah Ismail
ParticipantI hope the attached files, helps you Mohamed.
-
April 24, 2015 at 2:18 pm #20109
Abdallah Ismail
ParticipantThere is a problem with the attachments
-
April 24, 2015 at 2:19 pm #20111
Abdallah Ismail
Participantalso this one
-
April 24, 2015 at 2:21 pm #20112
Abdallah Ismail
Participantthat is it
-
April 24, 2015 at 2:22 pm #20113
Abdallah Ismail
Participanthope optiwave solve it
-
April 24, 2015 at 2:50 pm #20114
mohammed hosni
Participantthanks Abdallah Ismail
but how can i use this file with optisystem matlab component -
April 24, 2015 at 3:05 pm #20115
Abdallah Ismail
Participantit’ll not work without the function
-
April 24, 2015 at 3:06 pm #20116
Abdallah Ismail
Participanthope it attached for the last time
-
April 24, 2015 at 3:07 pm #20117
Abdallah Ismail
Participantcopy 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 #20118
Abdallah Ismail
ParticipantI 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 #20121
mohammed hosni
Participantif 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.