| 9 | 1/1 | 返回列表 |
| 查看: 4518 | 回復: 8 | |||||||||||
guolianshun銅蟲 (著名寫手)
|
[交流]
用Origin畫MS中電子云密度的技術。!超漂亮 已有8人參與
|
|
用Origin畫MS中電子云密度的技術!!超漂亮 有心人可以研究一下。! 版主請加精啊, 費了好大力氣的! 呵呵! 歡迎大家學習! MATLAB程序用來處理MS電荷密度的等高線做法,分四步: 1.找到電荷密度圖的存儲文件,一般是在自己所建的計算project中,隱藏的,castep模塊是.charg_frm格式的,可以用txt打開后另存為txt格式。(這個地方注意把文件單獨拷貝出來,以免破壞原始文件) 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是可以改的,分別對應a,b,c的重復單元。 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是可以改的,分別對應a,b,c的重復單元,即長寬高的值。 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(:,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); 一個牛人的總結(jié),轉(zhuǎn)自分子論壇。 |
第一性原理 | 物理前沿 | 第一性原理計算 | @個人收集@ |
Photochemistry | 計算化學相關 | 精華帖子 | 第一性 |
|
本帖內(nèi)容被屏蔽 |
金蟲 (小有名氣)

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

新蟲 (小有名氣)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|---|---|---|---|---|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +4 | Charlieyq 2026-03-19 | 4/200 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 8/400 |
|
|
[考研] 081700化工學碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊 2026-03-18 | 5/250 |
|
|
[考研] 0703化學調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 一志愿 西北大學 ,070300化學學碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 化學工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 0703化學調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|