 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4546 | 回復(fù): 8 | |||||||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||||||||||
guolianshun銅蟲 (著名寫手)
|
[交流]
用Origin畫MS中電子云密度的技術(shù)�。。〕� 已有8人參與
|
||||||||||
|
用Origin畫MS中電子云密度的技術(shù)�。�！超漂亮 有心人可以研究一下�。�！ 版主請加精啊， 費了好大力氣的！ 呵呵！ 歡迎大家學(xué)習(xí)！ 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的重復(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是可以改的，分別對應(yīng)a,b,c的重復(fù)單元，即長寬高的值。 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)自分子論壇。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一性原理 | 物理前沿 | 第一性原理計算 | @個人收集@ |
| Photochemistry | 計算化學(xué)相關(guān) | 精華帖子 | 第一性 |
» 本帖已獲得的紅花（最新10朵）
awmc2008» 猜你喜歡
309求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
求化學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分
已經(jīng)有5人回復(fù)
-
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600，專業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有16人回復(fù)
-
290求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- origin畫ms電荷密度圖
已經(jīng)有21人回復(fù)
- 【活動】Show 腳本，贏獎勵
已經(jīng)有23人回復(fù)
- 新手，請高手幫忙，計算小分子中某個碳原子上的電子云密度
已經(jīng)有15人回復(fù)
- 【討論】請教用origin畫一種圖
已經(jīng)有15人回復(fù)
- 一般用什么軟件計算離子的鍵能和電子云密度
已經(jīng)有11人回復(fù)
- 【求助】band.dat用origin畫能帶圖
已經(jīng)有5人回復(fù)
- 【求助】求電子云密度大小，滿意追加金幣
已經(jīng)有10人回復(fù)
- 用origin畫了個圖，似乎里面有Type 3 字體，生成pdf以后，提交的時候通不過。
已經(jīng)有22人回復(fù)
- 【求助】如何采用Matlab繪制電子云密度圖
已經(jīng)有8人回復(fù)
YUN杰2070
金蟲 (小有名氣)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 2376.3
- 紅花: 3
- 帖子: 119
- 在線: 95.6小時
- 蟲號: 1279373
- 注冊: 2011-04-27
- 性別: GG
- 專業(yè): 無機納米化學(xué)
3樓2011-10-28 21:20:32
★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
小木蟲(金幣+0.5):給個紅包，謝謝回帖
|
本帖內(nèi)容被屏蔽 |
2樓2011-10-11 15:40:54
風(fēng)之夢
木蟲 (小有名氣)
- 應(yīng)助: 37 (小學(xué)生)
- 金幣: 5863.2
- 紅花: 6
- 帖子: 278
- 在線: 352.4小時
- 蟲號: 211914
- 注冊: 2006-03-08
- 性別: GG
- 專業(yè): 無機非金屬材料
6樓2011-10-31 09:27:44
fanny0007
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 972.6
- 散金: 12
- 帖子: 220
- 在線: 120.5小時
- 蟲號: 533141
- 注冊: 2008-03-26
- 性別: MM
- 專業(yè): 電路與系統(tǒng)
7樓2011-10-31 10:24:55
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600，材料與化工321分求調(diào)劑 +10 | 大饞小子 2026-03-28 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
[考研] 材料與化工328分調(diào)劑 +7 | 。，。，。，。i 2026-03-23 | 7/350 |
|
|
|
[考研] 283求調(diào)劑（080500） +7 | A child 2026-03-27 | 7/350 |
|
|
|
[考研] 329求調(diào)劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
|
[碩博家園] 招收生物學(xué)/細胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +7 | 爭取九點睡 2026-03-28 | 8/400 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
|
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|