нормализованная гистограмма - PullRequest
0 голосов
/ 11 мая 2019

Здравствуйте, я пытаюсь создать гистограмму переменной Wkinet для различных времен j (это то, как распределение энергии некоторых частиц изменяется во времени), но я хочу 1) для каждого временного шага, не учитывающего ноль значения 2) Возьмите в результате все исторические данные для каждого значения времени в одной диаграмме, но только подгонку, а не столбцы. 3) Я также хотел бы, чтобы она нормализовалась по осям y

 nbins=100;
for j=1:ntM

 idxToKeep   = Wkinet(:,j ) ~= 0;
   Wkinet2     = Wkinet(idxToKeep,:);
   y=Wkinet2(:,j).*erg2eV;
  histfit(y,nbins)
  hold on
    delete(h(2))% deletes histogram
     end
      hold off
   end

1 Ответ

0 голосов
/ 05 июня 2019
        erg2eV = 6.2415e11;
        ntM=500;
          myTolerance=1e-15; % in erg units.
        nbins=1000;
    for j=1:ntM/10:ntM
        hold on
     myFit = fitdist((Wkinet_eV(Wkinet_eV(:,j)>myTolerance, j) ), 'Normal');
         maxEv = max(Wkinet_eV(Wkinet_eV(:,j)>myTolerance, j)); 
         minEv = min(Wkinet_eV(Wkinet_eV(:,j)>myTolerance, j)); %also

            myX   = minEv: (maxEv-minEv)/nbins: maxEv; 
            myY   = pdf(myFit, myX);
            loglog((myX),( myY)); %maybe just plot or %plot(log10(myX),myY)
                  ylim([10^-13,1])
         pause;
   end
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...