| 9 | 1/1 | 返回列表 |
| 查看: 4520 | 回復(fù): 8 | |||||||||||
guolianshun銅蟲(chóng) (著名寫手)
|
[交流]
用Origin畫MS中電子云密度的技術(shù)。!超漂亮 已有8人參與
|
|
用Origin畫MS中電子云密度的技術(shù)。!超漂亮 有心人可以研究一下。! 版主請(qǐng)加精啊, 費(fèi)了好大力氣的! 呵呵! 歡迎大家學(xué)習(xí)! MATLAB程序用來(lái)處理MS電荷密度的等高線做法,分四步: 1.找到電荷密度圖的存儲(chǔ)文件,一般是在自己所建的計(jì)算project中,隱藏的,castep模塊是.charg_frm格式的,可以用txt打開(kāi)后另存為txt格式。(這個(gè)地方注意把文件單獨(dú)拷貝出來(lái),以免破壞原始文件) 2.讀取數(shù)據(jù),在matlab命令里輸入一下命令: [path,fn]=uigetfile('*.txt','Open'); fp=fopen([fn,path],'r'); head=fscanf(fp,'%s',4); data=fscanf(fp,'%f',[4,18*18*144]); fclose(fp); 這是讀取剛才保存的txt數(shù)據(jù)的,其中18*18*144是可以改的,分別對(duì)應(yīng)a,b,c的重復(fù)單元。 3.作矩陣。 a=reshape(data(1,,[18,18,144]); b=reshape(data(2,,[18,18,144]); c=reshape(data(3,,[18,18,144]); d=reshape(data(4,,[18,18,144]); 其中18*18*144是可以改的,分別對(duì)應(yīng)a,b,c的重復(fù)單元,即長(zhǎng)寬高的值。 4.作等高線圖 layer=1; contour(reshape(a(:,layer,,[18,144]),reshape(c(:,layer,,[18,144]),reshape(d(:,layer,,[18,144]),8); layer=1(2,3,4....) 可以修改作圖的層數(shù),分別是從前面到后面。最后的8表示顯示線條數(shù)目,可以自己修改成9,10或其它。 hold on 表示累積作圖,hold off表示清空前面的 作圖 舉例;做表面層和第九層的等高線圖 >> layer=9; contour(reshape(a(:,layer,,[18,144]),reshape(c(:,layer,,[18,144]),reshape(d(:,layer,,[18,144]),8); >> hold on >> layer=1; contour(reshape(a(:,layer, ,[18,144]),reshape(c(:,layer, ,[18,144]),reshape(d(:,layer, ,[18,144]),8);這是沿側(cè)面做等高線(即xz平面),其它方向作圖類似。只要修改其中的參數(shù)就可以了,如把a(bǔ)(:,layer, ,[18,144]),reshape(c(:,layer, ,[18,144]),reshape(d(:,layer, ,[18,144]),8);改為a(layer,:, ,[18,144]),reshape(c(layer,l:, ,[18,144]),reshape(d(layer,:, ,[18,144]),8); 一個(gè)牛人的總結(jié),轉(zhuǎn)自分子論壇。 |
第一性原理 | 物理前沿 | 第一性原理計(jì)算 | @個(gè)人收集@ |
Photochemistry | 計(jì)算化學(xué)相關(guān) | 精華帖子 | 第一性 |
|
本帖內(nèi)容被屏蔽 |
金蟲(chóng) (小有名氣)

金蟲(chóng) (著名寫手)
送鮮花一朵
![]() ![]() |
木蟲(chóng) (小有名氣)
金蟲(chóng) (小有名氣)

新蟲(chóng) (小有名氣)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 324求調(diào)劑 +3 | lucky呀呀呀鴨 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +3 | 周周無(wú)極 2026-03-20 | 3/150 |
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 14/700 |
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變,給編輯又寫郵件了,沒(méi)回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +9 | 一鴨鴨喲 2026-03-14 | 11/550 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 化學(xué)求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
[考研] 【同濟(jì)軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|