RCWA在matlab里面的一维代码

电磁学在频率中的传播分析方法，RCWA是严格耦合波分析，可以分析各个衍射级的情况

% 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_1）*thickness_total 0 Lambda*number_of_plotted_period]）
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 0 Lambda*number_of_plotted_period]）
elseif plot_field_region_1==2 && plot_field_region_2==2 && plot_field_region_3==1
    set_view_grating_boundaries

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    .CA....     11251  2013-12-31 14:19  rcwa-1d\field_graphs.m

    .CA....      6914  2013-12-31 14:19  rcwa-1d\field_graphs_conical.m

    .CA....       187  2013-12-31 14:19  rcwa-1d\gratings\add_hom_layer.m

    .CA....      1700  2013-12-31 14:19  rcwa-1d\gratings\add_pol_layer.m

    .CA....      2882  2013-12-31 14:19  rcwa-1d\gratings\add_sin_layer.m

    .CA....      2992  2013-12-31 14:19  rcwa-1d\gratings\add_tri_layer.m

    .CA....      1700  2013-12-31 14:19  rcwa-1d\gratings\blazed_grating_1.m

    .CA....      1637  2013-12-31 14:19  rcwa-1d\gratings\blazed_grating_2.m

    .CA....      1233  2013-12-31 14:19  rcwa-1d\gratings\create_grating.m

    .CA....       149  2013-12-31 14:19  rcwa-1d\gratings\cSP_BK7_50nm_Au_water.m

    .CA....       133  2013-12-31 14:19  rcwa-1d\gratings\fill_NaN.m

    .CA....      1786  2013-12-31 14:19  rcwa-1d\gratings\gaussian_grating.m

    .CA....      2143  2013-12-31 14:19  rcwa-1d\gratings\general_triangular_grating.m

    .CA....      8049  2013-12-31 14:19  rcwa-1d\gratings\grating_example_37.m

    .CA....        59  2013-12-31 14:19  rcwa-1d\gratings\myfun.m

    .CA....     12351  2013-12-31 14:19  rcwa-1d\gratings\Nm20f50.m

    .CA....       217  2013-12-31 14:19  rcwa-1d\gratings\notch_filter.m

    .CA....      1619  2013-12-31 14:19  rcwa-1d\gratings\parabolic_grating.m

    .CA....      1293  2013-12-31 14:19  rcwa-1d\gratings\plot_grating_file.m

    .CA....     12105  2013-12-31 14:19  rcwa-1d\gratings\polyxpoly.m

    .CA....        86  2013-12-31 14:19  rcwa-1d\gratings\sample_grating.m

    .CA....      1585  2013-12-31 14:19  rcwa-1d\gratings\sinusoidal_grating_y.m

    .CA....      1086  2013-12-31 14:19  rcwa-1d\gratings\sinus_x.m

    .CA....       990  2013-12-31 14:19  rcwa-1d\gratings\test.m

    .CA....       100  2013-12-31 14:19  rcwa-1d\gratings\test_grating.m

    .CA....        62  2013-12-31 14:19  rcwa-1d\gratings\test_schuster.m

    .CA....      4450  2013-12-31 14:19  rcwa-1d\gratings\tetragonal_grating.m

    .CA....      2299  2013-12-31 14:19  rcwa-1d\gratings\trapezoidal_grating.m

    .CA....      1712  2013-12-31 14:19  rcwa-1d\gratings\triangular_grating.m

    .CA....     23089  2014-01-06 00:39  rcwa-1d\help_control_file.asv

............此处省略259个文件信息