processadados.m

% ##### Variaveis comuns:
T = 1;
omega = 1;
Kb = 1.380649e-23;
hbar = 1.0545718e-34;
beta = 1/(Kb*T);

arquivos = dir('*.csv'); % lista todos os arquivos terminados em .csv no...
% diretorio atual (resulta num vetor de structs, um elemento pra cada arquivo)

for i = 1:size(arquivos,1) % roda o codigo para cada arquivo
  xtemp = csvread(arquivos(i).folder); % pega a localizacao do arquivo e le...
  % os dados salvando numa matriz temporaria
  xtemp = xtemp(2:end,2:end); % tira a primeira linha e primera coluna

  y = sum(sum(xtemp)); % calcula o modulo/soma
  xtemp = xtemp/y; % divide pelo módulo

  for i = 1:size(xtemp,1)
    for j = 1:size(xtemp,2)
      elli = i - (size(xtemp,1)-1)/2;
      ellj = j - (size(xtemp,2)-1)/2;
      W(i,j) = (ellj-elli)*hbar*omega; // matriz com trabalho das transicoes
    endfor
  endfor

  ebW = exp(beta*W);

  val_med_ebW(i) = sum(sum(xtemp.*ebW)); %soma todos os elementos e multiplica cada termo do pll com o exp(bW)%
endfor

disp(val_med_ebW) % imprime os valores médios finais pros arquivo