因为科研需要,需要⾃⼰仿真产⽣⼀些晶粒⽅向,并对其进⾏极图和反极图展⽰,在⼩⽊⾍上到⼀个matlab⼯具包mtex,利⽤此⼯具包可以绘制极图和反极图,但是并没有现成的程序可供参考,于是⾃⼰摸索出可以利⽤欧拉⾓绘制极图和反极图,记录于此供⼤家交流学习。巩献田
网络安全事件管理此处⽤的欧拉⾓是内在旋转(Z-X-Z),Bunge约记符号,欧拉⾓相关知识参考博客
相应的matlab代码如下,⼯具包的其他功能暂时没⽤到,所以不是很清楚。
wuxiaworld提⽰:⼩⽊⾍的欧拉⾓数据只需要截取前三列数据即可。
clc,clear
fname=load('C:\Users\zhwdell\Desktop\');%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')%绘制反极图
绘制的极图和反极图结果如下:
极图
反极图
>弥撒音乐