查看: 116|回复: 3|关注: 0

[已答复] 如何用传输矩阵法求各项异性光子晶体的反射透射率

[复制链接]

新手

5 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
发表于 2020-5-23 21:46:30 | 显示全部楼层 |阅读模式
本帖最后由 TouAkira 于 2020-5-23 11:36 编辑
  1. c=3e8;ua=1;ub=1;u0=(4e-7)*pi;e0=8.85e-12;n_base=1.52;f=0.5;
  2. da=132e-9;db=25e-9;nb=2.12;eb=nb^2;f0=0.8213;
  3. wp=1.2983e16;r0=0.3420e15;w1=0.6330*wp;f1=0.4833;r1=0.9044e15;
  4. nc=2.12;ec=nc*nc;
  5. x=zeros(1,250);
  6. R0=zeros(1,250);
  7. T0=zeros(1,250); %要不要都不影响结果,要才严谨,而且运行更快
  8. for p=1:6
  9.     theta(p)=0+(p-1).*15;
  10.    the_0=theta(p)/180*pi;
  11. for n=1:250
  12.     lamda=(n+300)*(1e-9);
  13.      x(n)=lamda./(1e-9);
  14. w=2*pi*c./lamda;
  15. k0=w/c;
  16. ea0=1.0-f0*wp*wp./(w.^2+1i.*r0.*w)-f1.*wp.*wp./(w.^2-w1.*w1+1i.*r1.*w);
  17. eax=f.*ec+(1-f).*ea0;
  18. eaz=1./(f./ec+(1-f)./ea0);
  19. % ea=[eax 0 0;0 eax 0;0 0 eaz];
  20. na=sqrt(eax);%%%%%%%%%%%%%%%%%%%%%%%%%%怎样表示a的折射率,我这样是不正确的。
  21. the_a=asin(sin(the_0)./na);
  22. the_b=asin(sin(the_0)./nb);
  23. the_base=asin(sin(the_0)/n_base);
  24. k1z=w.*na.*cos(the_a)./c;
  25. k2z=w.*nb.*cos(the_b)./c;
  26. ita_a=na./cos(the_a);
  27. ita_b=nb./cos(the_b);
  28. % ita_a=w.*e0.*eax./k1z;
  29. % ita_b=w*e0*eb/k2z;
  30. M1=[cos(k1z*da),-1i.*sin(k1z*da)./ita_a;-1i.*ita_a.*sin(k1z*da),cos(k1z*da)];
  31. M2=[cos(k2z*db),-1i.*sin(k2z*db)./ita_b;-1i.*ita_b.*sin(k2z*db),cos(k2z*db)];
  32. Mt=(M1*M2)^3;
  33. A=Mt(1,1);
  34. B=Mt(1,2);
  35. C=Mt(2,1);
  36. D=Mt(2,2);
  37. ita_0=1/cos(the_0);
  38. ita_base=n_base/cos(the_base);
  39. r=(A.*ita_0+B.*ita_0.*ita_base-C-D.*ita_base)/(A.*ita_0+B.*ita_0.*ita_base+C+D.*ita_base);
  40. t=2.*ita_0/(A.*ita_0+B.*ita_0.*ita_base+C+D.*ita_base);
  41. R=abs(r).^2;
  42. T=abs(t).^2.*ita_base./ita_0;
  43. R0(n)=R;
  44. T0(n)=T;
  45. end
  46. plot(x,T0);
  47. hold on
  48. % plot(x,T0);
  49. end
复制代码





求大神指教该如何编写各项异性光子晶体的传输矩阵。

回复主题 已获打赏: 0 积分

举报

论坛优秀回答者

入门

188 麦片

财富积分


50500


0

主题

448

帖子

38

最佳答案
  • 关注者: 11
发表于 2020-5-24 09:37:23 | 显示全部楼层
你这样写大家肯定没办法给你解答,除非是专业的,才能明白异性光子晶体的反射透射率,最好是直接把求解的公式贴出来
回复此楼 已获打赏: 0 积分

举报

新手

5 麦片

财富积分


050


1

主题

2

