本程序为初学者使用,只考虑MT方向
下面的程序为matlab代码
只考虑MT方向
%This is a simple demo for Photonic Crystals simulation
%This is a simple demo for Photonic Crystals simulation
%This demo is for TE wave only, so only h wave is considered.
%for TM direction only,10 points is considered.
%---------------------------------------M
%| / |
%| / |
%| / |
%| --------------------|X
%| T |
%| |
%| |
%---------------------------------------
%equation :sum_{G',k}(K+G)(K+G')f(G-G')hz(k+G')=(omega/c)^2*hz(k+G)
%G' can considerd as the index of column, and G as index of rows
%[(K+G1)(K+G1)f(G1-G1) (K+G1)(K+G2)f(G1-G2) ][hz(G1)]=(omega/c)^2[hz(G1)]
%[(