求蒙特卡罗方法解一维ising模型的模拟程序,c或c++的,外场为零
来源:学生作业帮 编辑:大师作文网作业帮 分类:综合作业 时间:2024/11/13 14:33:19
求蒙特卡罗方法解一维ising模型的模拟程序,c或c++的,外场为零
我写过一个用mathematica的程序,比较简单.要的话可以给你.
再问: 好的,麻烦上传一下
再答: Clear[n, JT, s, M, Z, A, de]; n = 6; JT = .5; s = Table[1, {n}, {n}];(*选取初始configuration为6*6的自旋全部向上*) M = n^2.;(*M标记update过程中的总的磁化强度的绝对值,初始为36*) Z = 0;(*标记随机求和所进行的次数,初始为0*) data = Table[(A = 0; k = 10000; Do[ Module[{k1, k2, \[Xi]}, k1 = RandomInteger[{1, 6}]; k2 = RandomInteger[{1, 6}]; \[Xi] = .1 Random[]; de = 0; If[k1 > 1, de = de + s[[k1 - 1, k2]], de]; If[k1 < 6, de = de + s[[k1 + 1, k2]], de]; If[k2 > 1, de = de + s[[k1, k2 - 1]], de]; If[k2 < 6, de = de + s[[k1, k2 + 1]], de]; R = Exp[-2 JT de s[[k1, k2]]]; If[R > \[Xi], s[[k1, k2]] = -s[[k1, k2]], s[[k1, k2]]]; If[R > \[Xi], M = M + 2 s[[k1, k2]], M]; A = A + Abs[M]; ], {k}]; N[A/k]), {30}] 做图: Show[{Plot[Mean[data], {x, 0, 30}, PlotStyle -> Red], ListLinePlot[data]}, PlotRange -> {4, 7}, Frame -> True]
再问: 好的,麻烦上传一下
再答: Clear[n, JT, s, M, Z, A, de]; n = 6; JT = .5; s = Table[1, {n}, {n}];(*选取初始configuration为6*6的自旋全部向上*) M = n^2.;(*M标记update过程中的总的磁化强度的绝对值,初始为36*) Z = 0;(*标记随机求和所进行的次数,初始为0*) data = Table[(A = 0; k = 10000; Do[ Module[{k1, k2, \[Xi]}, k1 = RandomInteger[{1, 6}]; k2 = RandomInteger[{1, 6}]; \[Xi] = .1 Random[]; de = 0; If[k1 > 1, de = de + s[[k1 - 1, k2]], de]; If[k1 < 6, de = de + s[[k1 + 1, k2]], de]; If[k2 > 1, de = de + s[[k1, k2 - 1]], de]; If[k2 < 6, de = de + s[[k1, k2 + 1]], de]; R = Exp[-2 JT de s[[k1, k2]]]; If[R > \[Xi], s[[k1, k2]] = -s[[k1, k2]], s[[k1, k2]]]; If[R > \[Xi], M = M + 2 s[[k1, k2]], M]; A = A + Abs[M]; ], {k}]; N[A/k]), {30}] 做图: Show[{Plot[Mean[data], {x, 0, 30}, PlotStyle -> Red], ListLinePlot[data]}, PlotRange -> {4, 7}, Frame -> True]
求蒙特卡罗方法解一维ising模型的模拟程序,c或c++的,外场为零
阿尔法粒子散射的模拟实验c语言程序
C程程序填空题:输入一批正整数(以零或负数为结束标志),求其中的偶数和.
一道C语言题目.编写一个模拟两个人投掷骰子,比大小的的程序.
若 x≤b x>a,则(a-x)(x-b)的符号为 A 大于零 B 大于或等于零 C 小于零 D 小于或等于零
C程序中的一个变量可以被定义为两个或多个不同的类型.这句话对吗?为什么
用生成2维随机数的方法求pi的值.得出误差与n之间的关系曲线.c或c++程序
关于C语言,以下程序的运算结果为()
模拟电路试题6.理想二极管的正向电阻为( )A.零 B.无穷大 C.约几千欧 D.约几十欧
已知a,b,c为非零的有理数,...
已知:a、b、c分别为△ABC的三条边的长度,请用所学知识说明:(a-c)2-b2是正数、负数或零.
C语言,阅读下列程序,则在执行后,程序的运行结果为 ?