matlab 关于数值积分的问题(定积分但是数值很大且式子复杂)X=0.31;d=2*10^-6;E0=8.85*10^-12;T=77;e=1.6*10^-19;k=1.38*10^23;kb=8.625*10^-5;ND=5*10^20;NA=100*10^20;E=E0*(14*(1-X)+7.05*X)Eg=-0.295+1.87*X-0.28*X^2+(6-14*X+3*X^2)*10^-4*T
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/02 14:27:31
matlab 关于数值积分的问题(定积分但是数值很大且式子复杂)X=0.31;d=2*10^-6;E0=8.85*10^-12;T=77;e=1.6*10^-19;k=1.38*10^23;kb=8.625*10^-5;ND=5*10^20;NA=100*10^20;E=E0*(14*(1-X)+7.05*X)Eg=-0.295+1.87*X-0.28*X^2+(6-14*X+3*X^2)*10^-4*T
matlab 关于数值积分的问题(定积分但是数值很大且式子复杂)
X=0.31;
d=2*10^-6;
E0=8.85*10^-12;
T=77;
e=1.6*10^-19;
k=1.38*10^23;
kb=8.625*10^-5;
ND=5*10^20;
NA=100*10^20;
E=E0*(14*(1-X)+7.05*X)
Eg=-0.295+1.87*X-0.28*X^2+(6-14*X+3*X^2)*10^-4*T+0.35*X^4
nl=(1+3.25*kb*T/Eg)*9.56*10^14*Eg.^(3/2)*T.^(3/2)*10^6/(1+1.9*Eg.^(3/4)*exp(Eg/2/kb/T))
pn0=(-ND+(ND^2+4*nl^2)^(1/2))/2
nn0=(ND+(ND^2+4*nl^2)^(1/2))/2
pp0=(NA+(NA^2+4*nl^2)^(1/2))/2
np0=(NA+(NA^2+4*nl^2)^(1/2))/2
Vpn=kb*T*log(ND*NA/nl^2)/e
V0=(nn0+pn0-np0-pp0-ND*e*Vpn/k/T)*kb/e/(NA-ND)
E=-(2*kb*T*(-e*x*ND/kb/T+np0*exp(e*x/kb/T)+pp0*exp(e*x/kb/T)-nn0-pn0+ND*e*Vpn/k/T)/E).^(1/2);
对1/E求积分,上限是V0,下限是Vpn;
matlab 关于数值积分的问题(定积分但是数值很大且式子复杂)X=0.31;d=2*10^-6;E0=8.85*10^-12;T=77;e=1.6*10^-19;k=1.38*10^23;kb=8.625*10^-5;ND=5*10^20;NA=100*10^20;E=E0*(14*(1-X)+7.05*X)Eg=-0.295+1.87*X-0.28*X^2+(6-14*X+3*X^2)*10^-4*T
clear;clc
syms x
X=0.31;
d=2*10^-6;
E0=8.85*10^-12;
T=77;
e=1.6*10^-19;
k=1.38*10^23;
kb=8.625*10^-5;
ND=5*10^20;
NA=100*10^20;
E=E0*(14*(1-X)+7.05*X);
Eg=-0.295+1.87*X-0.28*X^2+(6-14*X+3*X^2)*10^-4*T+0.35*X^4;
nl=(1+3.25*kb*T/Eg)*9.56*10^14*Eg.^(3/2)*T.^(3/2)*10^6/(1+1.9*Eg.^(3/4)*exp(Eg/2/kb/T));
pn0=(-ND+(ND^2+4*nl^2)^(1/2))/2;
nn0=(ND+(ND^2+4*nl^2)^(1/2))/2;
pp0=(NA+(NA^2+4*nl^2)^(1/2))/2;
np0=(NA+(NA^2+4*nl^2)^(1/2))/2;
Vpn=kb*T*log(ND*NA/nl^2)/e;
V0=(nn0+pn0-np0-pp0-ND*e*Vpn/k/T)*kb/e/(NA-ND);
E=-(2*kb*T*(-e*x*ND/kb/T+np0*exp(e*x/kb/T)+pp0*exp(e*x/kb/T)-nn0-pn0+ND*e*Vpn/k/T)/E).^(1/2);
U=double(int(1/E,Vpn,V0))