f3和f4每次计算值不一样,有时0有时正有时负,或者帮我看看是我模型的问题,还是程序的问题?求大佬帮帮我。
Constant
q=1.6e-19,
mox=0.42*9.11e-31,
h=6.63e-34,
d=7.4e-7,
mm=0.1*9.11e-31,
F1=0.2,
ss=5.2,
s0=8.85E-14,
Tem=298,
K=1.38E-23,
pi=3.1415926,
Fb0=2.5,Ft=1.48,Nt=1.5e19,
E=1.5e6,
//t=2.2,
A=8*pi*sqrt(2*q*mox)/(3*h),
Fb=Fb0-sqrt((q*E)/(4*pi*ss*s0)),
Ct=(mm/mox)^(5/2)*(16*pi*q*F1^(3/2))/(3*h*sqrt(Ft-F1)),
C=(q*Ct*Nt)/E;
conststr
P1=exp((-A/E)*(t^(3/2)-Ft^(3/2))),
//P2_tri=exp((-A/E)*Ft^(3/2)),
//P2_tra=exp((-A/E)*(Ft^(3/2)-(t-E*d)^(3/2))),
fFD=1/(1+exp(q*(Fb-t)/(K*Tem)));
C1=1E160,C2=1E-150;
FUNCTION
//f0=A,
//f1=C2,
//f2=int(exp(x^(3/2)),x=1,2),
f3=(8*pi*sqrt(2*q*mox)/(3*h))*int(1/(1/(1/(1+exp(q*(Fb-t)/(K*Tem)))*exp((-A/E)*(t^(3/2)-Ft^(3/2))))+1/exp((-A/E)*(Ft^(3/2)-(t-E*d)^(3/2)))),t=Ft,2.5),
//f=1e200*Int((1/(fFD*P1)+1/P2_tra)^(-1),t=Ft,Ft+E*d);
f4=Int((1/((1/(1+exp(q*(Fb-t)/(K*Tem))))*(exp((-A/E)*(t^(3/2)-Ft^(3/2)))))+1/(exp((-A/E)*(Ft^(3/2)-(t-E*d)^(3/2)))))^(-1),t=Ft,Ft+E*d); |