[Matlab科学绘图] 由欧拉角绘制极图和反极图

因为科研需要,需要自己仿真产生一些晶粒方向,并对其进行极图和反极图展示,在小木虫上找到一个matlab工具包mtex,利用此工具包可以绘制极图和反极图,但是并没有现成的程序可供参考,于是自己摸索出可以利用欧拉角绘制极图和反极图,记录于此供大家交流学习。

Mtex工具包网址:https://mtex-toolbox.github.io/

欧拉角数据参考小木虫上已有的:http://muchong.com/bbs/viewthread.php?tid=12562139

此处用的欧拉角是内在旋转(Z-X-Z),Bunge约记符号,欧拉角相关知识参考博客欧拉角和坐标变换个人总结

相应的matlab代码如下,工具包的其他功能暂时没用到,所以不是很清楚。

提示:小木虫的欧拉角数据只需要截取前三列数据即可。

clc,clear
fname=load('C:\Users\zhwdell\Desktop\TEX_PH1.txt');%Bunge Euler angle convention
cs = crystalSymmetry('cubic');         %晶体对称性
ss = specimenSymmetry('triclinic'); %试样对称性
ori = orientation.byEuler(fname(:,1)*degree,fname(:,2)*degree,fname(:,3)*degree,cs,ss);
plotPDF(ori,Miller({1,0,0},{1,1,0},{1,1,1},cs),'all') %绘制极图
plotIPDF(ori,[vector3d.X,vector3d.Y,vector3d.Z],'all')%绘制反极图

绘制的极图和反极图结果如下:

极图
反极图

猜你喜欢

转载自blog.csdn.net/zhwzhaowei/article/details/88553892
今日推荐