green_fun_vector.m

来自「利用电磁场的源激发方法来计算光子晶体波导例如光子晶体光纤」· M 代码 · 共 27 行

function vGreen = green_fun_vector(strClass)
if lower(strClass(1)) == 'f'
        k = str2num(strClass(2));
        N = 6;
        vGreen = exp(j*2*pi*k/N*(0:N-1));
        return
end
switch lower(strClass)
        case 'p1'
                vGreen = ones(1,12);
        case 'p2'
                vGreen = (-1).^(0:11);
        case 'p3'
                vGreen = [1 -1 -1 1];
        case 'p4'
                vGreen = [1 1 -1 -1];
        case 'p5'
                vGreen = [1 1 1 1];
        case 'p6'
                vGreen = [1 -1 1 -1];
        case 'p7'
                vGreen = [1 -1 -1 1 1 -1 -1 1 1 -1 -1 1];
        case 'p8'
                vGreen = [1 1 -1 -1 1 1 -1 -1 1 1 -1 -1 ];
        otherwise
                error('green_fun_vector: unknown symmetry class')
end