Ce_YIG_Ce：YIG晶体_一维磁光晶体透射和法拉第旋转_法拉第_透射反射_

clc clear e0=8.854187817e-12;%F/m hbar=1.05457266e-34;%J*s ee=1.6e-19; c=3e8; lambda_PC=1550e-9; %石墨烯 dg=0.34e-9; vf=1e6;%m/s Ef=0.34*ee;%J F=1e7*ee;%J 散射时间 B=1;%T D=ee^2*abs(Ef)/hbar^2; wc=ee*B*vf^2/Ef; aaaa=wc^2; %A Ce:YIG e_xx1=4.884; n_A_n2=sqrt(4.884); e_xy1=0.009*1i; n_A_n=sqrt(e_xx1+1i*e_xy1);%正 n_A_p=sqrt(e_xx1-1i*e_xy1);%负 d_A=lambda_PC/4/n_A_n2; d_A2=lambda_PC/4/n_A_n; %B SiO2 n_B=sqrt(2.19); d_B=lambda_PC/4/n_B; for j=1:5000 lambda(1)=1546e-9;%30微米到130微米 lambda(j+1)=lambda(j)+8e-9/5000; w(j)=2*pi*c/lambda(j); sigma_xx(j)=(D/pi)*(F-1i*w(j))/(wc^2-(w(j)+1i*F)^2); sigma_xy(j)=(-D/pi)*wc/(wc^2-(w(j)+1i*F)^2); %石墨烯介电常数 e_xx(j)=1+1i*sigma_xx(j)/(w(j)*e0*dg); a(j)=1i*sigma_xx(j)/(w(j)*e0*dg); b(j)=w(j)*e0*dg; d(j)=1i*sigma_xx(j); ddd(j)=(F-1i*w(j)); aaa(j)=(w(j)+1i*F)^2; aa(j)=wc^2-(w(j)+1i*F)^2; e_xy(j)=1i*sigma_xy(j)/(w(j)*e0*dg); %石墨烯矩阵 G NG_positive(j)=sqrt(e_xx(j)+1i*e_xy(j));%可能不需要i NG_negative(j)=sqrt(e_xx(j)-1i*e_xy(j)); beta_G_p(j)=2*pi/lambda(j)*NG_positive(j)*dg; beta_G_n(j)=2*pi/lambda(j)*NG_negative(j)*dg; SG = zeros(4,4); SG(1,1)=cos(beta_G_p(j)); SG(1,2)=1i/NG_positive(j)*sin(beta_G_p(j)); SG(1,3)=0; SG(1,4)=0; SG(2,1)=1i*NG_positive(j)*sin(beta_G_p(j)); SG(2,2)=cos(beta_G_p(j)); SG(2,3)=0; SG(2,4)=0; SG(3,1)=0; SG(3,2)=0; SG(3,3)=cos(beta_G_n(j)); SG(3,4)=1i/NG_negative(j)*sin(beta_G_n(j)); SG(4,1)=0; SG(4,2)=0; SG(4,3)=1i*NG_negative(j)*sin(beta_G_n(j)); SG(4,4)=cos(beta_G_n(j)); %A Ce:YIG beta_A_n(j)=2*pi/lambda(j)*n_A_n*d_A; beta_A_p(j)=2*pi/lambda(j)*n_A_p*d_A; S_A = zeros(4,4); S_A(1,1)=cos(beta_A_n(j)); S_A(1,2)=1i/n_A_n*sin(beta_A_n(j)); S_A(1,3)=0; S_A(1,4)=0; S_A(2,1)=1i*n_A_n*sin(beta_A_n(j)); S_A(2,2)=cos(beta_A_n(j)); S_A(2,3)=0; S_A(2,4)=0; S_A(3,1)=0; S_A(3,2)=0; S_A(3,3)=cos(beta_A_p(j)); S_A(3,4)=1i/n_A_p*sin(beta_A_p(j)); S_A(4,1)=0; S_A(4,2)=0; S_A(4,3)=1i*n_A_p*sin(beta_A_p(j)); S_A(4,4)=cos(beta_A_p(j)); %B SiO2 beta_B(j)=2*pi/lambda(j)*n_B*d_B; S_B = zeros(4,4); S_B(1,1)=cos(beta_B(j)); S_B(1,2)=1i/n_B*sin(beta_B(j)); S_B(1,3)=0; S_B(1,4)=0; S_B(2,1)=1i*n_B*sin(beta_B(j)); S_B(2,2)=cos(beta_B(j)); S_B(2,3)=0; S_B(2,4)=0; S_B(3,1)=0; S_B(3,2)=0; S_B(3,3)=cos(beta_B(j)); S_B(3,4)=1i/n_B*sin(beta_B(j)); S_B(4,1)=0; S_B(4,2)=0; S_B(4,3)=1i*n_B*sin(beta_B(j)); S_B(4,4)=cos(beta_B(j)); %空气 D2=[1 1 0 0;1 -1 0 0;0 0 1 1;0 0 1 -1]; D1=inv(D2); N=6; M=D1*(S_A*S_B)^N*S_A*S_A*(S_B*S_A)^N*D2; % M=D1*(S_A*S_B)^N*S_A*SG*S_A*(S_B*S_A)^N*D2;%石墨烯嵌入中间 R_p(j)=(abs(M(2,1)/M(1,1)))^2; R_n(j)=(abs(M(4,3)/M(3,3)))^2; R(j)=(R_p(j)+R_n(j))/2; tp(j)=1/M(1,1); tn(j)=1/M(3,3); T_p(j)=(abs(1/M(1,1)))^2; T_n(j)=(abs(1/M(3,3)))^2; T(j)=(T_p(j)+T_n(j))/2; A_p(j)=1-R_p(j)-T_p(j); A_n(j)=1-R_n(j)-T_n(j); A(j)=1-R(j)-T(j); %法拉第旋转角 fai_t_p(j)=atan(imag(1/M(1,1))/real(1/M(1,1))); fai_t_n(j)=atan(imag(1/M(3,3))/real(1/M(3,3))); % fai_t_p(j)=angle(tp(j)); % fai_t_n(j)=angle(tn(j)); xita_t(j)=-1/2*(fai_t_p(j)- fai_t_n(j)); %椭圆率 kexi_t(j)=(1/M(3,3))/(1/M(1,1)); x_t(j)=atan((kexi_t(j)-1)/(kexi_t(j)+1)); end hold on figure(1) plot(lambda(1:5000),T(1:5000),'r') plot(lambda(1:5000),T(1:5000),'k',lambda(1:5000),1/2*T_n(1:5000),'r',lambda(1:5000),1/2*T_p(1:5000),'b') figure(2) plot(lambda(1:5000),xita_t(1:5000),'k',lambda(1:5000),fai_t_p(1:5000),'g',lambda(1:5000),fai_t_n(1:5000),'r',lambda(1:5000),x_t(1:5000),'b') % figure(1) % plot(lambda(1:5000),e_xx(1:5000),'r') % figure(2) % plot(lambda(1:5000),e_xy(1:5000),'r')