(57)【要約】　　　（修正有）【課題】太陽電池の出力（発電）を正確に評価するために、太陽の日射強度の変化に即応した（時間遅れのない）また正確な出力信号の出る日射強度の測定方法と日射計を提供する。【解決手段】日射強度の変化に対する出力信号の応答の早さから、太陽電池（セル）を利用し、その日射強度と太陽電池の光起電流の良好な応答特性を使い、日射強度を算出する。すなわち短絡電流Ｓ２２と太陽電池温度の測定値Ｓ２０および該太陽電池の特性値から太陽電池基本式を解くことにより光起電流を求め、この値から日射強度を求めるものである。特に太陽電池基本特性値の算出で種々の工夫をしている。【選択図】　　　　図２
2004077311.tif 000002
【特許請求の範囲】
【請求項１】
太陽電池基本式
Ｉ＝ＩＬ−ＩＯ＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／ｎＫｏＴ）−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ
ＩＯ＝ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）
ここに、
Ｉ　：出力電流［Ａ］　　　　　Ｃｏ　：飽和電流温度係数
Ｖ　：出力電圧［Ｖ］　　　　　　　　　　Ｅｇ　：エネルギーギャップ［ｅＶ］
ＩＬ　：光起電流［Ａ］　　　　　　　　　　Ｔ　：太陽電池（素子）絶対温度［Ｋ］
ＩＯ　：飽和電流［Ａ］　　　　　　　　　　Ｋｏ　：ボルツマン定数［Ｊ／Ｋ］
Ｒｓ　：直流抵抗［Ω］　　　　　　　　　　ｑ　：電子の電荷量［Ｃ］
Ｒｓｈ　：並列抵抗［Ω］
ｎ　：接合定数
における基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓ，Ｒｓｈ）のそれぞれについてａ℃，ｂ℃，ｃ℃の値を求め、それぞれのａ℃，ｂ℃，ｃ℃での値を使い温度の曲線補間をすることにより太陽電池の測定温度におけるそれぞれの基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓ，Ｒｓｈ）を求め、
これらの値に加え、測定時の短絡電流Ｉｓｊ（Ａ）　および太陽電池温度ｔｊ（℃）（絶対温度Ｔ（Ｋ）＝ｔｊ（℃）＋２７３）を日射強度Ｅ（ＫＷ／ｍ２）　の算出式
Ｅ＝（１／ＩＬ）＊［Ｉｓｊ＋Ｃｏ・Ｔ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑＲｓＩｓｊ／ｎＫｏＴ）−１｝＋Ｒｓ＊Ｉｓｊ／Ｒｓｈ］
に適用して日射強度Ｅを求めることを特徴とする光起電流による日射強度の測定方法。
【請求項２】
Ｒｓおよびその温度特性が与えられている場合、
太陽電池基本式
Ｉ＝ＩＬ−ＩＯ＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／ｎＫｏＴ）−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ
ＩＯ＝ＣｏＴ３　ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）
ここに、
Ｉ　：出力電流［Ａ］　　　　　Ｃｏ　：飽和電流温度係数
Ｖ　：出力電圧［Ｖ］　　　　　　　　　　Ｅｇ　：エネルギーギャップ［ｅＶ］
ＩＬ　：光起電流［Ａ］　　　　　　　　　　Ｔ　：太陽電池（素子）絶対温度［Ｋ］
ＩＯ　：飽和電流［Ａ］　　　　　　　　　　Ｋｏ　：ボルツマン定数［Ｊ／Ｋ］
Ｒｓ　：直流抵抗［Ω］　　　　　　　　　　ｑ　：電子の電荷量［Ｃ］
Ｒｓｈ　：並列抵抗［Ω］
ｎ　：接合定数
において、基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓｈ）のそれぞれについてａ℃，ｂ℃，ｃ℃の値を求め、それぞれのａ℃，ｂ℃，ｃ℃での値を使い温度の曲線補間をすることにより太陽電池の測定温度におけるそれぞれの基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓ，Ｒｓｈ）を求め、
これらの値とａ℃，ｂ℃，ｃ℃におけるＲｓの値に加え、測定時の短絡電流Ｉｓｊ（Ａ）および太陽電池温度ｔｊ（℃）（絶対温度Ｔ（Ｋ）＝ｔｊ（℃）＋２７３）を日射強度Ｅ（ＫＷ／ｍ２）の算出式
Ｅ＝（１／ＩＬ）＊［Ｉｓｊ＋Ｃｏ・Ｔ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑＲｓＩｓｊ／ｎＫｏＴ）−１｝＋Ｒｓ＊Ｉｓｊ／Ｒｓｈ］
に適用して日射強度Ｅを求めることを特徴とする光起電流による日射強度の測定方法。
【請求項３】
｛０１｝電圧Ｖ，電流Ｉ，日射強度１ＫＷ／ｍ２での光起電流ＩＬ，飽和電流温度係数Ｃｏ，接合定数ｎ，並列抵抗Ｒｓｈ　，直列抵抗Ｒｓ，太陽電池温度Ｔ（絶対温度）を含んだ
関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／（ｎＫｏＴ））−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ　を作成し、つぎに
接合定数ｎ，並列抵抗Ｒｓｈ　，直列抵抗Ｒｓ，太陽電池温度Ｔ（絶対温度）を含んだ
関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／（ｎＫ
ｏＴ））−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ　を作成し、つぎに
｛０２｝該関数Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）を変数Ｖで微分した関数：Ｄｉｆ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）を作成し、
｛０３｝太陽電池の基準状態（モジュール温度Ｔａ（２９８Ｋ（ｔａ＝２５℃））、日射強度Ｅａ（１ｋｗ／ｍ２））　での仕様値である、短絡電流Ｉｓｃａ、最適電流Ｉｏｐａ　−最適電圧Ｖｏｐａ　、開放電圧Ｖｏｃａの点Ｐ１（０，Ｉｓｃａ），Ｐ２（Ｖｏｐａ，Ｉｏｐａ），Ｐ３（Ｖｏｃａ，０）を選択し、
｛０４｝前記関数Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝　０のＴに基準状態の温度Ｔａ（２９８Ｋ），直列抵抗Ｒｓに基準温度での値Ｒｓａ，および前記Ｐ１，Ｐ２，Ｐ３の点の値を代入し、ＩＬ，Ｃｏ，ｎ，Ｒｓｈを未知数とする
関係式：Ｆｕｎｃ（０，Ｉｓｃａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，
関係式：Ｆｕｎｃ（Ｖｏｃａ，０，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，
関係式：Ｆｕｎｃ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０を作成し、
｛０５｝前記関数Ｄｉｆ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝　０　に、基準状態の温度Ｔａ（２９８Ｋ）、直列抵抗Ｒｓに基準温度での値Ｒｓａ　、および前記点Ｐ２の値（Ｖｏｐａ，Ｉｏｐａ）　を代入して、ＩＬ，Ｃｏ，ｎ，Ｒｓｈを未知数とする、
関係式：Ｄｉｆ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０を作成し、つぎに、
｛０６｝前記４つの関係式：Ｆｕｎｃ（０，Ｉｓｃａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，Ｆｕｎｃ（Ｖｏｃａ，０，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，Ｆｕｎｃ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，Ｄｉｆ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０　を満たす解Ａ（ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ）を、非線形解法のプログラムによって算出し、つぎに、
｛０７｝前記関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝　０のＩＬ，Ｃｏ，ｎ，Ｒｓｈに、前記解Ａ（ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ）　を代入して、さらにＴに基準状態の温度Ｔａ（２９８Ｋ）および直列抵抗Ｒｓに値Ｒｓａを代入し、変数Ｖ，Ｉの関係式：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ，Ｒｓａ，Ｔａ）＝　０を作成し、
｛０８｝この式を再び非線形解法のプログラムにより約４０〜５０点のＶに対するＩの解を求め、基準状態における電圧（Ｖ）−電流（Ｉ）および電圧（Ｖ）−電力（Ｐ）（電力（Ｐ）＝　電圧（Ｖ）×電流（Ｉ））　の関係の点（約４０〜５０点）を求め、これらの点を結んだ電圧Ｖ−電流Ｉカーブ（Ｉ−Ｖ　カーブ），電圧Ｖ−電力Ｐカーブ（Ｐ−Ｖ　カーブ）を作成し、つぎに、
｛０９｝前記太陽電池の日射強度Ｅｂ（ここでは１ｋｗ／ｍ２）、太陽電池モジュール温度Ｔｂ（絶対温度：Ｔｂ（Ｋ）＝ｔｂ（℃）＋２７３）　における電圧Ｖｂ，電流Ｉｂを求めるため、前記太陽電池の基準状態における日射強度Ｅａ（１ｋｗ／　ｍ２）、太陽電池温度ｔａ（２９８Ｋ：絶対温度：Ｔａ＝ｔａ＋２７３）における短絡電流Ｉｓｃａ，モジュール直列抵抗Ｒｓａ，温度が１℃変化したときの前記短絡電流Ｉｓｃａの変動値α、温度が１℃変化したときの前記開放電圧Ｖｏｃａ　の変動値β、曲線補正因子Ｋとする
変換式（Ｖａ，Ｉａ）→（Ｖｂ，Ｉｂ）：
Ｉｂ　＝　Ｉａ　＋　α＊（ｔｂ　−　ｔａ）
Ｖｂ　＝　Ｖａ　＋　β＊（ｔｂ　−　ｔａ）　−　Ｒｓａ＊（Ｉｂ　−　Ｉａ）　−　Ｋ＊Ｉｂ＊（ｔｂ　−　ｔａ）　を作成し、これに前記｛０８｝で作成した電圧−電流の各点またはそれを結んだＩ−Ｖ　カーブ上の各点を電流、電圧値すなわちＩａ，Ｖａ値として使用し、前式に適用し日射強度Ｅｂ（ここでは１ｋｗ／ｍ２）、モジュール温度Ｔｂ（Ｋ）　における各点（電圧Ｖｂ−電流Ｉｂ：約４０〜５０点）　を求め、これらの点を結んだＩ−Ｖ　カーブ，Ｐ−Ｖ　カーブを作成し、つぎに、
｛１０｝前記｛０９｝で作成したＩ−Ｖ　カーブ上から、近接しない任意の５点を選択し、これらの点（ＶＱ１，ＩＱ１），（ＶＱ２，ＩＱ２），（ＶＱ３，ＩＱ３），（ＶＱ４，ＩＱ４），（ＶＱ５，ＩＱ５）　を前記関係式：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｂ）＝　０の変数Ｖ，Ｉ、に代入して、ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数とする、
関係式：Ｆｕｎｃ（ＶＱ１，ＩＱ１，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｂ）＝　０，
関係式：Ｆｕｎｃ（ＶＱ２，ＩＱ２，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｂ）＝　０，
関係式：Ｆｕｎｃ（ＶＱ３，ＩＱ３，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｂ）＝　０，
関係式：Ｆｕｎｃ（ＶＱ４，ＩＱ４，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｂ）＝　０，
関係式：Ｆｕｎｃ（ＶＱ５，ＩＱ５，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｂ）＝　０を作成し、該５つの関係式の解Ｂ（ＩＬｂ，Ｃｏｂ，ｎｂ，Ｒｓｈｂ，Ｒｓｂ）を、非線形解法のプログラムによって算出し、つぎに、
｛１１｝前記太陽電池の日射強度Ｅｃ（ここでは１ｋｗ／ｍ２），太陽電池温度ｔｃ　（絶対温度：Ｔｃ（Ｋ）＝ｔｃ（℃）＋２７３）における電圧Ｖｃ，電流Ｉｃの関係についても前記｛０９｝と同様に、
変換式（Ｖａ，Ｉａ）→（Ｖｃ，Ｉｃ）：
Ｉｃ　＝　Ｉａ　＋　α＊（ｔｃ　−　ｔａ）
Ｖｃ　＝　Ｖａ　＋　β＊（ｔｃ　−　ｔａ）　−　Ｒｓａ＊（Ｉｃ　−　Ｉａ）　−　Ｋ＊Ｉｃ＊（ｔｃ　−　ｔａ）を使い　Ｉ−Ｖ　カーブ，Ｐ−Ｖ　カーブを作成し、
｛１２｝前記｛１１｝で作成したＩ−Ｖ　カーブ上から近接しない任意の５点を選択し、これらの値（ＶＲ１，ＩＲ１），（ＶＲ２，ＩＲ２），（ＶＲ３，ＩＲ３），（ＶＲ４，ＩＲ４），（ＶＲ５，ＩＲ５）　を前記関係式：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｃ）＝　０の変数Ｖ，Ｉに代入して、ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数とする、
関係式：Ｆｕｎｃ（ＶＲ１，ＩＲ１，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｃ）＝　０，
関係式：Ｆｕｎｃ（ＶＲ２，ＩＲ２，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｃ）＝　０，
関係式：Ｆｕｎｃ（ＶＲ３，ＩＲ３，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｃ）＝　０，
関係式：Ｆｕｎｃ（ＶＲ４，ＩＲ４，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｃ）＝　０，
関係式：Ｆｕｎｃ（ＶＲ５，ＩＲ５，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｃ）＝　０
を作成し、該５つの関係式の解Ｃ（ＩＬｃ，Ｃｏｃ，ｎｃ，Ｒｓｈｃ，Ｒｓｃ）を、非線形解法のプログラムによって、算出し、つぎに、
｛１３｝太陽電池の実測温度ｔｊ（摂氏：絶対温度はＴｊ＝ｔｊ＋２７３）を取り込み、
｛１４｝基準状態での前記温度ｔａ（摂氏２５℃：絶対温度Ｔａ（Ｋ）（＝ｔａ（℃）＋２７３）における前記｛０６｝の解Ａ（ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ）とＲｓａ　、前記温度ｔｂ（摂氏：Ｔｂ＝ｔｂ＋２７３）　における前記｛１０｝の解Ｂ（ＩＬｂ，Ｃｏｂ，ｎｂ，Ｒｓｈｂ，Ｒｓｂ）、前記温度ｔｃ（摂氏：Ｔｃ＝ｔｃ＋２７３）　における前記｛１２｝の解Ｃ（ＩＬｃ，Ｃｏｃ，ｎｃ，Ｒｓｈｃ，Ｒｓｃ）のそれぞれ（ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ）について曲線補間で、実測温度ｔｊ（摂氏：Ｔｊ＝ｔｊ＋２７３）　での特性値Ｍ（ＩＬｍ，Ｃｏｍ，ｎｍ，Ｒｓｈｍ，Ｒｓｍ）　を算出し、つぎに、
｛１５｝Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）
＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／（ｎＫｏＴ））−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ−Ｉ
を零とした式において、ＩＬのかわりにＩＬ＊Ｅ／１およびＶ＝０
とおき、Ｅ（ＫＷ／ｍ２）　について式を整理し、すなわち
Ｅ＝（１／ＩＬ）＊［Ｉ＋ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）｛ｅｘｐ（ｑＲｓＩ／ｎＫｏＴ）−１｝＋ＲｓＩ／Ｒｓｈ］
を作成し、
｛１６｝この式のＩ（電流）に短絡電流の測定値Ｉｓｊ，Ｔ（絶対温度）に太陽電池温度ｔｊの測定値を絶対温度に換算した値、およびＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，に｛１４｝で求めたＩＬｍ，Ｃｏｍ，ｎｍ，Ｒｓｈｍ，Ｒｓｍを代入して、日射強度Ｅを算出する
ことを特徴とする日射強度の測定方法。
【請求項４】
｛２０｝電圧Ｖ，電流Ｉ，日射強度１ＫＷ／ｍ２での光起電流ＩＬ，飽和電流温度係数Ｃｏ，接合定数ｎ，並列抵抗Ｒｓｈ　，直列抵抗Ｒｓ，太陽電池温度Ｔ（絶対温度）を含んだ
関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／（ｎＫｏＴ））−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ　を作成し、つぎに接合定数ｎ，並列抵抗Ｒｓｈ　，直列抵抗Ｒｓ，太陽電池温度Ｔ（絶対温度）を含んだ
関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／（ｎＫｏＴ））−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ　を作成し、つぎに
｛２１｝該関数Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）を変数Ｖで微分した関数：Ｄｉｆ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）を作成し、
｛２２｝太陽電池の基準状態（モジュール温度Ｔａ（２９８Ｋ（ｔａ＝２５℃））、日射強度Ｅａ（１ｋｗ／ｍ２））　での仕様値である、短絡電流Ｉｓｃａ、最適電流Ｉｏｐａ　−最適電圧Ｖｏｐａ　、開放電圧Ｖｏｃａの点Ｐ１（０，Ｉｓｃａ），Ｐ２（Ｖｏｐａ，Ｉｏｐａ），Ｐ３（Ｖｏｃａ，０）を選択し、
｛２３｝前記関数Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝　０のＴに基準状態の温度Ｔａ（２９８Ｋ），直列抵抗Ｒｓに基準温度での値Ｒｓａ，および前記Ｐ１，Ｐ２，Ｐ３の点の値を代入し、ＩＬ，Ｃｏ，ｎ，Ｒｓｈを未知数とする
関係式：Ｆｕｎｃ（０，Ｉｓｃａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，
関係式：Ｆｕｎｃ（Ｖｏｃａ，０，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，
関係式：Ｆｕｎｃ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０を作成し、
｛２４｝前記関数Ｄｉｆ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝　０　に、基準状態の温度Ｔａ（２９８Ｋ）、直列抵抗Ｒｓに基準温度での値Ｒｓａ　、および前記点Ｐ２の値（Ｖｏｐａ，Ｉｏｐａ）　を代入して、ＩＬ，Ｃｏ，ｎ，Ｒｓｈを未知数とする、
関係式：Ｄｉｆ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０を作成し、つぎに、
｛２５｝前記４つの関係式：Ｆｕｎｃ（０，Ｉｓｃａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，Ｆｕｎｃ（Ｖｏｃａ，０，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，Ｆｕｎｃ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，Ｄｉｆ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０　を満たす解Ａ（ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ）を、非線形解法のプログラムによって算出し、つぎに、
｛２６｝前記関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝　０のＩＬ，Ｃｏ，ｎ，Ｒｓｈに、前記解Ａ（ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ）　を代入して、さらにＴに基準状態の温度Ｔａ（２９８Ｋ）および直列抵抗Ｒｓに値Ｒｓａを代入し、変数Ｖ，Ｉの関係式：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ，Ｒｓａ，Ｔａ）＝　０を作成し、
｛２７｝この式を再び非線形解法のプログラムにより約４０〜５０点のＶに対するＩの解を求め、基準状態における電圧（Ｖ）−電流（Ｉ）および電圧（Ｖ）−電力（Ｐ）（電力（Ｐ）＝　電圧（Ｖ）×電流（Ｉ））　の関係の点（約４０〜５０点）を求め、これらの点を結んだ電圧Ｖ−電流Ｉカーブ（Ｉ−Ｖ　カーブ），電圧Ｖ−電力Ｐカーブ（Ｐ−Ｖ　カーブ）を作成し、つぎに、
｛２８｝前記太陽電池の日射強度Ｅｂ（ここでは１ｋｗ／ｍ２）、太陽電池温度Ｔｂ（絶対温度：Ｔｂ（Ｋ）＝ｔｂ（℃）＋２７３）　における電圧Ｖｂ，電流Ｉｂを求めるため、前記太陽電池の基準状態における日射強度Ｅａ（１ｋｗ／　ｍ２）、太陽電池温度ｔａ（２９８Ｋ：絶対温度：Ｔａ＝ｔａ＋２７３）における短絡電流Ｉｓｃａ，モジュール直列抵抗Ｒｓａ，温度が１℃変化したときの前記短絡電流Ｉｓｃａの変動値α、温度が１℃変化したときの前記開放電圧Ｖｏｃａ　の変動値β、曲線補正因子Ｋとする
変換式（Ｖａ，Ｉａ）→（Ｖｂ，Ｉｂ）：
Ｉｂ　＝　Ｉａ　＋　α＊（ｔｂ　−　ｔａ）
Ｖｂ　＝　Ｖａ　＋　β＊（ｔｂ　−　ｔａ）　−　Ｒｓａ＊（Ｉｂ　−　Ｉａ）　−　Ｋ＊Ｉｂ＊（ｔｂ　−　ｔａ）　を作成し、これに前記｛２７｝で作成した電圧−電流の各点またはそれを結んだＩ−Ｖ　カーブ上の各点を電流、電圧値すなわちＩａ，Ｖａ値として使用し、前式に適用し日射強度Ｅｂ（ここでは１ｋｗ／ｍ２）、太陽電池温度Ｔｂ（Ｋ）　における各点（電圧Ｖｂ−電流Ｉｂ：約４０〜５０点）　を求め、これらの点を結んだＩ−Ｖ　カーブ，Ｐ−Ｖ　カーブを作成し、つぎに、
｛２９｝前記｛２８｝で作成したＩ−Ｖ　カーブ上からＩｓｃｂ，Ｉｏｐｂ，Ｖｏｐｂ，Ｖｏｃｂを算出し、一方Ｒｓａをｂ℃に温度補正した値又は直接与えられた値であるＲｓｂを使い、そして前記｛２３｝〜｛２５｝の方法でＩｌｂ，Ｃｏｂ，ｎｂ，Ｒｓｈｂを算出し、つぎに
｛３０｝前記太陽電池の日射強度Ｅｃ（ここでは１ｋｗ／ｍ２）、太陽電池温度Ｔｃ（絶対温度：Ｔｃ（Ｋ）＝ｔｃ（℃）＋２７３）　における電圧Ｖｃ，電流Ｉｃを求めるため、前記太陽電池の基準状態における日射強度Ｅａ（１ｋｗ／　ｍ２）、太陽電池温度ｔａ（２９８Ｋ：絶対温度：Ｔａ＝ｔａ＋２７３）における短絡電流Ｉｓｃａ，モジュール直列抵抗Ｒｓａ，温度が１℃変化したときの前記短絡電流Ｉｓｃａの変動値α、温度が１℃変化したときの前記開放電圧Ｖｏｃａ　の変動値β、曲線補正因子Ｋとする
変換式（Ｖａ，Ｉａ）→（Ｖｃ，Ｉｃ）：
Ｉｃ　＝　Ｉａ　＋　α＊（ｔｃ　−　ｔａ）
Ｖｃ　＝　Ｖａ　＋　β＊（ｔｃ　−　ｔａ）　−　Ｒｓａ＊（Ｉｃ　−　Ｉａ）　−　Ｋ＊Ｉｃ＊（ｔｃ　−　ｔａ）　を作成し、これに前記｛２７｝で作成した電圧−電流の各点またはそれを結んだＩ−Ｖ　カーブ上の各点を電流、電圧値すなわちＩａ，Ｖａ値として使用し、前式に適用し日射強度Ｅｃ（ここでは１ｋｗ／ｍ２）、太陽電池温度Ｔｃ（Ｋ）　における各点（電圧Ｖｂ−電流Ｉｂ：約４０〜５０点）　を求め、これらの点を結んだＩ−Ｖ　カーブ，Ｐ−Ｖ　カーブを作成し、つぎに、
｛３１｝前記｛３０｝で作成したＩ−Ｖ　カーブのＩｓｃｃ，Ｉｏｐｃ，Ｖｏｐｃ，Ｖｏｃｃを算出し、一方Ｒｓａをｃ℃に温度補正した値又は直接与えられた値であるＲｓｃを使い、そして前記｛２３｝〜｛２５｝の方法でＩＬｃ，Ｃｏｃ，ｎｃ，Ｒｓｈｃを算出し、つぎに
｛３２｝太陽電池の実測温度ｔｊ（摂氏：絶対温度はＴｊ＝ｔｊ＋２７３）を取り込み、
｛３３｝基準状態での前記温度ｔａ（摂氏２５℃：絶対温度Ｔａ（Ｋ）（＝ｔａ（℃）＋２７３）における前記｛２５｝の解Ａ（ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ）とＲｓａ　、前記温度ｔｂ（摂氏：Ｔｂ＝ｔｂ＋２７３）　における前記｛２９｝の解Ｂ（ＩＬｂ，Ｃｏｂ，ｎｂ，Ｒｓｈｂ）とＲｓｂ、前記温度ｔｃ（摂氏：Ｔｃ＝ｔｃ＋２７３）　における前記｛３１｝の解Ｃ（ＩＬｃ，Ｃｏｃ，ｎｃ，Ｒｓｈｃ）とＲｓｃのそれぞれ（ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ）について曲線補間で、実測温度ｔｊ（摂氏：Ｔｊ＝ｔｊ＋２７３）　での特性値Ｍ（ＩＬｍ，Ｃｏｍ，ｎｍ，Ｒｓｈｍ，Ｒｓｍ）　を算出し、つぎに、
｛３４｝Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）
＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／（ｎＫｏＴ））−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ−Ｉ
を零とした式において、ＩＬのかわりにＩＬ＊Ｅ／１およびＶ＝０
とおき、Ｅ（ＫＷ／ｍ２）　について式を整理し、すなわち
Ｅ＝（１／ＩＬ）＊［Ｉ＋ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）｛ｅｘｐ（ｑＲｓＩ／ｎＫｏＴ）−１｝＋ＲｓＩ／Ｒｓｈ］
を作成し、
｛３５｝この式のＩ（電流）に短絡電流の測定値Ｉｓｊ，Ｔ（絶対温度）に太陽電池温度ｔｊの測定値を絶対温度に換算した値、およびＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，に｛３３｝で求めたＩＬｍ，Ｃｏｍ，ｎｍ，Ｒｓｈｍ，Ｒｓｍを代入して、日射強度Ｅを算出する
ことを特徴とする日射強度の測定方法。
【請求項５】
請求項３における前記｛１０｝、｛１２｝は作成したＩ−Ｖ　カーブ上から近接しない任意の５点を選択し、これらの値をＦｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝０の変数Ｖ、Ｉに代入してＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数５個の非線形連立方程式を作成して解いているが、これにかえ、作成したＩ−Ｖ　カーブ上の３点（（０，Ｉｓｃ），（Ｖｏｐ，Ｉｏｐ），（Ｖｏｃ，０）により前記｛０６｝のごとく、４つの関係式を作成し、これらの式からＲｓ，ｎの関係式を求め、次にこれらの点とは別にそのときのＩ−Ｖ　カーブ上の互いに近接しない任意の４点を選択し、これらの値をＦｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝０の変数Ｖ、Ｉに代入してＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数とする４つの式を作成し線形の要素であるＩＬ，Ｃｏ，Ｒｓｈ　を消去してＲｓ，ｎの関係式を求め、前記３点および４点により作成したＲｓ，ｎの関係式と連立させ、未知数２（Ｒｓ，ｎ）の非線形連立方程式をつくり、「連立のニュートン・ラプソン法（数値微分）」によりＲｓ，ｎを求め、そしてこの値を前記非線形の連立方程式に代入してＩＬ，Ｃｏ，Ｒｓｈ　を求める方法により、ｂ℃，ｃ℃におけるＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを算出し、前記請求項３の｛１３｝〜｛１６｝の方法で日射強度を算出することを特徴とする日射強度の測定方法。
【請求項６】
前記｛０６｝で作成した４つの非線形連立方程式に加え、｛０９｝で作成したＩ−Ｖカーブ上の、（Ｏ，Ｉｓｃ），（Ｖｏｐ，Ｉｏｐ），（Ｖｏｃ，Ｏ）とは異なり、これらの点に近接しない任意の点を選択し、該電圧−電流点を太陽電池特基本式に代入してもう１つの非線形方程式を作成し、これら５つの非線形連立方程式を「連立のニュートン・ラプソン法（数値微分）」により、ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数として算出し、次に｛０９｝〜｛１２｝によりｂ℃、Ｃ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを同様に算出し、前記｛１３｝〜｛１６｝の方法で日射強度を算出することを特徴とする日射強度の測定方法。
【請求項７】
前記請求項４の｛０６｝で作成した４つの非線形連立方程式に加え｛０９｝で作成したＩ−Ｖカーブ上の（Ｏ，Ｉｓｃ），（Ｖｏｐ，Ｉｏｐ），（Ｖｏｃ，Ｏ）とは異なり、これらの点に近接しない任意の点を選択し、該電圧−電流点を太陽電池基本式に代入してもう１つの非線形方程式を作成し、これら５つの非線形連立方程式を「連立のニュートン・ラプソン法（数値微分）」により、２５℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数として算出し、次に前記｛０９｝で作成したＩ−ＶカーブのＩｓｃｂ，Ｉｏｐｂ，Ｖｏｐｂ，Ｖｏｃｂを算出し、これらの点に近接しない任意の点を選択し、該電圧−電流点を太陽電池基本式に代入してもう１つの非線形方程式を作成し、これら５つの非線形連立方程式を「連立のニュートン・ラプソン法（数値微分）」により、ｂ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数として算出し、次にｃ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓについても前記｛１１｝で作成したＩ−Ｖカーブを使用して同様に算出し、前記｛１３｝〜｛１６｝の方法で日射強度を算出することを特徴とする日射強度の測定方法。
【請求項８】
前記請求項４の｛０６｝で作成した４つの非線形連立方程式において、ＩＬ＝ＩｓｃとしＣＯ，ｎ，Ｒｓｈ，Ｒｓを未知数として方程式を解いてＩＬ＝Ｉｓｃａ，Ｃｏａ，ｎａ，　Ｒｓｈａ，Ｒｓａを算出し、次にこれらの解と基準状態の温度Ｔａ（２８９Ｋ）を前記関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，ＣＯ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝０に代入し、変数Ｖ，Ｉの関係式を作成し、この式を再び非線形解法のプログラムによりＶに対するＩの解を４０〜５０ケース算出し、基準状態のＩ−Ｖカーブ、Ｐ−Ｖカーブを作成し、次に前記｛０９｝のように作成したＩ−ＶカーブのＩｓｃｂ，Ｉｏｐｂ，Ｖｏｐｂ，Ｖｏｃｂを算出し、前記のごとく３つの非線形方程式を作成し、さらにこれらの点に近接しない任意の点を選択し、該電圧−電流点を太陽電池基本式に代入してもう１つの非線形方程式を作成し、これら５つの非線形連立方程式を「連立のニュートン・ラプソン法（数値微分）」により、ｂ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数として算出し、次にｃ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓについても前記｛１１｝のように作成したＩ−Ｖカーブを使用して同様に算出し、前記｛１３｝〜｛１６｝の方法で日射強度を算出することを特徴とする日射強度の測定方法。
【請求項９】
前記請求項４の｛０６｝で作成した４つの非線形連立方程式に加え、与えられた基準状態の特性値（Ｉｓｃａ，Ｉｏｐａ，Ｖｏｐａ，Ｖｏｃａ）の（Ｏ，Ｉｓｃａ），（Ｖｏｐａ，Ｉｏｐａ），（Ｖｏｃａ，Ｏ）の３点を太陽電池特性にもとづく曲線により結び、Ｉ−Ｖカーブを作成し、このＩ−Ｖカーブ上の（Ｏ，Ｉｓｃａ），（Ｖｏｐ，Ｉｏｐａ），（Ｖｏｃａ，Ｏ）とは異なり、これらの点に近接しない任意の点を選択し、もう１つの非線形方程式を作成し、該５つの関係式ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数として解き、この解ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ，Ｒｓａおよび基準状態の温度Ｔａ（２８９Ｋ）を前記関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，ＣＯ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝０に代入して変数Ｖ，Ｉの関係式を作成し、この式を再び非線形解法のプログラムによりＶに対するＩの解を４０〜５０ケース算出し、基準状態のＩ−Ｖカーブ、Ｐ−Ｖカーブを作成し、次に前記｛０９｝で作成したＩ−ＶカーブのＩｓｃｂ，Ｉｏｐｂ，Ｖｏｐｂ，Ｖｏｃｂを算出し、前記のごとく３つの非線形方程式を作成し、さらにこれらの点に近接しない任意の点を選択し、該電圧−電流点を太陽電池基本式に代入してもう１つの非線形方程式を作成し、これら５つの非線形連立方程式を「連立のニュートン・ラプソン法（数値微分）」により、ｂ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数として算出し、次にｃ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓについても前記｛１１｝で作成したＩ−Ｖカーブを使用して同様に算出し、前記｛１３｝〜｛１６｝の方法で日射強度を算出することを特徴とする日射強度の測定方法。
【請求項１０】
前記｛０９｝で作成したＩ−Ｖカーブ上から近接しない任意の５点を選択し、これらの値をＦｕｎｃ（Ｖ，Ｉ，ＩＬ，ＣＯ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝０の変数Ｖ，Ｉに代入してＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数５個の非線形連立方程式を作成して２５℃の基本特性値を算出し、前記｛０９｝、｛１０｝によりｂ℃における｛１１｝、｛１２｝によってｃ℃における基本特性値ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを算出し、次に太陽電池の実測温度ｔｊ（摂氏：絶対温度はＴｊ＝ｔｊ＋２７３）を取り込み、実測温度ｔｊでの基本特性値Ｍ（Ｉｌｍ，Ｃｏｍ，ｎｍ，Ｒｓｈｍ，Ｒｓｍ）を算出し、前記｛１５｝、｛１６｝の方法で日射強度を算出することを特徴とする日射強度の測定方法。
【請求項１１】
太陽電池の実測された短絡電流Ｉｓｃｊ、太陽電池の実測された温度ｔｊを取り込むための実測値取り込みプログラムと、定数（Ｅｇ，Ｋｏ，ｑ）の入力プログラムと、太陽電池の基準状態での仕様値（Ｉｓｃ，Ｖｏｃ，Ｉｏｐ，Ｖｏｐ，α，β，Ｋ，Ｒｓ）を入力するプログラムと、請求項１，２，３，４，５，６，７，８，９，１０記載の日射強度Ｅを算出する処理プログラムとからなることを特徴とする日射強度算出プログラムを記録したコンピュータ読み取り可能なデータ記憶媒体。
【請求項１２】
前記太陽電池の短絡電流を実測する電流計と、前記太陽電池の温度を実測する温度計（熱電対）と、請求項１，２，３，４，５，６，７，８，９，１０記載の日射強度を算出する処理プログラムと請求項８の処理プログラムを記録したコンピュータ読み取り可能な記憶媒体が動作しうるコンピュータ、マイコンとからなることを特徴とする太陽電池の日射計。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
太陽電池は図４のごとく、太陽の光エネルギーを電気エネルギーに直接変換するものである。すなわち光電効果の一種である光起電力効果を応用し、起電力を発生させるものである。太陽電池中にエネルギーをもった光（光子）が入射すると、自由な電子と正孔が発生する。太陽電池半導体中のｐＭ接合近傍に達した電子と正孔は、それぞれｎ型半導体側、ｐ型半導体側に拡散し、両電極部に集まるので、電力が取り出せ、電圧および電流が発生するというわけである。
本発明は、前記のごとく、太陽の日射エネルギーにより発生する電流（ここでは短絡電流）の大きさ、太陽電池のモジュール温度および太陽電池の特性値（Ｉｓｃ，Ｉｏｐ，Ｖｏｐ，Ｖｏｃ，α，β，Ｒｓ，Ｋなど）を使い、太陽電池に入射する日射強度を光起電流を使い正確に算出する方法および本方法を使った日射計に関する。
【０００２】
【従来の技術】
太陽の日射強度を精密に測定する装置（精密日射計：気象管署で広く使われている日射計）としては熱電対を基本にしたサーモパイルで構成されている場合が多い（本発明ではこの方式の日射計を「全天日射計」と呼ぶ）。したがって原理上高価な装置であるとともに入射日射強度の大きさと出力電圧（出力信号）との間には１０〜３０秒程度の時間差が生じる。このため太陽電池の出力と日射計の出力が直接対応しにくく、また出力電圧には過渡的な変動部分も含まれている。なお、日射強度と太陽電池出力の時間的に厳密な対応が必要ない、通常の日射強度の測定においては、短時間における日射強度の変動はあまり意識する必要はなく、平均値的に扱った値でもほとんど問題はなかった。そのため上記日射計が一般に使用されていた。
【０００３】
一方、太陽電池素子を使い、その短絡電流の大きさにより日射強度を測定する方式の日射計も最近はよく見かける。この太陽電池素子を使った方式の日射計においては、一般に短絡電流の大きさと日射強度の関係が直線関係であるという前提のもとで日射強度を測定しているため、太陽電池素子温度や日射スペクトルを通常考慮しないので、その測定した精度に問題があった。すなわち、その適用している日射計の太陽電池の種類、諸特性、スペクトル特性などの差も考慮した精度と汎用性の高い日射計が望まれているのが現状である。
【０００４】
そして、最近太陽電池の出力（発電電力）を評価するために「Ｉ−Ｖカーブ測定器」が市販されている。この場合の日射測定においても前記２種類の日射計が使われているが、日射変動が激しいときには、データと処理結果にばらつきが多いため、測定精度上に問題があったり、それぞれの出力評価が十分できにくかった場合も多く見られる。
【０００５】
ここで、本発明のベースである太陽電池基本式と日射強度について簡単に説明を加える。一般に太陽電池基本式は下記のとおりである。
Ｉ＝ＩＬ−ＩＯ＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／ｎＫｏＴ）−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ　……（１）
ＩＯ＝ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）　　　　　　　　　……（２）
ここに各記号は次のとおりである。
Ｉ　：出力電流［Ａ］　　　　　Ｃｏ　：飽和電流温度係数
Ｖ　：出力電圧［Ｖ］　　　　　　　　　　Ｅｇ　：エネルギーギャップ［ｅＶ］
ＩＬ　：光起電流［Ａ］　　　　　　　　　　Ｔ　：太陽電池素子絶対温度［Ｋ］
ＩＯ　：飽和電流［Ａ］　　　　　　　　　　Ｋｏ　：ボルツマン定数［Ｊ／Ｋ］
Ｒｓ　：直流抵抗［Ω］　　　　　　　　　　ｑ　：電子の電荷量［Ｃ］
Ｒｓｈ　：並列抵抗［Ω］
ｎ　：接合定数
太陽電池の出力である電圧と電流は（１）　，（２）　式であらわされる。
そして、ＩＬ（光起電流）は太陽電池温度と日射強度によって変化するため、測定時の短絡電流と算出したＩＬとの比較により日射強度を算出している。
【０００６】
【発明が解決しようとする課題】
前記のような従来の技術に鑑み、本発明が解決しようとする課題を述べると次のとおりである。
第１に、従来から広く使われているサーモパイルを基本とした前記の全天日射計では測定精度については一般的に信頼されているものの、日射強度に対して出力信号は前記のとおり１０〜３０秒程度の時間遅れが生じ、したがって、日射強度の変動による出力信号の過渡的な変動も生じる。そして、時間遅れなどが許されない装置（例えば、前記太陽電池の電圧−電流カーブ（Ｉ−Ｖカーフ゛）　を測定し太陽電池の評価をする装置など…特願２００１−０７２０８１など参照）においては、根本的に適用上の問題点があった。また、この原理に基づく高精度の日射計は一般に高価（１５〜５０万円程度）でもあった。
そこで、太陽電池を使った日射計について検討することとした。前記のように従来にも太陽電池を使った日射計として前記のようにその短絡電流の大きさから日射強度を算出するものがあった。しかし、この方式のものでは太陽電池の温度による誤差への影響が大きいとともに、太陽電池の短絡電流と日射強度を直接的なものとしたため、太陽光のスペクトル分布の違いおよびその特定の太陽電池のスペクトル感度の違い、すなわちそれぞれの太陽電池の特性によって測定値が異なることが生じていた。
【０００７】
そこで、温度などの影響を少なくするため、温度による補正をかける方法を提案している例も一部に見られる
（電気学会論文２（伊賀：「実用的Ｉ−Ｖ　カーフ゛作成法を使った太陽電池日射計」，電学論Ｄ，１１７　巻１０号，（１９９７）））。
本発明ではより精度の高い日射強度を得るのに、太陽電池の短絡電流でなく、太陽電池の”光起電流”を使い、日射強度との間の関係を利用することに基づき、前記課題を解決しようとするものである。
すなわち、本発明は「光起電流は入射光のスペクトル分布の形が一定な（類似する）限り入射光の放射照度に一般に正比例する」ことを利用して、太陽電池により前記全天日射計で測定する日射強度を得ようとするものである。
そのため、太陽電池の基本式（次式）において、測定した短絡電流と太陽電池温度および測定時の太陽電池温度における特性値（ＩＬｍ，Ｃｏｍ，ｎｍ，Ｒｓｈｍ，Ｒｓｍ）　から太陽電池
基本式
Ｉ＝ＩＬ−ＩＯ＊｛ｅｘｐｑ（Ｖ＋ＲｓＩ）／ｎＫｏＴ｝−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ
ＩＯ＝ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）
を使い、
光起電流を求め、そして日射強度を求めるものである。
【０００８】
【課題を解決するための手段】
請求項１の光起電流による日射強度の測定方法は、
太陽電池基本式
Ｉ＝ＩＬ−ＩＯ＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／ｎＫｏＴ）−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ
ＩＯ＝ＣｏＴ３　ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）
ここに、
Ｉ　：出力電流［Ａ］　　　　　Ｃｏ　：飽和電流温度係数
Ｖ　：出力電圧［Ｖ］　　　　　　　　　　Ｅｇ　：エネルギーギャップ［ｅＶ］
ＩＬ　：光起電流［Ａ］　　　　　　　　　　Ｔ　：太陽電池（素子）絶対温度［Ｋ］
ＩＯ　：飽和電流［Ａ］　　　　　　　　　　Ｋｏ　：ボルツマン定数［Ｊ／Ｋ］
Ｒｓ　：直流抵抗［Ω］　　　　　　　　　　ｑ　：電子の電荷量［Ｃ］
Ｒｓｈ　：並列抵抗［Ω］
ｎ　：接合定数
における基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓ，Ｒｓｈ）のそれぞれについてａ℃，ｂ℃，ｃ℃の値を求め、それぞれのａ℃，ｂ℃，ｃ℃での値を使い温度の曲線補間をすることにより太陽電池の測定温度におけるそれぞれの基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓ，Ｒｓｈ）を求め、
これらの値に加え、測定時の短絡電流Ｉｓｊ（Ａ）　および太陽電池温度ｔｊ（℃）（絶対温度Ｔ（Ｋ）＝ｔｊ（℃）＋２７３）を日射強度Ｅ（ＫＷ／ｍ２）　の算出式
Ｅ＝（１／ＩＬ）＊［Ｉｓｊ＋Ｃｏ・Ｔ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑＲｓＩｓｊ／ｎＫｏＴ）−１｝＋Ｒｓ＊Ｉｓｊ／Ｒｓｈ］
に適用して日射強度Ｅを求めることを特徴とする。
【０００９】
請求項２の光起電流による日射強度の測定方法は、
Ｒｓおよびその温度特性が与えられている場合、
太陽電池基本式
Ｉ＝ＩＬ−ＩＯ＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／ｎＫｏＴ）−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ
ＩＯ＝ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）
ここに、
Ｉ　：出力電流［Ａ］　　　　　Ｃｏ　：飽和電流温度係数
Ｖ　：出力電圧［Ｖ］　　　　　　　　　　Ｅｇ　：エネルギーギャップ［ｅＶ］
ＩＬ　：光起電流［Ａ］　　　　　　　　　　Ｔ　：太陽電池（素子）絶対温度［Ｋ］
ＩＯ　：飽和電流［Ａ］　　　　　　　　　　Ｋｏ　：ボルツマン定数［Ｊ／Ｋ］
Ｒｓ　：直流抵抗［Ω］　　　　　　　　　　ｑ　：電子の電荷量［Ｃ］
Ｒｓｈ　：並列抵抗［Ω］
ｎ　：接合定数
において、基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓｈ）のそれぞれについてａ℃，ｂ℃，ｃ℃の値を求め、それぞれのａ℃，ｂ℃，ｃ℃での値を使い温度の曲線補間をすることにより太陽電池の測定温度におけるそれぞれの基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓ，Ｒｓｈ）を求め、
これらの値とａ℃，ｂ℃，ｃ℃におけるＲｓの値に加え、測定時の短絡電流Ｉｓｊ（Ａ）および太陽電池温度ｔｊ（℃）（絶対温度Ｔ（Ｋ）＝ｔｊ（℃）＋２７３）を日射強度Ｅ（ＫＷ／ｍ２）の算出式
Ｅ＝（１／ＩＬ）＊［Ｉｓｊ＋Ｃｏ・Ｔ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑＲｓＩｓｊ／ｎＫｏＴ）−１｝＋Ｒｓ＊Ｉｓｊ／Ｒｓｈ］
に適用して日射強度Ｅを求めることを特徴とする。
【００１０】
請求項３の日射強度の測定方法は、
｛０１｝電圧Ｖ，電流Ｉ，日射強度１ＫＷ／ｍ２での光起電流ＩＬ，飽和電流温度係数Ｃｏ，接合定数ｎ，並列抵抗Ｒｓｈ　，直列抵抗Ｒｓ，太陽電池温度Ｔ（絶対温度）を含んだ
関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／（ｎＫｏＴ））−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ　を作成し、つぎに
接合定数ｎ，並列抵抗Ｒｓｈ　，直列抵抗Ｒｓ，太陽電池温度Ｔ（絶対温度）を含んだ
関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／（ｎＫｏＴ））−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ　を作成し、つぎに
｛０２｝該関数Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）を変数Ｖで微分した関数：Ｄｉｆ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）を作成し、
｛０３｝太陽電池の基準状態（モジュール温度Ｔａ（２９８Ｋ（ｔａ＝２５℃））、日射強度Ｅａ（１ｋｗ／ｍ２））　での仕様値である、短絡電流Ｉｓｃａ、最適電流Ｉｏｐａ　−最適電圧Ｖｏｐａ　、開放電圧Ｖｏｃａの点Ｐ１（０，Ｉｓｃａ），Ｐ２（Ｖｏｐａ，Ｉｏｐａ），Ｐ３（Ｖｏｃａ，０）を選択し、
｛０４｝前記関数Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝　０のＴに基準状態の温度Ｔａ（２９８Ｋ），直列抵抗Ｒｓに基準温度での値Ｒｓａ，および前記Ｐ１，Ｐ２，Ｐ３の点の値を代入し、ＩＬ，Ｃｏ，ｎ，Ｒｓｈを未知数とする
関係式：Ｆｕｎｃ（０，Ｉｓｃａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，
関係式：Ｆｕｎｃ（Ｖｏｃａ，０，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，
関係式：Ｆｕｎｃ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０を作成し、
｛０５｝前記関数Ｄｉｆ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝　０　に、基準状態の温度Ｔａ（２９８Ｋ）、直列抵抗Ｒｓに基準温度での値Ｒｓａ　、および前記点Ｐ２の値（Ｖｏｐａ，Ｉｏｐａ）　を代入して、ＩＬ，Ｃｏ，ｎ，Ｒｓｈを未知数とする、
関係式：Ｄｉｆ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０を作成し、つぎに、
｛０６｝前記４つの関係式：Ｆｕｎｃ（０，Ｉｓｃａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，Ｆｕｎｃ（Ｖｏｃａ，０，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，Ｆｕｎｃ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，Ｄｉｆ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０　を満たす解Ａ（ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ）を、非線形解法のプログラムによって算出し、つぎに、
｛０７｝前記関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝　０のＩＬ，Ｃｏ，ｎ，Ｒｓｈに、前記解Ａ（ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ）　を代入して、さらにＴに基準状態の温度Ｔａ（２９８Ｋ）および直列抵抗Ｒｓに値Ｒｓａを代入し、変数Ｖ，Ｉの関係式：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ，Ｒｓａ，Ｔａ）＝　０を作成し、
｛０８｝この式を再び非線形解法のプログラムにより約４０〜５０点のＶに対するＩの解を求め、基準状態における電圧（Ｖ）−電流（Ｉ）および電圧（Ｖ）−電力（Ｐ）（電力（Ｐ）＝　電圧（Ｖ）×電流（Ｉ））　の関係の点（約４０〜５０点）を求め、これらの点を結んだ電圧Ｖ−電流Ｉカーブ（Ｉ−Ｖ　カーブ），電圧Ｖ−電力Ｐカーブ（Ｐ−Ｖ　カーブ）を作成し、つぎに、
｛０９｝前記太陽電池の日射強度Ｅｂ（ここでは１ｋｗ／ｍ２）、太陽電池モジュール温度Ｔｂ（絶対温度：Ｔｂ（Ｋ）＝ｔｂ（℃）＋２７３）　における電圧Ｖｂ，電流Ｉｂを求めるため、前記太陽電池の基準状態における日射強度Ｅａ（１ｋｗ／　ｍ２）、太陽電池温度ｔａ（２９８Ｋ：絶対温度：Ｔａ＝ｔａ＋２７３）における短絡電流Ｉｓｃａ，モジュール直列抵抗Ｒｓａ，温度が１℃変化したときの前記短絡電流Ｉｓｃａの変動値α、温度が１℃変化したときの前記開放電圧Ｖｏｃａ　の変動値β、曲線補正因子Ｋとする
変換式（Ｖａ，Ｉａ）→（Ｖｂ，Ｉｂ）：
Ｉｂ　＝　Ｉａ　＋　α＊（ｔｂ　−　ｔａ）
Ｖｂ　＝　Ｖａ　＋　β＊（ｔｂ　−　ｔａ）　−　Ｒｓａ＊（Ｉｂ　−　Ｉａ）　−　Ｋ＊Ｉｂ＊（ｔｂ　−　ｔａ）　を作成し、これに前記｛０８｝で作成した電圧−電流の各点またはそれを結んだＩ−Ｖ　カーブ上の各点を電流、電圧値すなわちＩａ，Ｖａ値として使用し、前式に適用し日射強度Ｅｂ（ここでは１ｋｗ／ｍ２）、モジュール温度Ｔｂ（Ｋ）　における各点（電圧Ｖｂ−電流Ｉｂ：約４０〜５０点）　を求め、これらの点を結んだＩ−Ｖ　カーブ，Ｐ−Ｖ　カーブを作成し、つぎに、
｛１０｝前記｛０９｝で作成したＩ−Ｖ　カーブ上から、近接しない任意の５点を選択し、これらの点（ＶＱ１，ＩＱ１），（ＶＱ２，ＩＱ２），（ＶＱ３，ＩＱ３），（ＶＱ４，ＩＱ４），（ＶＱ５，ＩＱ５）　を前記関係式：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｂ）＝　０の変数Ｖ，Ｉ、に代入して、ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数とする、
関係式：Ｆｕｎｃ（ＶＱ１，ＩＱ１，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｂ）＝　０，
関係式：Ｆｕｎｃ（ＶＱ２，ＩＱ２，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｂ）＝　０，
関係式：Ｆｕｎｃ（ＶＱ３，ＩＱ３，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｂ）＝　０，
関係式：Ｆｕｎｃ（ＶＱ４，ＩＱ４，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｂ）＝　０，
関係式：Ｆｕｎｃ（ＶＱ５，ＩＱ５，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｂ）＝　０を作成し、該５つの関係式の
解Ｂ（ＩＬｂ，Ｃｏｂ，ｎｂ，Ｒｓｈｂ，Ｒｓｂ）を、非線形解法のプログラムによって算出し、つぎに、
｛１１｝前記太陽電池の日射強度Ｅｃ（ここでは１ｋｗ／ｍ２），太陽電池温度ｔｃ　（絶対温度：Ｔｃ（Ｋ）＝ｔｃ（℃）＋２７３）における電圧Ｖｃ，電流Ｉｃの関係についても前記｛０９｝と同様に、
変換式（Ｖａ，Ｉａ）→（Ｖｃ，Ｉｃ）：
Ｉｃ　＝　Ｉａ　＋　α＊（ｔｃ　−　ｔａ）
Ｖｃ　＝　Ｖａ　＋　β＊（ｔｃ　−　ｔａ）　−　Ｒｓａ＊（Ｉｃ　−　Ｉａ）　−　Ｋ＊Ｉｃ＊（ｔｃ　−　ｔａ）を使い　Ｉ−Ｖ　カーブ，Ｐ−Ｖ　カーブを作成し、
｛１２｝前記｛１１｝で作成したＩ−Ｖ　カーブ上から近接しない任意の５点を選択し、これらの値（ＶＲ１，ＩＲ１），（ＶＲ２，ＩＲ２），（ＶＲ３，ＩＲ３），（ＶＲ４，ＩＲ４），（ＶＲ５，ＩＲ５）　を前記関係式：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｃ）＝　０の変数Ｖ，Ｉに代入して、ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数とする、
関係式：Ｆｕｎｃ（ＶＲ１，ＩＲ１，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｃ）＝　０，
関係式：Ｆｕｎｃ（ＶＲ２，ＩＲ２，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｃ）＝　０，
関係式：Ｆｕｎｃ（ＶＲ３，ＩＲ３，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｃ）＝　０，
関係式：Ｆｕｎｃ（ＶＲ４，ＩＲ４，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｃ）＝　０，
関係式：Ｆｕｎｃ（ＶＲ５，ＩＲ５，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔｃ）＝　０
を作成し、該５つの関係式の解Ｃ（ＩＬｃ，Ｃｏｃ，ｎｃ，Ｒｓｈｃ，Ｒｓｃ）を、非線形解法のプログラムによって、算出し、つぎに、
｛１３｝太陽電池の実測温度ｔｊ（摂氏：絶対温度はＴｊ＝ｔｊ＋２７３）を取り込み、
｛１４｝基準状態での前記温度ｔａ（摂氏２５℃：絶対温度Ｔａ（Ｋ）（＝ｔａ（℃）＋２７３）における前記｛０６｝の解Ａ（ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ）とＲｓａ　、前記温度ｔｂ（摂氏：Ｔｂ＝ｔｂ＋２７３）　における前記｛１０｝の解Ｂ（ＩＬｂ，Ｃｏｂ，ｎｂ，Ｒｓｈｂ，Ｒｓｂ）、前記温度ｔｃ（摂氏：Ｔｃ＝ｔｃ＋２７３）　における前記｛１２｝の解Ｃ（ＩＬｃ，Ｃｏｃ，ｎｃ，Ｒｓｈｃ，Ｒｓｃ）のそれぞれ（ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ）について曲線補間で、実測温度ｔｊ（摂氏：Ｔｊ＝ｔｊ＋２７３）　での特性値Ｍ（ＩＬｍ，Ｃｏｍ，ｎｍ，Ｒｓｈｍ，Ｒｓｍ）　を算出し、つぎに、
｛１５｝Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）
＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／（ｎＫｏＴ））−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ−Ｉ
を零とした式において、ＩＬのかわりにＩＬ＊Ｅ／１およびＶ＝０
とおき、Ｅ（ＫＷ／ｍ２）　について式を整理し、すなわち
Ｅ＝（１／ＩＬ）＊［Ｉ＋ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）｛ｅｘｐ（ｑＲｓＩ／ｎＫｏＴ）−１｝＋ＲｓＩ／Ｒｓｈ］
を作成し、
｛１６｝この式のＩ（電流）に短絡電流の測定値Ｉｓｊ，Ｔ（絶対温度）に太陽電池温度ｔｊの測定値を絶対温度に換算した値、およびＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，に｛１４｝で求めたＩＬｍ，Ｃｏｍ，ｎｍ，Ｒｓｈｍ，Ｒｓｍを代入して、日射強度Ｅを算出することを特徴とする。
【００１１】
請求項４の日射強度の測定方法は、
｛２０｝電圧Ｖ，電流Ｉ，日射強度１ＫＷ／ｍ２での光起電流ＩＬ，飽和電流温度係数Ｃｏ，接合定数ｎ，並列抵抗Ｒｓｈ　，直列抵抗Ｒｓ，太陽電池温度Ｔ（絶対温度）を含んだ
関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／（ｎＫｏＴ））−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ　を作成し、つぎに
接合定数ｎ，並列抵抗Ｒｓｈ　，直列抵抗Ｒｓ，太陽電池温度Ｔ（絶対温度）を含んだ
関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／（ｎＫｏＴ））−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ　を作成し、つぎに
｛２１｝該関数Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）を変数Ｖで微分した関数：Ｄｉｆ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）を作成し、
｛２２｝太陽電池の基準状態（モジュール温度Ｔａ（２９８Ｋ（ｔａ＝２５℃））、日射強度Ｅａ（１ｋｗ／ｍ２））　での仕様値である、短絡電流Ｉｓｃａ、最適電流Ｉｏｐａ　−最適電圧Ｖｏｐａ　、開放電圧Ｖｏｃａの点Ｐ１（０，Ｉｓｃａ），Ｐ２（Ｖｏｐａ，Ｉｏｐａ），Ｐ３（Ｖｏｃａ，０）を選択し、
｛２３｝前記関数Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝　０のＴに基準状態の温度Ｔａ（２９８Ｋ），直列抵抗Ｒｓに基準温度での値Ｒｓａ，および前記Ｐ１，Ｐ２，Ｐ３の点の値を代入し、ＩＬ，Ｃｏ，ｎ，Ｒｓｈを未知数とする
関係式：Ｆｕｎｃ（０，Ｉｓｃａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，
関係式：Ｆｕｎｃ（Ｖｏｃａ，０，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，
関係式：Ｆｕｎｃ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０を作成し、
｛２４｝前記関数Ｄｉｆ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝　０　に、基準状態の温度Ｔａ（２９８Ｋ）、直列抵抗Ｒｓに基準温度での値Ｒｓａ　、および前記点Ｐ２の値（Ｖｏｐａ，Ｉｏｐａ）　を代入して、ＩＬ，Ｃｏ，ｎ，Ｒｓｈを未知数とする、
関係式：Ｄｉｆ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０を作成し、つぎに、
｛２５｝Ｒｓａが与えられていないため、前記４つの関係式：Ｆｕｎｃ（０，Ｉｓｃａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，Ｆｕｎｃ（Ｖｏｃａ，０，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，Ｆｕｎｃ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０，Ｄｉｆ（Ｖｏｐａ，Ｉｏｐａ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓａ，Ｔａ）＝　０　を満たす解Ａ（ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ）を、非線形解法のプログラムによって算出し、つぎに、
｛２６｝前記関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝　０のＩＬ，Ｃｏ，ｎ，Ｒｓｈに、前記解Ａ（ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ）　を代入して、さらにＴに基準状態の温度Ｔａ（２９８Ｋ）および直列抵抗Ｒｓに値Ｒｓａを代入し、変数Ｖ，Ｉの関係式：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ，Ｒｓａ，Ｔａ）＝　０を作成し、
｛２７｝この式を再び非線形解法のプログラムにより約４０〜５０点のＶに対するＩの解を求め、基準状態における電圧（Ｖ）−電流（Ｉ）および電圧（Ｖ）−電力（Ｐ）（電力（Ｐ）＝　電圧（Ｖ）×電流（Ｉ））　の関係の点（約４０〜５０点）を求め、これらの点を結んだ電圧Ｖ−電流Ｉカーブ（Ｉ−Ｖ　カーブ），電圧Ｖ−電力Ｐカーブ（Ｐ−Ｖ　カーブ）を作成し、つぎに、
｛２８｝前記太陽電池の日射強度Ｅｂ（ここでは１ｋｗ／ｍ２）、太陽電池温度Ｔｂ（絶対温度：Ｔｂ（Ｋ）＝ｔｂ（℃）＋２７３）　における電圧Ｖｂ，電流Ｉｂを求めるため、前記太陽電池の基準状態における日射強度Ｅａ（１ｋｗ／　ｍ２）、太陽電池温度ｔａ（２９８Ｋ：絶対温度：Ｔａ＝ｔａ＋２７３）における短絡電流Ｉｓｃａ，モジュール直列抵抗Ｒｓａ，温度が１℃変化したときの前記短絡電流Ｉｓｃａの変動値α、温度が１℃変化したときの前記開放電圧Ｖｏｃａ　の変動値β、曲線補正因子Ｋとする
変換式（Ｖａ，Ｉａ）→（Ｖｂ，Ｉｂ）：
Ｉｂ　＝　Ｉａ　＋　α＊（ｔｂ　−　ｔａ）
Ｖｂ　＝　Ｖａ　＋　β＊（ｔｂ　−　ｔａ）　−　Ｒｓａ＊（Ｉｂ　−　Ｉａ）　−　Ｋ＊Ｉｂ＊（ｔｂ　−　ｔａ）　を作成し、これに前記｛２７｝で作成した電圧−電流の各点またはそれを結んだＩ−Ｖ　カーブ上の各点を電流、電圧値すなわちＩａ，Ｖａ値として使用し、前式に適用し日射強度Ｅｂ（ここでは１ｋｗ／ｍ２）、太陽電池温度Ｔｂ（Ｋ）　における各点（電圧Ｖｂ−電流Ｉｂ：約４０〜５０点）　を求め、これらの点を結んだＩ−Ｖ　カーブ，Ｐ−Ｖ　カーブを作成し、つぎに、
｛２９｝前記｛２８｝で作成したＩ−Ｖ　カーブ上からＩｓｃｂ，Ｉｏｐｂ，Ｖｏｐｂ，Ｖｏｃｂを算出し、一方Ｒｓａをｂ℃に温度補正した値又は直接与えられた値であるＲｓｂを使い、そして前記｛２３｝〜｛２５｝の方法でＩｌｂ，Ｃｏｂ，ｎｂ，Ｒｓｈｂを算出し、つぎに
｛３０｝前記太陽電池の日射強度Ｅｃ（ここでは１ｋｗ／ｍ２）、太陽電池温度Ｔｃ（絶対温度：Ｔｃ（Ｋ）＝ｔｃ（℃）＋２７３）　における電圧Ｖｃ，電流Ｉｃを求めるため、前記太陽電池の基準状態における日射強度Ｅａ（１ｋｗ／　ｍ２）、太陽電池温度ｔａ（２９８Ｋ：絶対温度：Ｔａ＝ｔａ＋２７３）における短絡電流Ｉｓｃａ，モジュール直列抵抗Ｒｓａ，温度が１℃変化したときの前記短絡電流Ｉｓｃａの変動値α、温度が１℃変化したときの前記開放電圧Ｖｏｃａ　の変動値β、曲線補正因子Ｋとする
変換式（Ｖａ，Ｉａ）→（Ｖｃ，Ｉｃ）：
Ｉｃ　＝　Ｉａ　＋　α＊（ｔｃ　−　ｔａ）
Ｖｃ　＝　Ｖａ　＋　β＊（ｔｃ　−　ｔａ）　−　Ｒｓａ＊（Ｉｃ　−　Ｉａ）　−　Ｋ＊Ｉｃ＊（ｔｃ　−　ｔａ）　を作成し、これに前記｛２７｝で作成した電圧−電流の各点またはそれを結んだＩ−Ｖ　カーブ上の各点を電流、電圧値すなわちＩａ，Ｖａ値として使用し、前式に適用し日射強度Ｅｃ（ここでは１ｋｗ／ｍ２）、太陽電池温度Ｔｃ（Ｋ）　における各点（電圧Ｖｂ−電流Ｉｂ：約４０〜５０点）　を求め、これらの点を結んだＩ−Ｖ　カーブ，Ｐ−Ｖ　カーブを作成し、つぎに、
｛３１｝前記｛３０｝で作成したＩ−Ｖ　カーブのＩｓｃｃ，Ｉｏｐｃ，Ｖｏｐｃ，Ｖｏｃｃを算出し、一方Ｒｓａをｃ℃に温度補正した値又は直接与えられた値であるＲｓｃを使い、そして前記｛２３｝〜｛２５｝の方法でＩＬｃ，Ｃｏｃ，ｎｃ，Ｒｓｈｃを算出し、つぎに
｛３２｝太陽電池の実測温度ｔｊ（摂氏：絶対温度はＴｊ＝ｔｊ＋２７３）を取り込み、
｛３３｝基準状態での前記温度ｔａ（摂氏２５℃：絶対温度Ｔａ（Ｋ）（＝ｔａ（℃）＋２７３）における前記｛２５｝の解Ａ（ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ）とＲｓａ　、前記温度ｔｂ（摂氏：Ｔｂ＝ｔｂ＋２７３）　における前記｛２９｝の解Ｂ（ＩＬｂ，Ｃｏｂ，ｎｂ，Ｒｓｈｂ）とＲｓｂ、前記温度ｔｃ（摂氏：Ｔｃ＝ｔｃ＋２７３）　における前記｛３１｝の解Ｃ（ＩＬｃ，Ｃｏｃ，ｎｃ，Ｒｓｈｃ）とＲｓｃのそれぞれ（ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ）について曲線補間で、実測温度ｔｊ（摂氏：Ｔｊ＝ｔｊ＋２７３）　での特性値Ｍ（ＩＬｍ，Ｃｏｍ，ｎｍ，Ｒｓｈｍ，Ｒｓｍ）　を算出し、つぎに、
｛３４｝Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）
＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／（ｎＫｏＴ））−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ−Ｉ
を零とした式において、ＩＬのかわりにＩＬ＊Ｅ／１およびＶ＝０
とおき、Ｅ（ＫＷ／ｍ２）　について式を整理し、すなわち
Ｅ＝（１／ＩＬ）＊［Ｉ＋ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）｛ｅｘｐ（ｑＲｓＩ／ｎＫｏＴ）−１｝＋ＲｓＩ／Ｒｓｈ］
を作成し、
｛３５｝この式のＩ（電流）に短絡電流の測定値Ｉｓｊ，Ｔ（絶対温度）に太陽電池温度ｔｊの測定値を絶対温度に換算した値、およびＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，に｛３３｝で求めたＩＬｍ，Ｃｏｍ，ｎｍ，Ｒｓｈｍ，Ｒｓｍを代入して、日射強度Ｅを算出することを特徴とする。
【００１２】
請求項５の日射強度の測定方法は、
請求項３における前記｛１０｝、｛１２｝は作成したＩ−Ｖ　カーブ上から近接しない任意の５点を選択し、これらの値をＦｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝０の変数Ｖ、Ｉに代入してＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数５個の非線形連立方程式を作成して解いているが、これにかえ、作成したＩ−Ｖ　カーブ上の３点（（０，Ｉｓｃ），（Ｖｏｐ，Ｉｏｐ），（Ｖｏｃ，０）により前記｛０６｝のごとく、４つの関係式を作成し、これらの式からＲｓ，ｎの関係式を求め、次にこれらの点とは別にそのときのＩ−Ｖ　カーブ上の互いに近接しない任意の４点を選択し、これらの値をＦｕｎｃ（Ｖ，Ｉ，ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝０の変数Ｖ、Ｉに代入してＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数とする４つの式を作成し線形の要素であるＩＬ，Ｃｏ，Ｒｓｈ　を消去してＲｓ，ｎの関係式を求め、前記３点および４点により作成したＲｓ，ｎの関係式と連立させ、未知数２（Ｒｓ，ｎ）の非線形連立方程式をつくり、「連立のニュートン・ラプソン法（数値微分）」によりＲｓ，ｎを求め、そしてこの値を前記非線形の連立方程式に代入してＩＬ，Ｃｏ，Ｒｓｈ　を求める方法により、ｂ℃，ｃ℃におけるＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを算出し、前記請求項３の｛１３｝〜｛１６｝の方法で日射強度を算出することを特徴とする。
【００１３】
請求項６の日射強度の測定方法は、
前記｛０６｝で作成した４つの非線形連立方程式に加え、｛０９｝で作成したＩ−Ｖカーブ上の、（Ｏ，Ｉｓｃ），（Ｖｏｐ，Ｉｏｐ），（Ｖｏｃ，Ｏ）とは異なり、これらの点に近接しない任意の点を選択し、該電圧−電流点を太陽電池特基本式に代入してもう１つの非線形方程式を作成し、これら５つの非線形連立方程式を「連立のニュートン・ラプソン法（数値微分）」により、ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数として算出し、次に｛０９｝〜｛１２｝によりｂ℃、Ｃ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを同様に算出し、前記｛１３｝〜｛１６｝の方法で日射強度を算出することを特徴とする。
【００１４】
請求項７の日射強度の測定方法は、
前記｛０６｝で作成した４つの非線形連立方程式に加え｛０９｝で作成したＩ−Ｖカーブ上の（Ｏ，Ｉｓｃ），（Ｖｏｐ，Ｉｏｐ），（Ｖｏｃ，Ｏ）とは異なり、これらの点に近接しない任意の点を選択し、該電圧−電流点を太陽電池基本式に代入してもう１つの非線形方程式を作成し、これら５つの非線形連立方程式を「連立のニュートン・ラプソン法（数値微分）」により、２５℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数として算出し、次に前記｛２８｝で作成したＩ−ＶカーブのＩｓｃｂ，Ｉｏｐｂ，Ｖｏｐｂ，Ｖｏｃｂを算出し、これらの点に近接しない任意の点を選択し、該電圧−電流点を太陽電池基本式に代入してもう１つの非線形方程式を作成し、これら５つの非線形連立方程式を「連立のニュートン・ラプソン法（数値微分）」により、ｂ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数として算出し、次にｃ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓについても前記｛３０｝で作成したＩ−Ｖカーブを使用して同様に算出し、前記｛１３｝〜｛１６｝の方法で日射強度を算出することを特徴とする。
【００１５】
請求項８の日射強度の測定方法は、
前記請求項４の｛０６｝で作成した４つの非線形連立方程式において、ＩＬ＝ＩｓｃとしＣＯ，ｎ，Ｒｓｈ，Ｒｓを未知数として方程式を解いてＩＬ＝Ｉｓｃａ，Ｃｏａ，ｎａ，　Ｒｓｈａ，Ｒｓａを算出し、次にこれらの解と基準状態の温度Ｔａ（２８９Ｋ）を前記関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，ＣＯ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝０に代入し、変数Ｖ，Ｉの関係式を作成し、この式を再び非線形解法のプログラムによるＶに対するＩの解を４０〜５０ケース算出し、基準状態のＩ−Ｖカーブ、Ｐ−Ｖカーブを作成し、次に前記｛０９｝のように作成したＩ−ＶカーブのＩｓｃｂ，Ｉｏｐｂ，Ｖｏｐｂ，Ｖｏｃｂを算出し、前記のごとく３つの非線形方程式を作成し、さらにこれらの点に近接しない任意の点を選択し、該電圧−電流点を太陽電池基本式に代入してもう１つの非線形方程式を作成し、これら５つの非線形連立方程式を「連立のニュートン・ラプソン法（数値微分）」により、ｂ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数として算出し、次にｃ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓについても前記｛１１｝のように作成したＩ−Ｖカーブを使用して同様に算出し、前記｛１３｝〜｛１６｝の方法で日射強度を算出することを特徴とする。
【００１６】
請求項９の日射強度の測定方法は、
前記請求項４の｛０６｝で作成した４つの非線形連立方程式に加え、与えられた基準状態の特性値（Ｉｓｃａ，Ｉｏｐａ，Ｖｏｐａ，Ｖｏｃａ）の（Ｏ，Ｉｓｃａ），（Ｖｏｐａ，Ｉｏｐａ），（Ｖｏｃａ，Ｏ）の３点を太陽電池特性にもとづく曲線により結び、Ｉ−Ｖカーブを作成し、このＩ−Ｖカーブ上の（Ｏ，Ｉｓｃａ），（Ｖｏｐ，Ｉｏｐａ），（Ｖｏｃａ，Ｏ）とは異なり、これらの点に近接しない任意の点を選択し、もう１つの非線形方程式を作成し、該５つの関係式ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数として解き、この解ＩＬａ，Ｃｏａ，ｎａ，Ｒｓｈａ，Ｒｓａおよび基準状態の温度Ｔａ（２８９Ｋ）を前記関数：Ｆｕｎｃ（Ｖ，Ｉ，ＩＬ，ＣＯ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝０に代入して変数Ｖ，Ｉの関係式を作成し、この式を再び非線形解法のプログラムによりＶに対するＩの解を４０〜５０ケース算出し、基準状態のＩ−Ｖカーブ、Ｐ−Ｖカーブを作成し、次に前記｛０９｝で作成したＩ−ＶカーブのＩｓｃｂ，Ｉｏｐｂ，Ｖｏｐｂ，Ｖｏｃｂを算出し、前記のごとく３つの非線形方程式を作成し、さらにこれらの点に近接しない任意の点を選択し、該電圧−電流点を太陽電池基本式に代入してもう１つの非線形方程式を作成し、これら５つの非線形連立方程式を「連立のニュートン・ラプソン法（数値微分）」により、ｂ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数として算出し、次にｃ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓについても前記｛１１｝で作成したＩ−Ｖカーブを使用して同様に算出し、前記｛１３｝〜｛１６｝の方法で日射強度を算出することを特徴とする。
【００１７】
請求項８の日射強度の測定方法は、
前記｛０９｝で作成したＩ−Ｖカーブ上から近接しない任意の５点を選択し、これらの値をＦｕｎｃ（Ｖ，Ｉ，ＩＬ，ＣＯ，ｎ，Ｒｓｈ，Ｒｓ，Ｔ）＝０の変数Ｖ，Ｉに代入してＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを未知数５個の非線形連立方程式を作成して２５℃の基本特性値を算出し、前記｛０９｝、｛１０｝によりｂ℃における｛１１｝、｛１２｝によってｃ℃における基本特性値ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを算出し、次に太陽電池の実測温度ｔｊ（摂氏：絶対温度はＴｊ＝ｔｊ＋２７３）を取り込み、実測温度ｔｊでの基本特性値Ｍ（Ｉｌｍ，Ｃｏｍ，ｎｍ，Ｒｓｈｍ，Ｒｓｍ）を算出し、前記｛１５｝、｛１６｝の方法で日射強度を算出することを特徴とする。
【００１８】
請求項９の日射強度算出プログラムを記録したコンピュータ読み取り可能なデータ記憶媒体は、
太陽電池の実測された短絡電流Ｉｓｃｊ、太陽電池の実測された温度ｔｊを取り込むための実測値取り込みプログラムと、定数（Ｅｇ，Ｋｏ，ｑ）の入力プログラムと、太陽電池の基準状態での仕様値（Ｉｓｃ，Ｖｏｃ，Ｉｏｐ，Ｖｏｐ，α，β，Ｋ，Ｒｓ）を入力するプログラムと、請求項１，２，３，４，５，６，７，８，９，１０記載の日射強度Ｅを算出する処理プログラムとからなることを特徴とする。
【００１９】
請求項１０の太陽電池日射計は、
前記太陽電池の短絡電流を実測する電流計と、前記太陽電池の温度を実測する温度計（熱電対）と、請求項１，２，３，４，５，６，７，８記載の日射強度を算出する処理プログラムと請求項８の処理プログラムを記録したコンピュータ読み取り可能な記憶媒体が動作しうるコンピュータ、マイコンとからなることを特徴とする。
【００２０】
本発明の請求項１，３，５，６，７，８，９，１０の特徴としては、太陽電池基本特性式中の基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓ，Ｒｓｈ）のそれぞれについて２５℃，ｂ℃（５５℃など），ｃ℃（４０℃など）の値について、Ｒｓを含んで太陽電池基本式を解くことにより基本特性値を求め、それぞれの２５℃，ｂ℃，ｃ℃の値を使った曲線補間により指定温度（ここでは太陽電池の測定温度）におけるＩＬ，Ｃｏ，ｎ，Ｒｓ，Ｒｓｈを求めている。従来の基本特性値を求める方法は、せいぜい実測値に基づき、特定の太陽電池の温度の関数としての基本特性式を作成した式を使う程度であった。また、一部個所では指定温度におけるＩＬ，Ｃｏ，ｎ，Ｒｓ，Ｒｓｈを求めるのに直線補間による方法（電気学会論文１（伊賀：「太陽電池の光照射状態での電圧−電流特性式を用いたＩ−Ｖカーフ゛　作成法とその活用」電学論Ｄ，１１６　巻１０号（１９９６）））も見られるが、本方法は、これらの方法に比べ曲線補間であるとともに各温度におけるＲｓも太陽電池基本式を解くことにより求める方法も適用しているため、さらに精度が高い基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓ，Ｒｓｈ）が求められる。
すなわち、次にこのようにして求めた太陽電池の測定温度における特性値（ＩＬｍ，Ｃｏｍ，ｎｍ，Ｒｓｍ，Ｒｓｈｍ）を、太陽電池基本式から導いた、日射強度を求める式（Ｅ＝（１／ＩＬ）＊［Ｉ＋ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）｛ｅｘｐ（ｑＲｓＩ／ｎＫｏＴ）−１｝＋Ｒｓ＊Ｉ／Ｒｓｈ］）に代入することにより日射強度求めている。なおこの式に適用することで従来の太陽電池日射計より優れているのは、従来は日射強度と太陽電池の短絡電流の間における直線関係（実際には太陽電池温度や日射スペクトルなどから厳密な直線関係には疑問があった）を利用する方法のかわりに、日射強度と光起電流の直線関係（この関係は「光起電流が入射光のスペクトル分布の形が一定な（類似する）限り、入射光の放射照度に正比例する」）特性を利用したものである。したがってより精度の高い日射強度が求まる。
【００２１】
請求項２，３，４，５の特徴は、各測定温度における１部又は全部のＲｓが与えられるため、非線形の連立方程式が確実に解け、基本特性値であるＩＬ，Ｃｏ，ｎ，Ｒｓｈの安定した値が得られる。
【００２２】
ここでさらに請求項５の特徴を具体的に述べると、太陽電池の基本特性値であるＩＬ，Ｃｏ，ｎ，Ｒｓｈ　の値を４つの非線形方程式を使い求める方法については、既にある程度公開しているが（電気学会論文１（伊賀：「太陽電池の光照射状態での電圧電流特性式を用いたＩ−Ｖ　カーブ作成法とその活用」、電学論Ｄ、１１６　巻１０号、１９９６））、しかしこれにＲｓを加えた５つの未知数について解くことは容易ではない。そこで前記のごとき方法でＩＬ，Ｃｏ，Ｒｓｈ　を消去した式（すなわち未知数ｎとＲｓよりなっている式）を作成し、次にこれらの点とは異なる４点を使い、前記太陽電池の基本特性式に代入してｎとＲｓを未知数とする別の式を作成して、前記式と合わせて２つの未知数（ｎ，Ｒｓ）に関する連立非線形方程式を作成し、この式を解いている。このような方法をとったのは基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ）のうちＩＬ，Ｃｏ，Ｒｓｈ　については容易に消去できるためである。またこの解法によると５つの未知数を含む、５つの非線形の方程式（すなわち非線形連立方程式）を直接解く方法に比べ、未知数が２つの連立方程式を解くことになるためコンピュータを使ったニュートン・ラプソン法による解法においても高速かつ、確実に解くことができる。またＩ−Ｖ　カーブ上の７〜８点を使った計算となるためより信頼性の高い解（ＩＬ，Ｃｏ，ｎ，Ｒｓ，Ｒｓｈ）が得られる。
【００２３】
請求項７，８，９の特徴は、Ｒｓの値も未知数として解くとともに、ａ℃（２５℃），　ｂ℃，　ｃ℃における算出方法は共通して基本３点とそれ以外の点による非線形の方程式により算出している。このため算出結果である基本特性値ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓが確実かつ精度よく求められることが期待できる。
【００２４】
請求項１０の特徴は、請求項５，６と同様に算出値がより確実かつ精度よく求められることが期待できる。
【００２５】
ここで本発明に関する用語などについて説明する。
○太陽電池温度は太陽電池の温度であり、通常は太陽電池のセルに埋込んだ熱電対により測定する。なお、太陽電池はその構成の段階によりセル→モジュール→アレイという名称で呼ばれている。
○太陽電池（モジュール）温度に使っている記号については、小文字で示したｔ（℃）　は摂氏を、大文字で示したＴ（Ｋ）　は絶対温度を示す。すなわちＴ（Ｋ）＝ｔ（℃）＋２７３　である。大文字Ｔは主として太陽電池基本式の中で使われている。
○太陽電池出力・発電量とは、太陽電池の発生電圧Ｖと発生電流の積を太陽電池出力（単位ＷまたはＫＷ）と言い、その時間積分値が発電量（単位ｗｈまたはｋｗｈ）　である。
○本発明では特性値については次のとおり使い分けている。
（１）　（太陽電池）基本特性値…ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ…主に太陽電池基本式の中で使用
（２）　（太陽電池）特性値…Ｉｓｃ，Ｉｏｐ，Ｖｏｐ，Ｖｏｃ，α，β，Ｒｓ，Ｋは「実用的Ｉ−Ｖ　カーブ変換式」などで使用。
なお、上記の通りＲｓは（１）　、（２）　の両方で使われている。
○太陽電池基本特性式は前記のように次の式である（前記（１），（２）式を１式で表現したもの）　。
Ｉ＝ＩＬ−Ｃｏ＊Ｔ３＊ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／ｎＫｏＴ）−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ
○温度、日射強度によるＩ−Ｖカーブの変換式は次のとおりである。
（１）　「実用的Ｉ−Ｖ　カーブ変換式」
Ｉ１＝Ｉ２＋Ｉｓｃ（Ｅ１／（Ｅ２）−１）＋α（ｔ１−ｔ２）
Ｖ１＝Ｖ２＋β（ｔ１−ｔ２）−Ｒｓ（Ｉ１−Ｉ２）−Ｋ・Ｉ１・（ｔ１−ｔ２）
（２）　「実用的Ｉ−Ｖ　カーブ変換式」の逆変換式
（（１）　の式のＶ２，Ｉ２について式を解いて整理したものである）
Ｉ２＝Ｉ１＋Ｉｓｃ（Ｅ２−Ｅ１）／Ｅ２＋α（ｔ２−ｔ１）
Ｖ２＝Ｖ１＋β（ｔ２−ｔ１）−Ｒｓ（Ｉ２−Ｉ１）−Ｋ・Ｉ１・（ｔ２−ｔ１）
なお添字の２は基準状態（日射強度１ｋｗ／ｍ２、太陽電池温度２５℃）をあらわし、添字１は測定時状態をあらわす。
また、α：温度が１℃変化したときのＩｓｃ　の変動値（Ａ／℃）（基準状態下）
β：温度が１℃変化したときのＶｏｃ　の変動値（Ｖ／℃）（基準状態下）
Ｒｓ：モジュールの直列抵抗（Ω）（基準状態下）
Ｋ：曲線補正因子（Ω）
Ｉｓｃ：短絡電流（基準状態下）である。
ここで、（１）　、（２）　は一般に知られているＪＩＳＣ８９１３，８９１４，８９１９の式と異なり、優れた特徴のある式である（図６参照）。
○太陽電池の出力特性はＩ−Ｖ　カーブ、Ｐ−Ｖ　カーブにより示される（図５参照）。
【００２６】
【発明の実施の形態】
つぎに、本発明の実施の形態を図面に基づき説明する。
図１は「光起電流による太陽電池日射計」の構成（イメージ）図である。本実施形態の太陽電池日射計は、日射強度の演算処理・出力の方式により２方式に分類できるので示す。すなわち第１は、太陽電池日射計本体に日射強度の演算処理とその信号出力機能を内蔵する方式、および第２は、前記太陽電池短絡電流と太陽電池温度を変換機で変換しＲＳ２３２Ｃケーブルなどを使いパソコンに入力してパソコン内で日射強度の演算をする方式である。
同図において、１は太陽電池日射計の中核部分である太陽電池セル（素子）を示している。
２は太陽電池諸特性値等設定用キーボードである。
３は前記第１の方式の演算処理・出力部を、４は第２の方式のための変換部を示している。またその変換信号の処理コンピュータを１７に示している。
太陽電池セル（素子）は、太陽放射スペクトルと太陽電池の分光感度特性が類似していることが望ましく、また安定した特性が望まれる等の理由からここでは単結晶太陽電池セルを利用する。そして太陽電池の温度は太陽電池に埋込まれた熱電対などにより測定される。
まず前記第１の方式に関して説明する。
太陽電池セルの短絡電流とセル温度のデータは請求項１〜８に述べた方式により内蔵した記憶媒体・演算処理機能を使い日射強度の算出処理をする。太陽電池の特性値であるＩｓｃ，Ｖｏｃ，Ｉｏｐ，α，β，Ｋ，Ｒｓなどや定数（Ｅｇ，Ｋｏ，ｑ）　は本体に設置されたキーボードよりキーインして入力する。もちろんこれらの特性値は、その使用する太陽セル電池の固有の特性値であるため、キーボードより入力するだけでなく、それぞれの太陽電池日射計の記憶媒体にあらかじめ入力されていてもよい。なお太陽電池の特性値（Ｉｓｃ，Ｖｏｃ，Ｉｏｐ，α，β，Ｒｓ，Ｋ　など）は、基準状態の値だけでなく他の温度（例えばｂ℃、ｃ℃）の値や基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓｈなど）が得られる場合には、それらの値を入力して請求項１〜８の演算処理の一部を省略して処理を進めることもできる。また、事前に別途にこれらの特性値を使い基本特性値（２５℃，ｂ℃，Ｃ℃のＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ）を算出し入力しておいてもよい。
演算処理・出力部３と設定用キーボード２は図１では本体の内部に組み込む構造としているが、実際にはこれらの部分を本体から切り離し、別置し、より良好な環境のもとに置くこともできる。そうすることにより、屋外の厳しい環境からマイコン等の電子機器を保護することにもなる。ただその間の接続線、信号線の抵抗などによる測定精度への影響は少ないよう考慮した設計が必要である。
次に前記第２の方式に関して説明する。
太陽電池セルの短絡電流と太陽電池温度についてはパソコン入力可能な信号に変換部で変換してＲＳ２３２Ｃケーブルなどでパソコンに接続する。変換部は前記のように本体内に設置するだけでなく本体から切り離して別置することもできる。この場合も前記第１の方式と同じであるが、その距離、設置場所については耐候性、精度なども考慮することが必要である。
次に、コンピュータ１７を説明する。
符号１７は、本発明の短絡電流と太陽電池温度のデータから日射強度を算出するプログラム１５Ｓ　を動作させるためのコンピュータである。符号１１〜１４はそれぞれコンピュータ本体、キーボード、マウス、モニタを示している。符号１５は、太陽電池の日射強度算出プログラム１５Ｓ　が記録されたコンピュータ読み取り可能なデータ記憶媒体であって、本体内部のディスクでもよいが、例えばＣＤ−ＲＯＭやフレキシブルディスクでもよい。
本コンピュータ本体は、ＣＰＵ　、メモリ、ハードディスク等のコンピュータ基本装置を備えている。このコンピュータ本体には、記憶媒体１５に格納された太陽電池の出力評価プログラム１５Ｓ　をインストールするためのＣＤ−ＲＯＭドライブやフレキシブルディスクドライブを備えている。もちろん、コンピュータ本体にはＷｉｎｄｏｗｓ（登録商標）ＮＴ、Ｗｉｎｄｏｗｓ（登録商標）２０００、Ｗｉｎｄｏｗｓ（登録商標）９８やＭａｃ（登録商標）等のＯＳはインストール済みである。なおコンピュータ１７は図ではノート型であるが、デスクトップ型、タワー型など特に限定されない。コンピュータ１７がノート型の場合、携帯できるので現場で使用でき、好適である。
【００２７】
図２は本発明の太陽電池日射計による日射強度の算出プログラム１５Ｓ　のを請求項３に例をとって示した概要ブロック図である。それぞれのブロックを中心に説明する。
まずＳ１１　は太陽電池の基準状態（日射強度１ｋｗ／ｍ２、太陽電池温度２５℃）の特性値を入力する部分である。Ｓ１３　では基準状態の基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓｈ）を算出し、Ｓ１４　でそのＩ−Ｖ　カーブを作成する。ここでは基準状態における特性値を入力しているが、例えば日射強度１ｋｗ／ｍ２、太陽電池温度５５℃の特性値がメーカーにより与えられ、得られる場合などでは、これらの値を使い直接５５℃の条件における基本特性値を算出しＳ１９　の算出入力データに使うこともできる。もちろん各温度における基本特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ）が得られている場合にはこれらを直接入力し使えばよい。また、Ｓ１５　、Ｓ１７　にあるｂ℃、ｃ℃は例えばそれぞれ５５℃、４０℃であるがそれぞれ７５℃、５０℃などでもよく、要は測定時の太陽電池温度におけるそれぞれの基本特性値を曲線補間により精度よく計算できるための温度であればよい。なお太陽電池基本特性値はＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓとしており、Ｅｇは定数のような取り扱いをしているがＥｇも特性値として扱いＩ−Ｖ　カーブから前記のごとき５個の未知数を６個にして求めることもできる。すなわち、エネルギーギャップＥｇは前記のとおり完全な定数とは扱い難い面もあるため、２５℃のエネルギーギャップＥｇはＳ１１　で与えておき、他の温度（ｂ℃、ｃ℃）のＥｇはＳ１６　、Ｓ１８　においてＩ−Ｖ　カーブ上の値の数を５点から６点に増加させることにより非線形連立方程式の数を１つ増加させることにより求めることもできる。
次にＳ１２　で入力する定数Ｅｇ（以下ここでは定数として説明する），Ｋｏ，ｅ（ｑ）はそれぞれ、エネルギーギャップ、ボルツマン定数および電子の電荷量であり、いずれも定数であるので、予め値を与えておけばよい。
Ｓ１３　では太陽電池の基準状態の特性値Ｉｓｃ，Ｖｏｃ，Ｉｏｐ，Ｖｏｐ，Ｒｓから次の太陽電池基本式
Ｉ＝ＩＬ−ＣｏＴ３ｅｘｐ（−ｑＥｇ／ｎＫｏＴ）＊｛ｅｘｐ（ｑ（Ｖ＋ＲｓＩ）／ｎＫｏＴ）−１｝−（Ｖ＋ＲｓＩ）／Ｒｓｈ
を使い基準状態の基本特性値ＩＬ，Ｃｏ，ｎ，Ｒｓｈを算出している。算出はＩ−Ｖカーブ上の３点（０，Ｉｓｃ），（Ｖｏｐ，Ｉｏｐ），（Ｖｏｃ，０）を上記太陽電池基本式に代入し、請求項３に記述するようにＩＬ，Ｃｏ，ｎ，Ｒｓｈを未知数とする非線形の４つの連立方程式をつくり、未知数のうちＩＬ，Ｃｏ，Ｒｓｈ　を消去の上ｎを未知数とする非線形方程式をつくり、ニュートン法によりｎをさらにＩＬ，Ｃｏ，Ｒｓを求めている。この部分の解法については前記の電気学会論文１（伊賀；「太陽電池の光照射状態での電圧−電流特性式を用いたＩ−Ｖ　カーブ作成法とその活用」、電学論、１１６　巻１０号、１９９６）に示してある。
Ｓ１４　では上記で求めたＩＬ，Ｃｏ，ｎ，Ｒｓｈ　とＲｓを太陽電池基本式に代入し、また定数であるＥｇ，Ｋｏ，ｅ　を入れ、Ｖの各値に対するＩの値をニュートン法によって解くことによって得られる電圧と電流の関係、約４０〜５０組の値を使い、Ｉ−Ｖ　カーブを作成するものである。
Ｓ１５　、Ｓ１７　はＳ１４　で得られたＩ−Ｖ　カーブ上の各点（約４０〜５０点）について、ｂ℃又はｃ℃の値に変換するため次式
Ｉ１　＝Ｉ２　＋α（ｔ１　−ｔ２　）
Ｖ１　＝Ｖ２　＋β（ｔ１　−ｔ２　）−Ｒｓ（Ｉ１　−Ｉ２　）−ＫＩ１　（ｔ１　−ｔ２　）
（ここで添字１のついているものはｂ℃又はｃ℃における値であり、添字２のついた値は２５℃における値である。）を使い、ｂ℃又はｃ℃におけるＩ−Ｖ　カーブを求めていることを示している。
Ｓ１６　、Ｓ１８　は上記Ｓ１５　、Ｓ１７　で求めたＩ−Ｖ　カーブ上の点（約５０点）のうち近接しない５点を任意に選択し、太陽電池基本特性式に代入して未知数５の非線形の連立方程式を作成して、これら非線形の連立方程式をニュートン法で解いてＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを求めるものである。このＳ１６　、Ｓ１８　の解法については、別の方法もあり、そのため別途請求項５および図３で説明しているので後述する。
Ｓ１９　は２５℃、ｂ℃、ｃ℃における太陽電池基本特性値ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓの各値から太陽電池モジュール温度ｔ℃における値を曲線補間により求める。すなわちｔ℃におけるＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを精度よく求めている。
なおＳ２０　は測定時の太陽電池温度の入力部分である。
Ｓ２１　はＳ１９　で求めた測定時太陽電池温度における太陽電池基本特性値ＩＬｍ，Ｃｏｍ，ｎｍ，Ｒｓｈｍ，Ｒｓｍ　および測定した太陽電池短絡電流Ｉｊ、太陽電池温度ｔｊ（絶対温度Ｔ＝ｔｊ＋２７３）　を太陽電池基本式から誘導した式に示した式に代入して日射強度Ｅを求めるものである。
【００２８】
図３は図２におけるＳ１６　、Ｓ１８　の別途の解法を説明したもので請求項５に相当する。すなわち太陽電池基本特性値ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓを効率よく確実に求める方法を示したものである。なおＳ１０１〜Ｓ１０５については既に前記論文１で述べた方法と一致している。すなわち、４つの未知数については既に解かれており、本発明では５つの未知数の非線形の連立方程式をその応用として解く方法を中心に示す。
ここでＳ１０１は太陽電池基本特性式の電圧−電流の値に短絡電流の点（０，Ｉｓｃ）　、最大出力の電圧−電流の点（Ｖｏｐ，Ｉｏｐ）　、および開放電圧の点（Ｖｏｃ，０）　を代入して３つの非線形方程式をつくり、またＳ１０２は太陽電池基本特性式を電圧Ｖで微分することを示し、Ｓ１０３はＳ１０２で微分した式に最大出力の電圧−電流の点を代入して非線形の方程式をつくっている。Ｓ１０４はＳ１０１とＳ１０３で作成した４つの非線形連立方程式を示しており、Ｓ１０５は４つの非線形連立方程式について、線形の未知数であるＩＬ，Ｃｏ，Ｒｓｈ　を消去してＲｓを含む未知数ｎのみの非線形方程式を作成することより解いていることを示している。この内容は既に前記論文Ｉに詳しく示している。
次にＳ１０６は、Ｉ−Ｖカーブ上の前記（０，Ｉｓｃ），（Ｖｏｐ，Ｉｏｐ），（Ｖｏｃ，０）　の３点に加え、これらの点以外でお互いに近接しないＩ−Ｖ　カーブ上の４点を任意に選択し、Ｓ１０７でこれらの点を太陽電池基本特性式に代入して４つの非線形連立方程式を作成することを示している。近接しない４点をとったのはできるだけ計算結果の算出精度をあげるためである。
Ｓ１０８はＳ１０７の４つの非線形連立方程式の線形の未知数ＩＬ，Ｃｏ，Ｒｓｈ　を消去してｎとＲｓの非線形方程式を求めることを示している。
Ｓ１０９はＳ１０５とＳ１０８で求めたＲｓ，ｎを未知数とする非線形方程式を連立させて、Ｓ１１０でニュートン・ラプソン法で解くことを示している。Ｒｓ，ｎが求まるとＳ１０５又はＳ１０８の非線形の連立方程式にもどして容易にＩＬ，Ｃｏ，Ｒｓｈ　を求めることができる（Ｓ１１２）。
【００２９】
図４は太陽電池の発電のしくみを示している。
【００３０】
図５は太陽電池の出力特性、すなわち太陽電池の出力端の電圧−電流特性（Ｉ−Ｖ　カーブ）、電圧−電力特性（Ｐ−Ｖ　カーブ）を示している。
【００３１】
図６は測定時の日射強度、太陽電池温度条件における電圧Ｖ−電流値Ｉを基準状態（日射強度１ｋｗ／ｍ２、太陽電池温度２５℃）の電圧値Ｖ−電流値Ｉに変換する式に関して発明者が開発した式を中心に比較したものを示した（電気学会論文２参照）。ここで示した変換式は、変形することにより基準状態の電圧Ｖ−電流値Ｉを任意の日射強度・モジュール温度条件の電圧Ｖ−電流値Ｉに変換する式にすることができる。この変換式は図２のＳ１５　、Ｓ１７　などで活用している。
【００３２】
【発明の効果】
請求項１〜８の方法により求めた日射強度は、前記のごとく短絡電流でなく太陽電池温度を考慮した光起電流にもとづくため、太陽電池の特性をより正確に反映し、精度よく日射強度を求めることができる。また一般に使用されているサーモパイルなどによる熱的な日射計に比べて応答が格段に早いため、太陽電池出力の評価（特願平２００１−０７２０８１参照）における必須の条件である時間的な応答の条件を満たしている日射計として最適である。
請求項１，３，５，６，７，８による太陽電池特性値（ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓ）を求める方法は、Ｉ−Ｖ　カーブをより忠実に反映した特性値の算出方法（すなわちＲｓも未知数とする算出方法）である。また確実に解ける方法であるので、信頼性と精度の高い日射強度の測定値が確実に得られる。
請求項９の日射強度演算プログラムを記録したコンピュータ読み取り可能なデータ記録媒体によれば、日射計と一体化した装置とすることも、またパソコンの中での処理もすることも可能である。
請求項１０のごとく前記方法を用いた日射計によると時間遅れのない日射強度が得られるため、正確かつ適正な太陽電池出力の評価ができる。
【００３３】
【図面の簡単な説明】
【図１】光起電流による太陽電池日射計の構成図である。
【図２】実施形態の日射強度算出のブロック図である。
【図３】太陽電池基本特性値ＩＬ，Ｃｏ，ｎ，Ｒｓｈ，Ｒｓの算出ブロック図である（請求項３関連）。
【図４】太陽電池の発電のしくみである。
【図５】太陽電池出力特性曲線（Ｉ−Ｖ　カーブ、Ｐ−Ｖ　カーブ）
【図６】測定時の日射強度・太陽電池温度条件における電圧Ｖ−電流Ｉ値の基準状態への変換式などである。
【符号の説明】
１　太陽電池セル（素子）
２　（諸特性値等）設定用キーボード
３　演算処理・出力部
４　変換部
１１　コンピュータ（パソコン）本体
１２　コンピュータキーボード
１３　コンピュータ用マウス
１４　コンピュータディスプレイ
１５　記憶装置
１７　コンピュータ（パソコン）
【図１】
2004077311.tif 000003
【図２】
2004077311.tif 000004
【図３】
2004077311.tif 000005
【図４】
2004077311.tif 000006
【図５】
2004077311.tif 000007
【図６】
2004077311.tif 000008