请问这个积分在matlab里用数值积分的方法应该如何做?
发布网友
发布时间：2023-07-12 16:04
共1个回答
热心网友
时间：2024-03-16 23:22
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);
