Question

请问这个积分在matlab里用数值积分的方法应该如何做？

积分下限是Eg，上限是10Eg，我用的数值积分方法，但是做不出来Eg=1.519-5.405*(10^(-4))*((T_0.^2)./(T_0+204));T_0=linspace(300,1500,1201);Tc=repmat(T_0,1e4,1);E0=repmat(Eg,1e4,1);k=1.38*10^(-23);C0=linspace(1,10,1e4);C=repmat(C0',1,1201);hv=E0.*C;R_matrix=hv.^2./(exp(hv./k./T)-1);R=sum(R_matrix)*10.0250719315895/1e4;plot(T_0,R);请问应该如何改进才能积出来？

Answer 1

T_0=linspace(300,1500,101);

%T_0=linspace(300,1500,1201);

Eg=1.519-5.405*(10^(-4))*((T_0.^2)./(T_0+204));

%T_0=linspace(300,1500,1201);

Tc=repmat(T_0,100,1);

%Tc=repmat(T_0,1e4,1);

E0=repmat(Eg,100,1);

%E0=repmat(Eg,1e4,1);

k=1.38*10^(-23);

C0=linspace(1,10,100);

%C0=linspace(1,10,1e4);

C=repmat(C0',1,101);

%C=repmat(C0',1,1201);

hv=E0.*C;

R_matrix=hv.^2./(exp(hv./k./T_0)-1);

%R_matrix=hv.^2./(exp(hv./k./T)-1);

R=sum(R_matrix)*10.0250719315895/100;

R1=sum(sum(R_matrix)*10.0250719315895/100)

plot(T_0,R);