normalized propagation constant vs v anyone please can u help me to correct the program

11 views (last 30 days)
oussama isaac on 28 Nov 2021
Answered: Jack Goodhew on 29 Nov 2021
clc;
clear all; % Given Parameters
lambda = 1550 * 10^-9; % in meter
n_2 = 1.444; % refractive index of the cladding
a = 8 * 10^-6 ; % Core radius % given core diameter = 16 micro meter
% n_eff = 1.451190;
n_eff = [(1.451190) (1.453455) (1.454277) (1.456903) (1.458948)]; % Effective refractive index for LP_01 mode
n_1 = 1.474; % refractive index of the core
%f=[12.7 10.8 8.3 7.5 5.3 3.7 1.9]; %delta n
% Calculation
k= (2*pi)/lambda;
%V=[2:2:10];
%V = linspace(0,10,5);
V=linspace(0,10,5);
% b =((n_eff).^2 -(n_2)^2)./((V *lambda)./(2*pi*a)).^2;
b = (n_eff-n_2)/(n_1-n_2);
%b = f*1e-3./(n_1-n_2);
plot (V,b);
xlabel('frequence normalise V')
ylabel('constante de propagation b')
%xlim([0 12])
%ylim([0 1])
0 CommentsShowHide -1 older comments

Sign in to comment.

Answers (1)

Jack Goodhew on 29 Nov 2021
clc;
clear all;
% Given Parameters
lambda = 1550*10^-9; % in meters
n_2 = 1.444; % refractive index of the cladding
a = 8*10^-6; % Core radius (given core diameter = 16 micro meters)
n_eff = [1.451190 1.453455 1.454277 1.456903 1.458948]; % Effective refractive index for LP_01 mode
n_1 = 1.474; % refractive index of the core
%f=[12.7 10.8 8.3 7.5 5.3 3.7 1.9]; % delta n
% Calculation
k = (2*pi)/lambda;
%V=[2:2:10];
%V = linspace(0,10,5);
V=linspace(0,10,5);
% b =((n_eff).^2 -(n_2)^2)./((V *lambda)./(2*pi*a)).^2;
b = (n_eff-n_2)/(n_1-n_2);
normalized_b = b./max(b);
%b = f*1e-3./(n_1-n_2);
plot (V,normalized_b);
xlabel('frequence normalise V')
ylabel('constante de propagation b')
%xlim([0 12])
%ylim([0 1])
0 CommentsShowHide -1 older comments

Sign in to comment.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!