(57)【要約】【課題】照明機器の発光光度を照度測定装置からの目標照度に制御するために、照明機器に対する照度測定装の距離的な位置関係を、常に、より正しく判断し得る照明システムを提供する。【解決手段】各照明機器の発光光度を、任意位置の照度測定装置からの測定照度および目標照度に基づき、少なくとも、当該照度測定装置の測定照度に影響を与える照明機器に対して、測定照度が目標照度に近づくように且つ電力計による消費電力も最適となるように制御する照明システムであって、上記測定照度に影響を与える照明機器を判断する指標として、所定の光度生成範囲内で変化されて照明機器に与えられる発光光度および照度測定装置にて測定される測定照度に基づき演算される回帰係数を用いるとともに、当該回帰係数を最小二乗法により求めるようにしたものである。【選択図】図１５
2008243389.tif 000002
【特許請求の範囲】
【請求項１】
複数の照明機器の影響下に位置する照度測定装置からの測定照度および目標照度に基づき、少なくとも、当該照度測定装置の測定照度に影響を与える照明機器に対して、測定照度が目標照度に近づくように、その発光光度を制御する照明制御方法であって、
上記測定照度に影響を与える照明機器を判断する指標として、所定の光度生成範囲内で変化されて照明機器に与えられる光度変化量および照度測定装置にて測定される照度変化量に基づき演算される回帰係数を用いることを特徴とする照明制御方法。
【請求項２】
回帰係数を、逐次型最小二乗法またはバッチ型最小二乗法を用いて求めることを特徴とする請求項１に記載の照明制御方法。
【請求項３】
複数の照明機器の発光光度およびこれら照明機器の影響下に位置する照度測定装置からの測定照度および当該照度測定装置に設定された目標照度並びに照明機器の消費電力に基づき、測定照度が目標照度の設定許容範囲内に入るように制御する照明制御装置であって、
予め設定された複数の光度生成範囲のうち、所定の光度生成範囲内で変化される発光光度の変化量およびそれに対応して得られる測定照度の変化量に基づき、各照明機器が照度測定装置の測定照度に与える影響度合いを判断するための回帰係数を演算する回帰係数演算手段と、
この回帰係数演算手段で求められた回帰係数および目標照度に対する測定照度の偏差に基づき、次の発光光度を生成するための上記所定の光度生成範囲を選択する光度生成範囲選択手段と、
この光度生成範囲選択手段で選択された光度生成範囲内で変化される光度値でもって発光光度を生成する光度生成手段と、
この光度生成手段にて生成された発光光度、この発光光度での発光による測定照度、回帰係数、および照明機器の消費電力に基づき、上記生成された発光光度の良否を評価するための評価値を演算する評価値演算手段と
を具備したことを特徴とする照明制御装置。
【請求項４】
照度測定装置が移動したことを検知し得る移動検知手段を具備したことを特徴とする請求項３に記載の照明制御装置。
【請求項５】
回帰係数演算手段において、回帰係数を逐次型最小二乗法またはバッチ型最小二乗法により求めるようにしたことを特徴とする請求項３または４に記載の照明制御装置。
【請求項６】
回帰係数演算手段は、
回帰係数に対する補正値を求める推定誤差共分散の逐次演算と、この補正値を用いて推定値である回帰係数の逐次演算とを有し、入力値である光度変化量を基準としたときの観測値である照度変化量の観測予測誤差を最小とする推定値である回帰係数を再帰的に推定する逐次型最小二乗法を用いるもので、
下記（１）式に基づき、或る時刻（演算ステップに相当する）（ｋ）における入力値ｘ（ｋ）、前の時刻（前回の演算ステップに相当する）（ｋ−１）における推定値の誤差共分散Ｐ（ｋ）および観測雑音ｗの共分散Ｒ（ｋ）を用いて、観測値の予測誤差共分散Ｓ（ｋ）を演算する観測予測誤差共分散演算部と、
下記（２）式に基づき、入力値ｘ（ｋ）、推定誤差共分散Ｐ（ｋ−１）および観測予測誤差共分散Ｓ（ｋ）を用いて、補正係数Ｗ（ｋ）を演算する補正係数演算部と、
下記（３）式に基づき、推定誤差共分散Ｐ（ｋ−１）、補正係数Ｗ（ｋ）、観測予測誤差共分散Ｓ（ｋ）を用いて、或る時刻（ｋ）における推定誤差共分散Ｐ（ｋ）を演算する推定誤差共分散演算部と、
下記（４）式に基づき、入力値ｘ（ｋ）および前の時刻における推定値β（ｋ−１）を用いて、或る時刻における観測値の予測値＊ｙ（ｋ）を演算する観測予測値演算部と（なお、記号「＊」は、以下に示す式中において、この後に続く文字にハット記号が付されていることを示す、以下同じ）、
下記（５）式に基づき、或る時刻（ｋ）における観測値ｙ（ｋ）および観測予測値＊ｙ（ｋ）を用いて、観測値の予測誤差ε（ｋ）を演算する観測予測誤差演算部と、
下記（６）式に基づき、前の時刻における回帰係数β（ｋ−１）、補正係数Ｗ（ｋ）および観測予測誤差ε（ｋ）を用いて、或る時刻における推定値＊β（ｋ）を演算する推定値演算部と
を有することを特徴とする請求項３または４に記載の照明制御装置。
Ｓ（ｋ）＝ｘ（ｋ）Ｐ（ｋ−１）Ｓ（ｋ）＋Ｒ（ｋ）  ・・・（１）
Ｗ（ｋ）＝Ｐ（ｋ−１）ｘ（ｋ）Ｓ−１（ｋ）＋Ｒ（ｋ）  ・・・（２）
Ｐ（ｋ）＝Ｐ（ｋ−１）−Ｗ（ｋ）Ｓ（ｋ）Ｗ（ｋ）  ・・・（３）
【数１】
2008243389.tif 000016
【数２】
2008243389.tif 000017
【数３】
2008243389.tif 000018
【請求項７】
複数の照明機器と、これら各照明機器に対応して配置されて各照明機器で発光させる発光光度を制御する複数の照明制御装置と、上記各照明機器にて消費される消費電力を計測する電力計と、上記照明機器による照度を測定し得る照度測定装置と、上記各照明制御装置、電力計および照度測定装置との間で、所定データの受け渡しを行い得るデータ通信手段とを具備する照明システムであって、
上記照度測定装置は、
受光により照度値を検出し得る光センサ部と、目標照度を設定し得る目標照度設定部と、上記光センサ部にて測定された測定照度、上記目標照度設定部で設定された目標照度、および当該照度測定装置の識別情報をデータ通信手段に出力するデータ出力部とから構成し、
上記照明制御装置は、
任意に変化される発光光度の変化量または所定の方法で変化される発光光度の変化量と測定照度の変化量とに基づき、各照明機器が照度測定装置の照度に与える影響度合いを判断するための回帰係数を演算する回帰係数演算手段と、
この回帰係数演算手段で求められた回帰係数および目標照度に対する測定照度の偏差に基づき、次に生成する光度範囲を選択する光度生成範囲選択手段と、
この光度生成範囲選択手段で選択された光度生成範囲内で任意に変化される光度値または所定の方法で変化される光度値でもって発光光度を生成する光度生成手段と、
この光度生成手段にて生成された発光光度、この発光光度での発光による測定照度、回帰係数、および照明機器の消費電力に基づき、当該生成された発光光度の良否を評価するための評価値を演算する評価値演算手段と
を有することを特徴とする照明システム。
【請求項８】
照明制御装置に、照度測定装置の移動を検知し得る移動検知手段を具備したことを特徴とする請求項７に記載の照明システム。
【請求項９】
回帰係数演算手段において、回帰係数を逐次型最小二乗法またはバッチ型最小二乗法により求めるようにしたことを特徴とする請求項７に記載の照明システム。
【請求項１０】
回帰係数演算手段は、
回帰係数に対する補正値を求める推定誤差共分散の逐次演算と、この補正値を用いて推定値である回帰係数の逐次演算とを有し、入力値である光度変化量を基準としたときの観測値である照度変化量の観測予測誤差を最小とする推定値である回帰係数を再帰的に推定する逐次型最小二乗法を用いるもので、
下記（１）式に基づき、或る時刻（演算ステップに相当する）（ｋ）における入力値ｘ（ｋ）、前の時刻（前回の演算ステップに相当する）（ｋ−１）における推定値の誤差共分散Ｐ（ｋ）および観測雑音ｗの共分散Ｒ（ｋ）を用いて、観測値の予測誤差共分散Ｓ（ｋ）を演算する観測予測誤差共分散演算部と、
下記（２）式に基づき、入力値ｘ（ｋ）、推定誤差共分散Ｐ（ｋ−１）および観測予測誤差共分散Ｓ（ｋ）を用いて、補正係数Ｗ（ｋ）を演算する補正係数演算部と、
下記（３）式に基づき、推定誤差共分散Ｐ（ｋ−１）、補正係数Ｗ（ｋ）、観測予測誤差共分散Ｓ（ｋ）を用いて、或る時刻（ｋ）における推定誤差共分散Ｐ（ｋ）を演算する推定誤差共分散演算部と、
下記（４）式に基づき、入力値ｘ（ｋ）および前の時刻における推定値β（ｋ−１）を用いて、或る時刻における観測値の予測値＊ｙ（ｋ）を演算する観測予測値演算部と（なお、記号「＊」は、以下に示す式中において、この後に続く文字にハット記号が付されていることを示す、以下同じ）、
下記（５）式に基づき、或る時刻（ｋ）における観測値ｙ（ｋ）および観測予測値＊ｙ（ｋ）を用いて、観測値の予測誤差ε（ｋ）を演算する観測予測誤差演算部と、
下記（６）式に基づき、前の時刻における回帰係数β（ｋ−１）、補正係数Ｗ（ｋ）および観測予測誤差ε（ｋ）を用いて、或る時刻における推定値＊β（ｋ）を演算する推定値演算部と
を有することを特徴とする請求項７または８に記載の照明システム。
Ｓ（ｋ）＝ｘ（ｋ）Ｐ（ｋ−１）Ｓ（ｋ）＋Ｒ（ｋ）  ・・・（１）
Ｗ（ｋ）＝Ｐ（ｋ−１）ｘ（ｋ）Ｓ−１（ｋ）＋Ｒ（ｋ）  ・・・（２）
Ｐ（ｋ）＝Ｐ（ｋ−１）−Ｗ（ｋ）Ｓ（ｋ）Ｗ（ｋ）  ・・・（３）
【数４】
2008243389.tif 000019
【数５】
2008243389.tif 000020
【数６】
2008243389.tif 000021
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、任意の場所に適切な照度を提供し得る照明制御方法および照明制御装置並びに照明システムに関するものである。
【背景技術】
【０００２】
例えば、展示場、図書館、会議場などの広い部屋（空間室）の天井には、非常に多くの照明機器が配置されており、通常、各照明機器の発光光度は一定値にされている。
しかし、展示場、図書館などにあっては、物品の展示位置または利用者の好みにより、要求される照度が異なる場合があり、従来、このような要求に応じるものとして、個々の照明機器または小範囲のエリア毎の照明機器の光度を、手動にて調節するものがあったが、その調節作業は、非常に面倒であるという欠点があった。
【０００３】
また、会議室などにあっては、照明機器の光度値が決められているが、経年変化などによりその光度が低下するため、やはり、定期的に照度を測定して調節する必要があり、やはり、その作業が非常に面倒であるという欠点があった。
【０００４】
さらに、部屋に窓があると、差し込む日光により手元の照度が大きく変化するとともに、照明機器による発光をそれ程必要としない場合には、電力が無駄に消費されるという欠点があった。
【０００５】
このような欠点に対処するものとして、照明機器からの発光光度による照度が目標照度に自動的に近づくように制御する照明システムが提案されている（例えば、特許文献１参照）。
【０００６】
この照明システムは、少なくとも、制御部を有する照明機器と、照明箇所に配置されてその照度を測定するとともに各照明機器の制御部にその照度および利用者が設定した目標照度を送信するようにした照度測定装置つまり照度センサとが具備されたものである。
【０００７】
例えば、図書館を例にして説明すると、利用者がその照度センサを持って或る場所に移動すると、その照度センサから目標照度および測定照度が制御部に送信され、この制御部にて、任意の照明機器に対する照度センサの位置関係を判断するとともに、照度センサからの目標照度に近づくように、その発光光度を制御するようにしたものである。
【０００８】
そして、この位置関係を判断する際に、照明機器における光度の変化量と照度センサで測定された照度の変化量との類似度を考慮して行われており、具体的には、所定の光度範囲内でランダムに変化される光度の時系列データと、照度センサにより取得された照度の時系列データとの相関係数により、互いの類似度が判断されていた。
【０００９】
すなわち、照度センサが複数ある場合、相関係数が１に近いほど、照明機器の光度による影響を最も受けており、言い換えれば、照明機器に近い位置に存在すると判断することにより、当該判断された照明機器について、その光度による照度が、目標照度に近づくように制御が行われるとともに、照明機器から遠いものについては、その光度を減らすように制御が行われている。
【特許文献１】特開２００６−３０２５１７
【発明の開示】
【発明が解決しようとする課題】
【００１０】
ところで、上述した照明システムの構成によると、照明機器に対する照度センサの位置関係を判断する際に相関係数が用いられているが、全ての判断が、精度良く行われるとは限らない。すなわち、相関係数は、２組のデータ列の傾向性が似ていることは判断できるが、照明機器と照度センサとの距離的な関係を常に正しく判断することができないという問題があった。
【００１１】
そこで、本発明は、照明機器の発光光度を照度測定装置からの目標照度に制御するために、照明機器に対する照度測定装の距離的な位置関係を、常に、より正しく判断し得る照明制御方法および照明制御装置並びに照明システムを提供することを目的とする。
【課題を解決するための手段】
【００１２】
上記課題を解決するため、本発明の照明制御方法は、複数の照明機器の影響下に位置する照度測定装置からの測定照度および目標照度に基づき、少なくとも、当該照度測定装置の測定照度に影響を与える照明機器に対して、測定照度が目標照度に近づくように、その発光光度を制御する照明制御方法であって、
上記測定照度に影響を与える照明機器を判断する指標として、所定の光度生成範囲内で変化されて照明機器に与えられる光度変化量および照度測定装置にて測定される照度変化量に基づき演算される回帰係数を用いる方法であり、
また上記回帰係数を、逐次型最小二乗法またはバッチ型最小二乗法を用いて求める方法である。
【００１３】
また、本発明に係る照明制御装置は、複数の照明機器の発光光度およびこれら照明機器の影響下に位置する照度測定装置からの測定照度および当該照度測定装置に設定された目標照度並びに照明機器の消費電力に基づき、測定照度が目標照度の設定許容範囲内に入るように制御する照明制御装置であって、
予め設定された複数の光度生成範囲のうち、所定の光度生成範囲内で変化される発光光度の変化量およびそれに対応して得られる測定照度の変化量に基づき、各照明機器が照度測定装置の測定照度に与える影響度合いを判断するための回帰係数を演算する回帰係数演算手段と、
この回帰係数演算手段で求められた回帰係数および目標照度に対する測定照度の偏差に基づき、次の発光光度を生成するための上記所定の光度生成範囲を選択する光度生成範囲選択手段と、
この光度生成範囲選択手段で選択された光度生成範囲内で変化される光度値でもって発光光度を生成する光度生成手段と、
この光度生成手段にて生成された発光光度、この発光光度での発光による測定照度、回帰係数、および照明機器の消費電力に基づき、上記生成された発光光度の良否を評価するための評価値を演算する評価値演算手段とを具備したものであり、
また照度測定装置が移動したことを検知し得る移動検知手段を具備したものであり、
さらに上記回帰係数演算手段において、回帰係数を逐次型最小二乗法またはバッチ型最小二乗法により求めるようにしたものである。
【００１４】
さらに、上記照明制御装置における回帰係数演算手段は、
回帰係数に対する補正値を求める推定誤差共分散の逐次演算と、この補正値を用いて推定値である回帰係数の逐次演算とを有し、入力値である光度変化量を基準としたときの観測値である照度変化量の観測予測誤差を最小とする推定値である回帰係数を再帰的に推定する逐次型最小二乗法を用いるもので、
下記（１）式に基づき、或る時刻（演算ステップに相当する）（ｋ）における入力値ｘ（ｋ）、前の時刻（前回の演算ステップに相当する）（ｋ−１）における推定値の誤差共分散Ｐ（ｋ）および観測雑音ｗの共分散Ｒ（ｋ）を用いて、観測値の予測誤差共分散Ｓ（ｋ）を演算する観測予測誤差共分散演算部と、
下記（２）式に基づき、入力値ｘ（ｋ）、推定誤差共分散Ｐ（ｋ−１）および観測予測誤差共分散Ｓ（ｋ）を用いて、補正係数Ｗ（ｋ）を演算する補正係数演算部と、
下記（３）式に基づき、推定誤差共分散Ｐ（ｋ−１）、補正係数Ｗ（ｋ）、観測予測誤差共分散Ｓ（ｋ）を用いて、或る時刻（ｋ）における推定誤差共分散Ｐ（ｋ）を演算する推定誤差共分散演算部と、
下記（４）式に基づき、入力値ｘ（ｋ）および前の時刻における推定値β（ｋ−１）を用いて、或る時刻における観測値の予測値＊ｙ（ｋ）を演算する観測予測値演算部と（なお、記号「＊」は、以下に示す式中において、この後に続く文字にハット記号が付されていることを示す、以下同じ）、
下記（５）式に基づき、或る時刻（ｋ）における観測値ｙ（ｋ）および観測予測値＊ｙ（ｋ）を用いて、観測値の予測誤差ε（ｋ）を演算する観測予測誤差演算部と、
下記（６）式に基づき、前の時刻における回帰係数β（ｋ−１）、補正係数Ｗ（ｋ）および観測予測誤差ε（ｋ）を用いて、或る時刻における推定値＊β（ｋ）を演算する推定値演算部とを有するものである。
【００１５】
Ｓ（ｋ）＝ｘ（ｋ）Ｐ（ｋ−１）Ｓ（ｋ）＋Ｒ（ｋ）  ・・・（１）
Ｗ（ｋ）＝Ｐ（ｋ−１）ｘ（ｋ）Ｓ−１（ｋ）＋Ｒ（ｋ）  ・・・（２）
Ｐ（ｋ）＝Ｐ（ｋ−１）−Ｗ（ｋ）Ｓ（ｋ）Ｗ（ｋ）  ・・・（３）
【００１６】
【数１】
2008243389.tif 000003
【００１７】
【数２】
2008243389.tif 000004
【００１８】
【数３】
2008243389.tif 000005
【００１９】
また、本発明に係る照明システムは、複数の照明機器と、これら各照明機器に対応して配置されて各照明機器で発光させる発光光度を制御する複数の照明制御装置と、上記各照明機器にて消費される消費電力を計測する電力計と、上記照明機器による照度を測定し得る照度測定装置と、上記各照明制御装置、電力計および照度測定装置との間で、所定データの受け渡しを行い得るデータ通信手段とを具備する照明システムであって、
上記照度測定装置は、
受光により照度値を検出し得る光センサ部と、目標照度を設定し得る目標照度設定部と、上記光センサ部にて測定された測定照度、上記目標照度設定部で設定された目標照度、および当該照度測定装置の識別情報をデータ通信手段に出力するデータ出力部とから構成し、
上記照明制御装置は、
任意に変化される発光光度の変化量または所定の方法で変化される発光光度の変化量と測定照度の変化量とに基づき、各照明機器が照度測定装置の照度に与える影響度合いを判断するための回帰係数を演算する回帰係数演算手段と、
この回帰係数演算手段で求められた回帰係数および目標照度に対する測定照度の偏差に基づき、次に生成する光度範囲を選択する光度生成範囲選択手段と、
この光度生成範囲選択手段で選択された光度生成範囲内で任意に変化される光度値または所定の方法で変化される光度値でもって発光光度を生成する光度生成手段と、
この光度生成手段にて生成された発光光度、この発光光度での発光による測定照度、回帰係数、および照明機器の消費電力に基づき、当該生成された発光光度の良否を評価するための評価値を演算する評価値演算手段とを有するものであり、
また上記照度測定装置の移動を検知し得る移動検知手段を具備したものであり、
また上記回帰係数演算手段において、回帰係数を逐次型最小二乗法またはバッチ型最小二乗法により求めるようにしたものである。
【００２０】
さらに、上記照明制御装置における回帰係数演算手段は、
回帰係数に対する補正値を求める推定誤差共分散の逐次演算と、この補正値を用いて推定値である回帰係数の逐次演算とを有し、入力値である光度変化量を基準としたときの観測値である照度変化量の観測予測誤差を最小とする推定値である回帰係数を再帰的に推定する逐次型最小二乗法を用いるもので、
下記（１）式に基づき、或る時刻（演算ステップに相当する）（ｋ）における入力値ｘ（ｋ）、前の時刻（前回の演算ステップに相当する）（ｋ−１）における推定値の誤差共分散Ｐ（ｋ）および観測雑音ｗの共分散Ｒ（ｋ）を用いて、観測値の予測誤差共分散Ｓ（ｋ）を演算する観測予測誤差共分散演算部と、
下記（２）式に基づき、入力値ｘ（ｋ）、推定誤差共分散Ｐ（ｋ−１）および観測予測誤差共分散Ｓ（ｋ）を用いて、補正係数Ｗ（ｋ）を演算する補正係数演算部と、
下記（３）式に基づき、推定誤差共分散Ｐ（ｋ−１）、補正係数Ｗ（ｋ）、観測予測誤差共分散Ｓ（ｋ）を用いて、或る時刻（ｋ）における推定誤差共分散Ｐ（ｋ）を演算する推定誤差共分散演算部と、
下記（４）式に基づき、入力値ｘ（ｋ）および前の時刻における推定値β（ｋ−１）を用いて、或る時刻における観測値の予測値＊ｙ（ｋ）を演算する観測予測値演算部と（なお、記号「＊」は、以下に示す式中において、この後に続く文字にハット記号が付されていることを示す、以下同じ）、
下記（５）式に基づき、或る時刻（ｋ）における観測値ｙ（ｋ）および観測予測値＊ｙ（ｋ）を用いて、観測値の予測誤差ε（ｋ）を演算する観測予測誤差演算部と、
下記（６）式に基づき、前の時刻における回帰係数β（ｋ−１）、補正係数Ｗ（ｋ）および観測予測誤差ε（ｋ）を用いて、或る時刻における推定値＊β（ｋ）を演算する推定値演算部とを有するものである。
【００２１】
Ｓ（ｋ）＝ｘ（ｋ）Ｐ（ｋ−１）Ｓ（ｋ）＋Ｒ（ｋ）  ・・・（１）
Ｗ（ｋ）＝Ｐ（ｋ−１）ｘ（ｋ）Ｓ−１（ｋ）＋Ｒ（ｋ）  ・・・（２）
Ｐ（ｋ）＝Ｐ（ｋ−１）−Ｗ（ｋ）Ｓ（ｋ）Ｗ（ｋ）  ・・・（３）
【００２２】
【数４】
2008243389.tif 000006
【００２３】
【数５】
2008243389.tif 000007
【００２４】
【数６】
2008243389.tif 000008
【発明の効果】
【００２５】
上記照明制御方法および照明制御装置並びに照明システムによると、複数配置された照明機器に対して、照度測定装置に設定された目標照度となるように、自動的に、その発光光度を制御するようにした自律分散型の照明制御アルゴリズムにおいて、照明機器に対する照度測定装置の測定照度に影響を与える影響度合いとして、回帰係数を用いるようにしたので、例えば相関係数を用いる場合に比べて、互いの位置関係を精度良く把握することができ、すなわち照明機器の光度変化量と照度測定装置での照度変化量とを回帰分析することにより得られる回帰係数は、照明機器が照度測定装置に影響するかどうかだけでなく、その距離感まで表すことができ、したがって相関係数を用いる場合に比べて、より正確な位置関係の把握が可能となり、各照明機器を、照度測定装置での目標照度を満たし且つ省電力な状態へと素早く収束させることができる。
【発明を実施するための最良の形態】
【００２６】
［実施の形態］
以下、本発明の実施の形態に係る照明制御方法および照明制御装置並びに照明システムおよび照明システムを、図１〜図２５に基づき説明する。
【００２７】
本実施の形態においては、例えば図書室などの広い部屋の天井（空間室の上方位置）に多数配置された蛍光灯などの照明機器を、図書室の利用者が携帯する照度測定装置により、任意の位置で利用者が希望する照度が得られるように制御する場合について説明する。
【００２８】
なお、以下の説明において、照明機器から発光される発光光度および照度測定装置により測定される測定照度については、単に、光度および照度と称したり、また必要に応じて、光度値および照度値と称する。また、電力計による測定される消費電力（使用電力でもある）についても、必要に応じて、消費電力値と称する。さらに、照明機器については、例えば蛍光灯（後述する）が用いられるが、以下においては、単に、「照明」と称して説明する場合があり、また照度測定装置についても、単に、「照度センサ」または「センサ」と称して説明する場合がある。
【００２９】
この照明システムは、図１に示すように、大きく分けて、部屋の天井に配置されるとともに電力供給源（例えば、商用電力である）１に電気配線２を介して接続された複数の照明機器３と、これら各照明機器３に一体的に設けられて（勿論、別置きでもよい）当該照明機器３にて発光される光度（カンデラＣｄで表される）を制御する照明制御装置（例えば、マイクロチップなどで構成される）４と、利用者が携帯（保持）して現在位置で測定された照度（ルクスｌｘで表される）を検出し得る移動式の照度測定装置５と、上記各照明機器３における全消費電力を計測する電力計６とから構成されており、さらに上記照明制御装置４、照度測定装置５および電力計６は、ネットワーク回線（データ通信手段の一例で、有線ＬＡＮ７ａおよびブルートゥースなどの無線方式の無線ＬＡＮ７ｂからなるものが用いられ、勿論、データ受信用の端末器７ｃが設けられている）７を介して接続されている。
【００３０】
上記照明機器３としては、発光源である蛍光灯３ａと、上記照明制御装置４から指示された光度値に基づき当該蛍光灯３ａへの供給電圧を制御してその発光光度を調節する光度調節部（例えば、インバータなどが用いられる）３ｂとから構成されている。勿論、発光源として蛍光灯以外のもの、例えば白熱灯などを用いてもよい。
【００３１】
上記照度測定装置５は、図２に示すように、例えば携帯し得るカード型のものであり、矩形状のカード本体部５ａの表面に配置されて受光するとともにその照度を検出し得る光センサ部５ｂと、目標照度および当該照度測定装置５の識別情報（所謂、ＩＤデータである）を設定し得る数字ボタンなどのデータ入力部５ｃと、これら入力された目標照度および識別情報を記憶するデータ記憶部５ｄと、上記光センサ部５ｂで得られた測定照度、設定された目標照度および識別情報を所定時間間隔（例えば、１秒間隔）でもってネットワーク回線７に出力するデータ出力部５ｅとから構成されている。勿論、データ出力部５ｅが無線ＬＡＮ７ｂに接続されている。上記カード本体部５ａには、データ入力部５ｃとしてのキーボードおよび目標照度、測定照度などを表示し得るデータ表示部５ｆ、表示データを切り替えるための切替ボタン５ｇなどが設けられている。なお、上記データ入力部５ｃとデータ記憶部５ｄとにより、目標照度設定部が構成される。
【００３２】
次に、照明制御装置４を図３に基づき説明する。
この照明制御装置４は、予め設定された複数の光度生成範囲（光度変化範囲ともいい、後述する）のうち、所定の光度生成範囲内で、照明機器３における光度をランダムに変化させて発光光度を生成するとともにこの変化された発光光度を順次出力し得る光度生成手段１１と、ネットワーク回線７を介して送られる識別情報、目標照度、測定照度、消費電力などの各種データを取り込むデータ取込手段１２と、このデータ取込手段１２にて取り込まれたデータを所定回数分（少なくとも、演算に必要とする回数分を超える回数）を記憶（格納）し得るデータ記憶手段１３と、このデータ記憶手段１３に記憶されたデータを取り込むとともにこのデータに基づき照度測定装置５が移動しているか否か、つまり利用者が移動したか否か（移動中であるか否かも含む）を検知する移動検知手段１４と、ランダムに変化される発光光度の変化量とそれに対応して得られた測定照度の変化量とに基づき、各照明機器３の照度測定装置５に対する照度の影響度合い（距離的位置関係ともいえる）を判断するための回帰係数を演算する回帰係数演算手段１５と、この回帰係数演算手段１５で求められた回帰係数および上記照度測定装置５により得られた照度状態（具体的には、目標照度に対する測定照度の偏差などである）に基づき、次に生成する発光光度（以下、次光度ともいう）の範囲を、予め設定されている複数の光度生成範囲から選択する光度生成範囲選択手段１６と、上記光度生成手段１１にて出力された次光度およびこの次光度に基づく発光による測定照度並びに照明機器３の消費電力に基づき、当該次光度の良否を評価するための評価値を評価式（目的関数ともいう）を用いて演算する評価値演算手段１７とから構成されている。
【００３３】
次に、回帰係数演算手段１５を図４に基づき説明する。
この回帰係数演算手段１５の機能を概略的に説明すると、照明機器３と照度測定装置５との因果関係を光度変化量［説明変数（入力値でもある）ｘ］と照度変化量（観測値ｙ）とで表される単回帰式を下記（１）式にて示すように定式化し、この単回帰式を説明する回帰係数（β）を求めることで、両者の影響度合いを得るようにしたものである。
【００３４】
ｙｉ＝εｉ＋α＋βｘｉ  ・・・（１）
この回帰係数演算手段１５は、回帰係数に対する補正値を求める推定誤差共分散の逐次演算と、この補正値を用いて推定値である回帰係数の逐次演算とを有し、入力値である光度変化量を基準としたときの観測値である照度変化量の観測予測誤差を最小とする推定値である回帰係数を、再帰的に推定する逐次型最小二乗法を用いるもので、
下記（２）式に基づき、或る時刻（演算ステップに相当する）（ｋ）における入力値ｘ（ｋ）、前の時刻（前回の演算ステップに相当する）（ｋ−１）における推定値の誤差共分散Ｐ（ｋ）および観測雑音ｗの共分散Ｒ（ｋ）を用いて、観測値の予測誤差共分散Ｓ（ｋ）を演算する観測予測誤差共分散演算部２１と、下記（３）式に基づき、入力値ｘ（ｋ）、推定誤差共分散Ｐ（ｋ−１）および観測予測誤差共分散Ｓ（ｋ）を用いて、補正値である補正係数Ｗ（ｋ）を演算する補正係数演算部２２と、下記（４）式に基づき、推定誤差共分散Ｐ（ｋ−１）、補正係数Ｗ（ｋ）および観測予測誤差共分散Ｓ（ｋ）を用いて、或る時刻（ｋ）における推定誤差共分散Ｐ（ｋ）を演算する推定誤差共分散演算部２３と、下記（５）式に基づき、入力値ｘ（ｋ）および前の時刻における推定値β（ｋ−１）を用いて、或る時刻における観測値の予測値＊ｙ（ｋ）を演算する観測予測値演算部２４と（なお、記号「＊」は、以下に示す式中において、この後に続く文字にハット記号が付されていることを示す、以下同じ）、下記（６）式に基づき、或る時刻（ｋ）における観測値ｙ（ｋ）および観測予測値＊ｙ（ｋ）を用いて、観測値の予測誤差ε（ｋ）を演算する観測予測誤差演算部２５と、下記（７）式に基づき、前の時刻における回帰係数β（ｋ−１）、補正係数Ｗ（ｋ）および観測予測誤差ε（ｋ）を用いて、或る時刻における推定値＊β（ｋ）を演算する推定値演算部２６とから構成されている。
【００３５】
Ｓ（ｋ）＝ｘ（ｋ）Ｐ（ｋ−１）Ｓ（ｋ）＋Ｒ（ｋ）  ・・・（２）
Ｗ（ｋ）＝Ｐ（ｋ−１）ｘ（ｋ）Ｓ−１（ｋ）＋Ｒ（ｋ）  ・・・（３）
Ｐ（ｋ）＝Ｐ（ｋ−１）−Ｗ（ｋ）Ｓ（ｋ）Ｗ（ｋ）  ・・・（４）
【００３６】
【数７】
2008243389.tif 000009
【００３７】
【数８】
2008243389.tif 000010
【００３８】
【数９】
2008243389.tif 000011
【００３９】
以下、この回帰係数演算手段１５での演算手順について詳しく説明する。
なお、以下の説明では、所定時間間隔おきに（つまり、所定周期毎に）データが取り込まれるとともに、これらデータおよび前回の回帰係数に基づき、次の回帰係数を推定値として求めるのであるが、この求める手順を「探索」と称して説明する。さらに、照度測定装置５を照度センサと称して説明する。
【００４０】
すなわち、逐次型最小二乗法では、現在の探索（時刻ｋ）において新たな観測値ｙ（ｋ）が得られると、まずその探索におけるシステムへの入力値ｘ（ｋ）、観測予測誤差共分散Ｓ（ｋ）、補正係数であるフィルタゲインＷ（ｋ）、および観測時に環境に生じる観測雑音Ｗの共分散であるＲ（ｋ）を用いて推定値の誤差共分散Ｐ（ｋ）を逐次に演算する。なお、時刻（ｋ）における推定誤差共分散Ｐは、時刻（ｋ）の入力値ｘ（ｋ）、時刻（ｋ−１）における推定誤差共分散Ｐ（ｋ−１）、およびＲ（ｋ）により算出される観測予測誤差共分散Ｓ（ｋ）と、その観測予測誤差共分散Ｓ（ｋ）、ｘ（ｋ）、およびＰ（ｋ−１）に基づき算出されるフィルタゲインＷ（ｋ）を用いて再帰的に推定する。推定値β（ｋ）の逐次演算では、時刻（ｋ−１）における推定値β（ｋ−１）およびｘ（ｋ）を用いて演算した観測予測値ｙ（ｋ）を用いて観測予測誤差ε（ｋ）を算出し、その観測予測誤差ε（ｋ）と既に求められている上記フィルタゲインＷ（ｋ）を用いて再帰的に推定する。なお、Ｒ（ｋ）は、観測値ｙ（ｋ）に重み付けを行う。具体的には、観測雑音ｗの分散が小さい時は、Ｒ（ｋ）も小さくなり観測値ｙ（ｋ）を重視するように働く。このため、Ｒ（ｋ）を含む観測予測誤差共分散Ｓ（ｋ）、およびＳ（ｋ）を含むフィルタゲインＷ（ｋ）も同様に重みの役割を担う。観測予測誤差共分散Ｓ（ｋ）とフィルタゲインＷ（ｋ）の間には、観測予測誤差ε（ｋ）が大きくなるとＳ（ｋ）は大きく、その影響を受けてＷ（ｋ）が小さくなるという関係にある。これらの重みが補正値（フィルタ）として用いることで、観測予測誤差ε（ｋ）を小さくし、そして探索毎に再帰的に推定値の補正を行うようにしているため、推定精度を向上させることができる。
【００４１】
各照明機器３は、上述したアルゴリズムを探索毎、全ての照度センサ５に対して行う。探索開始時を時刻（ｋ）とすると、探索開始時には、時刻（ｋ−１）の推定値β（ｋ−１）および推定誤差共分散Ｐ（ｋ−１）が必要となる。探索開始時、この推定値には適当な値を与えるが、その推定値は全くの未知であるため理論値との誤差は大きいと考えられる。そのため、推定値に重きを置かないように、推定値の誤差共分散Ｐ（ｋ−１）の初期値に十分に大きな値を与える。これにより、探索開始時における推定値の誤差に対処することができる。
【００４２】
ここで、推定した推定値βについて考えると、上述したように、これは回帰式の傾きであり、光度変化量と照度変化量の因果関係を説明する回帰係数である。すなわち、回帰係数βの大きさが照明の光度変化による各照度センサの照度変化量の大きさ（感度）を表しており、照明機器の照度センサに対する影響度合いを表していると考えることができる。なぜなら、照度センサの照度と照明機器の光度との間には、逆二乗の法則が成り立つからである。すなわち、理論的に照度センサの照度は照明機器の光度に比例し、距離の二乗に逆比例するという関係にある。つまり、推定した回帰係数βは照明機器と照度センサの距離が近いほど大きくなり、遠いほど小さくなる。なお、照明機器の照射方向に指向性がある場合、または照明を遮るものが存在する場合には、回帰係数の値が小さくなるが、基本的には、距離的位置関係を表しているものといえる。
【００４３】
ここで、回帰係数と、照明機器３および照度センサ５の位置関係とについて、図５に示す照明環境を用いて詳しく述べる。
図５に示される照明環境は、所定間隔おきに配置された３つの照明機器３（３Ａ，３Ｂ，３Ｃ）のうち、２つの照明機器３Ａ，３Ｂのそれぞれ下方に照度センサ５（５Ａ，５Ｂ）が配置されたものである。
【００４４】
図６に各照明機器３の各照度センサ５に対する回帰係数の探索結果の履歴を示す。図６の縦軸は回帰係数、横軸は探索回数である。図６（ａ）は各照度センサ５Ａ，５Ｂに対する照明機器３Ａの回帰係数を示したもので、照明機器３の真下にある照度センサ５Ａの回帰係数が高く、離れた位置にある照度センサ５Ｂの回帰係数は低くなっている。しかし、図６（ｂ）に示す照明機器３Ｂについては、照明機器３Ａとは逆に照度センサ５Ｂの回帰係数が高く、照度センサ５Ａの回帰係数が低くなっている。このように、照明機器３に、より近い照度センサ５の回帰係数の方が高くなっている。また、図６（ｃ）より、照明機器３Ｃでは、両照度センサ５Ａ，５Ｂに対する回帰係数は低いが、照明センサ５Ｂ、照明センサ５Ａの順で低くなっていることが判る。図５の照明環境から、これは照明機器３と各照度センサ５との距離の順番（近い順）になっていることが判る。このように、回帰係数は照明機器３の光度変化に対する各照度センサ５の感度として捉えることができる。したがって、照明機器３が各照度センサ５に対してそれぞれ推定した回帰係数を、それ自身に対す各照度センサの距離的位置関係を表す指標として用いることができる。
【００４５】
なお、回帰係数の演算には、光度変化量および照度変化量が用いられるが、勿論、各探索毎に、前回との発光光度の差、および測定照度の差が、それぞれ求められていることはいうまでもなく、図示しないが、これらの差を求める演算部が照明制御装置４に具備されている。
【００４６】
次に、上記光度生成範囲選択手段１６を図７〜図９に基づき説明する。
この光度生成範囲選択手段１６は、回帰係数を用いて、照明機器３と照度測定装置５との近傍度を検出するとともに、この近傍度と目標照度に対する測定照度の偏差の程度を示す照度状況とから得られる選択条件（以下、選択条件を決定することを近傍設計ともいう）により、予め、設定されている複数の光度生成範囲、例えば図８に示す７種類（Ａ〜Ｇ）から使用すべき、すなわち最適な次光度の生成範囲を選択するためのものである。
【００４７】
この光度生成範囲選択手段１６は、図７に示すように、求められた（推定された）回帰係数ＲＣを、下記の第１表に示すように、３つの閾値Ｔ（high,mid,low）（実験にて求められた値である）を用いて、４段階（ランク１〜３およびランク無し）にランク付けをするランク付け部３１と、このランク付け部３１で得られたランク値および照度測定装置５の照度状況［具体的には、目標照度、変動許容値（上限＋６％、下限−８％）に対する偏差量］に基づき選択条件を求める選択条件決定部３２とから構成されており、最終的には、ランク値と選択条件とに基づき、生成範囲（近傍）が決定される。つまり、近傍設計が行われる。
【００４８】
【表１】
2008243389.tif 000012
【００４９】
近傍設計を詳しく説明すると、照度測定装置５に対する照明機器３の近傍度を示すランクとその照度状況から当該照明機器３が必要とする次光度生成の傾向（増光、減光、調節など）を決定した後に、さらにこの傾向に対して光度を微調整することをいう。近傍設計で決定される生成範囲（近傍範囲）は、下記の第２表に示すように６種類設けられており、図９に示すツリー構造の分類規則により、それぞれ照明機器３に対し、６種類のいずれかが対応付けされる（分類される）。
【００５０】
【表２】
2008243389.tif 000013
【００５１】
ここで、上記光度生成範囲選択手段１６での選択手順、すなわち近傍設計について詳しく説明する。なお、ここでも、照度測定装置を照度センサと称して説明する。
すなわち、図８に示すように、７種類の光度生成範囲つまり近傍範囲Ａ〜Ｇが次光度生成に用いられる。近傍範囲Ａおよび近傍範囲Ｂは減光を重視した範囲であり、近傍範囲Ａは現在の光度から急激に光度を下げることを重視した範囲、近傍範囲Ｂは近傍範囲Ａよりは低速であるが減光を重視した範囲である。また、近傍範囲Ｃ、ＤおよびＥは光度の調節を行う範囲であり、近傍範囲Ｃは減光傾向、近傍範囲Ｅは増光傾向、近傍範囲Ｄは現在の光度を維持するように光度調整を行う範囲である。また、近傍範囲ＦおよびＧは近傍範囲ＡおよびＢとは逆に増光を重視した範囲であり、近傍範囲Ｆは低速、近傍範囲Ｇは急激に増光するようにした範囲である。なお、図８中の数値は、照明機器３の最大点灯光度を１００（％）とした際の光度変化量（％）であり、これらの値は実験的に求められたものである。数値的に説明すると、光度生成範囲は、許容可視変動範囲［人が知覚し得ない範囲であり、目標光度の１．０６（＋６％）〜０．９２（−８％）の範囲］が望ましいが、近傍範囲によっては、その範囲を少し超えるものもあるが、やはり、人の知覚では殆ど気が付くことがないため、支障は生じない。すなわち、図８にて示した光度の生成範囲を、拡大許容可視変動範囲ということができる。
【００５２】
近傍範囲を決定する場合、回帰係数を用いて照明機器と各照度センサとの回帰係数ＲＣを４段階（ランク１，２，３およびランク無し）でランク付けを行い、選択条件およびそのランクに従って７種類から適応的に１つ選択する。
【００５３】
以下、近傍範囲の決定手順について説明する。
（１）回帰係数を用いて各照度センサに対して、その照度に基づきランク付けを行う。
（２）各照度センサの照度状況に応じて選択条件を決定する。
（３）照明機器は照度センサのランクと選択条件から近傍範囲を一意に決定する。
【００５４】
まず、（１）項について説明すると、照明機器に対する照度センサのランク付けは、推定した回帰係数ＲＣを、閾値Ｔｈｉｇｈ、閾値Ｔｍｉｄ、および閾値Ｔｌｏｗの３つの値で比較して行う。これらの閾値を用いたランク付けの判断基準は第１表に示す通りである。なお、「ｉ」は照度センサの番号を示す。第１表より、最も高いランク１は、回帰係数ＲＣが閾値Ｔｈｉｇｈを超えた場合である。ランク２、ランク３およびランク無しについても、同様に、閾値Ｔｈｉｇｈ、閾値Ｔｍｉｄ、および閾値Ｔｌｏｗを用いて判断する。照度センサに影響が有る・無しの判断は、回帰係数ＲＣが閾値Ｔｌｏｗを下回るか否かである。なお、これらの閾値は、予め、実験より求められた値である。
【００５５】
次に、（２）項および（３）項について説明する。
（２）項において、各照明機器は各照度センサの照度状況、および（１）で得た自身のランクに基づき選択条件を決定する。その際には、図９に示すツリー構造が用いられる。図９の下部に、照度状況（分岐ルール）を示しておく。
【００５６】
選択条件１は急速な減光を行う範囲である。この選択条件１は、ランク無し、つまり照度センサに影響のない照明機器に用いられる。選択条件２は減光傾向を有する減光近傍群、選択条件３は急速な減光傾向を有する急速減光近傍群の照明機器に用いられる。これら選択条件２および３は、照度センサが目標照度を上回る場合に用いられ、上回る程度によって使い分けられる。また、選択条件４は急速な増光傾向を有する急速増光近傍群、選択条件５は増光傾向を有する増光近傍群の照明機器に用いられる。これら選択条件４および５は、照度センサが目標照度を下回る場合に用いられ、下回る程度によって使い分けられる。目標照度を上回るまたは下回る程度は、上述した許容可視変動範囲の上限である１．０６（＋６％）および下限である０．９２（−８％）を基準にして決定される。なお、目標照度を下回る場合には、下回る度合いが許容可視変動範囲の下限である０．９２に達する前に対処する必要がある。そこで、目標照度を下回る前に増光動作を行うように許容可視変動範囲の下限よりも高い閾値（０．９８）を用いて選択条件を決定する。また、目標照度を下回る場合には、人間が全く知覚できないほど微量に目標照度を下回ることが考えられる。この場合、その目標照度を下回る程度を誤差範囲と考えることができるが、選択条件４または５によって急速に増光されると略目標照度に収束している照度が荒れる危険性がある。そこで、選択条件６のような光度を調整するための近傍群が用意されている。選択条件６では、最も近いランク１の照明には、少しだけ増光傾向を持つ近傍範囲Ｄを与えて目標照度を満たすような動きをさせ、ランク２の照明には、調整する近傍範囲Ｂを与えて現状を維持するような動作をさせる。また、略目標照度を満たすような場合には、ランク３のような照明は、あまり必要でないと考えられるため、少し減光を行う近傍範囲Ｃを与え、電力の最小化を行うように動作させる。これにより急激な増光を避けることができ、目標照度付近での照度の安定を期待することができる。（３）項では、第２表に基づき（２）項で決定された選択条件からランクに応じて近傍範囲を一意に決定する。以下に、（２）項および（３）項の流れについて詳しく述べる。
【００５７】
まず、（２）においては、各照度センサの照度データとランクとを用いて、図９に示すツリー構造のルート１（ルートノード）から条件分岐を辿り、選択条件を決定する。条件分岐は全部で７種類ある。ノード１は、照明の影響のある照度センサが存在するかの条件分岐である。ここで影響のある照度センサが存在しなければ、選択条件１（ランク無し）を選択するが、存在する場合は条件分岐２に進む。ノード２は、自身が影響する照度センサが複数存在するかの条件分岐である。ここで、複数の照度センサに影響しなければ条件分岐３へ、影響すれば条件分岐４へと進む。ノード３は、影響する照度センサが目標照度を満たしているかの条件分岐であり、目標照度を満たしていれば条件分岐５へ、満たしていなければ条件分岐６に進む。ノード５は、現在照度が許容可視変動範囲内か否かを判断する条件分岐で、すなわち目標照度収束範囲内であれば選択条件２を選択し、収束範囲外の場合は選択条件３を選択する。また、ノード６は、目標照度を満たさない度合いが誤差と見なせる範囲であるか否かを判断する条件分岐であり、誤差範囲であれば選択条件６を選択し、誤差範囲外であれば条件分岐７に進む。なお、誤差範囲は、殆どの人が明るさの変動を感じない範囲［例えば、目標照度から２％（下限）］までとする。ノード７は、現在照度が許容可視変動範囲の下限より以下であるか否かを判断する条件分岐であり、目標照度収束範囲外であれば選択条件４を選択、範囲内であれば選択条件５を選択する。このように、ノード２から分岐する左側のツリーは、影響する照度センサが１個の場合に選択され、目標照度を満たすか否かで増光傾向または減光傾向を持つ選択条件が選択される。
【００５８】
次に、条件分岐４を進んだ場合について説明する。ノード４は、影響する全ての照度センサが目標照度を満たすか否かを判断する条件分岐である。ここで、全ての照度センサが目標照度を満たせば、条件分岐５へ進み、選択条件２または３が選択される。一方、影響する照度センサの中で１つでもその目標照度を満たさない場合は、条件分岐６へ進み、選択条件４，５または６が選択される。このように、ノード２から分岐する右側のツリーでは、影響する照度センサが複数の場合に選択され、その照度センサが目標照度を満たすか否かで増光傾向を持つまたは減光傾向を持つ選択条件が選択される。
【００５９】
以上のように、各照明機器（正確には、照明制御装置）は、取得したランクに基づき選択条件を決定し、自身に最も近い、すなわちランクが最も高い照度センサの照度状況に合わせ、第２表から近傍範囲を一意に決定する。但し、複数の照度センサに影響し、その照度センサのいずれかが目標照度を下回っている場合には、目標照度を満たさない照度センサの中で最もランクが高い照度センサのランクに応じて近傍範囲を一意に決定する。これにより、目標照度を満たす必要がある照度センサのための次光度生成が可能となる。なお、各照明機器（照明制御装置）が全ての照度センサに対して上述した手順を適用し、それぞれ適切な近傍範囲を決定する。
【００６０】
次に、上記評価値演算手段１７を図１０に基づき説明する。
この評価値演算手段１７は、得られた次光度が良好であるか不良であるかを判断するもので、つまりシステム全体の消費電力量と、目標照度に対する各照度測定装置５で得られる測定照度の偏差（値の大小）とを考慮して判断するためのもので、下記（８）式にて得られる評価値ｆに基づき判断される。
【００６１】
例えば、前回の回帰係数による評価値と、今回推定された回帰係数による評価値とを比較して、値が小さくなっていれば、システム全体の光度、言い換えれば各照明機器での光度が良好な方向に収束したことを意味しており、つまり、測定照度が目標照度に近づくとともに消費電力も少なくなるように制御されていることを示している。
【００６２】
【数１０】
2008243389.tif 000014
【００６３】
すなわち、上記（８）式から判るように、評価値ｆは、消費電力Ｐに、現在の測定照度Ｌｃと目標照度Ｌｔとの照度差である偏差に依存する照度依存量ｇｊを加算したものである。また、消費電力Ｐには、各照明機器３の発光光度Ｃｄｉの合計値を用いる。このように光度値を用いるのは、蛍光灯の電源電圧特性上、光度と電力とが線形関係にあるからである。但し、通常、消費電力Ｐとしては、電力計６にて測定された値が用いられる。発光光度Ｃｄｉの合計値は、電力計６による消費電力が得られない場合に用いられる。
【００６４】
そして、照度依存量ｇｊは、照度差が負である場合、または照度差が現在照度の６％以上の場合に加算されるもので、ペナルティとしての役割を持つ。ここで、現在照度に加算される６％は、人間が明るさ変動を知覚するかどうかの閾値（許容可視変動範囲の上限）である。なお、ここでは、許容可視変動範囲の上限（６％）とその下限（−８％）との間を目標照度収束状態とする。
【００６５】
このように、目標照度収束状態の上限値以上、または現在照度が目標照度を下回る場合は、ペナルティを与えることで評価値を上昇させ、急速に増光または減光させるようにしている。
【００６６】
また、照度依存量ｇｊには、照度差に対して回帰係数ＲＣｊが乗算される。なお、回帰係数ＲＣｊが閾値以下である場合には、「０（ゼロ）」が乗算される。すなわち、測定照度が目標照度を満たさない場合でも、その照度センサに対する影響度合いが低ければ、評価値にペナルティを与えない。これにより、影響度合いが高い照度センサ、すなわちその近くに位置する照度センサに最適化の対象を絞ることができ、目標照度を満たす精度を向上させることができる。また、照度依存量ｇｊには、重み係数ｗが乗算されており、この値を変化させることにより、目標照度を満足させることを優先するか、電力の最小化を優先するかを選択することができる。
【００６７】
したがって、上記上記評価値演算手段１７は、図１０に示すように、回帰係数ＲＣｊを入力して閾値Ｔと比較して評価すべき照度測定装置（評価対象）を絞るための選択指標（ここでは、回帰係数に相当する）Ｒｊを出力する選択指標決定部４１と、発光光度と測定照度との照度差に基づく照度依存量ｇｊを求める依存量演算部４２と、この依存量演算部４２で求められた値（照度測定装置が複数ある場合には、その合計値）を入力するとともにこの値に重み係数ｗを乗算したものに全ての照明機器の消費電力の合計値Ｐを加算する加算部４３と、この加算部４３にて得られた評価値ｆと前回の評価値（例えば、データ記憶手段に記憶されている）ｆとを比較して値が改善しているか否かを判断する、すなわち大きくなっているか否かを判断する評価値判断部４４とから構成されている。
【００６８】
そして、（８）式の但し書きの箇所で説明したように、選択指標Ｒｊは、得られた回帰係数ＲＣｊの値が閾値Ｔと等しいかまたは超える場合には、回帰係数ＲＣｊの値が用いられ、閾値未満である場合には「０（ゼロ）」が用いられる。
【００６９】
また、照度依存量ｇｊについては、偏差（偏差量）が負の場合、または偏差が現在照度の６％以上の場合には、偏差の二乗に、選択指標（つまり、回帰係数の値）Ｒｊが乗算されたものが用いられ、偏差が上記以外の範囲である場合には、「０（ゼロ）」が用いられる。
【００７０】
次に、照度測定装置５の移動の有無を検知する移動検知手段１４について、図１１および図１２に基づき説明する。
この移動検知手段１４は、図１１に示すように、取得した測定照度が、図１２に示す許容下限値以上で許容上限値以下の範囲（イ）、許容上限値超える範囲（ロ）、許容下限値未満の場合（ハ）の３つのいずれに属するか否かを判断する照度値判断部５１と、この照度値判断部５１にて上下限から外れていると［つまり（ロ）と（ハ）の場合］判断された場合には、その継続回数（継続探索回数）をカウントするカウント部５２と、このカウント部５２からのカウント数を入力して設定回数と比較して超える場合には移動と判断する移動判断部５３と、この移動判断部５３にて移動中であると判断された場合に、回帰係数および推定誤差共分散に外れ値を与えて回帰係数の推定演算をリセットするリセット部５４とから構成されている。
【００７１】
この移動検知手段１４での検知手順について、詳しく説明する。
すなわち、本発明に係る回帰係数を用いた照明制御アルゴリズム（自律分散型適応的近傍アルゴリズムともいう）においては、照度測定装置５の移動を照明機器３自ら検知する必要がある。ここで、実際の照明環境下における照度測定装置５の移動について考えると、照度測定装置５の移動発生時には、その照度が一時的に大きく変化する。そこで、照明機器３は、移動検知のために、探索過程において、常に全ての照度測定装置５の照度状態を監視している。なお、照度状態の監視は、図１２に示すように目標照度範囲内［照度範囲（イ）］および目標照度範囲を超えるおよび未満の２つの状態［照度範囲（ロ）および（ハ）］の３段階で行う。図１２の縦軸は、目標照度を基準とした照度変動率であり、その上下限範囲内は目標照度収束状態を示している。また、照度範囲の監視には、環境雑音（ノイズである）すなわち外的要因で生じる照度変化に注意が必要となる。例えば、照度測定装置５の光センサ部５ｂに人影が入るといった外的要因が考えられる。そこで、大きく照度変化が起こったことを感知した後、数回の探索を移動検知期間として設け、移動の有・無を確認する。具体的には、移動検知期間中に目標照度を満たさない探索回数をカウントし、連続して照度範囲（ロ）または（ハ）の状態が続くようであれば、照度測定装置５が移動したと見なす。これにより、外的要因による照度変化については、すぐに反応（検知）しないようにする。なお、移動を検知した際には、上述したように、推定誤差共分散に外れ値を与えることによって推定をリセットし、移動後の環境においても、正しい推定を行い得るようにする。
【００７２】
次に、目標照度を満たすのが、難しい環境への対処について説明する。
照明機器３が蛍光灯である場合、その構造により光の強さに指向性が存在するため、照度測定装置５の配置によっては、照明機器３と照度測定装置５との位置関係の把握が難しく、目標照度を下回る場合が生じる。目標照度を満たすのが難しい照明環境を、図１３に示す。なお、図１３の説明においては、照明機器を単に「照明」と称し、また照度測定装置を「センサ」と称し、またその後方の番号は、機器番号を示すものとする。
【００７３】
この照明環境下において、目標照度を満たすためには、影響が強い照明２および７が最大点灯し、次に影響が強い照明１，３，６および８の４つの蛍光灯中、２つの蛍光灯が８０％以上で点灯する必要がある。この図１３において、回帰係数を推定した結果を図１４に示す。この図１４より、照明２および７の蛍光灯は高く、照明１，３，６および８の４つの蛍光灯は低く、回帰係数を推定しており、その位置関係および距離は把握できている。しかし、照明１，３，６および８の４つの蛍光灯の回帰係数では、閾値Ｔｌｏｗを超えない場合がある。この場合、センサＡに対しランクを取得できず、近傍範囲Ａにより減光される。そのため、目標照度を満たすために必要な照明１，３，６および８の蛍光灯がセンサＡに対して光度を提供することができず、その他のランクを持つ照明だけでは目標照度を満たせない状態に陥る。
【００７４】
そこで、このような照明環境下においても、常にランクが得られるよう閾値Ｔｌｏｗを可変にする。これは、上述したように、センサの配置により適切な閾値が異なると推測されるためである。具体的には、上記と同様に、照度範囲を図１２の３段階で常に監視し、照度範囲（ハ）が連続して続く場合、閾値Ｔｌｏｗを下げる。また、閾値Ｔｌｏｗを下げた後に照度範囲（イ）で安定する場合は、閾値Ｔｌｏｗを元に戻す。これにより、必要な照明にランクを与えることができ、目標照度を満たすように増光し得る。
【００７５】
以下、上述した照明システムにおける自律分散型の適応的近傍アルゴリズムを用いた照明制御方法を、図１５のフローチャートに基づき簡単に説明する。なお、ここでの説明は全体的な流れについて簡単に行うものとし、詳細な手順については、上述した各構成手段にて説明しているため、省略する。
【００７６】
また、この照明制御方法は、上述したように、照明機器に対する照度測定装置での照度に与える影響度合いを把握するために、照明機器での光度変化量（時系列データ）と照度測定装置での照度変化量（時系列データ）との関係を示す単回帰式の回帰係数を用いるとともに、この回帰係数を求めるに際し、逐次型最小二乗法を用いたものであり、したがって或る時刻におけるデータによる演算時、すなわち或る探索ステップにて回帰係数を求めるのに前回の探索ステップでの回帰係数を必要とするため、システム起動時においては、必要とされる前回の値については、初期値が与えられている。また、以下の説明においても、照度測定装置を照度センサとして説明する。
【００７７】
システムが起動されると、初期光度、初期回帰係数などの初期パラメータが設定され（ステップ１）、そして各照明機器３は初期光度でもって点灯される（ステップ２）。
すなわち、このステップ２では、光度生成手段１１により、初期設定により設定された初期光度が生成されて照明機器３が点灯されるとともに、所定の光度生成範囲（例えば、光度生成範囲Ｄ）で、且つ所定時間間隔（例えば、１秒間隔）でもって、ランダムにその光度が変化される。具体的には、コンピュータ装置の演算処理部（マイクロチップ）で乱数（正確には、擬似乱数）を発生させることにより行われる。ここで発生された光度は、他の照明機器とは全く無関係に決定される。
【００７８】
なお、このとき、希望する目標照度を設定した移動式の照度測定装置（以下、照度センサと称す）５を携帯した利用者が照明機器（正確には、蛍光灯３ａである）３の下方の任意の場所に移動し、例えば図書室の或る場所に着席しているものとして説明する。
【００７９】
照明制御装置４は、ネットワーク回線７ａ，７ｂを介して、データ取込手段１２により、電力計６からの消費電力、照度センサ５からの目標照度および測定照度を、所定時間間隔（例えば、１秒間隔）でもって取得し、データ記憶手段１３に時系列データとして記憶する（ステップ３）。勿論、照明機器３に対する発光光度は、当該照明制御装置４にて把握されている。なお、少なくとも、回帰係数を演算する際に必要となる個数分のデータが記憶される。
【００８０】
次に、取得したデータに基づき、移動検知手段１４にて、照度センサ５が移動しているか否かが判断される（ステップ４）。
上記ステップ４で、照度センサ５の移動が検知されると、評価値を求めるために用いる回帰係数を、前回の回帰係数とするようにリセットされる（ステップ５）。
【００８１】
一方、ステップ４で、移動していないと判断された場合には、ステップ６に進む。
次に、評価値演算手段１７にて、取得した発光光度、照度センサの識別情報、目標照度、測定照度および消費電力を用いて評価値が演算される（ステップ６）。
【００８２】
次に、光度生成範囲選択手段１６にて、次に発光すべき光度生成範囲が選択されるとともに、この範囲内で、光度生成手段１１により次光度が生成される（ステップ７）。
次に、次光度に基づく測定照度、消費電力などが取得され（ステップ８）、そして回帰係数演算手段１５により新たに回帰係数が演算される（ステップ９）。
【００８３】
次に、評価値演算手段１７にて、再度、取得した発光光度、照度センサの識別情報、目標照度、測定照度および消費電力に基づき、次光度を評価するための評価値が演算される（ステップ１０）。
【００８４】
次に、ステップ１０で演算された評価値に基づき、前回の評価値より改善されているか否かが、すなわち次光度を受理するか否かが判断される（ステップ１１）。
そして、ステップ１１で改善されていると判断された場合には、演算終了判断部（例えば、演算を強制的に終了させるもので、無くてもよい）を介して（ステップ１３）、上記ステップ３に戻り、引き続き、上記各ステップ３〜１３が実行される。
【００８５】
一方、ステップ１１で改善されていないと判断された場合には、ステップ１２に進み前回の値、すなわち前回の回帰係数を採用した後、上記ステップ３に戻り、引き続き、上記各ステップ３〜１３が実行される。
【００８６】
ここで、上述した本発明に係る自律分散型で且つ回帰係数を用いた適応的近傍アルゴリズムを、実際に実験室に構築した照明システムに適用した場合について説明する。
図１６に示すように、部屋の天井に、照明（以下、１〜１５の番号を付して説明する）を１５個配置するとともに、３個の照度センサ５Ａ，５Ｂ，５Ｃが配置されたものであり、照度センサ５Ａの目標照度は７５０ルクスに、照度センサ５Ｂの目標照度は７００ルクスに、照度センサ５Ｃの目標照度は８００ルクスに設定されているものとする。このときのパラメータは、下記の第３表に示す通りである。なお、以下の説明では、１５個の照明機器３を区別するのに、単に、「照明」という語句の後に、その番号を付すとともに、３個の照度センサについては、「センサ」という語句の後に、Ａ〜Ｃの文字を付す。
【００８７】
【表３】
2008243389.tif 000015
【００８８】
実環境試験１（照明環境に変化がない場合）および実環境試験２（センサが移動した場合）を行った際の照度履歴を図１７に、消費電力履歴を図１８にそれぞれ示す。また、図１９に移動前のセンサＡに影響する照明７，８，９，１２，１３および１４の、図２０に移動後のセンサＡに影響する照明１，２，６および７の光度履歴および回帰係数をそれぞれ示す。図２１にセンサＢに影響する照明３および４の、図２２にセンサＣに影響する照明１３，１４および１５の、図２３に全てのセンサに影響しない照明５，１０および１１の光度履歴および回帰係数をそれぞれ示す。
【００８９】
図１７（ｂ）より、探索開始から各センサの照度は目標照度へと収束していき、探索回数３３回で各センサＡ，ＢおよびＣの照度は７４７，６７３および８４６ルクスになっている。例えば、相関係数を用いたものと比較すると、約１／２倍の探索回数で目標照度収束範囲へと到達していることが確認できる。また、相関係数を用いたものと比較して、照度の変動が少なく安定していることが確認できる。また、図１８より、探索が進むにつれて消費電力が減少し、照度履歴と同様、相関係数を用いたものより安定して電力削減を図り得ることが判る。また、光度履歴は、実環境であるため少々荒れるが、数値実験を行った場合と似た軌跡を取っていることが判る。すなわち、各照明は、自身に近い照度センサを回帰係数により判断し、その照度状況に応じて次光度を生成している。例えば、実環境試験１においては、照度センサに近い照明（照明３，４，８，１３および１４）は非常に高い光度で点灯しており、照度センサから遠い照明（照明１，２，５，６，１０および１１）は略最小点灯光度となっている。なお、相関係数を用いた場合の光度履歴と比較すると、光度の荒れが少なかった。すなわち、回帰係数を用いることにより、照度センサとの距離感を得ることができており、その距離感および照度状況に応じて次光度を与えているためである。これに対して、相関係数を用いる場合には、相関があるか否かおよび相関がある照度センサの照度状況だけで次光度を与えているため、距離感を無視した大雑把な次光度生成しか行うことができない。例えば、移動前のセンサＡおよびセンサＣに影響する照明１４では、真下に影響するセンサＣがあるにも拘わらず、少し離れた場所にあるセンサＡが目標照度を下回った場合には近傍範囲ＦやＧのように急速な増光を行う近傍範囲を選択して急激に増光を行うといったように、従来の制御方式では、自身が光度を提供するのが最も適していない場合でも、相関が有るか無いかで一律の行動を取ることがある。このため、光度が荒れてしまう。
【００９０】
また、図２３において、照明１０は、探索初期においてセンサＡに対する回帰係数が誤って高くなっていると予想できる。このため、探索初期において、本来必要のない照明１０が高い光度で点灯してしまう。しかし、探索回数２００を超えるとセンサＡに対する照明１０の回帰係数は低くなり、減光する様子が確認できる。このように、推定の初期段階では十分な逐次推定がなされていないため誤差を多く含むことがあるが、探索が進むにつれ正しく推定できることが確認できる。
【００９１】
次に、照度センサ（センサＡ）の移動時について説明する。
図１７（ｃ）より、探索回数１０００回目のセンサＡの移動により大幅に目標照度が下回るが、探索回数１０５４回でセンサＡ，ＢおよびＣの照度は、７４４，７４３および８０３ルクスとなり、相関係数を用いた場合の約１／２倍の探索回数で全ての照度センサが目標照度収束範囲に達することが確認できる。すなわち、センサＡの移動前に影響していた照明が即座に減光し、移動後に影響する照明が急速に増光している。これは、回帰係数の推定により、照度センサの照明機器に対する照度の影響度合い（距離的位置関係）が正しく把握できているためである。しかし、図２３（ａ）より、照明５の光度が一時的に増光していることが確認できる。これは、図２３（ｃ）より、センサＡの移動直後に行う推定値である回帰係数のリセットにより、一時的に誤ったランクを得てしまった結果である。勿論、探索を進めると、この誤りを訂正することができる。
【００９２】
最後に、センサＡの移動前の定常状態の点灯状況を図２４に、センサＡの移動後、十分に探索が進んだ定常状態における点灯状況を図２５にそれぞれ示す。図２４および図２５より、各照度センサは目標照度収束範囲内の照度であり、各照度センサの近くにある照明が明るく点灯し、照度センサに影響のない照明が最小点灯光度付近で点灯していることが確認できる。
【００９３】
上述した照明制御アルゴリズムによると、照明機器と照度センサの位置関係を、相関係数を用いた場合よりも、精度良く把握することができ、つまり、照明の光度変化量と照度センサの照度変化量とを回帰分析することによって得られる回帰係数は、照明が照度センサに影響するかどうかだけでなく、その距離感まで表すことができる。これにより相関係数より詳細な位置関係の把握が可能となり、目標照度を満たし省電力な状態へと素早く収束させることができる。
【００９４】
ここで、上記照明システムにおける照明制御方法をステップ形式にて（主なステップに基づき）記載した場合について、以下に示しておく。
すなわち、この照明制御方法は、複数の照明機器と、これら各照明機器に対応して配置されて各照明機器で発光させる発光光度を制御する複数の照明制御装置と、上記各照明機器にて消費される消費電力を計測する電力計と、上記照明機器による照度を測定し得る照度測定装置と、上記各照明制御装置、電力計および照度測定装置との間で、所定データの受け渡しを行い得るデータ通信手段とを具備し、
且つ上記複数の照明機器の影響下に位置する照度測定装置からの測定照度および目標照度に基づき、少なくとも、当該照度測定装置の測定照度に影響を与える照明機器に対して、測定照度が目標照度に近づくように、その発光光度を制御する際に、上記測定照度に影響を与える照明機器を判断する指標として、所定の光度生成範囲内で変化されて照明機器に与えられる発光光度の変化量および照度測定装置にて測定される測定照度の変化量に基づき演算される回帰係数を用いるようにした照明システムにおける照明制御方法であって、
上記各照明機器を、所定の光度生成範囲内で発光光度を変化させて点灯させる第１のステップと、上記照明機器の影響下にある照度測定装置から測定照度および識別情報並びに消費電力を取得する第２のステップと、取得した測定照度に基づき照度測定装置が移動しているか否かを判断する第３のステップと、上記第３のステップにて照度測定装置の移動を検知すると現在の回帰係数をリセットする第４のステップと、取得した発光光度、照度測定装置の識別情報、測定照度および消費電力を用いて現在の発光光度の評価を行うための評価値を演算する第５のステップと、予め設定された複数の光度生成範囲から、次に発光すべき光度生成範囲を選択する第６のステップと、上記第６のステップで選択された光度範囲内で変化される光度値、照度値および消費電力を取得する第７のステップと、上記第７のステップで取得されたデータに基づき、新たに回帰係数を演算する第８のステップと、上記第８のステップで求められた回帰係数、および上記第７のステップで取得された発光光度、照度測定装置の識別情報、測定照度および消費電力に基づき評価値を演算する第９のステップと、上記第９のステップで求められた評価値と上記第５のステップで求められた前回の評価値とを比較して当該評価値が改善されているか否かを判断する第１０のステップとを具備し、
且つ上記第１０のステップにて改善されていると判断された場合に最初のステップに戻るようにするとともに、改善されていないと判断された場合に、第６のステップで求められた発光光度を破棄して最初のステップに戻るようにした制御方法である。
【００９５】
ところで、上記実施の形態にて説明した各手段（例えば、光度生成手段、移動検知手段、回帰係数演算手段、光度生成範囲選択手段、評価値演算手段など）およびこれら各手段を構成する各構成部は、例えばプログラムにより実行されるものであり、それぞれ機能部と呼ぶこともでき、また場合によっては、それぞれの機能を発揮し得る回路部として構成することもできる。
【００９６】
さらに、上記実施の形態においては、回帰係数を求めるのに、逐次型最小二乗法を用いたが、例えばバッチ型最小二乗法を用いることもできる。このバッチ型最小二乗法は、発光光度の変化量および測定照度の変化量の時系列データを用いて、回帰係数をバッチ形式にて求めるものである。
【００９７】
また、上記実施の形態においては、回帰係数を求めるのに、照明機器の光度をランダム（任意に）に変化させたが、例えば各照明機器に対して互いに異なる所定の方法で変化させるようにしてもよい。すなわち、所定の変化パターンでしかも照明機器毎に異なる変化パターンでもって変化させるようにしてもよい。
【図面の簡単な説明】
【００９８】
【図１】本発明の実施の形態に係る照明システムの概略全体構成を示す図である。
【図２】同照明システムにおける照度測定装置のカード本体部の平面図である。
【図３】同照明制御装置の概略構成を示すブロック図である。
【図４】同照明制御装置における回帰係数演算手段の構成を示すブロック図である。
【図５】同回帰係数演算手段を説明するための照明環境を示す機器配置図である。
【図６】同照明環境での回帰係数の有効性を説明するグラフである。
【図７】同照明制御装置における光度生成範囲選択手段の構成を示すブロック図である。
【図８】同光度生成範囲選択手段における光度生成範囲を示すグラフ図である。
【図９】同光度生成範囲選択手段における近傍設計を説明するためのツリー構造を示す図である。
【図１０】同照明制御装置における評価値演算手段の構成を示すブロック図である。
【図１１】同照明制御装置における移動検知手段の構成を示すブロック図である。
【図１２】同移動検知手段における照度範囲を示す図である。
【図１３】同移動検知手段での移動検知を説明する照明機器の配置図である。
【図１４】同移動検知手段での移動検知時における光度の経過を示すグラフである。
【図１５】同照明システムにおける照明制御方法を説明するフローチャートである。
【図１６】同照明システムにおける実環境への適用実験を説明する機器配置図である。
【図１７】同適用実験における照度の収束状態の経過を示すグラフである。
【図１８】同適用実験における消費電力の経過を示すグラフである。
【図１９】同適用実験における照明センサＡの移動前における光度および回帰係数の経過を示すグラフである。
【図２０】同適用実験における照明センサＡの移動後における光度および回帰係数の経過を示すグラフである。
【図２１】同適用実験における照明センサＢにおける光度および回帰係数の経過を示すグラフである。
【図２２】同適用実験における照明センサＣにおける光度および回帰係数の経過を示すグラフである。
【図２３】同適用実験における影響を受けない照明センサにおける光度および回帰係数の経過を示すグラフである。
【図２４】同適用実験における定常状態での光度を示す機器配置図である。
【図２５】同適用実験における定常状態での光度を示す機器配置図である。
【符号の説明】
【００９９】
１      電力供給源
２      電気配線
３      照明機器
３ａ    蛍光灯
３ｂ    光度調節部
４      照明制御装置
５      照度測定装置
５ａ    カード本体部
５ｂ    光センサ部
５ｃ    データ入力部
５ｄ    データ記憶部
５ｅ    データ出力部
６      電力計
７      ネットワーク回線
１１      光度生成手段
１２      データ取込手段
１３      データ記憶手段
１４      移動検知手段
１５      回帰係数演算手段
１６      光度生成範囲選択手段
１７      評価値演算手段
２１      観測予測誤差共分散演算部
２２      フィルタゲイン演算部
２３      推定値誤差共分散演算部
２４      観測予測値演算部
２５      観測予測誤差演算部
２６      推定値演算部
３１      ランク付け部
３２      近傍選択部
４１      選択指標決定部
４２      依存値演算部
４３      加算部
４４      評価値判断部
５１      照度値判断部
５２      カウント部
５３      移動判断部
５４      リセット部
【図１】
2008243389.tif 000022
【図２】
2008243389.tif 000023
【図３】
2008243389.tif 000024
【図４】
2008243389.tif 000025
【図５】
2008243389.tif 000026
【図６】
2008243389.tif 000027
【図７】
2008243389.tif 000028
【図８】
2008243389.tif 000029
【図９】
2008243389.tif 000030
【図１０】
2008243389.tif 000031
【図１１】
2008243389.tif 000032
【図１２】
2008243389.tif 000033
【図１３】
2008243389.tif 000034
【図１４】
2008243389.tif 000035
【図１５】
2008243389.tif 000036
【図１６】
2008243389.tif 000037
【図１７】
2008243389.tif 000038
【図１８】
2008243389.tif 000039
【図１９】
2008243389.tif 000040
【図２０】
2008243389.tif 000041
【図２１】
2008243389.tif 000042
【図２２】
2008243389.tif 000043
【図２３】
2008243389.tif 000044
【図２４】
2008243389.tif 000045
【図２５】
2008243389.tif 000046