我想选择合适的函数以适合我的数据。这是我尝试过的:
myTolerance=1e-12; % in erg units.
nbins=100;
for j=1:ntM/10:ntM
data = Wkinet(abs(Wkinet(:,j))>myTolerance, j) * erg2eV;
pd = fitdist(data,'Normal');
x = linspace(min(data),max(data),1000);
y = pdf(pd,x);
loglog(x,y);
%%hold on;
pause;
end
这就是我得到的:
这就是我想要的样子:
在仿真开始时,分布是Maxwelian,最后是幂定律,如您所见。