为什么phi_d11=phi_s11啊？求大佬帮助

圣光之血

clear
clc
L=300e-9;
W=1e-6;
kox=11;
T      =300;
s=10e-9;
d=1.2e-9;
tox=10e-9;

c=0.511;
Ncnt=round(W/s);
q= 1.6e-19;
KB    = 1.38e-23;
h  = 6.626e-34;
acc=0.142e-9;
Vp   =3;
Cq    =4e-10;
lamda=130e-9;
gama=lamda./(lamda+L);
eps0    = 8.85e-12;
Cox    =2*pi*eps0*kox./log(2*s./d*sinh(2*pi*tox./s))
Rmetal=465.12;
E=[4.2e-20 7.5e-20;5.5e-20 8.5e-20;6.5e-20 1.06e-20];
vds0= 0.05;
vg = -0.71:0.01:0.7;
% Vcnt    =Cox./(Cox+Cq)*vg;
clear
clc
Vcnt    =-0.71:0.01:0.7;
q= 1.6e-19;
vds0=0.05;
c=0.511;
E=[4.2e-20 7.5e-20;5.5e-20 8.5e-20;6.5e-20 1.06e-20];
    if Vcnt > 0

        phi_s11     = Vcnt;
        phi_s12     = E(2)-(1-exp(c))*q*Vcnt;
        phi_s13     = E(3)-(1-exp(c))*q*Vcnt;
        phi_s21     = E(4)-(1-exp(c))*q*Vcnt;
        phi_s22     = E(5)-(1-exp(c))*q*Vcnt;
        phi_s23     = E(6)-(1-exp(c))*q*Vcnt;
  else
        phi_s11     = E(1)-q*Vcnt;
        phi_s12     = E(2)-q*Vcnt;
        phi_s13     = E(3)-q*Vcnt;
        phi_s21     = E(4)-q*Vcnt;
        phi_s22     = E(5)-q*Vcnt;
        phi_s23     = E(6)-q*Vcnt;
   end
   if Vcnt > vds0
        phi_d11     = E(1)-(1-exp(c))*q*Vcnt-exp(c)*q*vds0;
        phi_d12     = E(2)-(1-exp(c))*q*Vcnt-exp(c)*q*vds0;
        phi_d13     = E(3)-(1-exp(c))*q*Vcnt-exp(c)*q*vds0;   
        phi_d21     = E(4)-(1-exp(c))*q*Vcnt-exp(c)*q*vds0;
        phi_d22     = E(5)-(1-exp(c))*q*Vcnt-exp(c)*q*vds0;
        phi_d23     = E(6)-(1-exp(c))*q*Vcnt-exp(c)*q*vds0;
  else
        phi_d11     = Vcnt;
        phi_d12     = E(2)-q*Vcnt;
        phi_d13     = E(3)-q*Vcnt;
        phi_d21     = E(4)-q*Vcnt;
        phi_d22     = E(5)-q*Vcnt;
        phi_d23     = E(6)-q*Vcnt;

   end

  IE11     =4*q*KB*T./h*(log(1+exp(-phi_s11./KB./T))-log(1+exp(-(phi_d11+q*vds0)./KB./T)));
  IE12     =4*q*KB*T./h*(log(1+exp(-phi_s12./KB./T))-log(1+exp(-(phi_d12+q*vds0)./KB./T)));
  IE13     =4*q*KB*T./h*(log(1+exp(-phi_s13./KB./T))-log(1+exp(-(phi_d13+q*vds0)./KB./T)));
  IE21     =4*q*KB*T./h*(log(1+exp(-phi_s21./KB./T))-log(1+exp(-(phi_d11+q*vds0)./KB./T)));
  IE22     =4*q*KB*T./h*(log(1+exp(-phi_s22./KB./T))-log(1+exp(-(phi_d12+q*vds0)./KB./T)));
  IE23     =4*q*KB*T./h*(log(1+exp(-phi_s23./KB./T))-log(1+exp(-(phi_d13+q*vds0)./KB./T)));
%   K=[0.651 0.448;0.309 0.808;0.622 0.694];
%   A=[0.651*IE11 0.448*IE21;0.309*IE12 0.808*IE22;0.622*IE13 0.694*IE23];
%   I=K.*A;
I=0.651*IE11+0.448*IE21+0.309*IE12+0.808*IE22+0.622*IE13+0.694*IE23;
Id=Ncnt*gama*I+vds0./Rmetal;

figure(1);
plot(vg,Id,'linewidth',2);
set(gca,'linew',1,'fontsize',15)
ylabel('(log) I_d (A/\mum)','fontsize',20)
xlabel('V_g_s (V)','fontsize',20)