Info

This question is closed. Reopen it to edit or answer.

How could i calculate out Fi_v(0) and the raw values in numerical form from the code below..??

1 view (last 30 days)
%Inputs Wm = 120; Lembda0 = 1.55; Sigma = 0; nc = 1.456; nf = 1.46;
%Effective width calculation We = Wm+Lembda0/pi*((nc^2)/(nf^2))^Sigma*1/(sqrt(nf^2-nc^2));
%Beat length calculation Lpi = (4*nf*(We^2)/(3*Lembda0)); fprintf('Beat Length, Lpi = %.4f \n \n',Lpi)
%Propagation vector calculation k0=2*pi/Lembda0; fprintf('Free space wave number, k0 = %.4f \n \n',k0)
%Fundamental propagation constant calculation Beta0 = k0*nf- ((pi*Lembda0)/(4*nf*We^2)); fprintf('Beta0 = %.4f \n \n',Beta0)
%Propagation Constant Calculation & plotting X=[]; Vm = 13; for v=1:(Vm-1) Beta = Beta0 - ((v*(v+2)*pi)/(3*Lpi)); fprintf('Beta%d = %.4f \n \n',v , Beta) X(v) = Beta; end y=1:1:12; plot(y,X, '--rs','LineWidth',2,... 'MarkerEdgeColor','k',... 'MarkerFaceColor','g',... 'MarkerSize',5) xlabel('Mode Number'); ylabel('Beta_v');
%Effective index Calculation Q=[]; for v=0:(Vm-1) if v == 0 Neff = Beta0/k0; else Neff = X(v)/k0; end Q(v+1) = Neff; fprintf('Neff%d = %.4f \n \n',v , Neff) end q=0:1:(Vm-1); plot(q,Q, '--rs','LineWidth',2,... 'MarkerEdgeColor','k',... 'MarkerFaceColor','m',... 'MarkerSize',5) xlabel('Mode Number'); ylabel('N_effv');
%eigen mode calculation Fi_v(y) = sin(((v+1)*pi/We)*y); plot(Fi_v(y)) xlabel('Mode Number'); ylabel('Eigen mode');
I am trying to find out the eigen values. But in my graph it starts from 1, but i need to calculate the value for 0 also. Also it will be helpful if i get the values in numerical forms for each of the mode numbers.

Answers (0)

Community Treasure Hunt

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

Start Hunting!