版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

Gout

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 1351.4
紅花: 3
帖子: 345
在線: 227.4小時
蟲號: 1388887
注冊: 2011-09-03
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

[求助] 如何將取向矩陣中歸一化的{hkl}轉(zhuǎn)化成互質(zhì)化的{HKL}，求算法�。。。�

如何將取向矩陣中歸一化的{hkl}轉(zhuǎn)化成互質(zhì)化的整數(shù){HKL}，求程序算法！�。�！

具體就是將如何將三個小數(shù)化為互質(zhì)化的整數(shù)？？

回復(fù)此樓

仰望著那片藍(lán)天

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖置頂 ( 共有1個 )

zxhao168

捐助貴賓 (正式寫手)

應(yīng)助: 37 (小學(xué)生)
金幣: 51.7
散金: 1584
紅花: 18
帖子: 850
在線: 380.2小時
蟲號: 999345
注冊: 2010-04-17
性別: GG
專業(yè): 金屬材料的微觀結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
zhyq8767: 金幣+3, 感謝應(yīng)助 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]均非負(fù)
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    % 只有l(wèi)1 接近零
      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('互質(zhì)化后的指數(shù)[HKL]=')
disp([H K L])

贊一下(2人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

Gout

6樓2012-04-04 22:02:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

hljhrb_ccl

金蟲 (著名寫手)

應(yīng)助: 11 (小學(xué)生)
金幣: 4055.9
散金: 63
紅花: 5
帖子: 1893
在線: 395.9小時
蟲號: 859408
注冊: 2009-09-29
性別: GG
專業(yè): 金屬材料的微觀結(jié)構(gòu)

★
感謝參與，應(yīng)助指數(shù) +1
毛猴: 金幣+1, 感謝交流！ 2012-04-02 09:34:37
毛猴: 應(yīng)助指數(shù)-1, 感謝交流！ 2012-04-02 09:34:48

具體算法我不知道，使用textools這個軟件有一個功能就是可以進(jìn)行這樣的轉(zhuǎn)化的。我之前對鎂合金的織構(gòu)組分分析的時候，就利用這個軟件！

贊一下

回復(fù)此樓

讓知識轉(zhuǎn)化為能力！

2樓2012-04-02 07:51:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

huangbx3

鐵桿木蟲 (小有名氣)

應(yīng)助: 28 (小學(xué)生)
金幣: 8079
散金: 20
紅花: 5
帖子: 291
在線: 177.5小時
蟲號: 1354428
注冊: 2011-07-24
性別: GG
專業(yè): 金屬材料的合金相、相變及

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
毛猴: 金幣+1, 感謝交流！ 2012-04-02 09:34:56

個人感覺該問題應(yīng)該這樣解決：
先將三個小數(shù)整數(shù)化，即同乘以一個數(shù)后將三個小數(shù)都變?yōu)檎麛?shù)；
然后將整數(shù)化的三個數(shù)同除以公約數(shù)，即可得到三個互質(zhì)的數(shù)。

贊一下

回復(fù)此樓

Lovethepeopleyoumet

3樓2012-04-02 09:34:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Gout

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 1351.4
紅花: 3
帖子: 345
在線: 227.4小時
蟲號: 1388887
注冊: 2011-09-03
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

2樓: Originally posted by hljhrb_ccl at 2012-04-02 07:51:35:
具體算法我不知道，使用textools這個軟件有一個功能就是可以進(jìn)行這樣的轉(zhuǎn)化的。我之前對鎂合金的織構(gòu)組分分析的時候，就利用這個軟件！

嗯我看到楊平老師在他的那本書中介紹這個軟件時，提到textools這個軟件有這個功能�？墒牵椰F(xiàn)在授權(quán)使用的軟件只有channel 5 和labtex，而且我更想添加到一些小程序里面，看來還得好好想想。

贊一下

回復(fù)此樓

仰望著那片藍(lán)天

4樓2012-04-02 10:12:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Gout

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 1351.4
紅花: 3
帖子: 345
在線: 227.4小時
蟲號: 1388887
注冊: 2011-09-03
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

