| 5 | 1/1 | 返回列表 |
| 查看: 4533 | 回復: 8 | |||||||||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||||||||||
guolianshun銅蟲 (著名寫手)
|
[交流]
用Origin畫MS中電子云密度的技術(shù)。!超漂亮 已有8人參與
|
||||||||||
|
用Origin畫MS中電子云密度的技術(shù)!!!超漂亮 有心人可以研究一下!!! 版主請加精啊, 費了好大力氣的! 呵呵! 歡迎大家學習! 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是可以改的,分別對應(yīng)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是可以改的,分別對應(yīng)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 | 計算化學相關(guān) | 精華帖子 | 第一性 |
|
本帖內(nèi)容被屏蔽 |
金蟲 (小有名氣)

木蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +8 | 紙魚ly 2026-03-21 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 材料專業(yè)求調(diào)劑 +11 | hanamiko 2026-03-18 | 11/550 |
|
|
[考研] 333求調(diào)劑 +6 | 87639 2026-03-21 | 10/500 |
|
|
[考研] 352求調(diào)劑 +3 | 大米飯! 2026-03-22 | 3/150 |
|
|
[考研]
|
2117205181 2026-03-21 | 8/400 |
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
[考研] 一志愿深大,0703化學,總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
[考研] 一志愿南昌大學,327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 一志愿福大288有機化學,求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|