alpha.m

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

function ret=alpha(rebxn,pw,cw,m)
%rebxn=reverse permittivity in half cell
%pw=pixel width
%cw=cell width
%m=order index of Fourier coefficient
t=size(rebxn);
M=max(t);
t0=pi*m*pw/cw;
if t0==0
    t1=1;
else
    t1=sin(t0)/t0;
end
t2=1:M;
t3=rebxn.*cos(t0*2*(t2-1))*t1;
ret=2*pw/cw*sum(t3);