光纤光栅模拟

可实现不同结构，不同入射光条件下的光纤光栅光场分布求解。

clear;

for ns = [1.44];
    for theta = [5]

        clear -global;

        % function main
        yc_smf_28_1550;

        fiber.n = [nco; ncl; ns];
        fiber.r = [rco; rcl; rs];
        dr = 0.08;
        [r] = yc_gen_unif_mesh_fiber（fiberdr）;

        opts.pml.m = 2;
        opts.pml.R = 1e-12;
        opts.pml.n = round（10/dr）;
        opts.pml.kappa = 0;

        opts.sortmode = ‘lr‘;

        Dn.lam = 1.545/2/nco/cos（theta*pi/180）^2;
        Dn.chi = @（z） 5e-4;%*yc_gaussian（z05e3）;
        Dn.sigma =@（z） 2*Dn.chi（z）;
        Dn.dphi = @（z） zeros（size（z））;
        Dn.phi = @（z） zeros（size（z））;
        Dn.theta = theta *pi/180;
        Dn.L = 8e3;
        RI.n = @（r） nco;
        RI.P = @（r） 1;
        RI.rmin = 0;
        RI.rmax = rco;

        % z = linspace（-Dn.L/2Dn.L/21024）;
        % plot（zDn.chi（z））;

        grating.Dn = Dn;
        grating.RI = RI;
        grating.polarization = ‘ps‘;

        mode_solver.opts = opts;
        mode_solver.m1 = 1;
        mode_solver.nmodes1 = 5;
        mode_solver.sel_nmodes1 = 1;
        mode_solver.dir = ‘S0_FBG‘;
        mode_solver.m = 1:2:15;
        mode_solver.nmodes = 10;
        mode_solver.sel_nmodes = 10;
        mode_solver.r = r;
        mode_solver.savemodes = false;
        mode_solver.loadmodes = false;
        mode_solver.sol_modes = false;

        lambda = 1.536:0.00005:1.555;

        [RRTTRRRres] = yc_fbg（lambdafibermode_solvergrating）;

        name = [‘S1_‘num2str（theta）‘_‘...
            num2str（max（mode_solver.m））‘_‘...
            num2str（mode_solver.sel_nmodes）‘_‘...
            num2str（ns）‘.mat‘];
        save（name）;

        idx = 1:1:length（lambda）;
        figure
        plot（lambda（idx）10*log10（TT（idx:）））;
        hold on;
        % plot（lambda10*log10（RRR））;
        xlabel（num2str（theta））;
        title（num2str（ns））;
        % axis（[min（lambda） max（lambda）+eps -35 0]）;

        axis（[1.52 1.5575 -20 0]）;

    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1829  2008-12-14 23:20  yc_modesolver\examples\main_fiber_field_contour.m
     文件        1873  2014-03-09 07:18  yc_modesolver\examples\main_step_index_fiber.m
     文件         950  2014-03-09 07:18  yc_modesolver\examples\main_bragg_fiber_1.m
     文件        1137  2014-03-09 07:18  yc_modesolver\examples\main_bragg_fiber_2.m
     文件        1469  2014-03-09 07:18  yc_modesolver\examples\main_bragg_fiber.m
     文件        3038  2014-03-09 07:18  yc_modesolver\src\yc_circularc.m
     文件         395  2014-03-09 07:18  yc_modesolver\src\yc_gamma_fiber.m
     文件        1634  2014-03-09 07:18  yc_modesolver\src\yc_contour_mode_fiber.m
     文件        1011  2014-03-09 07:18  yc_modesolver\src\yc_sol_eig_fiber.m
     文件         618  2014-03-09 07:18  yc_modesolver\src\yc_cal_bandgap_fiber.m
     文件        1227  2014-03-09 07:18  yc_modesolver\src\yc_poynting_fiber.m
     文件        1466  2014-03-09 07:18  yc_modesolver\src\yc_gen_epsr_fiber.m
     文件        6762  2014-03-09 07:18  yc_modesolver\src\yc_fields_fiber.m
     文件        2447  2014-03-09 07:18  yc_modesolver\src\yc_fibermodes.m
     文件         358  2014-03-09 07:18  yc_modesolver\src\yc_gen_bragg_fiber.m
     文件         670  2014-03-09 07:18  yc_modesolver\src\yc_dispersion_fiber.m
     文件        5033  2014-03-09 07:18  yc_modesolver\src\yc_fiber_modes.m
     文件        1617  2014-03-09 07:18  yc_modesolver\src\yc_sort_modes_fiber.m
     文件         733  2014-03-09 07:18  yc_modesolver\src\yc_gen_unif_mesh_fiber.m
     文件        6339  2014-03-09 07:18  yc_modesolver\src\yc_build_mtx_fiber.m
     文件        1960  2014-03-09 12:46  yc_cmt\examples\main_fbg.m
     文件        1854  2014-03-09 12:44  yc_cmt\examples\main_lpf.m
     文件          52  2009-09-02 08:24  yc_cmt\examples\yc_smf_28_1550.m
     文件        1993  2008-10-16 17:18  yc_cmt\src\yc_patterson_test.m
     文件        3243  2014-03-09 12:37  yc_cmt\src\yc_cos_sin.m
     文件         645  2014-03-09 12:37  yc_cmt\src\yc_exp_phi.m
     文件        4000  2014-03-09 12:38  yc_cmt\src\yc_fbg.m
     文件        4761  2014-03-09 12:38  yc_cmt\src\yc_fbg_cmt.m
     文件         569  2014-03-09 12:38  yc_cmt\src\yc_fbg_r_t_cmt.m
     文件        2978  2008-10-16 15:58  yc_cmt\src\yc_patterson.m
     文件        4421  2014-03-09 12:38  yc_cmt\src\yc_kappa_chi.m
............此处省略9个文件信息