求文档: 用matlab编程求一维光子晶体的折射率和透射率

na=2.35 %A介质为硫化锌

nb=1.38 %B介质为二氟化镁

a=0.0000000357%中心波长为

600nm

b=0.0000001087

c=2.35

d=1.38

x=400:1:800

e=-2*pi*na*a./(x*0.000000001)

f=-2*pi*nb*b./(x*0.000000001)

num=length(e)

y=zeros(1,num)

for j=1:num

M=[cos(e(j)),i.*sin(e(j))./ci.*c.*sin(e(j)),cos(e(j))]

N=[cos(f(j)),i.*sin(f(j))./di.*d.*sin(f(j)),cos(f(j))]

H=M*N

T=H^10

F=[11]

E=T*F

B=E(1)

C=E(2)

rs=(B-C)/(B+C)

y(j)=1-(abs(rs))^2

end

plot(x,y,'k')

xlabel('波长')

ylabel(' 透射率')

耐心找

论坛里一大堆

还有你可以去维基百科搜FDTD

下面的相关链接里面有很多人自己编写的FDTD相关的网站

那里面很多都是开源的