matlab 等高线密度,CASTEP获得电荷密度等高线的Matlab作图法

资料来自小木虫论坛；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(:,layer,:),[18,144]),reshape(c(:,layer,:),[18,144]),reshape( d(:,layer,:),[18,144]),8);

layer=1(2，3，4....)可以修改作图的层数，分别是从前面到后面。最后的8表示显示线条数目，可以自己修改成9，10或其它。

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

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

>> layer=9;

contour(reshape(a(:,layer,:),[18,144]),reshape(c(:,layer,:),[18,144]),reshape(

matlab 等高线密度,CASTEP获得电荷密度等高线的Matlab作图法相关推荐

差分电荷密度 matlab,CASTEP 差分电荷密度操作流程
原标题:CASTEP 差分电荷密度操作流程本示例以 CO 吸附在 Pt(111) 表面上为例 1.对吸附的结构进行结构优化,如图1所示: 图1 「优化后的结构」 2.如果要研究吸附分子 CO 和表面 ...
CAD/CASS依据等高线等距或随机加高程点插件（等高线等距加高程点、等高线随机加高程点、等高线手动加高程点）
插件下载: QTools for AutoCADhttp://qtools.ysepan.com/ 使用方法: APPLOAD命令加载插件,dgx命令打开等高线工具(包含以下功能:修等高线.等高线升降 ...
差分电荷密度怎么画_MS电荷密度在Origin中作图方法
MS 电荷密度在 Origin 中作图方法很早就想把这个给大家分享,因为看到一些虫子还在为此发愁.我也用过 matlab ,但结果不甚理想,所以就想如何在 origin 中快速实现. 1 . Tot ...
matlab模拟薄膜振动,（数理方程）Matlab模拟琴弦振动发声并显示振动图像
1.物理模型建立由课题问题可转化为实际数理问题,即典型弦振动问题,一根均匀弦两端分别在x=0以及x=L处固定,设初始速度为零,初始时刻弦的形状为一抛物线,抛物线的顶点为(L/2,h),讨论其弦振动的 ...
MATLAB计算杨氏模量,四阶弹性模量Cijkl如何在matlab里表示啊？ - 计算模拟 - 小木虫 - 学术科研互动社区...
matlab 四元数运算计算包就可以了吧 Matlab 四元数操作函数 2012-06-03 21:02:55| 分类: MATLAB&Mathemati | 标签:四元数 quater ...
matlab系统的根轨迹,实验五利用MATLAB绘制系统根轨迹
<实验五利用MATLAB绘制系统根轨迹>由会员分享,可在线阅读,更多相关<实验五利用MATLAB绘制系统根轨迹(6页珍藏版)>请在人人文库网上搜索. 1.实验五利用MAT ...
在matlab中有几种获得帮助的途径,matlab经典习题及解答
精品文档 . 1欢迎下载第1章 MATLAB 概论 1.1 与其他计算机语言相比较,MATLAB 语言突出的特点是什么? MATLAB 具有功能强大.使用方便.输入简捷.库函数丰富.开放性强等特点. ...
【数字信号处理】卷积编程实现 ( Matlab 卷积和多项式乘法 conv 函数 | 使用 matlab 代码求卷积并绘图 )
文章目录一.Matlab 卷积和多项式乘法 conv 函数二.使用 matlab 代码求卷积并绘图一.Matlab 卷积和多项式乘法 conv 函数 Matlab 文档地址 : https:// ...
matlab利用作图法求圆周率
MATLAB本身,pi返回值最接近基本常数π的浮点数.要注意,MATLAB中定义的pi与基本常数π之间的微小差异会导致计算结果出现微小的误差. 例如,表达式sin(pi)(期望返回0)实际上在MATL ...
让书写的Matlab代码运行更快 Recipes for Faster Matlab Code
Matlab 在 Research 中用得非常多,确实也是非常方便实用,只是有一个问题就是写 Matlab 代码的时候经常需要用一些比较奇怪独特的方式来思考和处理问题,否则写出来的代码虽然同样能工作, ...

matlab 等高线密度,CASTEP获得电荷密度等高线的Matlab作图法

matlab 等高线密度,CASTEP获得电荷密度等高线的Matlab作图法相关推荐

最新文章

热门文章