CASTEP获得电荷密度等高线的Matlab作图法试题.docVIP

资料来自小木虫论坛；H.J.Zhao汇总 CASTEP获得电荷密度等高线的Matlab作图法 MATLAB程序用来处理MS电荷密度的等高线做法，分四步：1.找到电荷密度图的存储文件，一般是在自己所建的计算project中，隐藏的，castep模块是.charg_frm格式的，可以用txt打开后另存为txt格式。（这个地方注意把文件单独拷贝出来，以免破坏原始文件）2.读取数据，在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数据的，其中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(