基準状態の温度と日射強度における太陽電池の短絡電流、開放電圧、動作電圧、及び動作電流を示すデータの入力を受け、前記データを用いて前記太陽電池ごとに、前記基準状態における前記太陽電池の特性を決定するものであって前記太陽電池ごとにばらつきをもつパラメータを算出し、前記算出されたパラメータを用いて、前記開放電圧の温度特性と逆方向飽和電流の温度特性を算出し、前記開放電圧の温度特性と前記逆方向飽和電流の温度特性から、所定の温度における前記短絡電流、前記開放電圧、前記動作電圧、及び前記動作電流を算出する太陽電池の特性演算方法。請求項1に記載の太陽電池の特性演算方法において、前記パラメータのばらつき感度を算出し、前記ばらつき感度の高い順にフィッティングを行うことで前記パラメータの値を確定する太陽電池の特性演算方法。請求項1に記載の太陽電池の特性演算方法において、Iを前記太陽電池の出力電流、Isを前記太陽電池の逆方向飽和電流、Vを前記太陽電池の出力電圧、Iscを前記太陽電池の短絡電流、Tを前記太陽電池の温度、kをボルツマン定数、Rsを前記太陽電池の直列抵抗、qを電子の電荷量、Rshを前記太陽電池の並列抵抗、nをダイオード接合定数、pを日射強度、Ncellを前記太陽電池を構成する太陽電池セル数として、前記太陽電池の特性が下記の式によって表される場合、前記ばらつきをもつパラメータは、前記太陽電池の逆方向飽和電流、前記太陽電池の直列抵抗、前記太陽電池の並列抵抗、前記ダイオード接合定数である太陽電池の特性演算方法。IIscpIsexp//Rsh...請求項1に記載の太陽電池の特性演算方法において、前記所定の温度における前記短絡電流、前記開放電圧、前記動作電圧、及び前記動作電流を、前記太陽電池の特性を示す式に適用し、前記太陽電池の直列抵抗のフィッティングを行う太陽電池の特性演算方法。請求項1に記載の太陽電池の特性演算方法において、Isを前記太陽電池の逆方向飽和電流、Vocを前記太陽電池の開放電圧、nを前記太陽電池のダイオード接合定数、Tを前記太陽電池の温度、kをボルツマン定数、qを電子の電荷量、Egoを前記太陽電池のバンドギャップ電圧、Ncellを前記太陽電池を構成する太陽電池セル数とした場合、前記逆方向飽和電流の温度特性は、下記の式によって表され、前記開放電圧の温度特性は、下記の式によって表される太陽電池の特性演算方法。∂Is/∂TIs...∂Voc/∂T/q}Ncell...請求項5に記載の太陽電池の特性演算方法において、Tbを前記所定の温度、Isc[Tb]を前記所定の温度における短絡電流、Is[Tb]を前記所定の温度における逆方向飽和電流、Voc[Tb]を前記所定の温度における開放電圧、IopTbを前記所定の温度における動作電流、jを定数、VopTbを前記所定の温度における動作電圧とした場合、前記式と前記式を用いて、前記所定の温度における逆方向飽和電流と開放電圧を算出し、下記の式を用いて、前記所定の温度における短絡電流を算出し、下記の式を用いて、前記所定の温度における動作電流を算出し、下記の式を用いて、前記所定の温度における動作電圧を算出する太陽電池の特性演算方法。IscTb≒IsTbexpVocTb...IopTb≒jIscTb...VopTb/TTbVocTb...複数の太陽電池モジュールを直列に接続した太陽電池ストリングを複数有し、前記複数の太陽電池ストリングを並列接続することによって構成される太陽電池アレイと、前記太陽電池アレイに接続される負荷変動部と、前記太陽電池アレイの出力電流を検出する電流検出部と、前記太陽電池アレイの出力電圧を検出する電圧検出部と、前記出力電流及び前記出力電圧に基づいて、前記負荷変動部に制御信号を出力して前記太陽電池アレイの動作電圧を制御する制御部と、前記太陽電池アレイの特性を演算する特性演算部と、を有し、前記特性演算部は、基準状態の温度と日射強度における前記複数の太陽電池モジュールそれぞれの短絡電流、開放電圧、動作電圧、及び動作電流を示すデータの入力を受け、前記データを用いて前記太陽電池モジュールごとに、前記基準状態における前記太陽電池モジュールの特性を決定するものであって前記太陽電池モジュールごとにばらつきをもつパラメータを算出し、前記算出されたパラメータを用いて、前記開放電圧の温度特性と逆方向飽和電流の温度特性を算出し、前記開放電圧の温度特性と前記逆方向飽和電流の温度特性から、所定の温度における前記短絡電流、前記開放電圧、前記動作電圧、及び前記動作電流を算出し、前記複数の太陽電池モジュールごとの前記所定の温度における短絡電流、前記開放電圧、前記動作電圧、前記動作電流から、前記太陽電池アレイの特性を演算する太陽光発電システム。請求項7に記載の太陽光発電システムにおいて、前記特性演算部の算出結果から得られる最大電力と、前記出力電流と前記出力電圧から得られる電力とを比較する比較部を有し、前記制御部は、前記比較部の比較結果に基づいて前記制御信号を出力する太陽光発電システム。請求項7に記載の太陽光発電システムにおいて、前記特性演算部は、前記パラメータのばらつき感度を算出し、前記ばらつき感度の高い順にフィッティングを行うことで前記パラメータの値を確定する太陽光発電システム。請求項7に記載の太陽光発電システムにおいて、前記特性演算部は、前記所定の温度における前記短絡電流、前記開放電圧、前記動作電圧、及び前記動作電流を、前記太陽電池モジュールの特性を示す式に適用し、前記太陽電池モジュールの直列抵抗のフィッティングを行う太陽光発電システム。請求項7に記載の太陽光発電システムにおいて、前記開放電圧の温度特性と前記逆方向飽和電流の温度特性と前記出力電流と前記出力電圧から前記太陽電池アレイの温度を算出する環境データ算出部をさらに有する太陽光発電システム。請求項11に記載の太陽光発電システムにおいて、前記環境データ算出部は、前記太陽電池アレイの温度と前記開放電圧の温度特性と前記逆方向飽和電流の温度特性から、日射強度が前記基準状態であって前記太陽電池アレイの温度における短絡電流を算出し、前記出力電流に定数を乗じた値と前記算出された短絡電流との比をとることによって、前記太陽電池アレイにかかる日射強度を算出する太陽光発電システム。請求項12に記載の太陽光発電システムにおいて、前記定数は、基準状態の温度と日射強度における前記太陽電池モジュールの短絡電流と最大電力点における動作電流との比である太陽光発電システム。全量買取り制導入へ向けて、太陽光発電システムの発電量を高精度に把握する必要性が増している。太陽電池は、モジュール毎にデバイスばらつきをもち、さらに日射強度や温度によって特性が変動するため、太陽電池モジュールの特性を精度よく把握することが重要となる。ある発電サイトにおける太陽光発電システムを構築する時、図2に示すような太陽電池モジュールの検査シートが納入される。この検査シートには、日射強度1kWm2、常温298Kにおける各太陽電池モジュールの短絡電流Isc、開放電圧Voc、動作電圧Vop、動作電流Iopが記載されている。発電サイトにおける太陽電池モジュール毎の特性を高精度に把握するためには、この短絡電流Isc、開放電圧Voc、動作電圧Vop、動作電流Iopの情報から、デバイスばらつきや日射強度、温度特性を考慮して、高精度に太陽電池モジュールの特性を再現する手法が重要となる。特許文献1および2における方法では、常温298KにおけるIsc、Voc、Vop、IopとRs、α、β、Kもしくは、常温298Kと他の温度例えばTbにおけるIsc、Voc、Vop、Iopの値が必要となる。数万枚の太陽電池モジュールを必要とするメガソーラーや産業用発電サイトにおいては、太陽電池モジュールの検査シートの日射強度1kWm2、常温298Kにおける短絡電流Isc、開放電圧Voc、動作電圧Vop、動作電流Iopのみが与えられる。数万枚もの太陽電池モジュールに関して、他の温度におけるIsc、Voc、Vop、Iopの測定やRsの測定、α、β、Kの測定を行うことはコスト上においても困難であり、情報量として、常温以外のIsc、Voc、Vop、IopやRs、α、β、Kの値を設定することは現実的ではない。さらに、α、β、Kの値に関しては、太陽電池モジュールメーカーから入手できる温度特性から見積もることが可能である。しかしながら、入手可能なデータは、TYP値規格に位置する太陽電池モジュールに関する温度特性である。そのため、各太陽電池モジュールのα、β、Kのばらつきをモデルに反映することができず、発電量を推定する際の精度が低下するという欠点がある。このような事情に鑑み、本発明は、上述した太陽電池モジュールのパラメータを測定することなく、太陽電池モジュールメーカーから入手可能な検査シートに記載の日射強度、温度具体的には、日射強度1kWm2、常温298Kにおける短絡電流Isc、開放電圧Voc、動作電圧Vop、動作電流Iopを用いて、デバイスばらつきや日射強度や温度による特性変動に対しても精度の高い、太陽電池モジュールの特性の演算手法を提供することを目的とする。本願において開示される発明のうち、代表的なものについて簡単に説明すれば、下記の通りである。基準状態の温度と日射強度における太陽電池の短絡電流と開放電圧と動作電圧と動作電流を示すデータの入力を受け、データを用いて太陽電池ごとに基準状態における太陽電池の特性を決定するものであって太陽電池ごとにばらつきをもつパラメータを算出し、算出されたパラメータを用いて開放電圧の温度特性と逆方向飽和電流の温度特性を算出し、開放電圧の温度特性と逆方向飽和電流の温度特性から、所定の温度における短絡電流と開放電圧と動作電圧と動作電流を算出する太陽電池の特性演算方法である。本発明によれば、直列抵抗、短絡電流の温度係数、開放電圧の温度係数、曲線補正因子や他の温度における短絡電流、開放電圧、動作電圧、動作電流を太陽電池モジュールごとに測定することなく、検査シートに記載の日射強度、温度における短絡電流、開放電圧、動作電圧、動作電流を用いて、太陽電池モジュールの特性を精度よく演算することができる。図3は、本発明の実施例1に係るモデリング手法のフローチャートである。本実施例では、検査シートに記載される太陽電池モジュールごとの短絡電流、開放電圧、動作電圧、動作電流から、太陽電池モジュールの特性式1において太陽電池モジュールごとにばらつきをもつパラメータに関して、ばらつき演算を行ってそのパラメータを決定する。決定されたパラメータを用いて逆方向飽和電流、開放電圧の温度特性を算出し、所定の温度における短絡電流、開放電圧、動作電圧、動作電流を算出する。そして、算出結果を用いて直列抵抗Rsを決定することを特徴とする。これにより、直列抵抗Rs、短絡電流の温度係数α、開放電圧の温度係数β、曲線補正因子Kや他の温度における短絡電流、開放電圧、動作電圧、動作電流を太陽電池モジュールごとに測定することなく、検査シートに記載の短絡電流、開放電圧、動作電圧、動作電流から太陽電池モジュールの特性を演算できる。さらにはデバイスばらつき、温度特性を反映して演算するため、精度よく太陽電池モジュールの特性を演算することが可能となる。入力されるデータは、検査シートに記載のデータ動作電流、動作電圧、開放電圧、短絡電流、固定パラメータボルツマン定数k、電子の電荷量q、各太陽電池モジュールを構成する太陽電池セル数Ncell、環境パラメータである。まず、基準状態日射強度1kWm2、常温298Kにおけるばらつき感度の計算を行うS401。そのため、太陽電池の特性を示す式のパラメータを以下のように分類する。固定パラメータk:ボルツマン定数q:電子の電荷量Ncell各太陽電池モジュールを構成する太陽電池セル数環境パラメータここでは標準条件を代入p:日射強度1kWm2T:太陽電池素子絶対温度ばらつきパラメータIsc:短絡電流An:ダイオード接合定数Rsh:並列抵抗ΩRs:直列抵抗ΩIs:逆方向飽和電流AIscは検査シートに記載の値なので確定する。したがって、Rs:直列抵抗Ω、Rsh:並列抵抗Ω、n:ダイオード接合定数、Is:逆方向飽和電流に関して、ばらつきを考慮した演算を行う。モデリングしたい太陽電池モジュールの日射強度Ea1kWm2、常温Taにおける短絡電流Isc、開放電圧Voc、動作電圧Vop、動作電流Iopを入力し、基準状態演算S400において、Rs、Rsh、n、Isのばらつき演算S401とフィッティングS402を行う。これにより、基準状態における太陽電池モジュールのパラメータが確定する。次に、温度特性に注目する。温度特性の計算S500は、開放電圧と逆方向飽和電流の温度特性を算出するステップS501と、所定の温度における短絡電流、開放電圧、動作電圧、動作電流を算出するステップS502と、直列抵抗をフィッティングするステップS503からなる。ステップS501において、温度係数の大きい各太陽電池モジュールの開放電圧Vocの温度特性と逆方向飽和電流Isの温度特性を求めることによって高精度な温度依存性のモデリングを行うことが可能となる。次に、求めた開放電圧Vocの温度特性と逆方向飽和電流Isの温度特性を用い、発電条件である温度Tbにおける開放電圧VocTbと逆方向飽和電流IsTbを計算するS502a。その後、短絡電流Isc、動作電圧Vop、動作電流Iopの発電サイトにおける値、IscTb、VopTb、IopTbを算出するS502b。最後に、半導体デバイスに沿った演算が不可能である直列抵抗Rsの値をフィッティングするS503。これにより、発電条件での太陽電池の特性を示すパラメータが確定し、電流-電圧特性と電力-電圧特性といった太陽電池の特性を示すカーブを求めることが可能となるS504。図4に、基準状態演算S400の詳細なフローチャートを示す。上述に示したように太陽電池の特性を構成するパラメータを固定パラメータ、環境パラメータ、ばらつきパラメータに分類する。環境パラメータは、日射強度1kWm2、常温298Kで固定であり、短絡電流Iscは検査シートに記載の値なので確定する。したがって、Rs:直列抵抗Ω、Rsh:並列抵抗Ω、n:接合定数、Is:逆方向飽和電流に関して、ばらつき感度を計算する。ばらつき感度の計算をステップS401に示す。まず、其々のばらつきパラメータに対し、暫定値Rs0、Rsh0、n0、Is0を定める。Rsに関しては、発熱によるモジュール内の損失が所定値未満となるように設定する。例えば発熱によるモジュール内の損失が2未満、つまり、Iop×Rs2Wが定格の電力Vop×Iopの2未満となるように、Rs0001Ωと設定する。Rsh:並列抵抗Ω、n:接合定数に関しては、検査シートに記載されている日射強度1kWm2、常温298KのIsc、Voc、Vop、Iopを用いて、式と式から暫定値を決定することができる。式4は、式1のIに0を代入して変形することにより得られる。式5は、式1のIにIop、VにVopを代入して変形することにより得られる。n0q1ln...Rsh0IscIsexpIop...Is0に関しては、シリコン半導体の逆方向飽和電流の値として、168×10-5Aと定める。以上で求めた値を暫定値として、式で表される太陽電池モジュールのばらつき感度を求める。ばらつき感度は、Rs、Rsh、n、Isを其々変数xとする電流Iを用いて式によって算出することができる。SIxxI∂I∂x次に、ばらつき感度の高い順にフィッティング演算を行うS402。本実施例においては、ばらつき感度がn:ダイオード接合定数、Rs:直列抵抗Ω、Is:逆方向飽和電流A、Rsh:並列抵抗Ωの順に高いとして、その順にフィッティングを行う。n:ダイオード接合定数に対しては、開放電圧への感度が高いことに着目し、nを可変して式によって得られるVoc’と検査シートに記載されているVocを比較しながら、nの値をずらしていき、合わせこみを行う。図4に示すように、所定回数例えば10000回合わせこみを行い、その結果であるn0をnとする。Voc{qln}Ncell...Rs:直列抵抗Ω、Rsh:並列抵抗Ω、Is:逆方向飽和電流Aに対しては、動作電流、動作電圧への感度が高いことに着目し、検査シートから算出されるフィルファクタ:FFと式によって計算されるIop’を用いたフィルファクタ:FF’を比較しながら、各パラメータの値をずらしていき、合わせこみを行っていく。ステップS402aと同様に、所定回数合わせこみを行った結果であるRs0、Rsh0、Is0をそれぞれRs、Rsh、Isとする。これにより、基準状態日射強度1kWm2、常温298Kにおけるばらつきパラメータの設定が完了するS402e。以上のばらつき感度計算とフィッティングにより得られたばらつきパラメータを用いることにより、日射強度1kWm2、常温298Kの標準条件において、式を用いて得られる太陽電池モジュールの電流-電圧特性のモデルは、実測にて得られる電流-電圧特性に対して高精度に再現されることが可能となる。図5に、温度特性の計算S500の詳細なフローチャートを示す。温度特性は、ステップS501において、各太陽電池モジュールのVocの温度特性と逆方向飽和電流Isの温度特性を算出することによって高精度なモデリングをすることが可能となる。開放電圧Vocの温度特性は、式を温度で微分することにより、式のように表される。∂Voc∂TqlnNcellVocTNcell...ここで、バンドギャップ電圧Egoを用いると、IsKT3exp...∂Is∂T3KT2expkexpIs...Is:逆方向飽和電流の温度特性が求まる。式を式に代入すると、式に示すような開放電圧の温度特性は求まる。∂Voc∂Tq}Ncell...以上より、デバイス毎の開放電圧の温度係数βV°Cを求めることができる。この温度係数を用いることで、ステップS502において、発電条件である温度TbにおけるVocTb、IsTb、IscTb、VopTb、IopTbを算出する。式と式により、温度Tbが可変したときのIsTb、VocTbが求まり、式に代入することによって、温度Tbにおける短絡電流IscTbが求まる。IscTb≒IsTbexpVocTb...常温におけるIscとIscTbとの比率を求めることによって、デバイス毎の短絡電流の温度係数αA°Cを求めることができる。さらに、Vopに注目すると、式7と同様に、式13によりVopを求めることができる。VopqlnIsNcell...式と式よりIsを消去すると、式14が得られる。TqlnIscNcell...Tbに関しても同様に求めると、式15が得られる。TbqlnIscTbNcell...ここで、日射強度や温度といった環境が変動した場合においても、Iop≒jIscの関係が成り立つことが知られているので、式が成り立つ。qlnIscqlnIscTb...式と式から、Vopの温度変化の式が算出されることとなる。VopTbTTbVocTb...以上のように、ステップS502bにおける演算において、IscTbの算出には式を、VopTbの算出には式を、IopTbの算出には、Iop≒jIscで求まる定数jをIscTbにかけることを行う。これにより、所定の温度における短絡電流、開放電圧、動作電圧、動作電流を算出することができる。直列抵抗Rsのばらつきは、太陽電池モジュールの組立工程で現れるため、直列抵抗Rsは、半導体のデバイス特性を用いた演算が不可能である。そこで、温度Tbにおける開放電圧VocTb、短絡電流IscTb、動作電圧VopTb、動作電流IopTbを太陽電池モジュールの特性式に適用後、直列抵抗Rsのフィッティング演算をおいて行うS503。フィッティングの方法は、図4のステップS402dと同様に、動作電流、動作電圧への感度が高いことに着目し、検査シートから算出されるフィルファクタ:FFと式によって計算されるIop’を用いたフィルファクタ:FF’を比較しながら、各パラメータの値をずらしていき、合わせこみを行っていく。以上の温度特性の計算S500の演算によって得られたパラメータと式を用いて得られる太陽電池モジュールの電流-電圧特性のモデルは、ばらつきの感度の演算を行ってばらつきパラメータをフィッティングし、さらには短絡電流、開放電圧の温度特性を算出するため、発電サイトの実測にて得られる電流-電圧特性に対して高精度に再現することが可能となる。補足説明として、Rsh:並列抵抗Ωの温度特性について述べる。Rshは太陽電池のpn接合部の漏れ電流を示すものであり、一般的な半導体デバイスの法則に従うと温度が10°C増える毎に電流が2倍増えていく。しかし、現行市販されている太陽電池モジュールは、基準の漏れ電流が非常に小さく抑えられており、式によって求まるRshから算出すると数μAオーダーの漏れ電流と見積もることができ、実使用における温度範囲における漏れ電流の変動は、太陽電池の特性に影響を及ぼさない。したがって、Rshの温度変動はないものとした。図6は、本発明の実施例2に係るMPPT制御の監視機能を搭載した太陽光発電システム1のブロック図である。太陽光発電システム1は、太陽電池アレイ2、電流検出部3、電圧検出部4、負荷変動部5、電気的負荷6、アイソレーションアンプ7、AD変換器ADC18、ADC29、最大電力点追従部MPPT部10、パルス幅変調部PWM部11、特性演算部12、比較部13によって構成される。太陽電池アレイ2からの出力はDCであるため、DCACインバータ回路を介して、商用の系統電源に接続する。太陽光発電システムの中において、インバータ回路と商用の系統電源は、電気的負荷6の役割を担っているとみなすことができる。負荷変動部5は、昇圧チョッパ等によって実現され、スイッチング素子のスイッチング動作におけるオンとオフ比である通流率を可変することによって、負荷を可変し、太陽電池アレイ2の出力を制御する。電流検出部3により検出された電流値、電圧検出部4により検出された電圧値はそれぞれ、AD変換器ADC1、ADC2により、デジタル値に変換された後、MPPT部10に入力される。MPPT部10は、最大電力点を追従するために電圧値電流値の増加および減少を決定し、ある通流率を持った制御信号を発生させる。太陽光発電システムから高効率な電力を得るために一般的に山登り法と呼ばれている最大電力点追従制御法を用いることができる。この最大電力点追従制のことをMPPT制御MaximumPowerPointTrackingと呼ぶ。MPPT部10には、マイコン、DSPやCPUが記録装置に記録されている各種ソフトウェアプログラムを読み出して実行することにより実現される。あるいは、半導体集積回路により構築されたハードウェアなどで実現してもよい。PWM部11は、その制御信号を負荷変動部5に伝えるためのパルス幅変調を行う。PWM部11からの信号は、アイソレーションアンプ7などを介してレベル変換され、負荷変動部5内のスイッチング素子を駆動する。実施例1で説明したモデリング手法を用いることにより、MPPT制御によって太陽電池アレイの最大電力が得られているか否かを、判別することが可能である。特性演算部12は、検査シート14に記載されるデータと温度計15から取得される温度情報に基づいて、実施例1で示したように、基準状態演算S400、温度特性演算S500を行い、太陽電池モジュール毎のパラメータを算出する。その後、日射計16から取得される日射強度を用いてアレイ演算を行うことにより、太陽電池アレイとしての最大電力点動作電流、動作電圧あるいは最大電力を算出する。特性演算部12の機能の一部あるいは全部は、マイコン、DSPやCPUが記録装置に記録されている各種ソフトウェアプログラムを読み出して実行することにより実現される。あるいは、半導体集積回路により構築されたハードウェアなどで実現してもよい。太陽電池アレイは、太陽電池モジュールを複数枚直列に並べたストリングと呼ばれる単位が並列に並べられることによって構成される。また各太陽電池モジュールには、逆バイアスが掛かった時、逆方向電流が流れるのを防止するため、バイパスダイオードが取り付けられている。特性演算部12で実行されるアレイ演算は、ストリング解析とアレイ解析の組み合わせを行う。ストリング解析を行うときは、複数のモジュールに流れる電流は共通であるので、ある電流が流れている場合における各太陽電池のモジュール電圧を式から求めてその和を求める。式から電圧を計算する場合、逆関数となるが、ニュートン法などの繰り返し演算を適用することで簡単に求めることが可能である。モジュールに陰が掛かるなどの影響により、モジュール電圧が負になる場合は、バイパスダイオードが機能するため、モジュール電圧≒0として考える。アレイ解析を行うときは、複数のストリングに係る電圧は共通であるので、ある電圧が掛かっている場合における各ストリングから取り出される電流を式から求め、その和を求めればよい。以上により算出されるアレイ特性から得られる最大電力は、比較部13に入力される。比較部13が、特性演算部12により算出された最大電力と電流検出部3、電圧検出部4から得られるMPPT制御による電力の比較を行うことで、MPPTが正確に行われているかを監視することが可能となる。比較結果は、MPPT部に入力され、その比較結果に基づいてMPPT部10が制御信号を出力することで最大電力点の追従性を向上することができる。図7は、本発明の実施例3に係る太陽光発電システムのブロック図である。本実施例3では、実施例1で説明したモデリング手法を用い、電流検出部3、電圧検出部4から得られる電流値、電圧値に基づいて、太陽電池アレイ2の温度と日射強度を算出する環境データ算出部17を備えることを特徴とする。実施例1、2と同構成については同符号を付し、ここでの説明は割愛する。図8は、環境データ算出部17の処理の詳細を示したものである。環境データ算出部17の機能の一部あるいは全部は、マイコン、DSPやCPUが記録装置に記録されている各種ソフトウェアプログラムを読み出して実行することにより実現される。あるいは、半導体集積回路により構築されたハードウェアなどで実現してもよい。まず、環境データ算出部17は、検査シート14に記載されるデータに基づいて基準状態演算S400、開放電圧と逆方向飽和電流の温度特性算出S501を行う。その後、基準状態演算の結果に基づいて、アレイ演算行うS801。アレイ演算に関しては実施例2と同様の手法で行うが、本実施例では、基準状態日射強度1kWm2、常温298Kにおけるアレイの短絡電流、開放電圧、動作電圧、動作電流を算出するS802。ここで、発電サイトにて測定されて入力される動作電流I’opに定数jを掛けることにより、短絡電流I’scを算出する。その後、アレイとしての日射強度1kWm2、常温298Kにおける短絡電流で割ることにより、日射強度の暫定値p’0を算出する。ここで算出された日射強度は、温度の補正が入っていないため、あくまで暫定値である。次に、ステップS802で算出された基準状態のPVパラメータと、ステップS804で算出された日射強度暫定値p’0を用いて、アレイとしての暫定の日射強度p’0、常温298Kにおける、開放電圧Voc、動作電圧Vopを算出する。発電サイトにて測定されて入力される動作電圧V’opとステップS805で算出された開放電圧Voc、動作電圧Vopを用いて、発電サイトの温度T’を算出する。算出式は、式と式を変形した式である。Tz298qVocNcell298TTzqVocNcell...この演算により、発電サイトの温度が算出されるので、式と式より、アレイとしての日射強度1kWm2、発電サイトの温度T’における短絡電流を求める。動作電流I’opから求められる短絡電流I’scを、日射強度1kWm2、温度T’における短絡電流で割ることによって、実際の日射強度p’が求まる。以上、本発明者によってなされた発明を実施形態に基づいて具体的に説明したが、本発明はそれに限定されるものではなく、その要旨を逸脱しない範囲において種々変更可能であることは言うまでもない。