太赫兹和中红外频率下-石墨烯表面电导率的公式全套源码

%所有变量均取标准单位 e=1.602*10^(-19);%元电荷，C Ef=0.2*1.602*10^(-19);%费米能级，J h0=1.0546*10^(-34);%约化普朗克常数，Js kB=1.38*10^(-23);%玻尔兹曼常数，J/K T=300;%温度，K w=2*pi*36*10^12:2*pi*0.1*10^12:2*pi*60*10^12;%角频率，rad/s vf=10^6;%费米速度，m/s u=1;%载流子迁移率，m2/Vs t1=e*vf^2*(u*Ef)^(-1);%弛豫时间倒数，s-1 t=5*10^12;%同t1,s-1 sigma1=(2*1i*e^2*kB*T*(pi*h0^2)^(-1)*log(2*cosh(Ef/(2*kB*T))))./((w+1i*t));%石墨烯表面电导率，带内 sigma2=(e^2/(4*h0))*(0.5+(1/pi)*atan((w*h0-Ef)/(2*kB*T))-(1i/(2*pi))*log(((w*h0+Ef).^2)/((w*h0-Ef).^2+(2*kB*T)^2)));%石墨烯表面电导率，带间 sigma=sigma1+sigma2; tg=0.5*10^(-9);%石墨烯厚度，nm eps0=8.854*10^(-12);%真空介电常数，F/m eps=1+(1i*sigma)./(w.*eps0*tg); a=real(eps)'; b=imag(eps)'; c=(w/(2*pi*3*10^10))';