24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Gout

金虫 (正式写手)

应助: 2 (幼儿园)
金币: 1351.4
红花: 3
帖子: 345
在线: 227.4小时
虫号: 1388887
注册: 2011-09-03
性别: GG
专业: 凝聚态物性I:结构、力学和

[求助] 如何将取向矩阵中归一化的{hkl}转化成互质化的{HKL}，求算法！！！！

如何将取向矩阵中归一化的{hkl}转化成互质化的整数{HKL}，求程序算法！！！！

具体就是将如何将三个小数化为互质化的整数？？

回复此楼

仰望着那片蓝天

1楼 2012-04-01 23:25:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zxhao168

捐助贵宾 (正式写手)

应助: 37 (小学生)
金币: 51.7
散金: 1584
红花: 18
帖子: 850
在线: 380.2小时
虫号: 999345
注册: 2010-04-17
性别: GG
专业: 金属材料的微观结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhyq8767: 金币+3, 感谢应助 2012-04-04 22:05:10
Gout: 金币+10, ★★★★★最佳答案, 绝对的好人，非常好的程序，谢谢你啊~ 2012-04-05 09:10:41
Gout: 回帖置顶 2012-04-05 09:27:32

clear
clc
Indices=input('Please input the normalized indices [hkl]=');
h=Indices(1,1)/sqrt(dot(Indices,Indices));
k=Indices(1,2)/sqrt(dot(Indices,Indices));
l=Indices(1,3)/sqrt(dot(Indices,Indices));
h1=abs(h);
k1=abs(k);
l1=abs(l);
if h1>1e-3 && k1>1e-3 && l1>1e-3 % [hkl]均非负
if h1>=l1 && k1>=l1    % l1最小
      A=h1/l1;
      B=k1/l1;
      i=1;
      while (abs(A-round(A))>1e-2 || abs(B-round(B))>1e-2)
         A=A*(i+1)/i;
         B=B*(i+1)/i;
         i=i+1;
      end
      if sign(l)<0
         L=-i;
      else
         L=i;
      end
      if sign(h)<0
         H=-round(A);
      else
         H=round(A);
      end
      if sign(k)<0
         K=-round(B);
      else
         K=round(B);
      end
elseif h1>=k1 && l1>=k1 % k1最小
      A=h1/k1;
      B=l1/k1;
      i=1;
      while (abs(A-round(A))>1e-2 || abs(B-round(B))>1e-2)
         A=A*(i+1)/i;
         B=B*(i+1)/i;
         i=i+1;
      end
      if sign(k)<0
         K=-i;
      else
         K=i;
      end
      if sign(h)<0
         H=-round(A);
      else
         H=round(A);
      end
      if sign(l)<0
         L=-round(B);
      else
         L=round(B);
      end
elseif k1>=h1 && l1>=h1    % h1最小
      A=k1/h1;
      B=l1/h1;
      i=1;
      while (abs(A-round(A))>1e-2 || abs(B-round(B))>1e-2)
         A=A*(i+1)/i;
         B=B*(i+1)/i;
         i=i+1;
      end
      if sign(h)<0
         H=-i;
      else
         H=i;
      end
      if sign(k)<0
         K=-round(A);
      else
         K=round(A);
      end
      if sign(l)<0
         L=-round(B);
      else
         L=round(B);
      end
end
elseif h1<=1e-3          % h1 接近零
H=0;
if k1<=1e-3       % h1接近零的同时k1也接近零
      K=0;
      if sign(l)<0
         L=-round(l);
      else
         L=round(l);
      end
elseif l1<=1e-3 % h1接近零的同时，l1也接近零
      L=0;
      if sign(k)<0
         K=-round(k);
      else
         K=round(k);
      end
elseif k1>1e-3 && l1>1e-3    % 只有h1 接近零
      if k1>l1
         A=k1/l1;
         i=1;
         while(abs(A-round(A))>1e-3)
            A=A*(i+1)/i;
            i=i+1;
         end
         if sign(l)<0
            L=-i;
         else
            L=i;
         end
         if sign(k)<0
            K=-round(A);
         else
            K=round(A);
         end
      else
         A=l1/k1;
         i=1;
         while(abs(A-round(A))>1e-3)
            A=A*(i+1)/i;
            i=i+1;
         end
         if sign(k)<0
            K=-i;
         else
            K=i;
         end
         if sign(l)<0
            L=-round(A);
         else
            L=round(A);
         end
      end
end
elseif k1<=1e-3          % k1 接近零
K=0;
if h1<=1e-3       % k1接近零的同时h1也接近零
      H=0;
      if sign(l)<0
         L=-round(l);
      else
         L=round(l);
      end
elseif l1<=1e-3 % k1接近零的同时，l1也接近零
      L=0;
      if sign(h)<0
         H=-round(h);
      else
         H=round(h);
      end
elseif h1>1e-3 && l1>1e-3    % 只有k1 接近零
      if h1>l1
         A=h1/l1;
         i=1;
         while(abs(A-round(A))>1e-3)
            A=A*(i+1)/i;
            i=i+1;
         end
         if sign(l)<0
            L=-i;
         else
            L=i;
         end
         if sign(h)<0
            H=-round(A);
         else
            H=round(A);
         end
      else
         A=l1/k1;
         i=1;
         while(abs(A-round(A))>1e-3)
            A=A*(i+1)/i;
            i=i+1;
         end
         if sign(h)<0
            H=-i;
         else
            H=i;
         end
         if sign(l)<0
            L=-round(A);
         else
            L=round(A);
         end
      end
end
elseif l1<=1e-3          % l1 接近零
L=0;
if h1<=1e-3       % l1接近零的同时h1也接近零
      H=0;
      if sign(k)<0
         K=-round(k);
      else
         K=round(k);
      end
elseif k1<=1e-3 % l1接近零的同时，k1也接近零
      K=0;
      if sign(h)<0
         H=-round(h);
      else
         H=round(h);
      end
elseif h1>1e-3 && k1>1e-3    % 只有l1 接近零
      if h1>k1
         A=h1/k1;
         i=1;
         while(abs(A-round(A))>1e-3)
            A=A*(i+1)/i;
            i=i+1;
         end
         if sign(k)<0
            K=-i;
         else
            K=i;
         end
         if sign(h)<0
            H=-round(A);
         else
            H=round(A);
         end
      else
         A=k1/h1;
         i=1;
         while(abs(A-round(A))>1e-3)
            A=A*(i+1)/i;
            i=i+1;
         end
         if sign(h)<0
            H=-i;
         else
            H=i;
         end
         if sign(k)<0
            K=-round(A);
         else
            K=round(A);
         end
      end
end
end
disp('互质化后的指数[HKL]=')
disp([H K L])