ボルツマンマシン

メトロポリス法で，ボルツマンマシンからサンプリングするコードを書きました． とりあえず，状態xを与えたら，エネルギーを返してくれる関数を用意しておきます． using Plots using Random function energy(x,b,w) N = length(b) term_1_body = sum( x .* b ) term_2_body = sum( w[i,j] * x[i] * x[j] for i = 1:N for j = i:N) E = term_1_body + term_2_body return E end function delta_energy(x,…