温度と日射量が基準状態における太陽電池の短絡電流、開放電圧のデータを用いて、所定の逆方向飽和電流値および前記所定の逆方向飽和電流の式を構成するパラメータ値を基準として、前記太陽電池の性能を示す接合定数と前記太陽電池のキャリアの不純物濃度を算出することを特徴とする太陽電池。前記接合定数を用いて算出される第2の逆方向飽和電流と、前記接合定数を第2の逆方向飽和電流を用いて算出される第2の開放電圧と前記開放電圧を比較し、前記キャリアの不純物濃度の値を確定することを特徴とする請求項1に記載の太陽電池。前記所定の逆方向飽和電流をIs0とし、前記所定の逆方向飽和電流の式を構成するパラメータ値をK0:不純物濃度、Ego:バンドギャップeV、Ta:常温298K、k:ボルツマン定数JK、n0:接合定数とし、前記太陽電池について、Is1:逆方向飽和電流、K1:不純物濃度、Ego1:バンドギャップeV、n1:接合定数、Voc1:前記第2の開放電圧とし、前記太陽電池の検査仕様から得られる短絡電流をIscAとした場合、Is1:逆方向飽和電流とn1:接合定数は、下記の2式によって得られ、K1:不純物濃度は、前記太陽電池の検査仕様から得られるVoc:開放電圧と前記第2の開放電圧Voc1を比較しフィッティングを行うことを特徴とする請求項2に記載の太陽電池。Is1Is0expEgo1expEgoVoc1qln所定の温度におけるIs1:逆方向飽和電流A、K1:不純物濃度、n1:接合定数の値を前記2式に適用し、Rs:直流抵抗Ωのフィッティングを行うことを特徴とする請求項3に記載の太陽電池。所定の温度におけるIs1:逆方向飽和電流A、K1:不純物濃度、n1:接合定数の値を前記2式に適用し、Rs:直流抵抗Ωのフィッティングを行うことを特徴とする請求項3に記載の太陽電池。前記算出されたK1:不純物濃度、n1:接合定数を用いて、開放電圧と逆方向飽和電流の温度係数を算出し、それらによって得られた値から所定の温度における短絡電流、開放電圧、最大電力動作点を算出することを特徴とする請求項1に記載の太陽電池。前記算出されたK1:不純物濃度、n1:接合定数と、所定の温度におけるIs1:前記算出された短絡電流、開放電圧、最大電力動作点から所定の温度におけるRs:直流抵抗Ω、Rsh:シャント抵抗Ωを算出することによって得られる所定の温度、所定の日射強度における前記太陽電池の特性と計測により得られる特性の差分を求めることによって、前記太陽電池の光照射効果を定量化することを特徴とする請求項6に記載の太陽電池。複数の太陽電池モジュールを直列に接続した太陽電池ストリングを一つの単位として、複数の前記の太陽電池ストリングを並列接続することによって構成された太陽電池アレイと、前記太陽電池アレイの動作電圧を制御する電圧制御部と、前記電圧制御部での制御情報として前記太陽電池アレイの出力電圧を検出する電圧検出部と、前記太陽電池アレイの出力電流を検出する電流検出部と、を有し、請求項1に記載の算出を行うことを特徴とする太陽光発電システム。前記電圧制御部の制御によって前記電圧検出部と前記電流検出部から得られる電力と請求項1に記載の通り算出される最大電力を比較する比較手段を特徴とする請求項8に記載の太陽光発電システム。前記比較手段の結果をもとに、前記電圧制御部によって制御する電圧の値を可変させることを特徴とする請求項9に記載の太陽光発電システム。前記電流検出部から得られる電流と所定の温度におけるIs1:逆方向飽和電流Aから最大電力の電圧移動比を算出し、前記電圧検出部から得られる電圧比を比較する比較手段を特徴とする請求項8に記載の太陽光発電システム。前記請求項1に記載の通り算出される開放電圧と逆方向飽和電流の温度係数と前記電圧検出部と前記電流検出部から得られる動作電圧と動作電流の値から太陽電池の温度を算出することを特徴とする請求項8に記載の太陽光発電システム。前記算出された温度と開放電圧と逆方向飽和電流の温度係数から日射量が基準状態で前記算出された温度における短絡電流を算出し、前記短絡電流と前記動作電流に定数を掛けた値の比を取ることによって太陽電池に掛かっている日射量を算出することを特徴とする請求項8に記載の太陽光発電システム。前記定数は、温度と日射量が基準状態における太陽電池の短絡電流と最大電力点における動作電流の比であることを特徴とする請求項13に記載の太陽光発電システム。上記に挙げた特許文献1における方法では、常温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のばらつきをモデルに反映することができず、発電量を推定する際の精度が低下するという欠点がある。非特許文献1における方法では、日射強度1kWm2、常温298Kにおける短絡電流Isc、開放電圧Voc、動作電圧Vop、動作電流Iopのみを用いて、デバイスばらつきや日射量や温度による特性変動に対しても精度の高い太陽電池モジュールの特性を把握するモデリングの演算手法を開発することが可能である。非特許文献1に記載の方法は、太陽電池モジュール特性の分布が正規分布であることを利用した感度演算を行っているため、結晶欠陥が混じることによって正規分布を成さないような分布を示す太陽電池モジュールに関しては、モデリング精度が低下するという欠点がある。したがって、低コストで製造され、結晶欠陥の混じったシリコン結晶系太陽電池や製造上結晶欠陥を含むCIGS太陽電池においても日射強度1kWm2、常温298Kにおける短絡電流Isc、開放電圧Voc、動作電圧Vop、動作電流Iopのみを用いて、デバイス性能を示すパラメータを算出し、日射量や温度による特性変動に対しても精度の高い太陽電池モジュールの特性を把握するモデリングの演算手法を開発することが課題である。本願による課題を解決するための手段のうち、代表的なものを例示すれば、太陽電池であって、温度と日射量が基準状態における太陽電池の短絡電流、開放電圧のデータを用いて、所定の逆方向飽和電流値および前記所定の逆方向飽和電流の式を構成するパラメータ値を基準として、前記太陽電池の性能を示す接合定数と前記太陽電池のキャリアの不純物濃度を算出することを特徴とする。本願によれば、低コストで製造され、結晶欠陥の混じったシリコン結晶系太陽電池や製造上結晶欠陥を含むCIGS太陽電池においても日射強度1kWm2、常温298Kにおける短絡電流Isc、開放電圧Voc、動作電圧Vop、動作電流Iopのみを用いて、デバイス性能を示すパラメータを算出し、日射量や温度による特性変動に対しても精度の高い太陽電池モジュールの特性を把握するモデリングの演算手法を提供しうる。図3は、本発明の実施例1に係るモデリング手法のフローチャートである。式1を構成する逆方向飽和電流をIs1とすると、Is1は、K1:不純物濃度、Ego1:バンドギャップeV、T:太陽電池素子絶対温度K、k:ボルツマン定数JK、n1:接合定数のパラメータを用い、式によって表すことができる。Is1K1T3exp-...モデリングしたい太陽電池モジュールの日射強度Ea1kWm2、常温における短絡電流Isc、開放電圧Voc、動作電圧Vop、動作電流Iopを入力し、演算S400において、Is1、K1、n1の算出S401とフィッティングS402a、シャント抵抗の算出S402bを行う。これにより、基準状態における太陽電池モジュールのパラメータが確定する。次に、温度特性に注目する。温度特性の計算S500は、S501において、温度係数の大きい各太陽電池モジュールの開放電圧Vocの温度特性と逆方向飽和電流Is1の温度特性を求めることによって高精度な温度依存性のモデリングを行うことが可能となる。求めた開放電圧Vocの温度特性と逆方向飽和電流Is1の温度特性を用い、S502aにおいて暴露条件である温度Tbにおける開放電圧VocTbと逆方向飽和電流Is1Tbを計算する。その後、S502bにおいて、短絡電流Isc、動作電圧Vop、動作電流Iopの暴露サイトにおける値、IscTb、VopTb、IopTbを算出する。最後に、半導体デバイスに沿った演算が不可能である直列抵抗Rsの値をS503aにおいてフィッティングし、求まった直列抵抗Rsを用いて、S503bを用いてシャント抵抗Rshを算出することにより、暴露条件での太陽電池の特性を示すパラメータが確定し、電流‐電圧特性と電力‐電圧特性といった太陽電池の特性を示すカーブを求めることが可能となる。図4に基準状態演算S400の詳細なフローチャートを示す。上述に示したように式を用いS401において、Is1、K1、n1の算出を行う。リファレンスとなる正常な結晶シリコンの逆方向飽和電流:Is0を構成する各パラメータの値を基準値とし、求めたいモジュールの逆方向飽和電流:Is1を算出する。リファレンスとなる正常な結晶シリコンの逆方向飽和電流:Is0をK0:不純物濃度、Ego:バンドギャップeV、T:太陽電池素子絶対温度K、k:ボルツマン定数JK、n0:接合定数のパラメータを用い、式のように表す。Is0K0T3exp-...式と式からIs1を算出すると、Is1Is0expEgo1expEgo...ここで、リファレンスとして、Is0:168×10-5A、K0:01、Ego:1205eV、n0:18と定める。また、Ta298Kであり、k:13806503×1023である。Ego1は、モデリングの対象がシリコン結晶系の場合、1205eVを用い、CISの場合、115~120eVを用いる。まず、n1とK1の暫定値を定め、式を用いて、Is1を求める。次に、太陽電池モジュールの日射強度Ea1kWm2、常温における短絡電流Isc、開放電圧Vocと算出したIs1を式に代入して、n1を求める。n1q-1ln...式により求まったn1をさらに式に代入しIs1を求める。ここで算出されたn1とIs1を式に代入して開放電圧Voc1を求める。Voc1qln...S401の最後に、算出した開放電圧:Voc1と検査シートのVocを比較し、Voc1Vocであれば、K1K100001とし、それ以外であれば、K1K1-00001として、キャリアの不純物濃度を調整する。以上の式の計算からキャリアの不純物濃度を調整に至る動作を各パラメータ値が収束するまで繰り返し演算を行う。図4においては、1例として、繰り返し1000回とした。Rs:直流抵抗Ωに対しては、S402aに示すように、検査シートから算出されるフィルファクタ:FFと式によって計算されるIop’を用いたフィルファクタ:FF’を比較しながら、各パラメータの値をずらしていき、合わせこみを行っていく。Rsが求まったら、不明なパラメータはRshのみとなるので、式を用いて算出するS402b。以上の計算により求まったパラメータを用い、式から得られる日射強度1kWm2、常温条件下での太陽電池モジュールの電流‐電圧特性は、実測で得られる電流‐電圧特性と高い相関性を持つ。図5に温度特性演算部S500の詳細なフローチャートを示す。温度特性に注目する。温度特性は、S501において、各太陽電池モジュールのVocの温度特性と逆方向飽和電流Isの温度特性を算出することによって高精度なモデリングをすることが可能となる。開放電圧Vocの温度特性は式を温度で微分することにより、式のように表される。∂Voc∂Tqln-VocT-...ここで、式を変形して、∂Is1∂T3K1T2expexpIs1...より、Is1:逆方向飽和電流の温度特性が求まる。式を式に代入すると、式に示すような開放電圧の温度特性が求まる。∂Voc∂T-1Tq...以上より、デバイス毎の開放電圧の温度係数βV°Cを求めることができる。この温度係数を用いることで、S502において、暴露条件である温度TbにおけるVocTb、IsTb、IscTb、VopTb、IopTbを算出する。式と式により、温度Tbが可変したときのIs1Tb、VocTbが求まり、式に代入することによって、温度Tbにおける短絡電流IscTbが求まる。IscTb≒Is1TbexpVocTb...ここで、常温におけるIscとIscTbとの比率を求めることによって、デバイス毎の短絡電流の温度係数αA°Cを求めることができる。さらに、Vopに注目すると、VopqlnIs1...式と式よりIs1を消去すると、TqlnIsc...Tbに関しても同様に求めて、TbqlnIscTb...ここで非特許文献ComparisonofPhotovoltaicArrayMaximumPowerPointTrackingTechniquesEsramTrishan、IEEETransactionsonEnergyConversionVol11、No2、pp439449、2007の中において、日射量や温度といった環境が変動した場合においても、Iop≒jIscの関係が成り立つと記載されているので、qlnIscqlnIscTb...が成り立ち、Vopの温度変化の式が以下の通り算出されることとなる。VopTbTbVocb...つまり、S502における演算は、IscTbの算出には式を、VopTbの算出には式を、IopTbの算出には、Iop≒jIscで求まる定数jをIscTbにかけることを行う。直列抵抗:Rsは、半導体のデバイス特性を用いた演算が不可能なため、直列抵抗:Rsのフィッティング演算をS503において行う。フィッティングの方法は、図4のS402aと同様に、検査シートから算出されるフィルファクタ:FFと式によって計算されるIop’を用いたフィルファクタ:FF’を比較しながら、各パラメータの値をずらしていき、合わせこみを行っていく。最後に、S402bと同様に、S503において、式を用いてRshを算出する。以上のS500の演算によって得られたパラメータを用いて式を用いて得られる太陽電池モジュールの電流‐電圧特性は、暴露サイトの実測にて得られる電流‐電圧特性を高精度に再現することが可能となる。図6は、本発明の実施例2に係るMPPT制御の監視機能を搭載した太陽光発電システムのブロック図である。太陽光発電システム1は、電圧検出部3、電流検出部4、負荷変動部5、電気的負荷6によって構成される。負荷変動部5は、昇圧チョッパ等によって実現され、スイッチング素子のスイッチング動作におけるオンとオフ比である通流率を可変することによって、負荷を可変し、太陽電池の出力を制御する。ある通流率を持った制御信号は、MPPT部と記載した制御部において生成され、PWM化された後、アイソレーションアンプ7などを介してレベル変換され、負荷変動部5内のスイッチング素子を駆動する。太陽電池の出力電圧を検出する電圧検出部3、太陽電池の出力電流を検出する電流検出部4において検出された電圧情報と電流情報は、AD変換器ADC8とADC9よってデジタル値に変換された後、MPPT部に送信される。これにより、太陽電池パネルの出力をフィードバックしながら、通流率を制御することできる。太陽電池アレイからの出力はDCであるため、DC/ACインバータ回路を介して、商用の系統電源に接続する。太陽電池システムの中において、インバータ回路と商用の系統電源は、電気的負荷の役割を担っているとみなすことができる。太陽光発電システムから高効率な電力を得るために一般的に用いられている山登り法と呼ばれる最大電力点追従制御法が用いられている。この最大電力点追従制のことをMPPT制御MaximumPowerPointTrackingと呼ぶ。このMPPT制御によって太陽電池の最大電力が得られているか否かモデリング手法を用いることで判別することが可能である。実施例1で示したように、基準状態演算S400、温度特性演算S502を行い、太陽電池毎のパラメータを算出する。その後、アレイ演算を行うことにより、太陽電池アレイとしての最大電力点を算出する。太陽電池アレイ2は、太陽電池モジュールを複数枚直列に並べたストリングと呼ばれる単位が並列に並べられることによって構成される。また各太陽電池モジュールには、逆バイアスが掛かった時、逆方向電流が流れるのを防止するため、バイパスダイオードが取り付けられている。アレイ演算は、ストリング解析とアレイ解析の組み合わせを行う。ストリング解析を行うときは、複数のモジュールに流れる電流は共通であるので、ある電流が流れている場合における各太陽電池のモジュール電圧を式から求め、その和を求める。式から電圧を計算する場合、逆関数となるが、ニュートン法などの繰り返し演算を適用することで簡単に求めることが可能である。モジュールに陰が掛かるなどの影響により、モジュール電圧が負になる場合は、バイパスダイオードが機能するため、モジュール電圧≒0として考える。アレイ解析を行うときは、複数のストリングに係る電圧は共通であるので、ある電圧が掛かっている場合における各ストリングから取り出される電流を式から求め、その和を求めればよい。以上により算出されるアレイ特性から得られる最大電力と電圧検出部3、電流検出部4で得られるMPPT制御による電力の比較を行うことで、MPPT制御が正確に行われているかを監視することが可能となる。図7、図8は、本発明の実施例3に係るMPPT制御の監視機能を搭載した太陽光発電システムのブロック図である。図7の比較部13bについての算出フローを図8に示す。まず、温度計15により検出されたTbと式を用いて、Is1を求める。次に、あるタイミングで日射計16、電流検出部4によって検出された日射量をEbN、MPPにおける電流値をIopNとし、次に検出された日射量をEbN1、MPPにおける電流値をIopN1とすると、式は、式、式のように書き換えられ、理想的な電圧比は、式で表すことが出来る。VopNqlnIs1...VopN1qlnIs1...ln-lnlnlnここで、電圧検出部3によって検出されたVopN1とVopNの比と式を比較することによって、MPPT制御がきちんと行われているか監視することが可能となる。図9、図10は、本発明の実施例4に係る太陽電池を日射計と温度計として利用可能な太陽光発電システムのブロック図であり、図9の環境算出部17についての算出フローを図10に示す。まず、実施例2で示したアレイ演算までの手順を行うことによって、アレイとしての日射強度10kWm2、常温298Kにおける短絡電流、開放電圧、動作電圧、動作電流を算出しておく。ここで、暴露サイトにて測定される動作電流I’opに定数jを掛けることにより、短絡電流I’scを算出し、アレイとしての日射強度10kWm2、常温298Kにおける短絡電流で割ることにより、日射量の暫定値p’0を算出する。ここで算出された日射量は、温度の補正が入っていないため、あくまで暫定値である。次に、アレイとしての暫定の日射強度p’0kWm2、常温298Kにおける、開放電圧Voc、動作電圧Vopを算出しておく。ここで、暴露サイトにて測定される動作電圧V’opと開放電圧Voc、動作電圧Vopを用いて、暴露サイトの温度T’を算出する。算出式は、式と式を変形した式である。TzV’op-Voc298-q-VocNcell298T’Tzq-VocNcell...この演算により、暴露サイトの温度が判明するので、式と式より、アレイとしての日射強度10kWm2、温度T’における短絡電流が求まるので、この値で、短絡電流I’scを割ることによって、実際の日射量p’が求まる。図11は、本発明の実施例5に係るCIS太陽電池のモデリング結果を示したものである。S900はモデリングによる電流電圧特性を示しており、s901は暴露サイトで測定した電流電圧特性を示している。出荷時に測定した電流電圧特性は、s900と相関の高い電流電圧特性を示すが、CISの光照射効果によりs902の差分が発生する。これは、本発明を用いることにより、暴露サイトにおいて光照射効果の影響も定量化できることを意味する。