SHARE
TWEET

processadados.m

a guest Jul 19th, 2019 67 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. % ##### Variaveis comuns:
  2. T = 1;
  3. omega = 1;
  4. Kb = 1.380649e-23;
  5. hbar = 1.0545718e-34;
  6. beta = 1/(Kb*T);
  7.  
  8. arquivos = dir('*.csv'); % lista todos os arquivos terminados em .csv no...
  9. % diretorio atual (resulta num vetor de structs, um elemento pra cada arquivo)
  10.  
  11. for i = 1:size(arquivos,1) % roda o codigo para cada arquivo
  12.   xtemp = csvread(arquivos(i).folder); % pega a localizacao do arquivo e le...
  13.   % os dados salvando numa matriz temporaria
  14.   xtemp = xtemp(2:end,2:end); % tira a primeira linha e primera coluna
  15.  
  16.   y = sum(sum(xtemp)); % calcula o modulo/soma
  17.   xtemp = xtemp/y; % divide pelo módulo
  18.  
  19.   for i = 1:size(xtemp,1)
  20.     for j = 1:size(xtemp,2)
  21.       elli = i - (size(xtemp,1)-1)/2;
  22.       ellj = j - (size(xtemp,2)-1)/2;
  23.       W(i,j) = (ellj-elli)*hbar*omega; // matriz com trabalho das transicoes
  24.     endfor
  25.   endfor
  26.  
  27.   ebW = exp(beta*W);
  28.  
  29.   val_med_ebW(i) = sum(sum(xtemp.*ebW)); %soma todos os elementos e multiplica cada termo do pll com o exp(bW)%
  30. endfor
  31.  
  32. disp(val_med_ebW) % imprime os valores médios finais pros arquivo
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top