Ising

Z Skrypty dla studentów Ekonofizyki UPGOW

(Różnice między wersjami)
m (→Fast Ising (matlab))
m (→Fast Ising (matlab))
Linia 13: Linia 13:
[http://ocho.uwaterloo.ca/~pfieguth/Software/Ising/ising.html Paul Fieguth]
[http://ocho.uwaterloo.ca/~pfieguth/Software/Ising/ising.html Paul Fieguth]
 +
 +
Opis działania:
 +
 +
* załóżmy, że przeprowadzamy symulacje dla n=6
 +
 +
* definiowane są dwie tablice 8x8:
 +
 +
<source lang="matlab">
 +
t =
 +
 +
  0  0  0  0  0  0  0  0
 +
  0  1  0  1  0  1  0  0
 +
  0  0  1  0  1  0  1  0
 +
  0  1  0  1  0  1  0  0
 +
  0  0  1  0  1  0  1  0
 +
  0  1  0  1  0  1  0  0
 +
  0  0  1  0  1  0  1  0
 +
  0  0  0  0  0  0  0  0
 +
e = find(t);
 +
</source>
 +
<source lang="matlab">
 +
t =
 +
 +
  0  0  0  0  0  0  0  0
 +
  0  0  1  0  1  0  1  0
 +
  0  1  0  1  0  1  0  0
 +
  0  0  1  0  1  0  1  0
 +
  0  1  0  1  0  1  0  0
 +
  0  0  1  0  1  0  1  0
 +
  0  1  0  1  0  1  0  0
 +
  0  0  0  0  0  0  0  0
 +
o = find(t);
 +
</source>
 +
 +
dla których miejsca (indeksy) niezerowych elementów sa zapamiętanie w wektorach e i o.
===Problem czasu? ===
===Problem czasu? ===

Wersja z 10:15, 13 gru 2010

Model Isinga w ekonofizyce

Mamy Hamiltonian: \(H = -\sum_{i,j} J_{ij} s_i s_j - \mu B \sum_i s_i \, \)


Model Isinga dany jest przez ogólną

  • ferromagnetic coupling
  • anti-ferromagnetic coupling

Fast Ising (matlab)

Paul Fieguth

Opis działania:

  • załóżmy, że przeprowadzamy symulacje dla n=6
  • definiowane są dwie tablice 8x8:
t =
 
   0   0   0   0   0   0   0   0
   0   1   0   1   0   1   0   0
   0   0   1   0   1   0   1   0
   0   1   0   1   0   1   0   0
   0   0   1   0   1   0   1   0
   0   1   0   1   0   1   0   0
   0   0   1   0   1   0   1   0
   0   0   0   0   0   0   0   0
e = find(t);

t =
 
   0   0   0   0   0   0   0   0
   0   0   1   0   1   0   1   0
   0   1   0   1   0   1   0   0
   0   0   1   0   1   0   1   0
   0   1   0   1   0   1   0   0
   0   0   1   0   1   0   1   0
   0   1   0   1   0   1   0   0
   0   0   0   0   0   0   0   0
o = find(t);

dla których miejsca (indeksy) niezerowych elementów sa zapamiętanie w wektorach e i o.

Problem czasu?

Źródło „http://172.16.3.1/ekonofizyka/index.php/Ising”