作者Edwardkiller (台南金城武)
看板MATLAB
标题[问题] 手算没负值 Matlab却会算出负值
时间Mon Nov 7 15:43:16 2011
k=1.38066e-23; q=1.60218e-19; un=1500; up=450; E=11.9*8.854e-14;
Na=1e15; Nd=1e18; T=298; s=1e-7; G=1e13;
x=0:2e-6:2e-4
ni2=(9.38e19*(T/300)^2.*exp(-6884/T))^2;
Pn2=(ni2/Nd*(exp(q*(-5)/(k*T))-1)-G*s)*exp(-x/0.0011)+G*s+ni2/Nd;
问题在x=0时，算出负值，根本不合理啊
x=0时实际上只剩下ni2/Nd*exp(q*(-5)/(k*T))这一项
结果Matlab算出负的
请教各强者要如何做修正才会对
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 122.117.163.182
1F：推 alamabarry:我觉得 你要不要先无因次画一下 11/07 17:46
2F：→ alamabarry:还有正负号可能要看你是不是浮点运算的误差所致 11/07 17:47