level 1
代码原文
data = Import["Sigma2_realcm_eventbyevent_2M.txt", "Table"];
avQxy2alpha = {};
Filelength = 1000;
Do[avQxy2alpha = AppendTo[avQxy2alpha, data[[i, 2]]], {i, 1,
Filelength}];
avNalpha = Table[4, 1 + Filelength];
g = 1;
hbar = 197.3;
h = 197.3*2*Pi;
m = 4*938.6;
coe = g/h^3*4 Pi*(2 m)^1.5/2;
q1 = 2.0;
q2 = 8.0 - q1;
A = 1.44*4 Pi*hbar^2*q1*q2/2/m;
avN[V_, T_] :=
coe*V*T^1.5*
NIntegrate[x^0.5/(Exp[x + A/x/V/T^2] - 1), {x, 0, Infinity}];
avQxy2[V_, T_] := (2*m*T)^2*4/15*
NIntegrate[x^2.5/(Exp[x + A/x/V/T^2] - 1), {x, 0, Infinity}]/
NIntegrate[x^0.5/(Exp[x + A/x/V/T^2] - 1), {x, 0, Infinity}];
tab = {{20, 10}};
Do[tab = Append[tab,
Table[{V, T} /.
FindRoot[{avN[V, T] == avNalpha[[i]],
avQxy2[V, T] == avQxy2alpha[[i]]}, {V, 80}, {T, 10}]]], {i, 1,
Filelength}];
a = avNalpha/tab[[;; , 1]];
b = {};
Do[b = Append[b, {tab[[i + 1, 2]], a[[i + 1]]}], {i, 1, Filelength}];
filename = "outdata.txt"
Export["outdata.txt", b, "Table"];
2018年11月09日 01点11分
3
level 1
@xzcyr上面那个导入txt的帖子是你人工置顶的吗??我看到了,谢谢
2018年11月09日 02点11分
4
level 1
学习了,这次用这个语句导入数据,希望运算速度能有所加快
ReadList["ExampleData/numbers", {Number, Number}]
2018年11月09日 03点11分
5
吧务
level 11
用二进制保存 然后ReadBinaryList读取最快了
2018年11月11日 05点11分
8