版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

雁兒霏霏

木蟲(chóng) (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 2645.4
帖子: 883
在線: 160.1小時(shí)
蟲(chóng)號(hào): 934537

[交流] 【求助】緊急求教讀取CHGCAR的matlab腳本，謝謝

現(xiàn)在需要深入研究vasp的CHGCAR，希望用matlab讀取它，但對(duì)矩陣的讀取一直沒(méi)查到相關(guān)資料，求教哪位有讀能提供取CHGCAR的matlab小腳本，謝謝。

回復(fù)此樓

» 猜你喜歡

一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有6人回復(fù)
081700化工學(xué)碩調(diào)劑已經(jīng)有3人回復(fù)
材料學(xué)碩318求調(diào)劑已經(jīng)有4人回復(fù)
329求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有11人回復(fù)
一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一）已經(jīng)有3人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有14人回復(fù)
320求調(diào)劑0856 已經(jīng)有3人回復(fù)
307求調(diào)劑已經(jīng)有9人回復(fù)
288求調(diào)劑已經(jīng)有15人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

繼續(xù)求助Matlab繪圖！已經(jīng)有8人回復(fù)
弱弱問(wèn)一下：stereochemical activity 是什么意思？謝謝已經(jīng)有5人回復(fù)
緊急求助翻譯一篇德文文獻(xiàn)~翻譯其中合成的一部分~急求~謝謝~ 已經(jīng)有4人回復(fù)
matlab 計(jì)算過(guò)程中工作空間變量保存及讀取求解決內(nèi)存不足的方法已經(jīng)有5人回復(fù)
幫翻譯幾段文獻(xiàn)，請(qǐng)不要在線翻譯，謝謝已經(jīng)有2人回復(fù)
用matlab統(tǒng)計(jì)圖片中某種顏色成分占總圖片中所有顏色成分的比率已經(jīng)有4人回復(fù)
【求助】弱弱的問(wèn)一下Matlab數(shù)據(jù)提取的問(wèn)題已經(jīng)有5人回復(fù)
【求助】用matlab編程，需要將數(shù)據(jù)自檢分類（高手進(jìn)來(lái)指點(diǎn)下）已經(jīng)有5人回復(fù)
【求助】CHGCAR用什么來(lái)讀取已經(jīng)有10人回復(fù)
【求助】matlab如何只讀取圖像的一部分？已經(jīng)有11人回復(fù)
【求助】有哪位大俠做過(guò)次氯酸叔丁酯��？幫幫我啊！謝謝已經(jīng)有6人回復(fù)
【求助】matlab讀取fortran輸出的txt文件【已解決】已經(jīng)有3人回復(fù)
【求助】求教matlab解非線性方程組已經(jīng)有9人回復(fù)
【求助】Weickert的各向異性擴(kuò)散方程的濾波方法的matlab程序代碼已經(jīng)有8人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-01-15 10:33:48

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

09njtpw

鐵桿木蟲(chóng) (職業(yè)作家)

應(yīng)助: 1 (幼兒園)
金幣: 33684.9
帖子: 3231
在線: 329.2小時(shí)
蟲(chóng)號(hào): 860159

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖

引用回帖:

1樓: Originally posted by 雁兒霏霏 at 2011-01-15 10:33:48:
現(xiàn)在需要深入研究vasp的CHGCAR，希望用matlab讀取它，但對(duì)矩陣的讀取一直沒(méi)查到相關(guān)資料，求教哪位有讀能提供取CHGCAR的matlab小腳本，謝謝。

不懂樓主問(wèn)題解決了沒(méi)有？？方便傳一份腳本給我么？貼的那個(gè)有點(diǎn)亂

贊一下(1人)

回復(fù)此樓

3樓2011-12-30 22:59:32

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 3 個(gè)回答

liangab_234620

木蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1341.2
帖子: 131
在線: 145.5小時(shí)
蟲(chóng)號(hào): 461188

★ ★
sunyang1988(金幣+2): 謝謝幫助。這樣直接粘貼可能有格式錯(cuò)誤，如果方便請(qǐng)上傳到網(wǎng)盤 2011-01-16 14:10:22
雁兒霏霏(金幣+20): 謝謝大蝦，能給提供個(gè)程序文件嗎？粘貼有一些格式錯(cuò)誤。郵箱：shangyan2009@gmail.com ，謝謝！ 2011-01-16 19:37:24

% Read in CHG file to 3D matrix.

fid = fopen('CHG','r');

% Read POSCAR part of CHG file:
sysname = fgetl(fid);
lconst = fscanf(fid,'%f',1);
basisvec = fscanf(fid,'%f',[3,3])';
temp1 = fgetl(fid);
temp2 = fgetl(fid);
species = sscanf(temp2,'%i');
Natoms = sum(species);
temp3 = fgetl(fid);
atompos_f3 = fscanf(fid,'%f',[3,Natoms])';

volume = abs(dot(basisvec(1,

,cross(basisvec(2,

,basisvec(3,

)))*lconst^3;

% Read charge density data. (i,j,k) value in i+(j+k*Ny)*Nx, with coordinate
% system defined by basis vectors!:
gridsize = fscanf(fid,'%i',3);
Nx = gridsize(1);
Ny = gridsize(2);
Nz = gridsize(3);

temp4 = fscanf(fid,'%f',[Nx*Ny*Nz,1]);

density_sub_h1_surf=reshape(temp4,Nx,Ny,Nz)/volume;

% % Magnitization (if present):
% temp5 = fscanf(fid,'%i',3);
% temp6 = fscanf(fid,'%f',[Nx*Ny*Nz,1]);
% magn_lowacc_t2_N�蕆eshape(temp6,Nx,Ny,Nz)/volume;

% minval = min(min(min(density)));
% maxval = max(max(max(density)));

% For rectangular unit cell, unit axes:
% Plot isosurface:
% prect = patch(isosurface(density,7.5), 'FaceColor', 'red', 'EdgeColor', 'none');
% isonormals(density,prect);
% %axis tight;
% camlight; lighting phong;

% Non-rectangular (general) unit cell:
% Need to created position 3D matrices
xtemp = zeros(Nx*Ny*Nz,1);
ytemp = zeros(Nx*Ny*Nz,1);
ztemp = zeros(Nx*Ny*Nz,1);

for k = 0:Nz-1
for j = 0:Ny-1
      for i = 0:Nx-1
         xtemp(1+i+(j+k*Ny)*Nx) = i/Nx*basisvec(1,1)+j/Ny*basisvec(2,1)+k/Nz*basisvec(3,1);
         ytemp(1+i+(j+k*Ny)*Nx) = i/Nx*basisvec(1,2)+j/Ny*basisvec(2,2)+k/Nz*basisvec(3,2);
         ztemp(1+i+(j+k*Ny)*Nx) = i/Nx*basisvec(1,3)+j/Ny*basisvec(2,3)+k/Nz*basisvec(3,3);
      end
end
end

X = reshape(xtemp,Nx,Ny,Nz)*lconst;
Y = reshape(ytemp,Nx,Ny,Nz)*lconst;
Z = reshape(ztemp,Nx,Ny,Nz)*lconst;

z = squeeze(Z(1,1,

);
ave_z_h1 = squeeze(sum(sum(density_sub_h1_surf,1),2))/(Nx*Ny); %(sqrt(2)*Nx*Ny)*lconst^2;

figure
%subplot(211)
plot(z,ave_z_h1,'k')

% figure
% p = patch(isosurface(X,Y,Z,density,7.5), 'FaceColor', 'red', 'EdgeColor', 'none');
% isonormals(density,p);
% % maxx = max(basisvec(:,1))*lconst;
% % minx = min(basisvec(:,1))*lconst;
% % maxy = max(basisvec(:,2))*lconst;
% % miny = min(basisvec(:,2))*lconst;
% % maxz = max(basisvec(:,3))*lconst;
% % minz = min(basisvec(:,3))*lconst;
% % axis([minx maxx miny maxy minz maxz])
% camlight; lighting phong;

% Extract charge density along certain lines:

------------------------------

贊一下(1人)

回復(fù)此樓

2樓2011-01-16 12:51:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 3 個(gè)回答

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 材料學(xué)碩318求調(diào)劑 +4	February_Feb 2026-03-19	4/200	2026-03-19 23:29 by fmesaito
[考研] 321求調(diào)劑 +8	何潤(rùn)采123 2026-03-18	10/500	2026-03-19 16:46 by 何潤(rùn)采123
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4	脫穎而出 2026-03-16	12/600	2026-03-19 16:17 by 脫穎而出
[考研] 0817調(diào)劑 +3	沒(méi)有答案_ 2026-03-14	3/150	2026-03-19 09:51 by Xu de nuo
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +10	秋有木北 2026-03-14	10/500	2026-03-19 05:52 by anny19840123
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無(wú)懈可擊111
[考研] 化學(xué)工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 299求調(diào)劑 +5	△小透明* 2026-03-17	5/250	2026-03-18 11:49 by 盡舜堯1
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3	QGZDSYS 2026-03-13	5/250	2026-03-18 11:13 by QGZDSYS
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 301求調(diào)劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無(wú)際的草原
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無(wú)懈可擊111
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu