Info

この質問は閉じられています。 編集または回答するには再度開いてください。

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

1 回表示 (過去 30 日間)
Hossen_BD
Hossen_BD 2014 年 8 月 15 日
閉鎖済み: MATLAB Answer Bot 2021 年 8 月 20 日
%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.

回答 (0 件)

Community Treasure Hunt

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

Start Hunting!

Translated by