1d_photonic.m

一个用来计算一维光子晶体全息干涉分布的程序

>> theta=10*pi/180;lambda=441.6/1000;
>> [x,y]=meshgrid(0:0.05:15);
d=4;
a=size(x);b=size(y); a=a(1,2);b=b(1,2);a;b;
I2=zeros(a,b);
for m=1: a
for n=1: b
I2(m,n)=1.01+cos(2*pi*2*sin(theta)*y(m,n)./lambda);;
end
end
hold on;surf(x,y,I+I2)
view(0,90)
shading interp