帖子

0

最佳答案
 楼主| 发表于 3 天前 | 显示全部楼层
本帖最后由 ysn 于 2020-5-31 23:10 编辑
wjb986555360 发表于 2020-5-24 09:37
你这样写大家肯定没办法给你解答,除非是专业的,才能明白异性光子晶体的反射透射率,最好是直接把求解的公 ...

我知道哪里错了,现在我有一个新的问题。我已经得到了下面图的程序。我如果想得到上面的彩色投影图应该怎么改c=3e8;u0=4e-7*pi;e0=8.85e-12;n_base=1.52;f=0.5;
da=132e-9;db=25e-9;f0=0.8213;
nc=2.12;
ec=nc*nc;
wp=1.2983e16;
r0=0.3420e15;
w1=0.6330*wp;
f1=0.4833;r1=0.9044e15;
x=zeros(1,250);
R0=zeros(1,250);
T0=zeros(1,250);
theta=zeros(1,6);
for p=1:6
    theta(p)=0+(p-1).*15;%入射角
   the_0=theta(p)/180*pi;
for n=1:250
    lamda=(n+300)*(1e-9);
     x(n)=lamda./(1e-9);
     nb=2.12 ;
     eb=nb.^2;
w=2*pi*c./lamda;
k0=w/c;
ea0=1.0-f0*wp*wp./(w.^2+1i.*r0.*w)-f1.*wp.*wp./(w.^2-w1.*w1+1i.*r1.*w);
eax=f.*ec+(1-f).*ea0;
eaz=1./(f./ec+(1-f)./ea0);
the_b=asin(sin(the_0)./nb);
k1z=k0*sqrt(eax*(1-sin(the_0)^2./eaz));
k2z=k0*sqrt(eb*(1-sin(the_0)^2./eb));
the_base=asin(sin(the_0)./n_base);
ita_a=w.*e0.*eax./k1z;
ita_b=w.*e0.*eb./k2z;
M1=[cos(k1z.*da),-1i.*sin(k1z.*da)./ita_a;-1i.*ita_a.*sin(k1z.*da),cos(k1z.*da)];
M2=[cos(k2z.*db),-1i.*sin(k2z.*db)./ita_b;-1i.*ita_b.*sin(k2z.*db),cos(k2z.*db)];
Mt=(M1*M2)^3;
A=Mt(1,1);
B=Mt(1,2);
C=Mt(2,1);
D=Mt(2,2);
ita_0=sqrt(e0/u0)./cos(the_0);
ita_base=sqrt(e0/u0)./cos(the_base);
r=(A.*ita_0+B.*ita_0.*ita_base-C-D.*ita_base)./(A.*ita_0+B.*ita_0.*ita_base+C+D.*ita_base);
t=2.*ita_0./(A.*ita_0+B.*ita_0.*ita_base+C+D.*ita_base);
R=abs(r)^2;%反射率
T=abs(t)^2.*ita_base./ita_0;%透射率
R0(n)=R;
T0(n)=T;
end
subplot(2,1,2);
plot(x,T0);
hold on
legend('0','15','30','45','60','75');
set(gca,'XLim',[300,550]);
set(gca,'XTick',300:50:550);
set(gca,'YLim',[0,0.3]);
set(gca,'YTick',0:0.1:0.3);
subplot(2,1,1);
plot(x,R0);
hold on
legend('0','15','30','45','60','75');
set(gca,'XLim',[300,550]);
set(gca,'XTick',300:50:550);
set(gca,'YLim',[0,1]);
set(gca,'YTick',0:0.5:1);
hold on
end
微信图片_20200531225941.jpg
微信图片_20200531225930.jpg
回复此楼 已获打赏: 0 积分

举报

论坛优秀回答者

入门

188 麦片

财富积分


50500


0

主题

448

帖子

38

最佳答案
  • 关注者: 11
发表于 昨天 22:24 | 显示全部楼层
我看你程序画的图跟文献的不一样,而且我用你的程序参数画的图也与文献不太一样
untitled.jpg
回复此楼 已获打赏: 0 积分

举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表