1维严格耦合波分析matlab程序

1维严格耦合波分析matlab程序 适用于一维光栅 可计算多层结构

% Author: Pavel Kwiecien pavel.kwiecien@seznam.cz
% Czech Technical University in Prage Optical Physics Group Czech Republic

if new_figure==1
    figure;
end

if slice==1 || slice2==1
    display（‘slice is on‘）;
    if plot_field_region_1==1
        display（‘plot_field_region_1=2‘）;
    end
    
    if plot_field_region_3==1
        display（‘plot_field_region_1=2‘）;
    end
    
    
    plot_field_region_1=2;
    plot_field_region_2=1;
    plot_field_region_3=2;
end

set（gcf ‘Renderer‘ ‘zbuffer‘）;
axes（‘FontSize‘17）
if plot_field_region_1==1 && plot_field_region_2==1 && plot_field_region_3==1
    total_axis_z=[z_11 z_22 z_33];
    
    switch planar_field_component
        case 1
            if polarization==1
                total_field=view_field（[field_E_1 field_E_2 field_E_3]field_component_operation）;
            elseif polarization==2
                total_field=view_field（[field_H_1 field_H_2 field_H_3]field_component_operation）;
            end
        case 2
            if polarization==1
                total_field=view_field（[field_H_1 field_H_2 field_H_3]field_component_operation）;
            elseif polarization==2
                total_field=view_field（[field_E_1 field_E_2 field_E_3]field_component_operation）;
            end
        case 3
            if polarization==1
                total_field=view_field（[field_E_1 field_E_2 field_E_3]field_component_operation）;
            elseif polarization==2
                total_field=view_field（[field_H_1 field_H_2 field_H_3]field_component_operation）;
            end
    end
    change_scale_z;    
    surf（total_axis_z.*1E6x.*1E6total_field）;
    axis（[-size_region_1*thickness_total （1+size_region_3）*thickness_total x_min*1e6 x_max*1e6]）
elseif plot_field_region_1==2 && plot_field_region_2==1 && plot_field_region_3==1
    total_axis_z=[z_22 z_33];
    switch planar_field_component
        case 1
            if polarization==1
                total_field=view_field（[field_E_2 field_E_3]field_component_operation）;
            elseif polarization==2
                total_field=view_field（[field_H_2 field_H_3]field_component_operation）;
            end
        case 2
            if polarization==1
                total_field=view_field（[field_H_2 field_H_3]field_component_operation）;
            elseif polarization==2
                total_field=view_field（[field_E_2 field_E_3]field_component_operation）;
            end
        case 3
            if polarization==1
                total_field=view_field（[field_E_2 field_E_3]field_component_operation）;
            elseif polarization==2
                total_field=view_field（[field_H_2 field_H_3]field_component_operation）;
            end
    end
    change_scale_z;
    surf（total_axis_z.*1E6x.*1E6total_field）;
    axis（[0 （1+size_region_3）*thickness_total x_min*1e6 x_max*1e6]）
elseif plot_field_region_1==2 && plot_field_region_2==2 && plot_field_region_3==1
    set_view_grating_boundaries=2;
    total_axis_z=z_33;
 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-06-04 13:14  rcwa-1d\
     文件         199  2011-07-22 09:02  rcwa-1d\zzreplace.sh
     文件        6802  2014-03-05 13:45  rcwa-1d\field_graphs_conical.m
     文件       12909  2014-06-04 11:25  rcwa-1d\rcwa_procedure.m
     文件        1373  2011-08-30 10:30  rcwa-1d\load_saved_field.m
     文件         911  2014-06-04 12:47  rcwa-1d\setup_dispersion.m
     文件       15405  2009-10-29 08:14  rcwa-1d\License
     目录           0  2014-06-04 12:38  rcwa-1d\measurements\
     目录           0  2014-06-04 13:14  rcwa-1d\scripts\
     文件        1918  2011-07-22 11:06  rcwa-1d\scripts\volume_grating.m
     目录           0  2012-03-30 17:08  rcwa-1d\scripts\real_profile_grating\
     文件       12006  2011-07-22 11:07  rcwa-1d\scripts\real_profile_grating\create_real_profile.m
     文件       24639  2011-02-10 13:56  rcwa-1d\scripts\real_profile_grating\NM20f50_um.txt
     文件         217  2007-06-01 12:47  rcwa-1d\scripts\real_profile_grating\profile_2.txt
     文件         831  2011-02-10 14:39  rcwa-1d\scripts\real_profile_grating\plot_grating_file.m
     文件         216  2007-06-01 11:04  rcwa-1d\scripts\real_profile_grating\profile_3.txt
     文件        3535  2012-03-27 13:22  rcwa-1d\scripts\real_profile_grating\Nm20f50_output.m
     文件         231  2007-06-01 10:52  rcwa-1d\scripts\real_profile_grating\profile_1.txt
     文件        6186  2014-06-02 07:03  rcwa-1d\scripts\general_volume_grating.m
     文件       23464  2014-06-04 09:35  rcwa-1d\help_control_file.m
     目录           0  2014-06-04 12:40  rcwa-1d\lib\
     文件       16494  2014-03-05 13:55  rcwa-1d\lib\field2_5_local_ab.m
     文件        6022  2013-02-27 13:21  rcwa-1d\lib\S_matrix_W_S_grating_conical.m
     文件        2116  2011-07-22 09:05  rcwa-1d\lib\field2_correct.m
     文件        1350  2014-06-03 13:42  rcwa-1d\lib\grating_setup.m
     文件        3791  2013-01-01 09:52  rcwa-1d\lib\framework.m
     文件        2028  2012-12-30 15:54  rcwa-1d\lib\initialization2_2.m
     文件        4868  2011-06-01 07:57  rcwa-1d\lib\S_Li_final.m
     文件        1310  2014-03-06 10:16  rcwa-1d\lib\integrated_absorption.m
     文件        2132  2014-02-21 07:41  rcwa-1d\lib\field2_correct_local_ab.m
     文件        1251  2011-08-23 11:41  rcwa-1d\lib\S_matrix_initialization.m
............此处省略83个文件信息