Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- % ##### 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
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement