per_diff.m

利用电磁场的源激发方法来计算光子晶体波导例如光子晶体光纤

function out = per_diff(in)
% periodic diff
out = [in(2) in(3:end) in(2)] - [in(end-1) in(1:end-2) in(end-1)];