2024年5月1日发(作者：)

资料来自小木虫论坛；汇总

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(:,layer,:),[18,144]),reshape(c(:,layer,:),[18,144]),reshape(d(:,la

yer,:),[18,144]),8);

layer=1（2，3，4....） 可以修改作图的层数，分别是从前面到后面。最后的8表示

显示线条数目，可以自己修改成9，10或其它。

hold on 表示累积作图，hold off表示清空前面的 作图

举例；做表面层和第九层的等高线图

>> layer=9;