(57)【要約】 【課題】  月別温度係数が正確に算出できていないため正確な年間発電量を簡単に算出できなかった。また、太陽電池の基準状態（日射強度1kw/m2、太陽電池温度25℃）における特性値（Isc,Iop,Vop,Voc,α,β,Rs,K）から指定日射強度・太陽電池温度におけるI-V 、P-V カーブを描く一般的な方法が十分に確立していなかった。【解決手段】  月別温度係数を、基準状態のI-V カーブから指定した日射強度、太陽電池温度におけるI-V カーブを作成する方法を使い算出し、この係数を使い年間発電量を算出する。このI-V カーブを作成する方法として、（1）25℃、ｂ℃、ｃ℃の基本特性値（IL,Co,n,Rsh,Rs）を求め、指定温度における各基本特性値を曲線補間により求め、これらの値を使いニュートン法により各値を求めI-V カーブを作成する。又は（2）基準状態のI-V カーブの各点を「実用的I-V カーブ変換式」によって変換してI-V カーブを作成する２つの方法がある。【選択図】        図１
2005051014.tif 000002
【特許請求の範囲】
【請求項１】
対象地点の日射量・直達比率、太陽電池設置方位・傾斜角などから計算した受光面日射強度と、対象地点の外気温・風速・受光面日射から計算した太陽電池温度から月別時刻別発電量を算出して月間集計し、次に、
該月間発電量（ＫＷｈ）を月間受光面日射量（ＫＷｈ）と太陽電池定格容量で除して月別温度係数を予め算出しておき、次に、
該月別温度係数に月間受光面日射量（ＫＷｈ）と太陽電池定格容量を掛けて月間発電量を算出し、年間集計して年間太陽光発電量を計算する方法において、
太陽電池受光面の日射強度と太陽電池温度および太陽電池特性値から太陽電池電力を算出する過程で、
太陽電池基本式
I=IL-IO*{exp(q(V+RsI)/nKoT)-1}-(V+RsI)/Rsh   
IO=CoT3 exp(-qEg/nKoT)
ここに、
Ｉ  ：出力電流[Ａ]          Co  ：飽和電流温度係数
Ｖ  ：出力電圧[Ｖ]          Eg  ：エネルギーギャップ[eＶ] 
IL  ：光起電流[Ａ]          Ｔ  ：太陽電池素子絶対温度[Ｋ]
IO  ：飽和電流[Ａ]          Ko  ：ボルツマン定数[Ｊ/Ｋ] 
Rs  ：直流抵抗[Ω]          ｑ  ：電子の電荷量[Ｃ]
Rsh ：並列抵抗[Ω]
ｎ  ：接合定数
における基本特性値(IL,Co,n,Rsh,Rs)のそれぞれについて、ａ℃、ｂ℃、ｃ℃の値を求め、
それぞれのａ℃、ｂ℃、ｃ℃での値を補間することにより指定温度におけるそれぞれの基本特性(IL,Co,n,Rsh,Rs)を求め、再び太陽電池基本式を使い
指定した日射強度、太陽電池温度における電圧−電流カーブ(I-Vカーブ) を作成し、このI-Vカーブ 上のPmax点または指定電圧における電力を求める方法
を使うことを特徴とする年間太陽光発電量シミュレーション計算方法。
【請求項２】
対象地点の日射量・直達比率、太陽電池設置方位・傾斜角などから計算した受光面日射強度と、対象地点の外気温・風速・受光面日射から計算した太陽電池温度から月別時刻別発電量を算出して月間集計し、次に、
該月間発電量（ＫＷｈ）を月間受光面日射量（ＫＷｈ）と太陽電池定格容量で除して月別温度係数を予め算出しておき、次に、
該月別温度係数に月間受光面日射量（ＫＷｈ）と太陽電池定格容量を掛けて月間発電量を算出し、年間集計して年間太陽光発電量を計算する方法において、
太陽電池受光面の日射強度と太陽電池温度および太陽電池特性値から太陽電池電力を算出する過程で、
太陽電池の特性値(Isc,Iop,Vop,Voc) から基準状態の電圧−電流カーブを作成し、この電圧−電流カーブから日射・温度の変換式により指定日射強度、指定太陽電池温度における電圧−電流カーブ(I-Vカーブ) を作成し、
このI-Vカーブ 上のPmax点または指定電圧における電力を求める方法
を使うことを特徴とする年間太陽光発電量シミュレーション計算方法。
【請求項３】
対象地点の月平均１日合計日射量・月平均直達比率、太陽電池設置方位・傾斜角、経緯度、太陽赤緯から計算した月ごと時刻別受光面日射強度と、対象地点の外気温・風速・受光面日射から計算した月ごと時刻別太陽電池温度から月別時刻別発電量を算出して月間集計し、次に、
該月間発電量（ＫＷｈ）を、単位面積当り太陽電池受光面日射量（ＫＷｈ/ m2）と太陽電池定格容量（kw/kw/ m2）で除して月別温度係数を算出し、次に、
該月別温度係数に月間受光面日射量（ＫＷｈ/ m2）と太陽電池定格容量（kw/kw/ m2）を掛けて月間発電量を算出し、年間集計して年間太陽光発電量を計算する方法において、
太陽電池受光面の日射強度と太陽電池温度および太陽電池特性値から太陽電池電力を算出する過程で、
｛01｝電圧Ｖ、電流Ｉ、日射強度１kw/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛02｝該関数Func(V,I,IL,Co,n,Rsh,Rs,T)を変数Vで微分した
関数：Dif(V,I,IL,Co,n,Rsh,Rs,T)を作成し、
｛03｝太陽電池の基準状態（太陽電池温度Ta(298Ｋ(ta=25℃))、日射強度Ea（１kw/m2 ）) での仕様値である、短絡電流Isca、最適電流Iopa −最適電圧Vopa 、開放電圧Vocaの点P1(0,Isca),P2(Vopa,Iopa),P3(Voca,0)を選択し、
｛04｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0の太陽電池温度Tに基準状態の温度Ta(298Ｋ),直列抵抗Rsに基準温度での値Rsa,および前記P1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Isca,IL,Co,n,Rsh,Rsa,Ta)= 0,
関係式:Func(Voca,0,IL,Co,n,Rsh,Rsa,Ta)= 0,
関係式:Func(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0を作成し、
｛05｝前記関数Dif(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、基準状態の温度Ta(298Ｋ)、直列抵抗Rsに基準温度での値Rsa,および前記点P2の値(Vopa,Iopa) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Dif(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0を作成し、つぎに、
｛06｝前記４つの関係式:Func(0,Isca,IL,Co,n,Rsh,Rsa,Ta)= 0,Func(Voca,0,IL,Co,n,Rsh,Rsa,Ta)= 0,Func(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0,Dif(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0 を満たす解A(ILa,Coa,na,Rsha)を、非線形解法のプログラムによって算出し、つぎに、
｛07｝前記関数：Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のIL,Co,n,Rshに、前記解A(ILa,Coa,na,Rsha) を代入して、さらにＴに基準状態の温度Ta(298Ｋ)および直列抵抗Rsに値Rsaを代入し、変数V,Iの関係式:Func(V,I,ILa,Coa,na,Rsha,Rsa,Ta)= 0を作成し、
｛08｝この式を再び非線形解法のプログラムにより約40〜50点のVに対するIの解を求め、基準状態における電圧(V)−電流(I)および電圧(V)−電力(P)(電力(P)= 電圧(V)×電流(I))の関係の点(約40〜50点)を求め、これらの点を結んだ電圧V-電流Ｉカーブ(I-V カーブ),電圧V-電力Ｐカーブ(P-V カーブ)を作成し、つぎに、
｛09｝前記太陽電池の日射強度Eb(ここでは１kw/m2)、太陽電池温度Tb(絶対温度:Tb(Ｋ)=tb(℃)＋273)における電圧Vb，電流Ibを求めるため、前記太陽電池の基準状態における日射強度Ea(１kw/m2 )、太陽電池温度Ta(298Ｋ:絶対温度:Ta=ta＋273)における短絡電流Isca,直列抵抗Ｒｓa,温度が１℃変化したときの前記短絡電流Isc の変動値α、温度が１℃変化したときの前記開放電圧Voc の変動値β、曲線補正因子Ｋとする
変換式(Va,Ia)→(Vb,Ib)：
Ib = Ia + α*(tb - ta)
Vb = Va + β*(tb - ta) - Rsa*(Ib - Ia) - K*Ib*(tb - ta) を作成し、これに前記｛08｝で作成した電圧−電流の各点またはそれを結んだI-V カーブ上の各点をIa,Va値として使用し日射強度Eb(ここでは1kw/m2 )、太陽電池温度Tb(Ｋ)における各点(電圧Vb-電流Ib:約40〜50点)を求め、これらの点を結んだI-V カーブ，P-V カーブを作成し、つぎに、
｛10｝前記｛09｝で作成したI-V カーブ上から、近接しない任意の５点を選択し、これらの点(VQ1,IQ1),(VQ2,IQ2),(VQ3,IQ3),(VQ4,IQ4),(VQ5,IQ5) を前記関係式:Func(V,I,IL,Co,n,Rsh,Rs,Tb)= 0の変数V,I、に代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Func(VQ1,IQ1,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(VQ2,IQ2,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(VQ3,IQ3,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(VQ4,IQ4,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(VQ5,IQ5,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、該５つの関係式の
解B(ILb,Cob,nb,Rshb,Rsb)を、非線形解法のプログラムによって算出し、つぎに、
｛11｝前記太陽電池の日射強度Ec(ここでは１kw/ｍ2 ),太陽電池温度Tc (絶対温度:Tc(Ｋ)=tc(℃)＋273)における電圧Vc，電流Icの関係についても前記｛09｝と同様に、
変換式(Va,Ia)→(Vc,Ic)：
Ic = Ia + α*(tc - ta)
Vc = Va + β*(tc - ta) - Rsa*(Ic - Ia) - K*Ic*(tc - ta)を使い I-V カーブ，P-V カーブを作成し、
｛12｝前記｛11｝で作成したI-V カーブ上から近接しない任意の５点を選択し、これらの値(VR1,IR1),(VR2,IR2),(VR3,IR3),(VR4,IR4),(VR5,IR5) を前記関係式:Func(V,I,IL,Co,n,Rsh,Rs,Tc)= 0の変数V,Iに代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Func(VR1,IR1,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(VR2,IR2,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(VR3,IR3,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(VR4,IR4,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(VR5,IR5,IL,Co,n,Rsh,Rs,Tc)= 0
を作成し、該５つの関係式の解C(ILc,Coc,nc,Rshc,Rsc)を、非線形解法のプログラムによって、算出し、つぎに、
｛13｝指定した（すなわち月ごと各時刻における）日射強度Em, 太陽電池温度tmを取り込み、
｛14｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(Ｋ)=ta(℃)＋273) における前記｛06｝の解A(ILa,Coa,na,Rsha)、前記温度tb(℃:Tb(Ｋ)=tb(℃)＋273) における前記｛10｝の解B(ILb,Cob,nb,Rshb,Rsb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛12｝の解C(ILc,Coc,nc,Rshc,Rsc)および入力値Rsa のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間して、指定温度tmでの特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛15｝ILm を指定された日射強度EmによりIL′m=ILm×Em÷Ea (Ea=1(KW/m2))により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL′m,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL′m,Com,nm,Rshm,Rsm,Tm)=0 を作成し、電圧(V)-電流(I)の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)- 電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、
｛16｝このI-Vカーブ 上のPmax点の電力またはパラメータとして指定した電圧における電力を求める方法を
使うことを特徴とする年間太陽光発電量シミュレーション計算方法。
【請求項４】
対象地点の月平均１日合計日射量・月平均直達比率、太陽電池設置方位・傾斜角、経緯度太陽赤緯から計算した月ごと時刻別受光面日射強度と、対象地点の外気温・風速・受光面日射から計算した月ごと時刻別太陽電池温度から月別時刻別発電量を算出して月間集計し、次に、
該月間発電量（ＫＷｈ）を、単位面積当り太陽電池受光面日射量（ＫＷｈ/ m2）と太陽電池定格容量（kw/kw/ m2）で除して月別温度係数を算出し、次に、
該月別温度係数に月間受光面日射量（ＫＷｈ/ m2）と太陽電池定格容量（kw/kw/ m2）を掛けて月間発電量を算出し、年間集計して年間太陽光発電量を計算する方法において、
太陽電池受光面の日射強度と太陽電池温度および太陽電池特性値から太陽電池電力を算出する過程で、
｛17｝指定した（すなわち月ごと各時刻における）日射強度Em, 太陽電池温度tmを取り込み、
｛18｝前記｛01｝〜｛08｝により作成した基準状態(日射強度1kw/ｍ2 ,太陽電池温度25℃）における電圧Va−電流Ia値（約４０〜５０点）について、前記｛09｝のIsca,α,β,Rsa,kを使って
変換式(Va,Ia)→(Vk,Ik)：
Ik = Ia + Isca*(Em/Ea - 1) + α*(tm - ta)
Vk = Va + β*(tm - ta) - Rsa*(Ik - Ia) - K*Ik*(tm - ta)
によって変換し、これらの電圧V-電流Ｉ値または、それらを結んだI-V カーブ，P-V カーブを作成し、
｛19｝このI-Vカーブ 上のPmax点の電力またはパラメータとして指定した電圧における電力を求め、この値を各月の時刻ごとの電力値として、月・年間積算して月・年間発電量を求める方法を使う
ことを特徴とする年間太陽光発電量シミュレーション計算方法。
【請求項５】
請求項１、２、３、４記載の太陽光発電量シミュレーション計算方法の処理プログラムを記録したコンピュータ読み取り可能なデータ記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
太陽電池は、図８のごとく、太陽の光エネルギーを電気エネルギーに直接変換するものである。すなわち光電効果の一種である光起電力効果を応用し、起電力を発生させるものである。太陽電池中に適当なエネルギーを持った光（光子）が入射すると、自由な電子と正孔が発生する。太陽電池半導体中のｐｎ接合近傍に達した電子と正孔は、それぞれｎ型半導体側、ｐ型半導体側に拡散し、両電極部に集まるので、電力が取り出せ、電圧および電流が発生するというわけである。
本発明は、実際に設置した太陽電池の受光面に入射する月間日射エネルギー量を使い、その太陽電池の年間発電量を正確かつ簡単に予測計算（シミュレーション計算）する方法に関する。そのために、月別の温度係数を正確に算出する方法に関する。そのためには、日射強度と太陽電池温度および太陽電池特性値から正確にＩ−Ｖカーブを描く必要がある。したがって本発明は、太陽電池におけるI-Vカーブ 、P-Vカーブ を作成する方法およびこの方法を使った月別温度係数の算出方法と太陽光発電量シミュレーション計算方法に関する。
【背景技術】
【０００２】
太陽電池の受光面で受ける日射エネルギーを年間シミュレーション計算した例としては、新エネルギー総合開発機構の委託業務成果報告書（(財)日本気象協会：「太陽光発電システム実用化技術開発周辺技術の研究開発『発電量基礎調査』」（1987））がある。しかしこの報告書は受光面（傾斜面）の日射エネルギーの計算であって、太陽電池の発電量まで正確にシミュレーション計算したものではなく、このような事例はほとんど見当らない。
発電量までを正確にシミュレーション計算した例としては、本発明者が発表した論文、すなわち電気学会論文１（伊賀他：「I-Vカーブ 作成法を用いた太陽光発電量シミュレーション計算プログラムの開発」、電学論Ｄ115 巻６号、1995）に見られるのみである。このように発電量までの正確な予測計算がない理由は、指定した日射強度・太陽電池温度におけるI-Vカーブ を描く技術が十分確立していなかったことによる。本発明者はこのI-Vカーブ を描く技術として前記の電気学会論文１に示した「実用的I-Vカーブ 作成法」を発明し、これを適用した発電量のシミュレーション計算方法を既に開発している。
一方、一般に年間発電量を簡易に計算する方法としては、月間日射量（ＫＷｈ）に太陽電池定格出力と温度係数を掛けて、年間集計する方法が太陽電池メーカーなどでよく実施されている。
【発明の開示】
【発明が解決しようとする課題】
【０００３】
前記のような従来の技術の現状に鑑み、より汎用的で精度の高いI-Vカーブ の作成法を使い、より正確な月別温度係数算出しこれを使った年間太陽光発電量のシミュレーション計算をより簡単に実施しようとするものである。従来は月別温度係数として一般に、夏期に０.８、冬期に０．９、春・秋に０．８５を一律に使用することが多かった。しかし、この係数では発電量に大きい影響をあたえる地域差、個々の太陽電池特性値、設置架台などの影響が反映されていない。したがって、実態にあったこの係数を算出して適用しようとするものである。
そのために本発明では、この係数をより正確に算出するのに、指定した日射強度、太陽電池温度におけるI-Vカーブ の作成方法として、計算対象太陽電池の基準状態（日射強度1KW/m2，太陽電池温度２５℃）の特性値(Isc,Iop,Vop,Voc,Rs)を使い太陽電池基本特性式を解くことにより基本特性値(IL,Co,n,Rsh) を求め、そして基準状態のI-Vカーブ を作成する方法を適用する。
すなわち、この基準状態のI-Vカーブ をもとに、次の(1) 、(2) に示す２方法により指定した日射強度・太陽電池温度におけるI-Vカーブ を描くものである。
(1) 日射強度1kw/m2において、基準温度（ａ℃=25℃）におけるI-Vカーブから他の温度（ｂ℃，ｃ℃）におけるI-Vカーブ を作成し、それぞれの温度における基本特性値（IL,Co,n,RsR,Rs）を求め、指定温度におけるこの基本特性値を曲線補間によりそれぞれ求め、次にこれらの基本特性値から指定温度・日射強度のI-Vカーブ を作成する。
(2) 基準状態のI-Vカーブ 上の各点について、「実用的I-Vカーブ 変換式」によって指定温度・日射強度に変換する。
【課題を解決するための手段】
【０００４】
請求項１の年間太陽光発電量シミュレーション計算方法は、
対象地点の日射量・直達比率、太陽電池設置方位・傾斜角などから計算した受光面日射強度と、対象地点の外気温・風速・受光面日射から計算した太陽電池温度から月別時刻別発電量を算出して月間集計し、次に、
該月間発電量（ＫＷｈ）を月間受光面日射量（ＫＷｈ）と太陽電池定格容量で除して月別温度係数を予め算出しておき、次に、
該月別温度係数に月間受光面日射量（ＫＷｈ）と太陽電池定格容量を掛けて月間発電量を算出し、年間集計して年間太陽光発電量を計算する方法において、
太陽電池受光面の日射強度と太陽電池温度および太陽電池特性値から太陽電池電力を算出する過程で次の方法を使うことに特徴がある。
太陽電池基本式
I=IL-IO*{exp(q(V+RsI)/nKoT)-1}-(V+RsI)/Rsh   
IO=CoT3 exp(-qEg/nKoT)                 
ここに、
Ｉ  ：出力電流[Ａ]          Co  ：飽和電流温度係数
Ｖ  ：出力電圧[Ｖ]          Eg  ：エネルギーギャップ[eＶ] 
IL  ：光起電流[Ａ]          Ｔ  ：太陽電池素子絶対温度[Ｋ]
IO  ：飽和電流[Ａ]          Ko  ：ボルツマン定数[Ｊ/Ｋ] 
Rs  ：直流抵抗[Ω]          ｑ  ：電子の電荷量[Ｃ]
Rsh ：並列抵抗[Ω]
ｎ  ：接合定数
における基本特性値(IL,Co,n,Rsh,Rs)のそれぞれについて、ａ℃、ｂ℃、ｃ℃の値を求め、再び太陽電池基本式を使いそれぞれのａ℃、ｂ℃、ｃ℃での値を補間することにより指定温度におけるそれぞれの基本特性(IL,Co,n,Rsh,Rs)を求め、指定した日射強度、太陽電池温度における電圧−電流カーブ(I-Vカーブ) を作成し、このI-Vカーブ 上のPmax点または指定電圧における電力を求める方法、を使うことを特徴とする。
請求項２の太陽光発電量シミュレーション計算方法は、対象地点の日射量・直達比率、太陽電池設置方位・傾斜角などから計算した受光面日射強度と、対象地点の外気温・風速・受光面日射から計算した太陽電池温度から月別時刻別発電量を算出して月間集計し、次に、
該月間発電量（ＫＷｈ）を月間受光面日射量（ＫＷｈ）と太陽電池定格容量で除して月別温度係数を予め算出しておき、次に、
該月別温度係数に月間受光面日射量（ＫＷｈ）と太陽電池定格容量を掛けて月間発電量を算出し、年間集計して年間太陽光発電量を計算する方法において、
太陽電池受光面の日射強度と太陽電池温度および太陽電池特性値から太陽電池電力を算出する過程で次の方法を使うことに特徴がある。
太陽電池の特性値(Isc,Iop,Vop,Voc) から基準状態の電圧−電流カーブを作成し、この電圧−電流カーブから日射・温度の変換式により指定日射強度、指定太陽電池温度における電圧−電流カーブ(I-Vカーブ) を作成し、
このI-Vカーブ 上のPmax点または指定電圧における電力を求める方法
を使うことを特徴とする。
請求項３の太陽光発電量シミュレーション計算方法は
対象地点の月平均１日合計日射量・月平均直達比率、太陽電池設置方位・傾斜角、経緯度太陽赤緯から計算した月ごと時刻別受光面日射強度と、対象地点の外気温・風速・受光面日射から計算した月ごと時刻別太陽電池温度から月別時刻別発電量を算出して月間集計し、次に、
該月間発電量（ＫＷｈ）を、単位面積当り太陽電池受光面日射量（ＫＷｈ/ m2）と太陽電池定格容量（kw/kw/ m2）で除して月別温度係数を予め算出しておき、次に、
該月別温度係数に月間受光面日射量（ＫＷｈ/ m2）と太陽電池定格容量（kw/kw/ m2）を掛けて月間発電量を算出し、年間集計して年間太陽光発電量を計算する方法において、
太陽電池受光面の日射強度と太陽電池温度および太陽電池特性値から太陽電池電力を算出する過程で、次の方法を使うことに特徴がある。
｛01｝電圧Ｖ、電流Ｉ、日射強度１kw/m2 での光起電流IL、飽和電流温度係数Co、接合定数n、並列抵抗Rsh、直列抵抗Rs、太陽電池温度T （絶対温度）を含んだ
関数：Func(V,I,IL,Co,n,Rsh,Rs,T) = IL - CoT3 exp(-qEg/nk0T)*(exp( q*(V+Rs*I)/(n*k0*T) )-1) - (V+Rs*I)/ Rsh - I  を作成し、つぎに、
｛02｝該関数Func(V,I,IL,Co,n,Rsh,Rs,T)を変数Vで微分した
関数：Dif(V,I,IL,Co,n,Rsh,Rs,T)を作成し、
｛03｝太陽電池の基準状態（太陽電池温度Ta(298Ｋ(ta=25℃))、日射強度Ea（１kw/m2 ）)での仕様値である、短絡電流Isca、最適電流Iopa−最適電圧Vopa 、開放電圧Vocaの点P1(0,Isca),P2(Vopa,Iopa),P3(Voca,0)を選択し、
｛04｝前記関数Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のTに基準状態の温度Ta(298Ｋ),直列抵抗Rsに基準温度での値Rsa,および前記P1,P2,P3の点の値を代入し、IL,Co,n,Rshを未知数とする
関係式:Func(0,Isca,IL,Co,n,Rsh,Rsa,Ta)= 0,
関係式:Func(Voca,0,IL,Co,n,Rsh,Rsa,Ta)= 0,
関係式:Func(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0を作成し、
｛05｝前記関数Dif(V,I,IL,Co,n,Rsh,Rs,T)= 0 に、基準状態の温度Ta(298Ｋ)、直列抵抗Rsに基準温度での値Rsa,および前記点P2の値(Vopa,Iopa) を代入して、IL,Co,n,Rshを未知数とする、
関係式:Dif(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0を作成し、つぎに、
｛06｝前記４つの関係式:Func(0,Isca,IL,Co,n,Rsh,Rsa,Ta)= 0,Func(Voca,0,IL,Co,n,Rsh,Rsa,Ta)= 0,Func(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0,Dif(Vopa,Iopa,IL,Co,n,Rsh,Rsa,Ta)= 0 を満たす解A(ILa,Coa,na,Rsha)を、非線形解法のプログラムによって算出し、つぎに、
｛07｝前記関数：Func(V,I,IL,Co,n,Rsh,Rs,T)= 0のIL,Co,n,Rshに、前記解A(ILa,Coa,na,Rsha) を代入して、さらにＴに基準状態の温度Ta(298Ｋ)および直列抵抗Rsに値Rsaを代入し、変数V,Iの関係式:Func(V,I,ILa,Coa,na,Rsha,Rsa,Ta)= 0を作成し、
｛08｝この式を再び非線形解法のプログラムにより約40〜50点のVに対するIの解を求め、基準状態における電圧(V)−電流(I)および電圧(V)−電力(P)(電力(P)= 電圧(V)×電流(I))の関係の点(約40〜50点)を求め、これらの点を結んだ電圧V-電流Ｉカーブ(I-V カーブ),電圧V-電力Ｐカーブ(P-V カーブ)を作成し、つぎに、
｛09｝前記太陽電池の日射強度Eb (ここでは１kw/m2 )、太陽電池温度Tb(絶対温度:TbＫ=tb(℃)＋273)における電圧Vb，電流Ibを求めるため、前記太陽電池の基準状態における日射強度Ea(１kw/m2 )、太陽電池温度Ta(298Ｋ:絶対温度:Ta=ta＋273)における短絡電流Isca,直列抵抗Ｒｓa,温度が１℃変化したときの前記短絡電流Isc の変動値α、温度が１℃変化したときの前記開放電圧Voc の変動値β、曲線補正因子Ｋとする
変換式(Va,Ia)→(Vb,Ib)：
Ib = Ia + α*(tb - ta)
Vb = Va + β*(tb - ta) - Rsa*(Ib - Ia) - K*Ib*(tb - ta) を作成し、これに前記｛08｝で作成した電圧−電流の各点またはそれを結んだI-V カーブ上の各点をIa,Va値として使用し日射強度Eb(1kw/m2 )、太陽電池温度Tb(Ｋ)における各点(電圧Vb-電流Ib:約40〜50点)を求め、これらの点を結んだI-V カーブ，P-V カーブを作成し、つぎに、
｛10｝前記｛09｝で作成したI-V カーブ上から、近接しない任意の５点を選択し、これらの点(VQ1,IQ1),(VQ2,IQ2),(VQ3,IQ3),(VQ4,IQ4),(VQ5,IQ5) を前記関係式:Func(V,I,IL,Co,n,Rsh,Rs,Tb)= 0の変数V,I、に代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Func(VQ1,IQ1,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(VQ2,IQ2,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(VQ3,IQ3,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(VQ4,IQ4,IL,Co,n,Rsh,Rs,Tb)= 0,
関係式:Func(VQ5,IQ5,IL,Co,n,Rsh,Rs,Tb)= 0を作成し、該５つの関係式の
解B(ILb,Cob,nb,Rshb,Rsb)を、非線形解法のプログラムによって算出し、つぎに、
｛11｝前記太陽電池の日射強度Ec(ここでは１kw/ｍ2 ),太陽電池温度Tc (絶対温度:Tc(Ｋ)=tc(℃)＋273)における電圧Vc，電流Icの関係についても前記｛09｝と同様に、
変換式(Va,Ia)→(Vc,Ic)：
Ic = Ia + α*(tc - ta)
Vc = Va + β*(tc - ta) - Rsa*(Ic - Ia) - K*Ic*(tc - ta)を使い I-V カーブ，P-V カーブを作成し、
｛12｝前記｛11｝で作成したI-V カーブ上から近接しない任意の５点を選択し、これらの値(VR1,IR1),(VR2,IR2),(VR3,IR3),(VR4,IR4),(VR5,IR5) を前記関係式:Func(V,I,IL,Co,n,Rsh,Rs,Tc)= 0の変数V,Iに代入して、IL,Co,n,Rsh,Rsを未知数とする、
関係式:Func(VR1,IR1,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(VR2,IR2,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(VR3,IR3,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(VR4,IR4,IL,Co,n,Rsh,Rs,Tc)= 0,
関係式:Func(VR5,IR5,IL,Co,n,Rsh,Rs,Tc)= 0
を作成し、該５つの関係式の解C(ILc,Coc,nc,Rshc,Rsc)を、非線形解法のプログラムによって、算出し、つぎに、
｛13｝指定した（すなわち月ごと各時刻における）日射強度Em, 太陽電池温度tmを取り込み、
｛14｝基準状態での前記温度ta(摂氏25℃:絶対温度Ta(Ｋ)=ta(℃)＋273) における前記｛06｝の解A(ILa,Coa,na,Rsha)、前記温度tb(摂氏:Tb=tb＋273) における前記｛10｝の解B(ILb,Cob,nb,Rshb,Rsb)、前記温度tc(摂氏:Tc=tc＋273) における前記｛12｝の解C(ILc,Coc,nc,Rshc,Rsc)および入力値Rsa のそれぞれ(IL,Co,n,Rsh,Rs)に関して３点について曲線補間して、指定温度tmでの特性値M(ILm,Com,nm,Rshm,Rsm) を算出し、つぎに、
｛15｝ILmを実測された日射強度EmによりIL′m=ILm×Em÷Ea(Ea=1(KW/m2 )により補正した後、前記関係式：Func(V,I,IL,Co,n,Rsh,Rs,T)=0にIL′m,Com,nm,Rshm,Rsmを代入して、Func(V,I,IL′m,Com,nm,Rshm,Rsm,Tm)=0 を作成し、電圧(V)-電流(I) の関係(約40〜50点)を非線形解法のプログラムによって求め、電圧(V)-電流(I)の関係またはそれを結んだI-V カーブ,P-V カーブを作成し、
｛16｝このI-Vカーブ 上のPmax点の電力またはパラメータとして指定した電圧における電力を求め、｛17｝この値を各月時刻ごとの電力値として、月・年間積算して月・年間発電量を求めることを特徴とする。
請求項４の太陽光発電量シミュレーション計算方法は、
対象地点の月平均１日合計日射量・月平均直達比率、太陽電池設置方位・傾斜角、経緯度太陽赤緯から計算した月ごと時刻別受光面日射強度と、対象地点の外気温・風速・受光面日射から計算した月ごと時刻別太陽電池温度から月別時刻別発電量を算出して月間集計し、次に、
該月間発電量（ＫＷｈ）を、単位面積当り太陽電池受光面日射量（ＫＷｈ/ m2）と太陽電池定格容量（kw/kw/ m2）で除して月別温度係数を予め算出しておき、次に、
該月別温度係数に月間受光面日射量（ＫＷｈ/ m2）と太陽電池定格容量（kw/kw/ m2）を掛けて月間発電量を算出し、年間集計して年間太陽光発電量を計算する方法において、
太陽電池受光面の日射強度と太陽電池温度および太陽電池特性値から太陽電池電力を算出する過程で、次の方法を使うことに特徴がある。
｛18｝指定した（すなわち月ごと各時刻における）日射強度Em, 太陽電池温度tmを取り込み、
｛19｝前記｛01｝〜｛08｝により作成した基準状態(日射強度1kw/ｍ2 ,太陽電池温度25℃）における電圧Va−電流Ia値（約４０〜５０点）について、前記｛09｝のIsca,α,β,Rsa,kを使って
変換式(Va,Ia)→(Vk,Ik)：
Ik = Ia + Isca*(Em/Ea - 1) + α*(tm - ta)
Vk = Va + β*(tm - ta) - Rsa*(Ik - Ia) - K*Ik*(tm - ta)
によって変換し、これらの電圧V-電流Ｉ値または、それらを結んだI-V カーブ，P-V カーブを作成し、
｛20｝このI-Vカーブ 上のPmax点の電力またはパラメータとして指定した電圧における電力を求め、この値を各月の時刻ごとの電力値として、月・年間積算して月・年間発電量を求めることを特徴とする。
請求項５の太陽光発電量シミュレーション計算方法の処理プログラムを記録したコンピュータ読み取り可能なデータ記憶媒体は、請求項１、２、３、４記載の発明において、これらの方法が実施可能な処理プログラムを含むことを特徴とする。
【０００５】
請求項１，２，３，４の特徴としては、月別温度係数を使うことにより簡単かつ正確に年間発電量を算出できるところに特徴がある。
また、請求項１，２，３，４の特徴としては、特定の単結晶太陽電池に限定されず他の単結晶、多結晶、アモルファス太陽電池などについても、指定した日射強度・太陽電池温度のI-Vカーブ が汎用的で精度よく作成でき、そして太陽電池の発電量が正確にシミュレーション計算できることにある。
また、請求項１，３の特徴としては、請求項２，４と同様に汎用的で精度よく太陽電池の発電量がシミュレーション計算できることの他に、日射強度は1KW/m2で他の温度における特性値（例えば５５℃などにおけるIsc,Iop,Vop,Voc など）が与えられている場合には、これらの値を使うことにより、α，β，Ｒｓ，Ｋが入手できなくても指定日射強度・太陽電池温度のI-Vカーブ が作成でき、発電量のシミュレーション計算ができることにある。
請求項５の特徴としては、請求項１，２，３，４による日射量および発電量のシミュレーション計算プログラムを記録した記憶媒体よりコンピュータ読み取り可能なため、容易に任意の条件における太陽光発電量のシミュレーション計算ができ、太陽光発電システムの設計、運用などに幅広い用途に直接役立つ。
【０００６】
ここで本発明に関する事項、用語などについて説明する。
まず年間太陽光発電量の計算の全体を理解するために、本発明者らが既に開発している太陽光発電量シミュレーション計算プログラムの概要を説明する（電気学会論文１参照）。図３は太陽光発電量シミュレーション計算プログラムのブロック図を示すものであり、３つのサブプログラムすなわち、(1) （太陽電池）受光面日射エネルギー算出サブプログラム、(2) 太陽電池（モジュール）温度算出サブプログラム、(3) 太陽電池出力算出サブプログラムより構成されている。
【０００７】
(1) 受光面の日射エネルギー算出サブプログラムの特徴は、
・任意の地点、傾斜・方位角に設置した太陽電池や、種々の追尾方式の太陽電池の指定期間の日射量が精度よく計算できるよう、計算時間のベースは後述の太陽電池算出計算を含め、３０分としている。
・太陽電池受光面日射量を求める方法としては図４で示すごとく、水平面全天日射量を直達光成分と散乱光成分に分離し、それぞれについて計算精度の高い方法で受光面日射強度を算出し、地面反射日射量も含めて総合して太陽の受光面日射量を求める方法としている。
・月平均−日合計日射量から月平均時刻別日射量を算出するため、次の手順（図５参照）で示すごとく計算精度の高い方法で求めている。
ａ．各月１５日の日出・日入時刻を求める。
ｂ．月平均−日合計日射量が基本サインカーブの面積（図５の横方向ハッチングの部分）になるようにy=a＊sin(π/to＊t)のａを求める。ただしｔ0は(日入−日出)の時間である。
ｃ．月平均−日合計日射量が周期２倍のサインカーブの面積（図５の縦方向ハッチングの部分）になるようにy=b＊sin((2π/to)＊t-π/2)+b のｂを求める
ｄ．各時刻のｙの値を上記基本サインカーブ（上記ｂ）、周期２倍のサインカーブ（上記ｃ）の式から求め、その平均値を各時刻の日射強度とする。
【０００８】
(2) 太陽電池モジュール温度算出サブプログラムの特徴は、
・上記(1) で算出した月平均時刻ごと（３０分ごと：以下同様）の受光面日射強度、月平均最高外気温・最低外気温から算出した月平均時刻ごとの外気温、および風速の入力値を使い、重回帰式により太陽電池（モジュール）温度を求めている。
【０００９】
(3) 太陽電池出力算出サブプログラムの特徴は、前記のごとく、従来の技術ではせいぜい電気学会論文１のごとく特定の太陽のI-Vカーブ をベースに、各月・各時刻におけるI-Vカーブ を作成している。本発明では特定の種類の太陽電池に限定されない汎用的I-Vカーブ 作成法とそれを使った太陽電池発電量シミュレーション計算方法を提供するものである。
【００１０】
つぎに、用語などについて説明する。
○太陽電池の温度は一般には、「太陽電池モジュール温度」と呼ばれているが、本発明では短縮して「太陽電池温度」とも呼んでいる。そしてこの温度は通常太陽電池モジュールに埋込んだ熱電対により測定する。なお、太陽電池はその構成の段階によりセル→モジュール→アレイという名称で呼ばれている。
【００１１】
○太陽電池温度に使っている記号については、小文字で示したt(℃) は摂氏を、大文字で示したT(Ｋ) は絶対温度を示す。すなわちT(Ｋ)=t(℃)+273 である。大文字Ｔは主として太陽電池基本式の中で使われており、その他のところではｔが使われている。
【００１２】
○太陽電池出力・発電量については、太陽電池の瞬間の発生電圧Ｖと瞬間の発生電流Ｉの積を太陽電池の出力（単位ＷまたはKW）といい、その時間積算値が発電量（単位whまたはkwh） である。
【００１３】
○太陽電池基本特性式は次の式である。
I=IL-Co*T3*exp(-qEg/nKoT)*(exp(q(V+RsI)/nKoT)-1)-(V+RsI)/Rsh 
ここに各記号は次のとおりである。
Ｉ  ：出力電流[Ａ]          Co  ：飽和電流温度係数
Ｖ  ：出力電圧[Ｖ]          Eg  ：エネルギーギャップ[eＶ] 
IL  ：光起電流[Ａ]          Ｔ  ：太陽電池素子絶対温度[Ｋ]
IO  ：飽和電流[Ａ]          Ko  ：ボルツマン定数[Ｊ/Ｋ] 
Rs  ：直流抵抗[Ω]          ｑ  ：電子の電荷量[Ｃ]
Rsh ：並列抵抗[Ω]
ｎ  ：接合定数
上記式は半導体の基本式に基づく式で理論的式である。この式からI-Vカーブ 上の各点を求めるには、この式の左辺のＩを右辺に移項し
Func(V,I,IL,Co,n,Rsh,Rs,T)
=IL-Co*T3*exp(-qEg/nKoT)*(exp q(V+RsI)/nKoT-1)-（V+RsI）/Rsh-I
の関数をつくり基本特性値(IL,Co,n,Rsh,Rs)の値を代入のうえ、V,I の関係を非線形解法プログラムで解いている。電気学会論文２（伊賀：「太陽電池の光照射状態での電圧−電流特性式を用いたI-Vカーブ 作成法とその活用」、電学論Ｄ116 巻１０号、1996）を参照
【００１４】
○本発明では特性値については次のとおり使い分けている。
基本特性値は太陽電池基本式中にある基本特性である。
(1) （太陽電池）基本特性値…IL,Co,n,Rsh,Rs
(2) （太陽電池）特性値…Isc,Iop,Vop,Voc,α,β,Rs,K
なお、上記の通りRsは(1) 、(2) の両方で使われている。
【００１５】
○変換式については次のとおりである。
(1) 「実用的I-V カーブ変換式」
I1=I2+Isc(E1/(E2)-1)+α(t1-t2)
V1=V2+β(t1-t2)-Rs(I1-I2)-K＊I1＊(t1-t2)
(2) 「実用的I-V カーブ変換式」の逆変換式（逆の適用）
（(1) の式のV2,I2について式を解いて変形したものである）
I2=I1+Isc(E2-E1)/E2+α(t2-t1)
V2=V1+β(t2-t1)-Rs(I2-I1)-K＊I1＊(t2-t1)
ここで、(1) 、(2) は一般に知られているJIS8913,8914,8919 の式と異なり、新しい優れた式である。(1) の式については特願平６−２６２６および前記電気学会論文１で使っている（図９参照）。
また、これらの式で使っている記号は基準状態での電圧値、電流値、日射強度、太陽電池温度をそれぞれＶ２，Ｉ２，Ｅ２，ｔ２としている。また指定状態での電圧値、電流値、日射強度、太陽電池温度をＶ１,Ｉ１,Ｅ１,ｔ１とする。
また、α：温度が１℃変化した時のIsc の変動値（A/℃）
β：温度が１℃変化した時のVoc の変動値（V/℃）
Rs：モジュールの直列抵抗（Ω）
Ｋ：曲線補正因子（Ω）
Isc：短絡電流  
Top：最適電流
Vop：最適電圧
Voc：解放電圧  である。
基準状態の電圧−電流値と測定時日射強度・太陽電池温度条件における電圧−電流値の変換式の一覧を図９に示す。
そして、図９の下欄(1)、(2)式が「実用的I-V カーブ変換式」に、また図９の右欄の(3)、(4)式が「実用的I-V カーブ変換式」の逆変換式に相当する。   
この図は電気学会論文３（伊賀：「実用的I-V カーブ作成法を使った太陽電池日射計」，電学論Ｄ，117 巻１０号，1997）に掲載されたもので右端の欄および下欄が著者の式で下欄は特願平６−２６２６の中で論文発表の前に特許申請がなされている。一般には、まだ左端又は中央の欄の式が使われることもある。
太陽光を受光している状態での太陽電池の電圧Ｖ−電流値Ｉの関係を図７に示す。
図７に示すように電圧−電流値の関係曲線をI-V 曲線（又はI-V カーブ）といい電圧−電力値の関係曲線をP-V 曲線（又はP-V カーブ）という。そして最大電力（Pmax）のときの電圧（Vop ）を最適電圧と、またそのときの電流（Iop ）を最適電流という。
【発明の効果】
【００１６】
請求項１，２，３，４の方法により、地域・太陽電池特性値・太陽電池架台など別にあらかじめ月別温度係数を計算しておくことにより、正確な年間発電量が簡単に算出できる。
また請求項１，２，３，４のI-V カーブ作成方法により、指定の日射強度・太陽電池温度におけるI-V カーブを作成し、太陽光発電量をシミュレーション計算する方法で次の効果が期待できる。
(1) 特定の種類や特定の太陽電池だけでなく一般的に、太陽電池の特性値（Isc,Iop,Vop,Voc,Rs,α,β,K）を入力し、その太陽電池を使ったシステムの月・年間発電量が正確に算出できる。すなわち適用する太陽電池の種類・設置個所・太陽電池架台等により、その発電量が正確に計算でき、効果的な太陽光発電システムの設計・運用などに役立つ。
(2) 前記「太陽光発電量シミュレーション計算プログラム（電気学会論文１）」と本発明を組合せてシステム化することにより、より精度の高い的確な発電量の計算が可能となるため、太陽光発電量シミュレーション計算がより効果的なものとなる。
(3) 太陽光発電量のシミュレーション計算において、従来からほとんどの個所で実施されてきた、受光面日射量の計算値に変換効率を掛けて発電量を求める方法に比べ、きめの細かい正確な計算ができているため、より深く細かい解析が可能である。すなわち該パラメーターによる効果の詳細な評価も可能となる。
(4) 請求項１，３の指定日射強度・太陽電池温度におけるI-V カーブを作成し、太陽光発電量をシミュレーション計算する方法では上記(1) ，(2) ，(3) で述べた効果の他に、太陽電池の特性値（Isc,Iop,Vop,Voc,α,β,Rs,K）がすべて入手できなくても、すなわち基準温度（25℃）以外の温度（例えば55℃）におけるI-V カーブやIsc,Iop,Vop,Voc の値が得られれば計算が可能となる場合がある。
請求項５では請求項１，２，３，４によるI-V カーブの作成法を含んだ太陽光発電量シミュレーション計算方法によるプログラムを記録媒体に入れ、必要時にコンピュータにより読み取り計算することが容易である。したがって太陽光発電システムの設計・工事・運用に広く役立つ。
【発明を実施するための最良の形態】
【００１７】
つぎに、本発明の実施形態を図面に基づき説明する。
図１は月別温度係数を算出するフロー図を示している。ここで、月平均時刻別発電量計算プログラムでは、本発明にかかわるＩ−Ｖカーブ作成法だけでなく本発明者らが既に開発している方法を適用してもよい。また図の中の破線でかこった部分については一般に発表されているデータがある場合はそのデータを活用してもよい。
図２は本実施形態の太陽光発電量シミュレーション計算方法の中核部分の一つである「指定の日射強度・太陽電池温度のI-V カーブの作成方法」を説明している。
S11 で太陽電池の特性値（Isc,Iop,Vop,Voc,α,β,Rs,K）を入力し、S12 で定数（Eg,Ko,e ）を入力し、S13 で太陽電池の基準状態（日射強度1kw/m2、太陽電池温度25℃）における基本特性値（IL,Co,n,Rsh ）を太陽電池基本特性式
I=IL-Co*T3exp(-qEg/nKoT)*(exp(q(V+RsI)/nKoT)-1)-(V+RSI)/Rsh
解くことにより求める。すなわち上記基本特性式に与えられた特性値（Isc,Iop,Vop,Voc ）を代入して非線形の４つの連立方程式をつくり解く。
求めた基本特性値をILa,Coa,na,Rsha とすると、S14 において上記太陽電池基本特性式に再びこれらの基本特性値を代入し、電圧Ｖと電流Ｉの関係式を得る。この関係式において電圧Ｖに対する電流Ｉをこの特性式を解くことにより求め、I-V の関係を求め、基準状態のI-V カーブを作成する。
【００１８】
S15 はS14 で作成したI-V カーブ上の各点について「実用的I-V カーブ変換式」を適用することにより、指定の日射強度、指定太陽電池温度におけるI-V カーブ（S19 ）を求める。
【００１９】
一方S16 においては日射強度は1kw/m2のままで、25℃以外の温度ｂ℃，ｃ℃におけるI-V カーブを「実用的I-V カーブ変換式」によって求める。
S20 はｂ℃におけるI-V カーブ上の５点〜８点の値を使い、前記基本特性式を解くことにより基本特性値IL,Co,n,Rsh,Rsを求める。またS21 では同様にｃ℃における基本特性値を求めている。S22 は25℃,ｂ℃,ｃ℃におけるそれぞれの基本特性値から指定太陽電池温度における基本特性値を曲線補間により求める。
【００２０】
ここで求めた指定電圧における基本特性値ILは日射強度補正を加える。S23 ではこのようにして求めた基本特性値を、前記太陽電池基本特性式に適用して電圧Ｖと電流Ｉの関係すなわちI-V カーブを求める。
図３は、本発明者が開発した「太陽光発電量シミュレーション計算プログラム」のブロック図で既に電気学会論文１で示されている。前述のとおり本発明の中核部分である「指定の日射強度・太陽電池温度におけるI-V カーブの作成方法」は、「太陽電池出力算出サブプログラム」の中に組み込まれる。
【００２１】
図４、図５は前記「太陽光発電量シミュレーション計算プログラム」の中に含まれている特徴的な機能を説明したものである。
【００２２】
図６は本発明を含んだ「太陽光発電量シミュレーション計算プログラム」を記録したデータ記憶媒体とコンピュータの図である。
【図面の簡単な説明】
【００２３】
【図１】月別温度係数算出のフロー図である。
【図２】指定の日射強度・太陽電池温度のI-V カーブ作成方法のブロック図である。
【図３】太陽光発電量シミュレーション計算プログラム（既開発）のブロック図である。
【図４】太陽電池受光面日射量の算出フロー図である。
【図５】月平均時刻別日射量カーブ作成図である。
【図６】太陽光発電量シミュレーション計算プログラムを記録したデータ記憶媒体とコンピュータの図である。
【図７】太陽電池出力特性曲線（I-V カーブ、P-V カーブ）である。
【図８】太陽光発電のしくみの図である。
【図９】測定時の日射強度・太陽電池温度における電圧Ｖ−電流値Ｉの基準状態への変換式である（電気学会論文３：「実用的I-V カーブ作成法を使った太陽電池日射計、電学論Ｄ、117 巻１０号、1997より）
【符号の説明】
【００２４】
１    コンピュータ本体
２    キーボード
３    マウス
４    モニタ
５    記憶媒体
１７    コンピュータ
【図１】
2005051014.tif 000003
【図２】
2005051014.tif 000004
【図３】
2005051014.tif 000005
【図４】
2005051014.tif 000006
【図５】
2005051014.tif 000007
【図６】
2005051014.tif 000008
【図７】
2005051014.tif 000009
【図８】
2005051014.tif 000010
【図９】
2005051014.tif 000011