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/24 05:18:10
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
xTMkA+=d>v7L&^Þt vi*ě-EQ{X ^ГJdS/|w&Ŧ zvv ײ̞_L݋W(;=2ĀO{ ;ώ^X̵r-,=ZҞȶTdlV6Vsq0t1[J~SB,Oy# ,Ra(*S+U+8ȝ*Ae gEܮ[7K.8"{ ӔhkY"Dұ+4*WW 3G? [IM\ C5V[@7`]rn_pPcq8Eހ;|e&K[-"X1s!iCi~͐N|Co\xdB`f (4q!a%G@xAѬaD#<~KoO _1CIG(Ћ,g6fow-L_-4u6SЏޖ}n)ZIIU4-W\ArYt-B/)Z߼lC+pQK?oi

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))