引用回帖:

3樓: Originally posted by huangbx3 at 2012-04-02 09:34:13:
個人感覺該問題應(yīng)該這樣解決：
先將三個小數(shù)整數(shù)化，即同乘以一個數(shù)后將三個小數(shù)都變?yōu)檎麛?shù)；
然后將整數(shù)化的三個數(shù)同除以公約數(shù)，即可得到三個互質(zhì)的數(shù)。

我也曾這樣想過，但是有一個問題。得到的取向矩陣每個元素都是小數(shù)點(diǎn)后面好多位的小數(shù)，所以取整數(shù)再除以公約數(shù)會有很大的困難，而且計(jì)算得可能是得到，例如1000001:2000001:3000001這種答案等等。

贊一下

回復(fù)此樓

仰望著那片藍(lán)天

5樓2012-04-02 10:19:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Gout

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 1351.4
紅花: 3
帖子: 345
在線: 227.4小時
蟲號: 1388887
注冊: 2011-09-03
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和

送鮮花一朵

引用回帖:

6樓: Originally posted by zxhao168 at 2012-04-04 22:02:47:
clear
clc
Indices=input('Please input the normalized indices =');
h=Indices(1,1)/sqrt(dot(Indices,Indices));
k=Indices(1,2)/sqrt(dot(Indices,Indices));
l=Indices(1,3)/sqrt(dot(Indices,Indices) ...

很棒的程序，非常感謝�。�！

贊一下

回復(fù)此樓

仰望著那片藍(lán)天

7樓2012-04-05 09:14:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

夢中人sss

鐵蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 217.5
散金: 50
帖子: 105
在線: 159.4小時
蟲號: 2086649
注冊: 2012-10-25
性別: GG
專業(yè): 金屬材料的合金相、相變及

引用回帖:

7樓: Originally posted by Gout at 2012-04-05 09:14:19
很棒的程序，非常感謝�。�！...

哥們，你好，請問下這個程序如何執(zhí)行嗎，是以矩陣形式輸入一個行向量嗎，如輸入【1，0，0】這種嗎，感覺運(yùn)行會出錯，不是這么輸入的

贊一下

回復(fù)此樓

Hello,ladiesandgentleman

8樓2014-09-05 23:17:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jeffchenhust

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 62.4
散金: 76
帖子: 184
在線: 78.5小時
蟲號: 1666191
注冊: 2012-03-05
專業(yè): 集成電路制造與封裝

Textools有試用版的。http://www.sciencesoftware.com.c ... detail12.asp?id=688
在右側(cè)那邊可以下載試用版

贊一下

回復(fù)此樓

9樓2015-01-25 13:20:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lwp8989

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 2.3
散金: 12
帖子: 6
在線: 5.8小時
蟲號: 1490261
注冊: 2011-11-14
性別: GG
專業(yè): 機(jī)械結(jié)構(gòu)強(qiáng)度學(xué)

引用回帖:

9樓: Originally posted by jeffchenhust at 2015-01-25 13:20:34
Textools有試用版的。http://www.sciencesoftware.com.cn/share/search_soft_detail12.asp?id=688
在右側(cè)那邊可以下載試用版

Textools這個軟件還能找到嗎？我也需要用！感謝！

贊一下

回復(fù)此樓

10樓2018-05-29 10:39:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 Gout 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 274求調(diào)劑 +5	S.H1 2026-03-18	5/250	2026-03-18 21:27 by guosr9609
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 化學(xué)工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 268求調(diào)劑 +7	好運(yùn)連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考博] 26申博 +4	八旬速覽 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[基金申請] 國自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 機(jī)械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz

24小時熱門版塊排行榜

[求助] 如何將取向矩陣中歸一化的{hkl}轉(zhuǎn)化成互質(zhì)化的{HKL}，求算法�。。。�

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

» 本帖已獲得的紅花（最新10朵）

【答案】應(yīng)助回帖

[求助] 如何將取向矩陣中歸一化的{hkl}轉(zhuǎn)化成互質(zhì)化的{HKL}，求算法�。。。�

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: