JIS K 7020:1998 ガラス強化熱硬化性プラスチック（ＧＲＰ）管及び継手―回帰分析法及びその使用 | ページ 6

                                                                                             25
                                                                    K 7020 : 1998 (ISO 10928 : 1997)
                            附属書A(規定) 数学的手順
A.1 行列の演算
  3.2.3(解法)に示す式(27a),(27b)及び(27c)を用いるとき,行列(マトリックス)による解法を記述する。
これらの式に基づき,行列(A : 1)をつくる。
                                         2
                           n      xi   xi    c       yi
                                   2     3
                            xi   xi    xi    d      xiyi   (A : 1)
                             2     3     4            2
                           xi    xi    xi    e      xi yi
  それぞれの変数に対し3.2.6に示す数値を代入すれば,行列(A : 1)は次の式(A : 2)となる。
                            15   21.671  62.989   c     56.728
                          21.671 62.989  180.623  d     80.932   (A : 2)
                         62.989 180.623 584.233    e    235.175
  数学の公式に従い,式(A : 2)の左辺第一項の3×3の行列の逆行列をつくれば次の (A: 3) のようになる。
                          .0132 7   .0040 8  .00017
                           .0040 8 .0152 5   .0042 7  (A : 3)
                           .00017   .0042 7 .00151
  行列(A : 2)と(A : 3)の両辺を乗じて,次の新しい式(A : 4)を得る。
                         1  0  0   c     56.728    .0132 7   .0040 8  .00017
                         0  1  0   d     80.932     .0040 8 .0152 5   .0042 7  (A : 4)
                         0  0  1   e     235.175    .00017   .0042 7 .00151
  この行列(A : 4)を解けば,結果は式(A : 5)となる。
                         c     .3828 8
                         d      .0026 2 (A : 5)
                         e      .0002 2
  このベクトル行列から3.2.3及び3.2.6に示す変数,c,d,eの解を得る。
A.2 代入による計算
  同一セットの式(3.2.3参照)及び3.2.6に示すデータについて上記の解法に代替する解き方を次に示す。
    備考 次の部分は数値を丸めてあるので異なる計算機で計算すると結果は,わずかながら相違する可
          能性がある。
  式(27a)から,
                        56.728=15.000×c+21.671×d+62.989×e  (A : 6)
  式(27b)から,
                        80.932=21.671×c+62.989×d+180.623×e  (A : 7)
  式(27c)から,
                        235.175=62.989×c+180.623×d+584.233×e  (A : 8)
  式(A : 8)を解いて,式(A : 9)を得る。

――――― [JIS K 7020 pdf 26] ―――――

26
K 7020 : 1998 (ISO 10928 : 1997)
                        e=0.402 5−0.107 8×c−0.309 2×d  (A : 9)
  式(A : 9)を式(A : 6)及び(A : 7)に代入して式(A : 10)及び式(A : 11)を得る。
                        31.375=8 210×c+2.195×d  (A : 10)
                        8.231=2.200×c+7.140×d  (A : 11)
  式(A : 11)を解いて式(A : 12)を得る。
                        d=1.152 8−0.308 1×c  (A : 12)
  式(A : 12)のdを式 (A : 10)に代入すれば,式(A : 13)を得る。
                        7.533 7×c=28.844 6  (A : 13)
  したがって,
                        c=3.828 8  (A : 14)
  式(A : 14)の結果を式(A : 13)のcに代入すれば,式(A : 15)によってdを得る。
                        d=−0.026 2  (A : 15)
  式(A : 14)のc,式(A : 15)のdの値をそれぞれ式(A: 9)に代入すれば,eの値を得る。
                        e=−0.002 2  (A : 16)
  このようにして,両法はいずれも同じ結果を与える(ただし,A.2の備考を参照)。
                               JIS K 7020 原案作成委員会 構成表
                                氏名                               所属
 (委員長)              宮 入 裕 夫         東京医科歯科大学医用器材研究所
                         宗 宮     詮        慶應義塾大学理工学部機械工学科
                         安 宅 信 行         昭和女子大学大学院生活機構研究科
                         西 出 徹 雄         通商産業省基礎産業局化学課
                         (増  田     優)
                         大 嶋 清 治         通商産業省工業技術院標準部材料規格課
                         橋 本 繁 晴         財団法人日本規格協会技術部国際規格整合室
                         青 木     茂        青木技術事務所
                         磯 貝 智 彦         東京電力株式会社原子力技術部
                         (山 本 晋 児)
                         西 本 直 樹         三井化学株式会社生産技術本部エンジニアリング部
                         鳥 山     義        興和化成株式会社化成品本部
                         中 井 邦 彦         富士化工株式会社技術部
                         野間口 兼 政        日立化成工業株式会社山崎工場
                         河 内 秀 二         栗本化成株式会社湖東工場技術部
                         北 村 達 人         社団法人強化プラスチック協会
 (事務局)              田 村 正 勝         日本プラスチック工業連盟
                         備考  ( ) 内は,前任者