"対象の挙動をコンピュータ上に表現したモデルを用いて複数の要素を有する物理量である分布量に関するシミュレーションの結果得られる数値の品質を、その目的とする実機体系を模擬した実験の結果を用いて自動的に推定する品質評価装置において、前記実験の結果得られる分布量の測定値QEを記憶する実験体系分布量測定値記憶部と、前記実験体系について前記モデルを用いたシミュレーションの結果得られた前記分布量の計算値QCを記憶する実験体系分布量演算値記憶部と、前記実機体系について前記モデルを用いたシミュレーションの結果得られた前記分布量の計算値QCを記憶する実機体系分布量演算値記憶部と、前記実験についての前記モデルを用いてシミュレーションした結果得られた前記分布量の計算値QCの、前記分布量の当該実験における測定値QEに対する前記分布量相対差diffを、前記各要素間の差に基づいて計算し記憶する分布量相対差演算記憶部と、実験の体系において前記分布量に関して前記モデルを用いてシミュレーションした結果得られた値QEの、前記モデル中の複数のパラメータの各入力値単位変化に対する変化量を表す感度係数ベクトルSEを計算して記憶する実験体系感度係数ベクトル演算記憶部と、目的とする実機体系において前記分布量に関して前記モデルを用いてシミュレーションした結果得られた値QCの、モデル中の複数のパラメータの各入力値単位変化に対する変化量を表す感度係数ベクトルSRを計算して記憶する実機体系感度係数ベクトル演算記憶部と、シミュレーションに用いる前記モデルへの前記パラメータの各入力値の不確かさの割合の関係を示す共分散誤差行列Wを演算し記憶する共分散誤差行列演算記憶部と、RF=SETWSR/{SETWSE1/2SRTWSR1/2}として実験体系と目的とする実機体系との模擬性を模擬性評価因子として定義する該RF値により評価する模擬性評価因子演算記憶部と、前記分布量相対差演算記憶部からの出力である分布量相対差が第1の規定値以下か否かを判定しその規定値を越えれば対象の挙動をコンピュータ上に表現したモデルを用いて複数の要素を有する物理量である分布量に関するシミュレーションの手法を変更する分布量相対差判定部と、前記模擬性評価因子演算記憶部からの出力である模擬性評価因子と1との差が第2の規定値以下か否かを判定しその規定値を超えれば目的とする実機体系を模擬した実験の結果を変更すべき旨を出力する模擬性評価因子判定部と、を有することを特徴とする品質評価装置。前記分布量相対差diffの計算は、計算で得られた分布量をQC=q1c、q2c、qncとし、測定で得られた分布量をQE=q1E、q2E、qmEとし、まず比較に際してQCとQEの成分数であるnとmとが一致していることを確認した後、QCとQEの大きさを以下の式で一致させ、その後に、分布量相対差diffをにより計算することを特徴とする請求項1に記載の品質評価装置。前記分布量相対差diffの計算は、計算で得られた分布量をQC=q1c、q2c、qncとし、測定で得られた分布量をQE=q1E、q2E、qnEとし、まず比較に際してQCとQEの成分数であるnとmとが一致していることを確認した後、QCとQEの大きさを以下の式で一致させ、その後に、分布量相対差diffをdiff=QC-QE/QE=QC-QEQE/QCQE=QCQE/QEQE-1ここでは内積を示すにより計算することを特徴とする請求項1に記載の品質評価装置。前記分布量相対差diffの計算は、計算で得られた分布量をQC=q1c、q2c、qncとし、測定で得られた分布量をQE=q1E、q2E、qnEとし、まず比較に際してQCとQEの成分数であるnとmとが一致していることを確認した後、QCとQEの大きさを以下の式で一致させ、その後に、分布量相対差diffをで計算することを特徴とする請求項1に記載の品質評価装置。前記分布量相対差diffの計算は、計算で得られた分布量をQC=q1c、q2c、qncとし、測定で得られた分布量をQE=q1E、q2E、qnEとし、まず比較に際してQCとQEの成分数であるnとmとが一致していることを確認した後、QCとQEの大きさを以下の式で一致させ、その後に、分布量相対差diffをdiff=QC-QE/QE=QC-QEQE/QCQE=QCQE/QEQE-1ここでは内積を示すにより計算することを特徴とする請求項1に記載の品質評価装置。前記感度係数ベクトルSEおよびSRの計算は、計算で得られた分布量をQC=q1c、q2c、qncとし、測定で得られた分布量をQE=q1E、q2E、qnEとし、まず比較に際してQCとQEの成分数であるnとmとが一致していることを確認した後、QCとQEの大きさを以下の式で一致していることを確認した後、計算上の入力パラメータで着目するものをpiとして、感度係数ベクトルSの次元Lはパラメータの種類の数として、Sのi番目の成分sii=1、2、...、Lをで定義する際、k番目の成分qkk=1,2,...,mそれぞれについて特定のpii=1,2,...,Lに関しての微分を考え、それを成分とするm次元の小ベクトルsii=1、2、...、Lを考え、siを次元mのベクトルとしてsiの第k番目の成分ssikk=1,2,...,mをと定義し、実験体系の感度係数ベクトルをSE、目的とする体系実機体系の感度係数ベクトルをSRとした場合、SE=ss11、、ss1mE、ssi1、、ssimESR=ss11、、ss1mR、ssi1、、ssimRとして、更にSEとSRの成分SuEとSvRどうしの乗算が生じた場合は、SuE×SvR=ssu1、ssu2、...、ssumEssv1、ssv2、...、ssvmRによりベクトルの内積の演算方式で計算すること、を特徴とする請求項1に記載の品質評価装置。前記感度係数ベクトルSEおよびSRの計算は、計算で得られた分布量をQC=q1c、q2c、qncとし、測定で得られた分布量をQE=q1E、q2E、qnEとし、まず比較に際してQCとQEの成分数であるnとmとが一致していることを確認した後、QCとQEの大きさを以下の式で一致していることを確認した後、計算上の入力パラメータで着目するものをpiとして、感度係数ベクトルSの次元Lはパラメータの種類数として、Sのi番目の成分sii=1、2、...、Lをで定義する際、の式の分子部分をQrefはパラメータ変化前の基準系での物理量ベクトルでQPertはパラメータ変化後の摂動系での物理量ベクトルとして、基準系におけるQrefも摂動系におけるQPertもベクトルの大きさを規格化した後、下記の式ここでQrefQrefはQrefどうしの内積、|Qref|はベクトルQrefの大きさ、cosθはベクトル量Qrefとベクトル量QPertが形成する角度で計算すること、を特徴とする請求項1に記載の品質評価装置。対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果得られる分布量の数値の品質を、その対象を模擬した実験の結果を用いて推定する品質評価装置において、実験についての前記モデルを用いたシミュレーションで得られたある分布量Rの計算値の当該実験で測定された当該分布量Rの測定値に対する分布量相対差diff値を計算し記憶する相対差記憶部と、実験の体系について前記モデルを用いてシミュレーションした結果のそのモデルへの入力値の単位変化に対する変化量を表す感度係数ベクトルSEを計算して記憶する実験体系感度係数ベクトル記憶部と、目的とする実機体系について前記モデルを用いてシミュレーションした結果のそのモデルへの入力値の単位変化に対する変化量を表す感度係数ベクトルSRを計算して記憶する実機体系感度係数ベクトル記憶部と、シミュレーションに用いる前記モデルへの入力値の不確かさの割合の関係を示す共分散誤差行列として単位行列Eを記憶する共分散誤差行列演算記憶部と、RF=SETESR/{SETESE1/2SRTESR1/2}として目的とする体系実機体系と実験体系についての模擬性RFを評価する模擬性評価演算部を有することを特徴とする請求項1記載の品質評価装置。対象の挙動をコンピュータ上に表現したモデルを用いて複数の要素を有する物理量である分布量Qに関するシミュレーションの結果得られる数値の品質を、その目的とする実機体系を模擬した実験の結果を用いて自動的に推定する品質評価方法であって、実験体系分布量測定値記憶部でコンピュータがする処理が、前記実験の結果得られる分布量の測定値QEを記憶するステップであり、実験体系分布量演算値記憶部でコンピュータがする処理が、前記実験体系についての前記モデルを用いたシミュレーションの結果得られた前記分布量の計算値QCを記憶するステップであり、実機体系分布量演算値記憶部においてコンピュータがする処理が、前記実機体系についての前記モデルを用いたシミュレーションの結果得られた前記分布量の計算値QCを記憶するステップであり、分布量相対差演算記憶部においてコンピュータがする処理が、前記実験についての前記モデルを用いてシミュレーションした結果得られた前記分布量Qの計算値QCの、前記分布量の当該実験における測定値QEに対する分布量相対差diffを、において前記各要素間の差に基づいて計算し記憶するステップであり、分布量相対差判定部においてコンピュータがする処理が、前記分布量相対差演算記憶部からの出力である分布量相対差が第1の規定値以下か否かを判定し、その規定値を越えればシミュレーション手法を変更するべき旨を出力し計算モデルや手法を変更するステップであり、実験体系感度係数ベクトル演算記憶部においてコンピュータがする処理が、実験の体系において前記分布量Qに関して前記モデルを用いてシミュレーションした結果得られた値QEの、モデル中の複数のパラメータの各入力値単位変化に対する変化量を表す感度係数ベクトルSEを計算して記憶するステップであり、実機体系感度係数ベクトル演算記憶部においてコンピュータがする処理が、目的とする実機体系において前記分布量Qに関して前記モデルを用いてシミュレーションした結果得られた値QCの、モデル中の複数のパラメータの各入力値単位変化に対する変化量を表す感度係数ベクトルSRを計算して記憶するステップであり、共分散誤差行列演算記憶部においてコンピュータがする処理が、模擬性評価因子演算記憶部においてコンピュータがする処理が、シミュレーションに用いる前記モデルへの前記パラメータの各入力値の不確かさの割合の関係を示す共分散誤差行列Wを演算し記憶するステップであり、模擬性評価因子演算記憶部においてコンピュータがする処理が、RF=SETWSR/{SETWSE1/2SRTWSR1/2}として実験体系と目的とする実機体系との模擬性を、模擬性評価因子として定義する該RF値により評価するステップであり、模擬性評価因子判定部においてコンピュータがする処理が、前記模擬性評価因子演算記憶部からの出力である模擬性評価因子と1との差が第2の規定値以下か否かを、判定しその規定値以下でなければ別の模擬実験体系の選定を行うべき旨を出力するステップである、ことを特徴とする品質評価方法。臨界実験の測定値と計算値は、臨界量臨界固有値及び燃料棒や特定の空間位置での核物質の核分裂率分布出力分布に関して多くの比較が実施されてきた。一方ここまでに述べた内容について、少なくとも補正因子バイアス法や断面積アジャストメントに利用されてきた物理量は臨界量臨界固有値などスカラー量である。分布量ベクトル量である核分裂率分布は、単に測定値と計算値をお互い規格化して比較可能な状態にしてそれぞれの空間点で空間点のその差を把握することが行われるだけで終わってしまい、例えば補正因子バイアス法や断面積アジャストメントを始めとする計算コードシステムへのフィードバックには利用されてこなかった。その理由の一つに計算で得られた分布量と実験で得られた分布量とを比較する際に、ベクトル量である分布量どうしの差、あるいは誤差について定義されていないことがある。分布量どうしの一致の程度が定量的に評価できなければ計算手法、計算コード、計算に用いるパラメータの品質について定量的な議論は行えない。加えて断面積アジャストメントに分布量が適用されてこなかった理由の一つとしてベクトル量とベクトル量の差はベクトルとなるので、特定の核種の、特定の反応種類の、特定のエネルギーの断面積についてどのように変化させれば差のベクトルが小さくできるかという感度係数が定義されなかったことがある。この事情は原子力分野に限るものではなくて、科学技術分野で計算に用いるパラメータの品質を高めることを目的に測定値と計算値とを比較して、数理的に計算に用いるパラメータを調整するフィードバックを実施する手法が確立しているのはスカラー量での比較の際だけである。核分裂率分布出力分布等の分布量に対して定量的、数学的に詳しい検討がなされなかったため、実験を行う際に、目的とする体系の核分裂率分布すなわち出力分布を明らかにするには実験においてどの項目を重要視すべきかの判断ができないので、使用できる実験資源を有効に利用した最適な実験を行ったかどうか判断できないという問題がある。さらに、核分裂率測定による出力分に関して複数の実験を行った際、出力分布のような単一のスカラー量ではない物理量に注目した結果に関して、どの実験が実際に目的とする実機体系の計算値の品質確認のために優れていたかを、定量的に比較する方法が提供されていない。その他、過去に実施した実験の情報を有効利用する明確な手法、実験に必要な資源人的、時間的、経済的資源の最適化を行う有効な手法がないのが現状である。そこで、本発明は、対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果得られる分布系の物理量の信頼性を、その対象を模擬した実験の結果を用いて定量的に評価し判定できるようにすることを目的とする。上述の目的を達成するため、本発明は、対象の挙動をコンピュータ上に表現したモデルを用いて複数の要素を有する物理量である分布量に関するシミュレーションの結果得られる数値の品質を、その目的とする実機体系を模擬した実験の結果を用いて自動的に推定する品質評価装置において、前記実験の結果得られる分布量の測定値QEを記憶する実験体系分布量測定値記憶部と、前記実験体系について前記モデルを用いたシミュレーションの結果得られた前記分布量の計算値QCを記憶する実験体系分布量演算値記憶部と、前記実機体系について前記モデルを用いたシミュレーションの結果得られた前記分布量の計算値QCを記憶する実機体系分布量演算値記憶部と、前記実験についての前記モデルを用いてシミュレーションした結果得られた前記分布量の計算値QCの、前記分布量の当該実験における測定値QEに対する前記分布量相対差diffを、前記各要素間の差に基づいて計算し記憶する分布量相対差演算記憶部と、実験の体系において前記分布量に関して前記モデルを用いてシミュレーションした結果得られた値QEの、前記モデル中の複数のパラメータの各入力値単位変化に対する変化量を表す感度係数ベクトルSEを計算して記憶する実験体系感度係数ベクトル演算記憶部と、目的とする実機体系において前記分布量に関して前記モデルを用いてシミュレーションした結果得られた値QCの、モデル中の複数のパラメータの各入力値単位変化に対する変化量を表す感度係数ベクトルSRを計算して記憶する実機体系感度係数ベクトル演算記憶部と、シミュレーションに用いる前記モデルへの前記パラメータの各入力値の不確かさの割合の関係を示す共分散誤差行列Wを演算し記憶する共分散誤差行列演算記憶部と、RF=SETWSR/{SETWSE1/2SRTWSR1/2}として実験体系と目的とする実機体系との模擬性を模擬性評価因子として定義する該RF値により評価する模擬性評価因子演算記憶部と、前記分布量相対差演算記憶部からの出力である分布量相対差が第1の規定値以下か否かを判定しその規定値を越えれば対象の挙動をコンピュータ上に表現したモデルを用いて複数の要素を有する物理量である分布量に関するシミュレーションの手法を変更する分布量相対差判定部と、前記模擬性評価因子演算記憶部からの出力である模擬性評価因子と1との差が第2の規定値以下か否かを判定しその規定値を超えれば目的とする実機体系を模擬した実験の結果を変更すべき旨を出力する模擬性評価因子判定部と、を有することを特徴とする。また、本発明に係る品質評価方法は、対象の挙動をコンピュータ上に表現したモデルを用いて複数の要素を有する物理量である分布量Qに関するシミュレーションの結果得られる数値の品質を、その目的とする実機体系を模擬した実験の結果を用いて自動的に推定する品質評価方法であって、実験体系分布量測定値記憶部でコンピュータがする処理が、前記実験の結果得られる分布量の測定値QEを記憶するステップであり、実験体系分布量演算値記憶部でコンピュータがする処理が、前記実験体系についての前記モデルを用いたシミュレーションの結果得られた前記分布量の計算値QCを記憶するステップであり、実機体系分布量演算値記憶部においてコンピュータがする処理が、前記実機体系についての前記モデルを用いたシミュレーションの結果得られた前記分布量の計算値QCを記憶するステップであり、分布量相対差演算記憶部においてコンピュータがする処理が、前記実験についての前記モデルを用いてシミュレーションした結果得られた前記分布量Qの計算値QCの、前記分布量の当該実験における測定値QEに対する分布量相対差diffを、において前記各要素間の差に基づいて計算し記憶するステップであり、分布量相対差判定部においてコンピュータがする処理が、前記分布量相対差演算記憶部からの出力である分布量相対差が第1の規定値以下か否かを判定し、その規定値を越えればシミュレーション手法を変更するべき旨を出力し計算モデルや手法を変更するステップであり、実験体系感度係数ベクトル演算記憶部においてコンピュータがする処理が、実験の体系において前記分布量Qに関して前記モデルを用いてシミュレーションした結果得られた値QEの、モデル中の複数のパラメータの各入力値単位変化に対する変化量を表す感度係数ベクトルSEを計算して記憶するステップであり、実機体系感度係数ベクトル演算記憶部においてコンピュータがする処理が、目的とする実機体系において前記分布量Qに関して前記モデルを用いてシミュレーションした結果得られた値QCの、モデル中の複数のパラメータの各入力値単位変化に対する変化量を表す感度係数ベクトルSRを計算して記憶するステップであり、共分散誤差行列演算記憶部においてコンピュータがする処理が、模擬性評価因子演算記憶部においてコンピュータがする処理が、シミュレーションに用いる前記モデルへの前記パラメータの各入力値の不確かさの割合の関係を示す共分散誤差行列Wを演算し記憶するステップであり、模擬性評価因子演算記憶部においてコンピュータがする処理が、RF=SETWSR/{SETWSE1/2SRTWSR1/2}として実験体系と目的とする実機体系との模擬性を、模擬性評価因子として定義する該RF値により評価するステップであり、模擬性評価因子判定部においてコンピュータがする処理が、前記模擬性評価因子演算記憶部からの出力である模擬性評価因子と1との差が第2の規定値以下か否かを、判定しその規定値以下でなければ別の模擬実験体系の選定を行うべき旨を出力するステップである、ことを特徴とする。本発明によれば、対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果得られる分布系の物理量の品質信頼性を、その対象を模擬した実験の結果を用いて定量的に評価判定できる。本発明の実施形態は、原子力分野等で使用される計算コードとそこに含まれる物理量のパラメータ断面積ライブラリなどを用いて得られた、目的とする実機体系の分布量の計算値の品質を定量的に評価確認するための装置、方法である。確認には目的とする体系を模擬した実験体系で分布量として得られた測定値とその計算値を用いる。以下の実施形態は、次に述べる第1の評価と第2の評価を行うことにより、最終的に目的とする実機体系における分布量の計算値品質の定量的評価を可能とするものである。第1の評価は、核分裂率分布出力分布等に代表される分布量ベクトル量について、計算値と実験で得られた測定値との一致の程度を数値化し定量的な比較を可能とするものである。第2の評価は、分布量について計算に用いるパラメータに関しての感度係数を分布系に拡張定義した上でこれを求め、得られた感度係数ベクトルを用いて実験体系で測定した分布量が目的とする実機体系の分布量をどれだけ模擬しているかを模擬性評価因子として数値化して、最終的に目的とする実機体系の分布量の計算値の品質の定量的な評価を可能とするものである。すなわち着目する物理量分布量に関して、実験体系ができる限り目的とする体系を模擬した状態であって、その実験体系について計算した計算値と実験で得られた測定値との一致が良いほど、計算に用いた手段計算コードと数値パラメータの品質が良いとの保証が得られる。本発明の実施形態は、そのような判断を定量化された数値に基づいて行うことを可能とする装置を構築することである。第1の対応として対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果の品質を、その対象を模擬した実験の結果を用いて判定する計算値品質判定部において、核分裂率分布出力分布等に代表される分布量ベクトル量について計算値と実験で得られた測定値との一致の程度を定量化する。いま計算で得られた分布量をQC=q1c、q2c、qncとし、測定で得られた分布量をQE=q1E、q2E、qmEとする。まず比較に際してQCとQEの次元成分の数が一致していること、すなわち、n=mが成立することが必要である。次にQCとQEの大きさを一致させる。この場合、大きさについては分布量の物理的意味によって以下の2通りの評価の何れかで大きさを定義してその大きさを一致させ、規格化する。ここでconstとは一定の数値、例えば1を指す。分布量の物理的意味から上記の1式、2式の方式の何れかを選びQCとQEの大きさを一致させたのち、次に以下の何れかの式を用いたQCとQEの差として分布量相対差diffを定義する。diff=QC-QE/QE=QC-QEQE/QCQE=QCQE/QEQE-143式は、統計学で用いられるRootMeanSquarederror差の二乗平均平方根と呼ばれる式と同じ式であり、4式は、2式を適用していた場合、ベクトル量QCとQEがなす角度をθとすれば、cosθ-1と等しくなる。なお3式も4式もQCとQEが一致するときには、分布量相対差diffはゼロとなる。これで分布量QCとQEの差が定義され、目的とする実機体系における試験の代替実験として実施した実験の測定値と、設計コードによる計算値が、どの程度一致しているかについての第1の評価の指標が得られる。分布量相対差diffができる限りゼロに近く、模擬性評価因子RFができる限り1に近いほど、実機体系について計算で求めた分布量の品質が高いと判断される。次に、第2の評価として、実験体系で測定した分布量が目的とする体系の分布量をどれだけ模擬しているかの指標について説明する。まず、物理量が分布系の場合についての感度係数ベクトルを定義する。続いて、得られた感度係数ベクトルを用いて、実験体系で測定した分布量が目的とする体系の分布量をどれだけ模擬しているかを定量化した指標である模擬性評価因子を定義する。最後に、実験体系の分布量の計算値の品質を確認することで目的とする体系の分布量の計算値の品質を確認できるか、定量的に評価する。模擬性評価因子RFRepresentativityfactorは、感度係数ベクトルSと誤差行列Wを用いて以下の式で定義されるRF=SETWSR/{SETWSE1/2SRTWSR1/2}5ここで実験体系も目的とする実機体系も同じ種類の物理量Qに注目するとして、SEはその物理量に関する実験体系の感度係数ベクトル、SRはその物理量に関する目的とする実機体系の感度係数ベクトルである。また、Wは着目するパラメータに関しての共分散誤差行列であり、実験体系と、目的とする実機体系とに共通するパラメータ入力値の不確かさを表す。共分散誤差行列の成分wijは、実験体系のパラメータのi番目の要素の相対誤差ΔpEi/pEiと、目的とする実機体系のパラメータのj番目の要素の相対誤差ΔpRj/pRjとの積である。なお、パラメータの次元をLとして、i、jとも、1からLの範囲である。一方、感度係数ベクトルSは、その次元をLとして、Sのi番目の成分sii=1、2、、Lは、で表される。ここでpiは着目するパラメータである。Sの次元Lはパラメータの種類の数に等しい。上式は物理量Qがスカラー量であれば一意的に定義されるが、物理量Qが分布量すなわちベクトル量と扱える量であればベクトルの微分は自動的には定義されず、上式の扱いを拡張することが必要である。ここで、実験体系と目的とする実機体系とで、注目する物理量である分布量の次元が同一の場合と異なる場合の2通りの評価方法でベクトル量Qの感度係数としての感度係数ベクトルSを定義する。同一次元の場合の評価方法まず、実験体系と目的とする実機体系とで、注目する物理量である分布量の次元が同一の場合について説明する。数学的に自然な方法で模擬性評価因子の適用範囲をベクトル量に拡張する。着目するパラメータpii=1、2、、Lに関しての微分から導かれる感度係数ベクトルであるベクトルSの次元をLとして、Sのi番目の成分をsjj=1、2、、Lとすることは同じである。ただし成分sjj=1、2、、Lをスカラーとはせず、ベクトルに拡張する。Qはベクトル量成分の数、すなわち次元がmなのでk番目の成分qKk=1、2、、mそれぞれについて特定のpii=1、2、、Lに関しての微分を考え、それを成分とする小ベクトルsjj=1、2、、Lを考える。この場合、この小ベクトルsjの成分の数、すなわち次元はmとなる。すなわちここでsjを次元mのベクトルとして、sjの第k番目の成分ssikk=1、2、、mを、と定義する。形式上の表現ではSQ=s1,,sL=ss11,,ss1m,,ssL1,,ssLm8と表現できる。また、SEやSRは、SE=ssE11,,ssE1m,,ssEL1,,ssELmSR=ssR11,,ssR1m,,ssRL1,,ssRLm9と考える。また、なお7式の評価上、sEuとsRvとの掛け算が生じた場合は、でベクトルの内積と設定する。ただし内積が定義できるためには当然ながらsEuとsRvの次元がどちらもmで等しいことが必要となる。この場合、共分散誤差行列W=wuvu、v=1、2、...、Lを介した感度係数ベクトルの2次形式は、と表現できる。ここで、同一次元の場合の評価方法の適用可能範囲を再び確認する。5式の分子の演算においてまず分母の計算については小ベクトル自身の二段階の積が表れる。数学的に自然な定義として内積を適用した。形式上、1式の分母の計算方法は問題が無いが、5式の分子SETWSRの計算を、これも二段階の内積とするならば、sEuの小ベクトルsuEとSRの小ベクトルとsRvとの掛け算が生じる。この際、成分数、すなわち実験体系の物理量であるベクトル量Qの次元と目的とする体系実機の物理量であるベクトル量Qの次元が一致する必要がある。加えて小ベクトルの成分すなわち感度係数ベクトルが表す物理的意味において互いに掛け合わせることに意味をもつことが必要である。翻って、以下のようなケースについては演算に問題が生じ、適用性は保証できない。ケース1物理量を出力分布としたとき、出力分布について注目している空間点数がある。例えば注目している燃料棒本数あるいは位置がある。実験体系と目的とする体系実機の出力分布の注目点数が異なる場合がある。例えば目的とする体系がBWR燃料集合体で9×9型すなわち81点であり、実験体系がBWR燃料集合体で8×8型すなわち64点のような場合である。また、注目する燃料棒本数が異なる場合、SEの小ベクトルsuEとSRの小ベクトルsvRとの掛け算において次元、すなわち燃料棒本数が一致しない場合がある。このような場合には、内積の計算が成立しない。また、例えばPWR17×17燃料集合体の出力分布を実験体系で13×13の空間で模擬した場合などもある。ケース2実験体系と目的とする体系実機の出力分布の注目点数が同じ場合でも、燃料棒の種類が全く異なる場合は内積の物理的意味が定まらずに疑問点が生ずる。例えば、実験体系と注目する体系のガドリニウム棒の挿入位置が異なる場合、内積ではウラン燃料棒の感度係数ベクトルとガドリニウム棒の感度係数ベクトルの掛け算が生じて妥当性が十分とは言えなくなる。これらの場合の解決方法としては更に以下に述べる異なる次元の場合の評価方法で説明する。異なる次元の場合の評価方法次に、実験体系と目的とする実機体系とで、注目する物理量である分布量の次元が同一の場合と異なる場合、すなわち、同一次元の評価方法では適切でなかった前記ケース1およびケース2においても感度係数ベクトルを求めて模擬性評価因子を計算できる評価方法を示す。以下に、評価方法2の評価手法を説明する。感度係数Sはベクトルであり、この次元をLとして、Sのi番目の成分sii=1、2、、nは、であった。の式の分母はスカラーであり、計算について何の問題もない。そこで分子部分を以下の操作でスカラー化する。QRefはパラメータ変化前、これを基準系と呼ぶが、そのときの物理量であり、ベクトル量である。QPertはパラメータ変化後、これを摂動系と呼ぶが、そのときの物理量であり、ベクトル量である。ここでQRefQRefはQRefどうしの内積を表し、|QRef|はベクトルQRefの大きさ、cosθはベクトル量QRefとベクトル量QPertが形成する角度を表す。なお、式13の最も左に示された値を計算する際は、同式2行目右側に示された式ではなく、同式2行目左側の式で計算するのが数値を扱う面では容易である。またその場合、ベクトル量Qは出力分布であれば、計算の前に基準系QRefも摂動系QPertも出力の値の合計が同じかあるいはベクトルの大きさが同じになるような規格化を行う必要がある。なぜなら規格化を行わなければ上式の最右辺に示された式では式の上で値が定まらないからである。この演算を繰り返す事で、計算に必要な感度係数ベクトルを得ることができる。以上により、異なる次元の場合の評価方法、すなわち、物理量が分布系の場合について拡張して定義した感度係数ベクトルを用いて、拡張定義した模擬性評価因子RFを計算する方法により、実験体系で測定した分布量が、目的とする実機体系の分布量をどれだけ模擬しているかを評価することができる。通常RFが0.9以上の場合、比較した2つの体系の模擬性、類似性が高いとされている。以上に述べた方法に基づけば、計算で得られた分布量のベクトル量と実験で得られた分布量の一致の程度を定量的に評価し、加えて着目するパラメータに関して、実験で得られた分布量のベクトル量の目的とする実機体系の分布量に対する模擬性を模擬性評価因子RFを用いて定量評価することによって実機の分布量を計算で評価する際の妥当性の定量的比較が可能となる。ただし分布量相対差diffが0に近い数値としてどの値で品質の良し悪しを判定するか、RFの値がどの数値よりも小さいときには実験体系を計算値の品質評価には用いないかは当該工学分野での技術者や設計責任者の判断により、評価する対象によって決定されるものであり、評価を実施する前の段階であらかじめ設定するものとする。以下、図面を参照して本発明に係る品質評価装置の、以上述べた評価方法に基づく実施形態について説明する。ここで、同一または類似の部分には、共通の符号を付して、重複説明は省略する。なお、実施形態は例示であり、発明の範囲は実施形態に限定されない。ここで、本発明に係る品質評価装置の実施形態の構成とその処理の流れを説明する。図1は、本発明に係る品質評価装置の実施形態におけるハードの全体構成を示すブロック図である。コンピュータ30では、CPUCentralProcessingUnit31、RAMRandomAccessMemory32、ROMReadOnlyMemory33、入力制御部35、表示制御部37、ハードディスク制御部39などが、バス40を介して接続されている。入力制御部35には、キーボードやマウスなどの入力装置34が接続されている。表示制御部37には、液晶ディスプレイなどの表示装置36が接続されている。ハードディスク制御部39には、ハードディスク38が接続されている。CPU31は、計算誤差を推定する誤差推定処理、およびその他の様々な演算処理や制御処理を行う。RAM32は、CPU31が処理を行う際に一時的にデータを記憶する。ROM33あるいはハードディスク38は、測定値を介して得られる相対誤差、誤差行列、計算して得られた感度係数ベクトルなどをCPU31が処理するために必要なデータやプログラムを長期的に記憶する。演算開始の指示などコンピュータ30への入力は、入力装置34を介して行われる。推定された誤差などの必要な情報は、表示装置36に表示される。図2は、本発明に係る品質評価装置の実施形態の機能的な全体構成を示すブロック図である。品質評価装置10は、実験体系分布量測定値記憶部11、実験体系分布量演算部12、実験体系分布量演算値記憶部13、実機体系分布量演算部14、実機体系分布量演算値記憶部15、計算値測定値分布量相対差演算記憶部16、実験体系感度係数ベクトル演算記憶部17、実機体系感度係数ベクトル演算記憶部18、共分散誤差行列演算記憶部19、模擬性評価因子演算記憶部20、分布量相対差判定部21、模擬性評価因子判定部22、を有している。この品質評価装置10は、コンピュータ30上に構築することができる。図3は、本発明に係る品質評価装置の実施形態における誤差推定方法のフローチャートである。図3に記載したS01~S12は処理ステップを示す。以下、図2の機能的構成を踏まえ、図3にしたがって、処理の流れを説明する。まずシミュレーション計算を開始するに当たり、判断のための第1の規定値、第2の規定値が設定されるS01。次に、模擬実験体系が選定されるS02。また、模擬実験体系での分布量の測定値が、実験体系分布量測定値記憶部11に記憶されるS03。更に模擬更に分布量を計算する手法シミュレーション方法が選定されるS04。分布量計算のためのパラメータ入力の選定と設定が行われるS05。次に、共分散誤差行列演算記憶部19において、共分散誤差行列が導出され、結果が記憶されるS06。続いて、ステップS04において選定された計算手法とステップS05によって入力されたパラメータに基づき、実験体系分布量演算部12において実験体系の分布量のシミュレーション計算が実施されS07、得られた実験体系についての分布量の計算結果は、実験体系分布量演算値記憶部13に格納、記憶されるS08。ステップS03で実験体系分布量測定値記憶部11に記憶された模擬実験体系での分布量の測定値と、ステップS08で実験体系分布量演算値記憶部13に記憶された実験体系についての分布量の計算結果は、計算値測定値分布量相対差演算記憶部16に出力され、分布量相対差diffの値が計算され、結果が格納、記憶されるS09。ステップS09の結果は、分布量相対差判定部21に出力され、前記のdiff値の絶対値が第1の規定値以下であるかどうかの判定がなされるS10。diff値が第1の規定値を越える場合は、シミュレーションに用いた手法の見直しが行われるS04。ステップS10でdiff値の絶対値が第1の規定値以下であると判定された場合には、引き続いて、ステップS04において選定された計算手法とステップS05によって入力されたパラメータに基づき、実機体系分布量演算部14において目的とする実機体系の分布量のシミュレーション計算が実行され、実機体系についての分布量について得られた計算結果は、実機体系分布量演算値記憶部15に格納、記憶されるS11。その後、実験体系感度係数ベクトル演算記憶部17において模擬実験の分布量の感度係数が導出され、記憶されるS12。また、併せて、実機体系感度係数ベクトル演算記憶部18において目的体系での分布量の感度係数が導出され、記憶されるS13。ステップS06で共分散誤差行列演算記憶部19において記憶された共分散誤差行列、ステップS12で実験体系感度係数ベクトル演算記憶部17に記憶された模擬実験の分布量の感度係数、ステップS13で実機体系感度係数ベクトル演算記憶部18に記憶された実機体系での分布量の感度係数に基づき、模擬性評価因子演算記憶部20において模擬性評価因子RFが計算、記憶されるS14。模擬性評価因子判定部22において、ステップS14で計算された模擬性評価因子RFと1.0との差の絶対値が第2の規定値以下であるか否かの判定がなされるS15。第2の規定値を越える場合は、模擬実験体系は目的体系を十分模擬できていないとの判定がなされ、模擬実験体系の見直しが行われるS15、S02。模擬性評価因子RFと1.0との差の絶対値が第2の規定値以下である場合は、すでに実験体系に関する分布量の計算結果は、実験結果による分布量と良好に一致しており、計算手法と選定された入力パラメータは妥当であり、実機体系の計算値の品質は十分であり合格、採用できるとの判断がなされS16、以上の自動プロセスを終了する。なお共分散誤差行列は原理的にはシミュレーションに用いる入力の全てについて定義されるものであるが、その評価が非常に困難である場合は共分散誤差行列を単位行列としてもよい。本実施の形態によれば、対象の挙動をコンピュータ上に表現したモデルを用いてシミュレーションした結果得られる目的とする体系の分布系の物理量の値の品質や信頼性を、その対象を模擬した実験の結果を用いて定量的に評価判定できる。以上、実施形態について説明したが、以上の説明は単なる例示であり、本発明は上述の実施形態に限定されず、様々な形態で実施することができる。たとえば、実験体系における分布量の測定値と、その計算値との分布量相対差の演算にあたり、両者の大きさを合わせる方法は、実験体系の性格により、装置に式1を組み込み使用するか、式2を組み込み使用するか、あるいは、式1、2のいずれも組み込み使用時に選択できるようにするか、の方法、装置とすることが可能である。また、実験体系と目的とする実機体系における注目する物理量たる分布量の次元が同一の場合と異なる場合の装置への反映についても、実験体系と実機体系との対応で、いずれかの場合のみを装置に反映するか、いずれも装置に反映し、選択できるようにするかの装置、方法も可能である。さらに、以上のケースを組み合わせた装置、方法が可能である。"