Drude模型计算介电常数matlab（可直接运行）

clear
clc
lamda=1:0.001:5;
c=3*10^8;
q=1.60218910^(-19);%单个电子的电荷量
ee=5.5;%背景折射率
me1=9.1095310^(-31);%单个电子的质量
me=0.21*me1;
n=3.710^26;%载流子浓度
u=359*10^(4);%载流子迁移率
e0=8.85*10^(-12);%真空的介电常数
gama=q/(u*me);%阻尼
wp=((n*q^2)/(me*e0));%等离子体频率
w=(2*pi*c)./(lamda*10^(-6));

jiedian=ee-((wp)^2)./(w.*(gama*1i+w));%矩阵需要用.*

plot(lamda,real(jiedian))
figure
plot(lamda,imag(jiedian))