[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')%绘制反极图

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

极图

 

 

反极图