squarepixel_ebxn_1d_fc.m

来自「通过计算晶体结构来确定其光谱分布」· M 代码 · 共 18 行

function ret=SquarePixel_ebxn_1D_FC(ebxn,pw,m)
%return Fourier Coefficients of 1D ebxn with Square pixels
%ebxn=the first quadrant in permittivity 
%pw=pixel width
%m=orders index of Fourier coefficient along z axis
t=size(ebxn);
M=max(t);
cw=2*M*pw;
t0=2*pw/cw;
mt=pi*m*pw/cw;
if mt==0
    sm=1;
else
    sm=sin(mt)/mt;
end
ret=t0*sum(ebxn.*cos(mt*(2*(1:M)-1)))*sm;