0

我编写了以下代码来使用 24 个晶体对称运算符计算大型数据集中两点之间的定向障碍。尽管最终结果不正确，但代码似乎工作正常。我有一个巨大的数据集，其中包含每个点的欧拉角。我发现第 1 点和第 2 点之间以及第 2 点和第 3 点之间的方向错误，依此类推。

我希望最终结果包含每两个数据点对应的错误定向角度。以下是您更好地理解我的要求的完整代码。 

LA=[phi1 phi phi2];

function gi=get_gi(pvec)
%pvec is a 1x3 vector of [phi1 phi phi2]
g_11=((cosd(pvec(1)).*cosd(pvec(3)))-(sind(pvec(1)).*sind(pvec(3)).*cosd(pvec(2))));
g_12=((sind(pvec(1)).*cosd(pvec(3)))+(cosd(pvec(1)).*sind(pvec(3)).*cosd(pvec(2))));
g_13= (sind(pvec(3)).*sind(pvec(2)));
g_21 =((-cosd(pvec(1)).*sind(pvec(3)))-(sind(pvec(1)).*cos(pvec(3)).*cos(pvec(2)))); 
g_22 = ((-sin(pvec(1)).*sind(pvec(3)))+(cosd(pvec(1)).*cosd(pvec(3)).*cosd(pvec(2))));
g_23 = (cosd(pvec(3)).*sind(pvec(2)));
g_31 = (sind(pvec(1)).* sind(pvec(2)));
g_32 = -cosd(pvec(1)).* sind(pvec(2));
g_33 = cosd(pvec(2)); 

gi =[g_11 g_12 g_13;g_21 g_22 g_23;g_31 g_32 g_33];


f = [1 1 1 -1 1 -1 -1 -1 1 1 -1 -1];
l= [1 1 1];
for i=1:3:10    
l1= [f(i) 0 0;0 f(i+1) 0;0 0 f(i+2)];
l2= [0 f(i) 0;0 0 f(i+1);f(i+2) 0 0];
l3= [0 0 f(i);f(i+1) 0 0;0 f(i+2) 0];
l4= -[0 0 f(i);0 f(i+1) 0;f(i+2) 0 0];
l5= -[0 f(i) 0;f(i+1) 0 0;0 0 f(i+2)];
l6= -[f(i) 0 0;0 0 f(i+1);0 f(i+2) 0];
l=[l;l1;l2;l3;l4;l5;l6];
end
k=1;
t=1;


for m=1:(length(a)-1)
for i=2:3:71
    for j=2:3:71
        g_r= (get_gi(LA(m,:)*l(i:i+2,1:3))*(inv(get_gi(LA(m+1,:)))*inv(l(j:j+2,1:3))));
        tr(k)= g_r(1,1) + g_r(2,2) +g_r(3,3);
        angle(k) = real(acosd((tr(k)-1)/2));
        k=k+1;
    end
    angle(angle==0)=360;
    del_theta=min(angle)
    del(t)=del_theta;
    t=t+1;
end
4

0 回答 0