




版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、 原創(chuàng) 用 Monte-Carlo 方法研究二維 Ising 模型的相變問題, 求出臨界溫度 T,比熱 C 及磁化率 圖片:圖片:/*程序用途:用 Monte-Carlo 方法研究二維 Ising 模型的相變問題,求出臨界溫度T,比熱 C及磁化率*編譯環(huán)境: Visual C+ 6.0*作 者:lxg* 單 位:*完成時(shí)間: 2005 年 12 月 28 日*/#include#include#include#include#define L 150#define preN 100000#define lateN 5100000#define deltaN (lateN-preN)/*L L
2、個(gè)格點(diǎn)*/* 舍棄前 preN 個(gè) E, M 值*/*總共計(jì)算 lateN 個(gè) E, M 值*/*Monte-Carlo 方法抽樣 E ,M 個(gè)數(shù) */double recordEnergydeltaN; int recordMagnetdeltaN;int latticeLL;int tempLL;double energy;int magnet;double T;struct CXdouble C;double X;/* 存儲(chǔ)抽樣能量值/* 存儲(chǔ)抽樣磁矩量/*格點(diǎn),存儲(chǔ)格點(diǎn)狀態(tài)/*保留初始化格點(diǎn)狀態(tài)/*記錄初試化及試探后系統(tǒng)能量,磁矩*/*/*/* 并賦給 recordEnergy 和 r
3、ecordMagnet /* 溫度*/* 比熱 C 和磁化率 */*/*/*/* *函數(shù)原型*/ struct CX Caculate_C_and_X();void voidvoidvoid void voidInitialLattice();SaveInitial();GetInitial();CaculateFirstMagnet();CaculateFirstEnergy();GenerateNext();/* *主函數(shù)*/ int main()int i;struct CX cx200; FILE *ft,*fx,*fc; ft=fopen(T.txt,w); fx=fopen(X.t
4、xt,w); fc=fopen(C.txt,w);srand(time(NULL);T = 1.0;for(i=0;T4.0;i+)InitialLattice();cx =Caculate_C_and_X(); fprintf(ft,%fn,T); fprintf(fx,%fn,cx.X); fprintf(fc,%fn,cx.C);printf(T=%f ,T); printf( X=%f ,cx.X); printf( C=%fn,cx.C);T = T+0.05;fclose(ft); fclose(fx); fclose(fc); return 0; /* 保存第一次初始化狀態(tài) */
5、 /* 若為冷初試化,可以不 */ /*調(diào)用該函數(shù)*/* 提取第一次初始化狀態(tài) */ /* 若為冷初試化,可以不 */ /*調(diào)用該函數(shù)*/*函數(shù)實(shí)現(xiàn) */ void InitialLattice()/*初始化格點(diǎn)狀態(tài),可以 */* 通過熱初始化和冷初始化 */int i,j;for(i=0;iL;i+)for(j=0;j5) latticej=1;else latticej=1;/*1 冷初始化, -1 熱初試化 */void SaveInitial()int k,m;for(k=0;kL;k+)for(m=0;mL;m+) tempkm=latticekm;void GetInitial()
6、int k,m;for(k=0;kL;k+)for(m=0;mL;m+)latticekm=tempkm;void CaculateFirstMagnet()/* 計(jì)算初試化后系統(tǒng)磁矩 */int i,j;int iMagnet=0;for(i=0;iL;i+) for(j=0;jL;j+) iMagnet += latticej;magnet = iMagnet;void CaculateFirstEnergy()/* 計(jì)算初試化后系統(tǒng)能量 */double dEnergy = 0.0;int up,down,right,left,px,py;for(px=0;pxL;px+)for(py=
7、0;py0 )/*判斷是否接受試探 */if(rand()*1.0/RAND_MAX) = exp(-1.0)*dE/T)latticepxpy = -latticepxpy; magnet += 2*latticepxpy; energy += dE;/* 更新磁矩大小 */* 更新能量大小 */elselatticepxpy = -latticepxpy; magnet += 2*latticepxpy; energy += dE;/*更新磁矩大小 */ /* 更新能量大小 */struct CX Caculate_C_and_X()struct CX cx;int i;double Av
8、gEnergy=0.0;double AvgSquareEnergy=0.0;double AvgMagnet=0.0;double AvgSquareMagnet=0.0;/*計(jì)算比熱 C及磁化率 */*初始化磁矩 */*初始化能量 */CaculateFirstMagnet();CaculateFirstEnergy();for(i=0; i=preN)recordEnergyi-preN/*舍棄前preN 個(gè) */= energy;/* 記錄能量值 */recordMagneti-preN= magnet;/* 記錄磁矩值 */for(i=0;ideltaN;i+)/* 計(jì)算物理量的平均值 */AvgEnergy += (double)recordEnergy/deltaN;AvgSquareEnergy += (double)recordEnergy*recordEnergy/deltaN;AvgMagnet += (double)recordMagnet/deltaN;AvgSquareMagnet += (double)recordMagnet*recordMagnet/deltaN;cx.
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 精通營(yíng)養(yǎng)師考試的關(guān)鍵試題及答案
- 2012年1月國(guó)家開放大學(xué)漢語言文學(xué)本科《心理學(xué)》期末紙質(zhì)考試試題及答案
- 縣城任教招考題及答案大全
- 物理中考題及答案江蘇
- 吉林省延邊州安圖縣聯(lián)考2023-2024學(xué)年中考數(shù)學(xué)押題試卷含解析
- 融資面試試題及答案
- 江蘇旅游職業(yè)學(xué)院《城鄉(xiāng)社會(huì)綜合調(diào)查與研究》2023-2024學(xué)年第二學(xué)期期末試卷
- 伊通滿族自治縣2024-2025學(xué)年三下數(shù)學(xué)期末考試模擬試題含解析
- 四川省廣元天立國(guó)際學(xué)校2024-2025學(xué)年高三下學(xué)期第一次月考文綜試卷含解析
- 蘇州百年職業(yè)學(xué)院《醫(yī)學(xué)細(xì)胞及分子生物學(xué)》2023-2024學(xué)年第二學(xué)期期末試卷
- 血液透析個(gè)案護(hù)理
- 2025年電工(技師)考試100題及答案
- 《遙感原理與應(yīng)用》全冊(cè)配套完整課件
- GB/T 4706.28-2024家用和類似用途電器的安全第28部分:吸油煙機(jī)及其他烹飪煙氣吸排裝置的特殊要求
- 部隊(duì)食品采購(gòu)合同范本
- 2024江蘇蘇州市滸墅關(guān)經(jīng)開區(qū)招聘村(社區(qū))人員筆試歷年典型考題及考點(diǎn)剖析附答案帶詳解
- 華東政法大學(xué)2017-2018學(xué)年期末測(cè)試《國(guó)際經(jīng)濟(jì)法概論》試卷
- 文化創(chuàng)意與策劃實(shí)訓(xùn)報(bào)告總結(jié)
- 江城鎮(zhèn)矛盾糾紛排查調(diào)處考核獎(jiǎng)懲辦法(試行)
- 2024年江西農(nóng)業(yè)工程職業(yè)學(xué)院?jiǎn)握新殬I(yè)適應(yīng)性測(cè)試題庫(kù)含答案
- 工程文化概論機(jī)械工程與文化
評(píng)論
0/150
提交評(píng)論