Сбор статистики из большого набора файлов

Question

Есть ли более быстрый способ сделать это

function [A,B]=powerplot(basepath,N) 
    A=zeros(N,1); 
    B=zeros(N,1); 
    for k=0:N-1 
     filename=sprintf('%s-%d.txt',basepath,k); 
     data=load(filename); 
     [a,b]=powerfit(data(:,1),data(:,2)); 
     A(k+1)=a; 
     B(k+1)=b; 
    end 
end 

где

function [a,b]=powerfit(X,Y) 
    p=polyfit(log(X),log(Y),1); 
    a=exp(p(2)); 
    b=p(1); 
end 

Файлы содержат измерение на ансамбле, каждый файл захваченного в разное время, содержащем разное количество строк. Могу ли я лучше упорядочить свои данные? Чтобы ускорить его, нужно ли делать этот анализ в собственном коде с помощью библиотеки ROOT?

Answer 1

Проблема была load. Заменив load на dlmread, код стал намного быстрее.

function [A,B]=powerplot(basepath,N) 
    A=zeros(N,1); 
    B=zeros(N,1); 
    for k=0:N-1 
     filename=sprintf('%s-%d.txt',basepath,k); 
     data=dlmread(filename); 
     [a,b]=powerfit(data(:,1),data(:,2)); 
     A(k+1)=a; 
     B(k+1)=b; 
    end 
end