(57)【要約】  【課題】  各種プラント、装置の最適制御等に用いられる大域的最適化装置の操作を高速化し、速やかに最適値を得る。【解決手段】大域的最適化装置の操作方法であって、ポテンシャル関数2006059237.tif 000035の真の最小値2006059237.tif 000036の取り得る領域Ω全体をサーベイ探索する前に、2006059237.tif 000037を力学システムの座標と見做し、その量子力学的トンネル効果を使うことにより真の最小値2006059237.tif 000038を計算する大域最適化操作手順において、該力学システムのプランク定数に相当する部分に置換されるべき設計定数ＨＴと、力学システムの時間発展を規定する時間微分方程式における該力学システムの散免の変化率を規定する摩擦係数γとを具備することを特徴とする大域的最適化装置の操作方法。【選択図】  図１
2006059237.tif 000002
【特許請求の範囲】
【請求項１】
大域的最適化装置の操作方法であって、ポテンシャル関数
2006059237.tif 000027
の真の最小値
2006059237.tif 000028
の取り得る領域Ω全体をサーベイ探索する前に、
2006059237.tif 000029
を力学システムの座標と見做し、その量子力学的トンネル効果を使うことにより真の最小値
2006059237.tif 000030
を計算する大域最適化操作手順において、該力学システムのプランク定数に相当する部分に置換されるべき設計定数ＨＴと、力学システムの時間発展を規定する時間微分方程式における該力学システムの散免の変化率を規定する摩擦係数γとを具備することを特徴とする大域的最適化装置の操作方法。
【請求項２】
請求項１の設計定数ＨＴが、該力学システムに対する量子的揺らぎ力の作用がその2乗ＨＴ２に比例して計算算出されることを特徴とする大域的最適化装置の操作方法。
【請求項３】
請求項１の設計定数ＨＴが、ポテンシャル障壁を量子学トンネル効果で透過する割合が、ｅｘｐ（−定数／ＨＴ）に比例するような量として規定されていることを特徴とする大域的最適化装置の操作方法。
【請求項４】
請求項１の摩擦係数γが、該力学システムのエネルギ散免率の比例定数として規定されていることを特徴とする大域的最適化装置の操作方法。
【請求項５】
大域的最適化装置の操作方法であって、ポテンシャル関数
2006059237.tif 000031
の真の最小値
2006059237.tif 000032
の取り得る領域Ω全体をサーベイ探索する前に、
2006059237.tif 000033
を力学システムの座標と見做し、その量子力学的トンネル効果を使うことにより真の最小値
2006059237.tif 000034
を計算する大域最適化操作手順において、探索空間Ωに拘束条件が付帯している際はこの拘束条件を、該力学システムのプランク定数に相当する部分に置換されるべき設計定数ＨＴと、力学システムの時間発展を規定する時間微分方程式における該力学システムの散免の変化率を規定する摩擦係数γとともに入力し、この入力式に応じて該装置で交換関係の代数とそれに基づくハミルトニアン演算子が構成され波動方程式が作成され、これに基づき大域最適化演算が行われることを特徴とする大域的最適化装置の操作方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、各種プラント、装置の最適制御、シミュレーション、数値解析等の各種理工学的課題の解決に適用される大域的最適化技術に関するものである。
【背景技術】
【０００２】
局所的最少点を複数個有するような誤差関数に対する大域的最小点を量子トンネル効果を利用して検索するための計算手順が発表されている。ここでは遺伝アルゴリズム等の従来の最適化手法においては、局所的最小点にトラップされ易いことを避けるために量子力学的なトンネル効果を利用している。
【０００３】
すなわち、局所的最小点にトラップされた場合でも、物理システムはいずれは真のエネルギ最小点に遷移するので、時間さえ掛ければ大域的最小点にトンネル効果を利用して遷移することができる。このアルゴリズムを用いれば遺伝アルゴリズムのような従来の最適化では真の最小点を得ることが不可能な場合でも原理的には大域解を求めることができ、その応用範囲は飛躍的に拡がる。
【０００４】
すなわち、図２に示すような居所的最少点を複数個有するような誤差関数を対象とし、任意点
2006059237.tif 000003
から出発した力学的システムがシュレディンガ方程式に従って時間発展することにより、図２に示す
2006059237.tif 000004
のような局所的最小点にトラップされることなく量子学的トンネル効果によって原理的にはいずれ大域的最小点
2006059237.tif 000005
に到達する。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかし、この従来技術においては量子トンネル効果をシステムに生起させる必要があるためシュレディンガ方程式を解かねばならず、どの程度の時間、該方程式を積分すれば真の大域的最小点に到達できるか、が自明ではない。このため、例えばリアルタイムでの大域最適点を必要とされる場合等ではおのずと最適化のために費やされる時間は制限されているので、その時間だけシュレディンガ方程式を積分して到達した解が大域解になっているという保証がない。要するに従来技術は遺伝アルゴリズム等を凌駕する良い性質を持っているが、その長所を生かすためにアルゴリズム設計者がチューニングできるようなパラメータを有しないという第一の課題がある。
【０００６】
また、上記従来技術においては探索空間Ωに拘束条件が付帯している場合、シュレディンガ方程式を設定できないので、量子トンネル効果を使えないという第二の課題がある。
【課題を解決するための手段】
【０００７】
上記課題を解決するため、本願で特許請求される発明は下記のとおりである。
（１）大域的最適化装置の操作方法であって、ポテンシャル関数
2006059237.tif 000006
の真の最小値
2006059237.tif 000007
の取り得る領域Ω全体をサーベイ探索する前に、
2006059237.tif 000008
を力学システムの座標と見做し、その量子力学的トンネル効果を使うことにより真の最小値
2006059237.tif 000009
を計算する大域最適化操作手順において、該力学システムのプランク定数に相当する部分に置換されるべき設計定数ＨＴと、力学システムの時間発展を規定する時間微分方程式における該力学システムの散免の変化率を規定する摩擦係数γとを具備することを特徴とする大域的最適化装置の操作方法。
（２）（１）の設計定数ＨＴが、該力学システムに対する量子的揺らぎ力の作用がその2乗ＨＴ２に比例して計算算出されることを特徴とする大域的最適化装置の操作方法。
（３）（１）の設計定数ＨＴが、ポテンシャル障壁を量子学トンネル効果で透過する割合が、ｅｘｐ（−定数／ＨＴ）に比例するような量として規定されていることを特徴とする大域的最適化装置の操作方法。
（４）（１）の摩擦係数γが、該力学システムのエネルギ散免率の比例定数として規定されていることを特徴とする大域的最適化装置の操作方法。
【０００８】
（５）大域的最適化装置の操作方法であって、ポテンシャル関数
2006059237.tif 000010
の真の最小値
2006059237.tif 000011
の取り得る領域Ω全体をサーベイ探索する前に、
2006059237.tif 000012
を力学システムの座標と見做し、その量子力学的トンネル効果を使うことにより真の最小値
2006059237.tif 000013
を計算する大域最適化操作手順において、探索空間Ωに拘束条件が付帯している際はこの拘束条件を、該力学システムのプランク定数に相当する部分に置換されるべき設計定数ＨＴと、力学システムの時間発展を規定する時間微分方程式における該力学システムの散免の変化率を規定する摩擦係数γとともに入力し、この入力式に応じて該装置で交換関係の代数とそれに基づくハミルトニアン演算子が構成され波動方程式が作成され、これに基づき大域最適化演算が行われることを特徴とする大域的最適化装置の操作方法。
【発明の効果】
【０００９】
本願の請求項１〜４の発明によれば、操作者が誤差関数の特性を予め評価した上でＨＲ，γを設定する自由度を持ち、また、このパラメータ設定によって可能な限り速やかに最適値に到達することを可能にする。
【００１０】
また、本願の請求項５の発明によれば、探索空間Ωの上に付帯条件（例えばＸ１２＋Ｘ２２＋・・・・ＸＮ２＝γ２等の等式で表現される付帯条件）が付せられている場合であっても、探索空間の上の点
2006059237.tif 000014
が、あたか自由に動き得るものとして取り扱うことができ、したがって計算量が低減される。
また、探索空間上の付帯条件を満足しつつ揺らぎを利用するので、真の最適解への到達能力が向上する。
【発明を実施するための最良の形態】
【００１１】
図１は、上記第１〜４の発明における実施例図である。図１において、
2006059237.tif 000015
は複数個の局所的最小値を有するポテンシャル関数である。
2006059237.tif 000016
の真の最小値
2006059237.tif 000017
の取り得る領域Ω全体をサーベイ探索する前に、最適化アルゴリズムにおける力学的システムのプランク定数に相当する部分に置換されるべき設計定数ＨＴと、力学システムの時間発展を規定する時間微分方程式における力学システムの散免の変化率を規定する摩擦係数γと、初期状態
2006059237.tif 000018
を入力する。
【００１２】
なお、ＨＴは該力学システムに対する量子的揺らぎ力の作用がその2乗ＨＴ２に比例して計算算出される。またＨＴはポテンシャル障壁を量子学トンネル効果で透過する割合がｅｘｐ（−定数／ＨＴ）に比例するような量として規定されているものとする。
【００１３】
また、摩擦係数γは該力学システムのエネルギ散免率の比例定数として規定されているものである。最適化演算では一定の時間Ｃの間に摩擦つき量子トンネル効果によって点
2006059237.tif 000019
に到達する。これが、すなわち
2006059237.tif 000020
の最小値を与えるか否かを判定し、もし最小値であると判定されれば本アルゴリズムは解
2006059237.tif 000021
を出力する。もしそうでないと判断されれば、得られた
2006059237.tif 000022
を 
2006059237.tif 000023
とし、ＨＴ、γをいずれも増量して再度摩擦つき量子トンネルによる新の最小値サーベイを続行する。ここで真の最小値か否かの判断は、
2006059237.tif 000024
よりできるだけ遠距離数点をその個数ｌも含め乱数を用いて選出し、これを
2006059237.tif 000025
であればＯＫと判断し、
2006059237.tif 000026
ならＮＯと判断するものである。
【００１４】
次に、上記課題２を解決するための第５の発明の実施例を図２および図３により説明する。図２のような拘束条件が探索空間Ωの上に付せられている。このような場合の演算手順は図3に示す通りである。まず操作者は付帯条件を等式拘束条件として入力する。この入力式に応じて装置内では変換関係の代数とそれに基づくハミルトニアン演算子が構成され波動方程式が作られる。この波動方程式にはＨＴ２に比例する揺らぎが含まれているので、対象とする評価関数の最小点を図４のようになトンネル効果によって探索することができる。
【図面の簡単な説明】
【００１５】
【図１】本願請求項１〜４の発明における演算装置の構成図。
【図２】本願請求項５の発明において、探索空間Ωに付せられた付帯条件の説明図。
【図３】本願請求項５の発明における演算手順の説明図。
【図４】大域最適解の探索手順の構成図。
【図１】
2006059237.tif 000039
【図２】
2006059237.tif 000040
【図３】
2006059237.tif 000041
【図４】
2006059237.tif 000042