基于matlab实现光电导混频器中，交叉电极的光电导宽度对THz功率的影响.rar

clear Le=0.2e-6; Lg=1.8e-6; S1=400e-12; E=12.8; E0=8.85e-12; V=25; P=25e-3; e0=1.6e-19; tau=0.27e-12; mu=200e-4; hv=2.2784e-019; R=72; w=1e11:0.01e12:1e13; k=(1-(tan(pi*Le/(4*(Le+Lg))))^4)^0.5; C1=pi*(1+E)*E0*S1 C2=2*(Le+Lg)*log(2*(1+k^0.5)/(1-k^0.5)); C=C1/C2; S2=400e-12; G0=(P*e0*tau*mu)/(hv*S2); P1=(V*G0)^2*R; P2=2*(1+(w.*tau).^2).*(1+(w.*R*C).^2); y=P1./P2; y1=1e9*y; x=w./(2*pi*1e12); Le=0.5e-6; Lg=1.8e-6; S1=400e-12; E=12.8; E0=8.85e-12; V=25; P=25e-3; e0=1.6e-19; tau=0.27e-12; mu=200e-4; hv=2.2784e-019; R=72; w=1e11:0.01e12:1e13; k=(1-(tan(pi*Le/(4*(Le+Lg))))^4)^0.5; C1=pi*(1+E)*E0*S1 C2=2*(Le+Lg)*log(2*(1+k^0.5)/(1-k^0.5)); C=C1/C2; S2=400e-12; G0=(P*e0*tau*mu)/(hv*S2); P1=(V*G0)^2*R; P2=2*(1+(w.*tau).^2).*(1+(w.*R*C).^2); y=P1./P2; y2=1e9*y; x=w./(2*pi*1e12); Le=1.0e-6; Lg=1.8e-6; S1=400e-12; E=12.8; E0=8.85e-12; V=25; P=25e-3; e0=1.6e-19; tau=0.27e-12; mu=200e-4; hv=2.2784e-019; R=72; w=1e11:0.01e12:1e13; k=(1-(tan(pi*Le/(4*(Le+Lg))))^4)^0.5; C1=pi*(1+E)*E0*S1 C2=2*(Le+Lg)*log(2*(1+k^0.5)/(1-k^0.5)); C=C1/C2; S2=400e-12; G0=(P*e0*tau*mu)/(hv*S2); P1=(V*G0)^2*R; P2=2*(1+(w.*tau).^2).*(1+(w.*R*C).^2); y=P1./P2; y3=1e9*y; x=w./(2*pi*1e12); loglog(x,y1,'r',x,y2,'g',x,y3,'k') set(gca,'linewidth',1); axis([0.06 6 0.01 5]) legend('\fontsize{12}\rm le=0.2','\fontsize{12}\rm le=0.5','\fontsize{12}\rm le=1.0') xlabel('\fontsize{12}\rm frequentcy (THz) ') ylabel('\fontsize{12}\rm output power(nw)')