ネットワーク内での伝達過程に関するシミュレーションを行なおうとしている振動或いは信号の発生源であり、前記ネットワークの一ノードであるスタートノードが、前記ネットワーク内のノード間の繋がり具合を示すネットワークの結合状態情報を、前記ネットワークから取得するネットワーク情報取得手段と、前記結合状態情報に基づいて構築された計算モデルに従って、前記スタートノードの演算装置が前記シミュレーションを実行することにより、前記ネットワーク内の一のノードに接続されたn本nは自然数のリンクを、そのノードに先に前記振動或いは信号を伝達したリンクから順に第1着リンク、第2着リンク、、第n着リンクと呼ぶとき、少なくとも第1着リンクを前記ネットワーク内の各ノードについて求めるネットワーク状態分析手段と、前記スタートノード以外の前記ネットワーク内の一ノードであるターゲットノードから前記スタートノードに向かい、途中の各ノードの第1着リンクを辿ってなる経路を、前記スタートノードと前記ターゲットノードとの間の所望の経路として前記スタートノードが選択するネットワーク経路選択手段とを備えることを特徴とする経路選択装置。前記計算モデルは、前記結合状態情報に基づいて定められた方程式に従って振動し、ノード間のリンクを介して互いに振動を伝達する振動子として、前記ネットワークの各ノードを見立てた計算モデルであり、前記スタートノードを振動の発生源として運動する振動子として計算モデル化し、前記振動子それぞれの物理的位相をシミュレートすることにより、前記スタートノードの振動子から他の振動子への振動の伝達をシミュレートすることを特徴とする請求項1に記載の経路選択装置。前記振動子を非線形振動子として計算モデル化し、前記スタートノードの振動子の周波数の初期値は、当該他のノードの振動子の初期値よりも高い周波数を有し、前記スタートノードの振動子から当該他の振動子への振動の伝達を、周波数の引き込み現象による前記スタートノードの振動子との周波数の一致に基づいて検出することを特徴とする請求項2に記載の経路選択装置。ノードの振幅の変化に基づいてそのノードへの位相の伝達を検出することを特徴とする請求項2に記載の経路選択装置。リンクに沿った一の方向への振動の伝達を検出するため、ノードに対応しない仮想的な振動子の運動を前記シミュレーションに含めることを特徴とする請求項2に記載の経路選択装置。ネットワーク内のすべてのリンクをそのリンクに関する物理量に応じて定められる数の区間に分割し、前記区間のひとつと単位時間とを相関させて、ネットワーク内の信号の伝搬状態をシミュレートする計算モデルであることを特徴とすることを特徴とする請求項1に記載の経路選択装置。前記計算モデルは、一のリンクに関する単位時間当たりの物理量である単位物理量を示す値の整数倍として、各リンクの物理量を近似し、前記単位時間が経過する毎に前記物理量が前記単位物理量刻みに離散的に変化するとして、前記信号の移動をモデル化し、対応するリンクに信号が入った時刻からカウントを開始し、単位時間経過毎にカウントを1ずつ進め、対応するリンクの物理量に対応する値を上限としてカウントするカウンタを、ネットワークの各リンクに対応させて設けることを特徴とする請求項6に記載の経路選択装置。前記リンクの物理量は当該リンクの距離であり、前記単位物理量はスタートノードを発生源とする信号が単位時間にリンク内を移動する距離であることを特徴とする請求項7に記載の経路選択装置。前記リンクの距離として、そのリンクの伝送速度の逆数に比例した係数をかけた値を用いることを特徴とする請求項8に記載の経路選択装置。ネットワークの各リンクに対応させて、リンクに沿った一の方向に対応するカウンタと、当該他の方向に対応するカウンタの動作を前記シミュレーションに含めることを特徴とする請求項6乃至9のいずれかに記載の経路選択装置。ネットワーク内での伝達過程に関するシミュレーションを行なおうとしている振動或いは信号の発生源であり、前記ネットワークの一ノードであるスタートノードが、前記ネットワーク内のノード間の繋がり具合を示すネットワークの結合状態情報を、前記ネットワークから取得するネットワーク情報取得手段、前記結合状態情報に基づいて構築された計算モデルに従って、前記スタートノードの演算装置が前記シミュレーションを実行することにより、前記ネットワーク内の一のノードに接続されたn本nは自然数のリンクを、そのノードに先に前記振動或いは信号を伝達したリンクから順に第1着リンク、第2着リンク、、第n着リンクと呼ぶとき、少なくとも第1着リンクを前記ネットワーク内の各ノードについて求めるネットワーク状態分析手段、及び、前記スタートノード以外の前記ネットワーク内の一ノードであるターゲットノードから前記スタートノードに向かい、途中の各ノードの第1着リンクを辿ってなる経路を、前記スタートノードと前記ターゲットノードとの間の所望の経路として前記スタートノードが選択するネットワーク経路選択手段としてコンピュータを機能させるためのプログラム。前記計算モデルは、前記結合状態情報に基づいて定められた方程式に従って振動し、ノード間のリンクを介して互いに振動を伝達する振動子として、前記ネットワークの各ノードを見立てた計算モデルであり、前記スタートノードを振動の発生源として運動する振動子として計算モデル化し、前記振動子それぞれの物理的位相をシミュレートすることにより、前記スタートノードの振動子から他の振動子への振動の伝達をシミュレートすることを特徴とする請求項11に記載のプログラム。前記振動子を非線形振動子として計算モデル化し、前記スタートノードの振動子の周波数の初期値は、当該他のノードの振動子の初期値よりも高い周波数を有し、前記スタートノードの振動子から当該他の振動子への振動の伝達を、周波数の引き込み現象による前記スタートノードの振動子との周波数の一致に基づいて検出することを特徴とする請求項12に記載のプログラム。ノードの振幅の変化に基づいてそのノードへの位相の伝達を検出することを特徴とする請求項12に記載のプログラム。リンクに沿った一の方向への振動の伝達を検出するため、ノードに対応しない仮想的な振動子の運動を前記シミュレーションに含めることを特徴とする請求項12に記載のプログラム。ネットワーク内のすべてのリンクをそのリンクに関する物理量に応じて定められる数の区間に分割し、前記区間のひとつと単位時間とを相関させて、ネットワーク内の信号の伝搬状態をシミュレートする計算モデルであることを特徴とする請求項11に記載のプログラム。前記計算モデルは、一のリンクに関する単位時間当たりの物理量である単位物理量を示す値の整数倍として、各リンクの物理量を近似し、前記単位時間が経過する毎に前記物理量が前記単位物理量刻みに離散的に変化するとして、前記信号の移動をモデル化し、対応するリンクに信号が入った時刻からカウントを開始し、単位時間経過毎にカウントを1ずつ進め、対応するリンクの物理量に対応する値を上限としてカウントするカウンタを、ネットワークの各リンクに対応させて設けることを特徴とする請求項16に記載のプログラム。前記リンクの物理量は当該リンクの距離であり、前記単位物理量はスタートノードを発生源とする信号が単位時間にリンク内を移動する距離であることを特徴とする請求項17に記載のプログラム。前記リンクの距離として、そのリンクの伝送速度の逆数に比例した係数をかけた値を用いることを特徴とする請求項18に記載のプログラム。ネットワークの各リンクに対応させて、リンクに沿った一の方向に対応するカウンタと、当該他の方向に対応するカウンタの動作を前記シミュレーションに含めることを特徴とする請求項16乃至19のいずれかに記載のプログラム。請求項1乃至請求項10のいずれかに記載の経路選択装置、及び、請求項11乃至請求項20のいずれかに記載のプログラムを格納した記録媒体を備え、前記プログラムに従って動作するコンピュータのいずれかをノードとし、複数のノードを備えるネットワークシステム。前記ノードはそれぞれ、当該他のノードに情報を伝達しようとする度毎に、当該ノードを前記スタートノードとして前記スタートノードと前記ターゲットノードとの間の所望の経路を求めることを特徴とする請求項21に記載のネットワークシステム。ネットワーク内での伝達過程に関するシミュレーションを行なおうとしている振動或いは信号の発生源であり、前記ネットワークの一ノードであるスタートノードが、前記ネットワーク内のノード間の繋がり具合を示すネットワークの結合状態情報を、前記ネットワークから取得する段階と、前記結合状態情報に基づいて構築された計算モデルに従って、前記スタートノードの演算装置が前記シミュレーションを実行することにより、前記ネットワーク内の一のノードに接続されたn本nは自然数のリンクを、そのノードに先に前記振動或いは信号を伝達したリンクから順に第1着リンク、第2着リンク、、第n着リンクと呼ぶとき、少なくとも第1着リンクを前記ネットワーク内の各ノードについて求める段階と、前記スタートノード以外の前記ネットワーク内の一ノードであるターゲットノードから前記スタートノードに向かい、途中の各ノードの第1着リンクを辿ってなる経路を、前記スタートノードと前記ターゲットノードとの間の所望の経路として前記スタートノードが選択する段階とを含むことを特徴とする経路選択方法。前記ネットワークのノードのいずれかが、当該他のノードに情報を伝達しようとする度毎に、その伝達元となるノードを前記スタートノードとして請求項23に記載の経路選択方法を実行することを特徴とする経路選択方法。前記計算モデルは、前記結合状態情報に基づいて定められた方程式に従って振動し、ノード間のリンクを介して互いに振動を伝達する振動子として、前記ネットワークの各ノードを見立てた計算モデルであり、前記スタートノードを振動の発生源として運動する振動子として計算モデル化し、前記振動子それぞれの物理的位相をシミュレートすることにより、前記スタートノードの振動子から他の振動子への振動の伝達をシミュレートすることを特徴とする請求項23及び請求項24のいずれかに記載の経路選択方法。前記振動子を非線形振動子として計算モデル化し、前記スタートノードの振動子の周波数の初期値は、当該他のノードの振動子の初期値よりも高い周波数を有し、前記スタートノードの振動子から当該他の振動子への振動の伝達を、周波数の引き込み現象による前記スタートノードの振動子との周波数の一致に基づいて検出することを特徴とする請求項25に記載の経路選択方法。ノードの振幅の変化に基づいてそのノードへの位相の伝達を検出することを特徴とする請求項25に記載の経路選択方法。リンクに沿った一の方向への振動の伝達を検出するため、ノードに対応しない仮想的な振動子の運動を前記シミュレーションに含めることを特徴とする請求項25に記載の経路選択方法。ネットワーク内のすべてのリンクをそのリンクに関する物理量に応じて定められる数の区間に分割し、前記区間のひとつと単位時間とを相関させて、ネットワーク内の信号の伝搬状態をシミュレートする計算モデルであることを特徴とすることを特徴とする請求項23及び請求項24のいずれかに記載の経路選択方法。前記計算モデルは、一のリンクに関する単位時間当たりの物理量である単位物理量を示す値の整数倍として、各リンクの物理量を近似し、前記単位時間が経過する毎に前記物理量が前記単位物理量刻みに離散的に変化するとして、前記信号の移動をモデル化し、対応するリンクに信号が入った時刻からカウントを開始し、単位時間経過毎にカウントを1ずつ進め、対応するリンクの物理量に対応する値を上限としてカウントするカウンタを、ネットワークの各リンクに対応させて設けることを特徴とする請求項29に記載の経路選択方法。前記リンクの物理量は当該リンクの距離であり、前記単位物理量はスタートノードを発生源とする信号が単位時間にリンク内を移動する距離であることを特徴とする請求項30に記載の経路選択方法。前記リンクの距離として、そのリンクの伝送速度の逆数に比例した係数をかけた値を用いることを特徴とする請求項31に記載の経路選択方法。ネットワークの各リンクに対応させて、リンクに沿った一の方向に対応するカウンタと、当該他の方向に対応するカウンタの動作を前記シミュレーションに含めることを特徴とする請求項29乃至32のいずれかに記載の経路選択方法。以上見てきたように、従来の最短経路検索アルゴリズムは計算時間がかかる、あるいは計算リソースを膨大に消費するという問題を抱えている。このため、ルータのような計算資源が限られた機器では特に、リアルタイムに最短経路を検索し、経路を選択することができなかった。これは、ダイナミックにネットワークの状態が変わる情報ネットワークではとりわけ重要な問題で、情報伝達の遅延により、ユーザーが満足のいくネットワークレスポンスを供給できないという問題を生じさせている。近年、情報ネットワークは増大の一途である。これは最短経路検索がさらに難しくなることを意味しており、このままでは年々ユーザーの利便性が損なわれることが予見される。本発明は、そのような状況を鑑みて考案されたものであり、従来の最短経路探索の計算時間、計算リソースを大幅に減じる技術を提供するものである。上述の課題を解決するため、本発明は、その一態様として、ネットワーク内での伝達過程に関するシミュレーションを行なおうとしている振動或いは信号の発生源であり、ネットワークの一ノードであるスタートノードが、ネットワーク内のノード間の繋がり具合を示すネットワークの結合状態情報を、ネットワークから取得するネットワーク情報取得手段と、結合状態情報に基づいて構築された計算モデルに従って、スタートノードの演算装置がシミュレーションを実行することにより、ネットワーク内の一のノードに接続されたn本nは自然数のリンクを、そのノードに先に振動或いは信号を伝達したリンクから順に第1着リンク、第2着リンク、、第n着リンクと呼ぶとき、少なくとも第1着リンクをネットワーク内の各ノードについて求めるネットワーク状態分析手段と、スタートノード以外のネットワーク内の一ノードであるターゲットノードからスタートノードに向かい、途中の各ノードの第1着リンクを辿ってなる経路を、スタートノードとターゲットノードとの間の所望の経路としてスタートノードが選択するネットワーク経路選択手段とを備えることを特徴とする経路選択装置を提供する。また、本発明は、他の一態様として、ネットワーク内での伝達過程に関するシミュレーションを行なおうとしている振動或いは信号の発生源であり、ネットワークの一ノードであるスタートノードが、ネットワーク内のノード間の繋がり具合を示すネットワークの結合状態情報を、ネットワークから取得するネットワーク情報取得手段、結合状態情報に基づいて構築された計算モデルに従って、スタートノードの演算装置がシミュレーションを実行することにより、ネットワーク内の一のノードに接続されたn本nは自然数のリンクを、そのノードに先に振動或いは信号を伝達したリンクから順に第1着リンク、第2着リンク、、第n着リンクと呼ぶとき、少なくとも第1着リンクをネットワーク内の各ノードについて求めるネットワーク状態分析手段、及び、スタートノード以外のネットワーク内の一ノードであるターゲットノードからスタートノードに向かい、途中の各ノードの第1着リンクを辿ってなる経路を、スタートノードとターゲットノードとの間の所望の経路としてスタートノードが選択するネットワーク経路選択手段としてコンピュータを機能させるためのプログラムを提供する。また、本発明は、他の一態様として、ネットワーク内での伝達過程に関するシミュレーションを行なおうとしている振動或いは信号の発生源であり、ネットワークの一ノードであるスタートノードが、ネットワーク内のノード間の繋がり具合を示すネットワークの結合状態情報を、ネットワークから取得する段階と、結合状態情報に基づいて構築された計算モデルに従って、スタートノードの演算装置がシミュレーションを実行することにより、ネットワーク内の一のノードに接続されたn本nは自然数のリンクを、そのノードに先に振動或いは信号を伝達したリンクから順に第1着リンク、第2着リンク、、第n着リンクと呼ぶとき、少なくとも第1着リンクをネットワーク内の各ノードについて求める段階と、スタートノード以外のネットワーク内の一ノードであるターゲットノードからスタートノードに向かい、途中の各ノードの第1着リンクを辿ってなる経路を、スタートノードとターゲットノードとの間の所望の経路としてスタートノードが選択する段階とを含むことを特徴とする経路選択方法を提供する。本発明によれば、必要なリソースが少なく、計算速度が速い経路選択方法及び装置を提供することができる。また、これにより、大自由度ノードやリンクの数が膨大のネットワークにおける経路選択を行なうことができる。また、経路選択方法を高い頻度で実行できるようになるので、ダイナミックに状況が変わるようなネットワークであっても適切な経路選択を実行することができるようになる。本発明の経路選択方法では、ネットワーク内の各経路の長さを比較して最短の経路を求めるものではなく、各経路の長さを直接的に計算する必要はない。本発明の経路選択方法は、ネットワーク内のスタートノードを振動子或いは信号発生源と見立て、スタートノードから他のノードへの振動/信号の伝達過程を計算機上でシミュレーションし、そのノードに振動/信号を最初に伝達したリンクを求める。本発明のネットワークシステムではネットワークの各ノードが経路選択を行なう。また、経路選択に先立って、各ノードは、ネットワークの全体像に関する情報をネットワークからある程度正確に取得してあるという前提で機能する。この前提は上述したOSPFと同じ前提である。ネットワークの全体像情報は、具体的には、ネットワークのトポロジー構造情報と、ネットワークの結合状態情報を含む。以下、それぞれ単に全体像情報、トポロジー情報、結合状態情報とも呼ぶものとする。ネットワークの結合状態情報とは、ネットワークの結合状態、或いは、それぞれのノードの繋がり具合を示す情報であり、例えば、あるノード間のリンクの距離、あるノード間のリンクを流れる情報流の上限容量、あるノード間のリンクを流れる情報流の速度である。本実施の形態では、結合状態情報としてノード間の距離を用いるものを中心に説明するが、ネットワークの結合状態に関する物理量全般を用いることができることは当業者には容易に理解できるだろう。また、ここでの説明は、情報ネットワークを想定して説明するが、交通ネットワーク等の他のネットワークでも同じことである。ここである程度と書いたのは、ネットワークが大きくなればなるほど、遠いノード間で最新の情報を共有することが難しくなるため、ある時刻でネットワークの全体像を正確に知っていることは難しいからである。また、情報ネットワークのように、時々刻々とリンクの状況繋がり具合が変化する場合、やはり正確なネットワークの全体像を知っていることは難しいからである。本発明によれば、ネットワークの全体像情報が必ずしも完全ではなくとも有効な経路探索を行うことができる。その理由については口述する。第1の実施の形態本発明の第1の実施の形態である経路選択方法についてネットワークシステム1に基づいて説明する。これは本願発明の最も基本的な形態である。ネットワークシステム1はN個Nは自然数のノード1、ノード2、、ノードNを順に直列に接続したネットワークである。ノード1、ノード2、、ノードNはそれぞれ演算処理装置を備える機器であり、具体的にはインターネットアクセスプロバイダのアクセスサーバやルータ等の機器である。これらネットワークのノードそれぞれが以下に説明するような経路選択法を実行する。先に記載したように、各ノードはネットワークの全体像をある程度正確に知っており、その上で、図1に示すような計算モデルを構築し、経路計算を行う。本実施の形態では、ノード間の距離を結合状態情報として用いるので、最短経路の経路計算の例となる。結合状態情報を情報流の流れやすさの逆数などにすれば、最速情報伝達経路の経路計算となる。図1は、説明を簡略化するために、複数ノードが1直線に並んでいるネットワークの場合に本願発明を適用した図となっている。各ノードは隣のノードとリンク1本のみで繋がっている。各ノードに対応する位置に非線形振動子が設置されている。本実施の形態では、非線形振動子としてファンデルポール振動子を想定し、各ノードを繋ぐリンクの距離を振動の拡散係数Dで置き換えている。振動子は、ファンデルポール振動子に限らず、いろいろな振動子で代用できる。結合状態情報は拡散係数Dと所定の関係を有し、本実施の形態ではリンクの距離長さを拡散係数Dに何らかの関数で変換する必要がある。本実施の形態では、熱力学からの類推の元、振動子の振動の拡散を拡散方程式で表す。振動子の振動の拡散を表す方程式は、振動の拡散を記述可能であれば他の方程式でも良い。結果として、図1の非線形振動子のネットワークは以下の方程式で記述できる。ここで、Dの添え字のijはi番目のノードとj番目のノード間であることを示している。Dijがついた項は、振動子iと振動子jの間の関係を表す。ε、αは振動の特徴を決定するパラメータである。上記方程式において、xは振動子の振幅であり、yは物理的にはばねの影響を表す項で、ばねを意味し、rは振動子の位置である。tは時刻である。こうした数1のような、振動の拡散を記述可能な方程式に従って、ネットワークのノード1、ノード2、...、ノードNに対応する振動子1、振動子2、...、振動子Nが運動すると考えて、全振動子の運動をシミュレーションする。出発点のノードに対応する振動子の初期値を変更すると振動をはじめる。出発点以外のノードに対応する振動子の運動もそれぞれ上記方程式に支配されているので、出発点ノードの振動子での振動は、Dijの項で表される振動の伝達状態により、他の振動子に次々と伝達されていく。その際、ノードに最初に振動を伝達するリンクをノード毎に判定し、これらの判定結果に基づいて、スタートノードから各ノードに至る最短経路を求める。上記説明したような計算モデルを設定し、その後、スタートノード計算を行う自分自身に相当するノードの非線形振動子にある外力または初期値を設定すると、その非線形振動子が振動を始める。ε、αを適当に設定しておけば、その振動は外力、初期値によらず、リミットサイクルという一定の振動に収斂する。スタートノードで生じた振動は、拡散係数Dのリンクを通して、周辺のノードに伝達され、やがて、スタートノード周辺の非線形振動子も振動を始める。周辺の非線形振動子が振動を始める時刻は、スタートノードの振動の位相が到達した時刻に等しいため、スタートノードから近いノードほど振動を開始する時刻が早くなる。図2に直列に繋いだ5つのノードNode1、Node2、、Node5のNode1にて発生した振動の位相が伝達する様子を示す。時刻0においてNode2~5は振動していない。スタートノードであるNode1の振動が、Node2、Node3と次々に伝達している様子がわかる。図2では、スタートノード以外は最初振動していない例を示したが、非線形振動子の場合、高い周波数の振動に低い周波数の振動が収斂するという引き込みという現象が存在する。このため、スタートノードを最も高い周波数で振動させ、他のノードをスタートノードより低い周波数で振動させておくと、スタートノードの位相が到達すると、他のノードの周波数が次々とスタートノードの高い周波数に収斂する。この周波数の伝達を位相伝達の検出に使用することもできる。線形振動子ではこの便利な引き込みという現象が起きないので、位相伝達の検出のバリエーションが減るが、単純な振動の伝達図2のようなケースであるならば、線形振動子でも問題ない。説明の簡略化のため、図1では1次元の直線ネットワークを用いて説明したので、上記現象でネットワークの経路選択ができるのかイメージしにくいが、上記現象を利用すれば、ネットワークの経路選択が可能となる。ひとつのノードに複数のリンクが接続されている場合でも、位相の伝達は最も早く位相が伝わる経路、つまり最短経路を通って伝達される。したがって、あるノードに注目した場合、そのノードに最も早く位相を伝達したリンクが、そのノードへの最短経路を形成すると考えればよい。ここで、あるノードにm個mは自然数のリンクが接続されているとき、そのノードに位相を1番目、2番目、、m番目に伝達したリンクを順に第1着リンク、第2着リンク、第m着リンクと呼ぶものとする。ノードに接続されたあるリンクを介して位相が伝達され、第1着リンクが決定した時点で、そのノードの第2着以下のリンクは最短経路を形成する資格を失う。このため、最短経路検索だけを行うのであれば、第1着リンクのみを記録し、他のリンクは無視してもよい。スタートノードからそのノードへの最短経路を求めるには、そのノードの第1着リンクから順に途中の各ノードの第1着リンクを経てスタートノードに戻るような経路を求める。同様にしてネットワーク内のスタートノード以外の各ノードからスタートノードに戻る経路を求める。これにより、スタートノードを頂点として、他のノードへの最短経路がツリー状に示されるのである。これを用いれば、所望のターゲットノードへ情報を伝達するには、自分はどのリンクに情報を伝達すれば良いかがわかる。図3左側のネットワーク概略図は、各ノードにおける第1着リンクを太線で示し、第2着以下のリンクを細線で示した例である。このようなネットワーク概略図から第1着リンクのみを抽出すると図3右側のツリー構造を得ることができる。例えば、図中左側のネットワーク概略図の出発点、スタートノードであるノードN1からノードN7に至る最短経路は、図中右側のツリー構造によるとノードN1、N2、N5、N6、N7であることが分かる。最短経路だけでなく、2番目、3番目にあるノードに到達するリンクを調べたければ、あるノードに最も早く位相を伝達したリンクだけを記録するのではなく、各リンクに対して、位相を伝達した時刻を記憶し、リンク間で伝達時刻を比較することによりリンクの着順を判定すればよい。ノードに位相が伝達されたことを検出する方法は幾つか考えられるので、図4に典型的なものを2つまとめて示しておく。一つ目はそのノードに対応する振動子のシミュレーションの振幅変化を検出するものである。ある振動子から位相が伝達されると、振動していなかった振動子が振動を開始する。これを利用して、そのノードに対応する振動子である程度の振幅が検出されたら、その振動子、即ちノードに位相が伝達されたと判断する方法である。図4aに示すように、位相伝達検出用スライスレベルを定め、そのノードに対応する振動子の運動を、例えば数1等の方程式に従ってシミュレーションする。シミュレーションした振動子の振幅が位相伝達検出用スライスレベルを超えた時点でそのノードに位相が伝達されたと判定する。2つ目は振動の周期を検出するものである。ある周波数で振動する振動子から振動が伝達された振動子は、伝達元の振動子と同じ周期で振動するようになる。これを利用して、数1等の方程式に従ってシミュレーションした振動子の周期をモニターし、スタートノードの振動子の周期に相当する周期用スライスレベルと同じになったとき、スタートノードの振動子からそのノードの振動子に位相が伝達されたと判定する。上記に位相の検出方法を示したが、ひとつのノードに幾つかのリンクが結合している場合、どのリンクから位相が到達したかを判断するのはやや面倒になるときがある。というのは、各リンクからの位相伝達の差が微妙で、上記の検出方法で検出される前に、複数のリンクから位相が実際は伝達されてしまうような場合、どちらのノードからの位相が早く伝達したかを確認するには、リンク内の振動の様子、一つ前のノードの振動の様子などを分析するアルゴリズムを使用しなければならないからである。そこで、本願発明者は、図5に示すような改良型の経路選択方法も考案した。説明の簡略化のため、図1と同様、1次元の直線ネットワークを用いて説明している。図1との相違は、各リンクにどちらのノードから位相が伝達されたかを検出する位相検出用の振動子を仮想的に設置した点である。既に説明したように、各ノードに対応する振動子は、数1のような振動の拡散を記述可能な方程式に従って運動するものとしてシミュレーションするが、位相検出用振動子の運動についても同じ方程式に従うものとする。拡散係数Dについても同様である。便宜上、この2つの振動子を往路、復路と呼ぶことにする。本実施の形態では、この位相検出用の振動子も、ファンデルポール振動子で形成している。この位相検出用の振動子は、各リンクに往路用復路用の2つのものが設置され、振動子の反対側に位置するノードの振動のみを伝達させ、自身に振動が伝わった時刻に、振動子の反対側に位置するノードから振動が伝わったと判断する。このようにすることで、明確に往路側と復路側の振動の伝達が切り分けられ、どちら側からより早く位相が伝わったかを判断することができる。図5を参照すると、ネットワーク内で互いに隣接する2つのノードがあり、これらのノードに対応する振動子として振動子1、振動子2がある。振動子1と振動子2の間はリンクにて接続されている。リンクの振動子1側の端部付近に、ネットワーク内のノードに対応しない仮想的な振動子として、位相検出用振動子復路を仮想的に設置する。同様に、リンクの振動子2側の端部付近に位相検出用振動子往路を仮想的に設置する。位相検出用振動子復路には、振動子2からの振動のみを伝達させ、振動子1からの振動は伝達させない。また、位相検出用振動子往路には、振動子1からの振動のみを伝達させ、振動子2からの振動は伝達させない。ここで、5つのノードNode1、Node2、...、Node5が図6のように接続されたネットワークを考える。図中の丸数字は対応するリンクを示すリンク番号であり、例えば、Node1とNode2の間のリンクはリンク番号1を有するものとする。図6にはLink1、Link2、...、Link5が存在する。リンク番号の横に記載したDの値は、そのリンクにおけるノードの繋がり具合を示す値であり、上記の数1における拡散係数Dに相当する。Link1、2、...、5にはそれぞれ、往路及び復路の位相検出用振動子不図示が設置されているものとする。このとき、Node1~5に対応する振動子、Link1~5における位相検出用振動子往路、及び、Link1~5における位相検出用振動子復路の振動をシミュレートした結果を図7~9に示す。特に、図8と図9とを比較することにより、往路と復路の振動が微妙にずれていることがわかり、それにより、最短経路はどの経路かが明確にわかるようになっている。本実施の形態は、計算の並列化が容易なため、GPGPUなどの並列計算が得意なハードウェアとの相性が良い。また、後述する第2の実施の形態とは異なり、ノード間の距離を連続的な値として扱うことができる。第2の実施の形態次に、本発明の第2の実施の形態である経路選択方法について説明する。第1の実施の形態では、ノードを振動子と見立て、リンクを振動の伝達媒体と見立てて、物理的な位相をノード間で伝達するときの各ノードの運動を計算機上でシミュレーションした。これに対して、本実施の形態では、ネットワーク内のすべてのリンクをそのリンクに関する物理量例えばそのリンクの距離、伝送速度などに応じて定められる個数の区間に分割し、分割した区間のひとつと単位時間とを相関させることにより、ネットワーク内の信号の伝搬状態、即ち、スタートノードが発信した信号が他のノードに伝搬していく過程をシミュレートするような計算モデルを構築する。ネットワーク或いはリンクの結合状態を示す物理量の単位量を単位時間に関連づけて定めた上で、そのリンクの物理量が単位時間毎に離散的に変化させることにより、リンクにおける信号の伝搬状態をシミュレーションするようなモデルを構築した。例えば、ネットワークの結合状態を示す物理量としてリンクの距離を用いる場合、リンクの一端から他端までの間のどこまで信号が伝達しているのかを、単位時間における信号の移動距離、即ち、単位距離の整数倍で表したものが、リンクにおける信号の伝達状態となる。ここでは、リンクの物理量としてリンクの距離を用いる場合を考える。スタートノードから発信された信号がリンク内を移動し、他のノードそれぞれに到達するまでの間のリンク内における信号の移動を仮想的にシミュレーションする。信号はリンク内を一定速度で伝達されるものとし、信号があるリンクからノードに入力されると、信号は入力時のリンクを含む全てのリンクから分岐して出力されるものとする。また、ここではネットワーク内のどのリンクも等しい伝送速度を有するとする。第1の実施形態と同様に、各ノードは、経路選択に先立って、ネットワークの全体像情報、即ち、ネットワークのトポロジー情報と結合状態情報をネットワークから取得する。ここでは結合状態情報としてノード間の距離つまりリンクの距離を取得するものとして説明する。単位時間に信号がリンク内を伝達される距離を単位距離と呼ぶものとする。単位距離は最大でもネットワーク内で最短のリンクの距離であり、必要な精度に応じて定める。より小さい単位距離を用いると計算精度は上がるが計算負荷は大きくなる。ネットワーク内のリンク間の距離の比を全て整数で表すことができる場合は、これら整数の最大公約数を単位距離として用いることとすれば、高い計算精度と低い計算負荷を両立することができる。当然ながら、信号が2単位距離のリンクを通過するには2単位時間を要し、3単位距離のリンクを通過するには3単位時間を要することになる。このように、ネットワーク内の全リンクを単位距離の整数倍で近似すると、信号は、ネットワークのどのリンクであっても単位時間の整数倍で通過し、あるリンクに信号が入った時刻から経過した単位時間は、そのままそのリンク内での信号の現在位置を示すことになる。そこで、リンクの一端にあるノードに信号が入ると起動してその時刻からカウントを開始し、単位時間経過毎にカウントを1ずつ進め、対応するリンクの距離を単位距離の整数倍で近似した整数値と等しい値を上限としてカウントするカウンタを、ネットワークの各リンクに対応させて設ける。このようなカウンタにおいて、カウンタの値が上限に達することは、対応するリンクの一端のノードから入った信号が、そのリンクの他端のノードに到着することを意味する。このようなカウンタをネットワーク内のリンク毎に設け、単位時間毎に、ネットワーク内で現在起動中の全カウンタの値を1ずつ進めることにより、スタートノードから発信した信号がネットワークの他ノードに伝達していく過程をシミュレートする。第1の実施の形態において妥当な結果を得るためには、1周期を数千程度のオーダーの単位時間に分割し、単位時間が経過する度毎に、その時刻について数1のような方程式を解く演算処理が発生する。あるリンクの一端から他端に振動が伝達される過程をシミュレートするとき、仮に100周期分の振動をするのであれば、その間に数十万程度の単位時間が経過することになり、従って上述の演算処理を数十万回程度実行することになる。また、この演算処理は差分計算等の科学演算を含む演算処理であり、浮動小数点演算である。このため、シミュレーションに必要な計算負荷の全体を比較すると、第1の実施の形態はダイクストラ法ベースの手法に比べてそれほど有利ではないが、第1の実施の形態の手法は並列化が容易なため、並列計算を行なうことによりダイクストラ法ベースの手法よりも短時間で処理を実行することができる。尚、一般にダイクストラ法ベースの手法は並列化が難しい。このような第1の実施の形態に対して第2の実施の形態は単位時間を長くし、必要な計算負荷を下げることができる。例えば、本実施の形態では、必要な計算精度や、ネットワーク内のリンクの距離のばらつき具合にもよるが、ネットワーク内で最短のリンクの距離を単位距離とすることも可能であり、この場合、最短リンクの信号通過のシミュレーションで経過する単位時間は1である。また、第2の実施の形態において、単位時間が経過する毎に実行すべき処理は、基本的にはリンク毎にカウンタが起動しているか否か判定し、起動中のカウンタの値を1ずつ加算する処理である。この処理は差分計算等の科学計算の演算処理と比較して負荷が小さい。つまり、第2の実施の形態では、実行すべき演算処理の回数が少なく、かつ、各回で実行すべき演算処理の負荷が小さい。第1の実施形態では、ノードに振動を伝達したリンクを特定するため、図5のような位相検出用振動子を設けた。信号の伝達方向を特定するため、本実施の形態では、リンク毎に往路用カウンタと復路用カウンタの2つのカウンタを設けることとしてもよい。一般にひとつのリンクの両端にはそれぞれノードが接続されている。往路用カウンタはそのうちの一方のノードを担当ノードとし、復路用カウンタは他方のノードを担当ノードとする。往路用カウンタの担当ノードが信号の伝達を検出すると、往路用カウンタはこれに応じてカウントを開始する。復路用カウンタも同様にしてカウントを開始する。次に、本実施の形態にて行なうシミュレーション処理について説明する。この処理ではスタートノードから発信された信号がネットワーク内のリンクを移動する過程をシミュレーションする。スタートノードから信号を発信するのと同時に、スタートノードを担当ノードとするカウンタを起動する。スタートノードに複数のリンクが接続されているのであれば、それらリンクに対応する複数のカウンタ全てを起動する。以後、スタートノードから信号を発信した時刻をt=0とし、1、2、...、n単位時間経過後の時刻をそれぞれt=1、2、...nと表すとする。1単位時間経過毎に、起動中の全カウンタの値を1だけ進める。t=1においてスタートノードを担当ノードとするカウンタの値はどれも1になっており、スタートノードに接続されたリンクの、スタートノードから1単位距離離れた位置を信号が進んでいることを示している。スタートノードに接続されたリンクのひとつが2単位距離であるとすると、t=2でそのリンクのカウンタはカウントを完了する。このようにカウンタがカウントを完了したということは、そのリンクに接続された2つのノードのうち、一方のノードで発信或いは入力された信号が、他方のノードに到着したことを意味する。このように信号がノードに到着するのと同時に、そのリンク、そのノード及び時刻を互いに関連づけて記録し、そのノードを担当ノードとするカウンタを起動する。スタートノードでのカウンタの起動と同様、そのノードを担当ノードとするカウンタが複数ある場合は全てのカウンタを起動する。このようにして起動したカウンタを、以後1単位時間経過毎に、既に起動中のカウンタと共に1ずつカウントアップしていく。起動中のカウンタいずれかでのカウントが完了して信号がノードに到着すると、そのノードを担当ノードとする全カウンタを起動する、といった動作を繰り返していく。カウント完了時に信号が到着したノードを担当ノードとするカウンタが存在しない場合は新しいカウンタを起動しない。このようなシミュレーションの結果、ネットワーク内の各ノードについて、ノード、そのノードに信号を伝達させたリンク、及び、その到着時刻が関連づけて記録する。最短経路のみを求める場合には各ノードの第1着リンクのみを記録すればよい点では第1の実施の形態と同様であり、この場合、既知のネットワークのトポロジーと各リンクの第1着リンクに基づいて、図3と同様のツリー構造を生成することができる。本実施の形態では、計算ステップ時刻とリンクの距離とを、単位時間と単位距離として相関させるため、各リンクの距離を連続に表現することは苦手で、ある程度離散的な値として近似する必要がでてくるが、物理的な位相を計算しないため、現在の直列処理を得意とするCPUであっても、計算時間を大幅に減らすことができる。これは本実施の形態の大きなメリットである。第1の実施の形態と同様に計算の並列化が容易なため、GPGPUを用いて計算することは更に好ましい。上述の説明ではネットワーク内のどのリンクも等しい伝送速度を有するとしたが、異なる伝送速度のリンクが混在する場合であっても本実施の形態を用いることができる。例えば、距離を2とし、伝送速度を1とするようなリンクL1と、距離を4とし、伝送速度を2とするようなリンクL2とを考えると、リンクL2の距離はリンクL1の倍だが、リンクL2の伝送速度もリンクL1の倍のため、リンクL1とリンクL2を信号が通過する際の所要時間はどちらも等しいと考えられる。これを考慮して、各リンクに対してそのリンクの伝送速度の逆数に比例した係数をかけたものを実質距離と呼ぶものとすると、単位距離を定める前に各リンクの実質距離を求めて、実質距離に基づいて単位距離を定めることとすればよい。このようにして求めた単位距離を実質単位距離と呼ぶものとする。リンクの実質距離を実質単位距離の整数倍で近似した値をそのリンクのカウンタの上限値とする。このようにすれば、伝送速度が異なるリンクが混在するネットワークであっても本実施の形態を適用することができる。上述の説明においては、カウンタは単位時間毎に1ずつカウント値を加算するが、減算することとしてもよい。即ち、ネットワーク内の各リンクにそのリンクの距離に比例するカウンタ値を予めセットしておき、そのリンクの一端のノードへの信号到達に応じてカウンタを起動する。起動されたカウンタは1単位時間経過する毎にカウント値を1ずつ減算する。カウント値がゼロになるとそのリンクの他端のノードに信号が到達したと考えるのである。また、今回、仮想的な位相伝達情報を伝えるのにカウンタを用いたが、位相が物理的に伝わるのと相関がとれる位相が伝わったという情報を使用するなら、本願発明と同様のことが実施できる。次に本発明の実施例1である経路選択装置10について図10を参照して説明する。経路選択装置10はネットワーク情報取得部11、ネットワーク状態分析部12、伝搬方向識別部13、ネットワーク経路選択部14を備える。経路選択装置10及び各部11~14は、記憶装置に格納したプログラムに従って動作するコンピュータとして実現され、より具体的には例えばインターネットアクセスプロバイダのサーバや、ルータの一部を構成する。ネットワーク情報取得部11は、ネットワークインタフェース装置を備え、ネットワークの全体像情報、即ち、ネットワークのトポロジー情報と、ネットワークの結合状態情報をネットワークから取得する。本願発明を最も特徴付けるのは、ネットワーク状態分析部12である。ネットワーク状態分析部12はその内部に更に伝播方向識別部13を持つ場合があり、本実施例では、伝播方向識別部13が搭載されている場合を示している。ネットワーク状態分析部12の役目は、ネットワーク情報取得部11が獲得したネットワークの結合状態に関する情報から、ネットワーク経路選択部14が経路選択を実施できるまでの情報に加工することにある。ネットワーク状態分析部12で行う処理こそが、第1の実施の形態において述べた経路選択方法によって経路選択処理を実行する装置である。つまり、振動の拡散を記述可能な方程式に従って運動し、ノードの繋がり具合を示す物理量に応じて振動を伝達するリンクにて互いに接続された振動子をノードとするネットワークにおいて、スタートノードに対応する振動子での振動が他の振動子に伝達される過程をシミュレートすることにより、各ノードの第1着リンクを求め、予め用意した例えば図3右側のネットワーク概略図のようなネットワークのトポロジーと、各ノードの第1着リンクから例えば図3右側のようなツリー構造を求める処理を実行する。スタートノードは当該経路選択装置10を備えるサーバ、ルータ等である。本実施例では、非線形振動子の位相伝達を用いた。非線形振動子はファンデルポール振動子を用いた。伝播方向識別部13は、第図5及び6を参照して説明したように、各リンクに設置される往路用、復路用の仮想リンクの状況をモニターする部位である。ネットワーク状態分析部12は、伝播方向識別部13による仮想リンクのモニター情報を用いて、図3に示したようなツリー構造を作成していく。本実施例では、最短経路探索を行った例を示すが、本願発明は発明の効果のところで述べたように、各リンクからの位相伝達時刻を記憶しておけば、様々な経路選択に応用できる。ネットワーク経路選択部14は、ネットワーク状態分析部12が求めたツリー構造に基づいて、当該経路選択装置10を備えるサーバ、ルータ等、スタートノードとなるノード装置から、外部のプログラムや入力装置から指定された、ターゲットノードとなるノード装置に至る最短経路を選択し、出力する。図11に本実施例で最短経路検索を行ったネットワークの例を示す。ノード数をNとして、立方体型のネットワークを構築、使用した。また、本実施例ではリンクの距離拡散係数Dはすべて同じとした。本実施例では、複数の非線形振動子の物理的位相を計算する例であるので、通常の単一のCPUでは計算が遅くなってしまう。そこで、本実施例では、GPGPUを用いて、並列計算を行った。GPGPUGeneralPurposeGraphicProcessingUnitとは、汎用目的で使用できるグラフィック用のプロセッサで、超並列計算を安価なボードで実現できる近年注目されているプロセッサである。なぜGPGPUが注目されているかと言うと、CPU単体の計算性能はすでに頭打ちであり、単純な処理を超並列で行うしかトータルの計算性能を上げる術がなくなってきたからである。本願発明のように、ひとつひとつの処理各非線形振動子の計算は独立して行えるような並列計算に適したアルゴリズムは、GPGPUなどの力を借りて今後主流になると思われる。従来用いられてきたダイクストラ法などのアルゴリズムは並列計算に適していない。上記の条件で、本実施例を実施したところ、任意のスタートノードから他の全ノードまでの最短経路検索を行ったところ、妥当なツリー構造の出力を得ることができた。図12にダイクストラ法と本実施例の計算速度の結果をノード数の依存性として示す。計算速度は、1000ノードのときの本実施例の計算時間で規格化して算出している。本実施例は従来のダイクストラ法よりもノードの数が大きくなるにつれ、速い計算結果が得られている。このことから本願発明が非常に有効な経路選択方法であることがわかった。本発明の実施例2の経路選択装置について説明する。本実施例は第2の実施の形態に対応する。本実施例の経路選択装置も実施例1の経路選択装置と同様の構成を有するので、図10を参照して説明する。本実施例では、ネットワーク状態分析部12は、非線形振動子の位相伝達を計算するのではなく、第2の実施の形態に相当するシミュレーションを行なう。即ち、スタートノードから信号を発信した時刻を起点として、単位時間毎のリンク上での信号の位置を計算することにより、各ノードの第1着リンク等を求める。第1の実施の形態ではシミュレーションに必要な計算を並列化しやすいため、実施例1ではGPGPUを用いてシミュレーションを実行したが、第2の実施の形態では、計算の並列化が容易なだけではなく、シミュレーションに要する演算負荷が低いために並列化されていないCPUであっても計算速度が向上する。この効果をより明らかにするため、一般的なCPUを用いた結果をグラフにして示す。尚、第2の実施の形態も並列化が容易な点では変わりなく、計算処理を並列化してGPGPUで実装すれば更なる高速化を実現することができる。使用したネットワークは実施例1と同じ図11に示したような立方体型のネットワークである。また、リンクの距離はすべて同じとして、1時刻後に隣接ノードに情報が到達するという条件とした。このような条件で本実施例を実施し、任意のスタートノードから他の全ノードまでの最短経路検索を行ったところ、妥当なツリー構造の出力を得ることができた。図13にダイクストラ法と本実施例の計算速度の結果をノード数の依存性として示す。計算速度は、1000ノードのときの本実施例の計算時間で規格化して算出している。本実施例は従来のダイクストラ法よりも速い計算結果が得られている。本発明の方式は、ネットワークの全体像をある程度正確に知っているという前提で動作することは、第1及び第2の実施の形態の説明において述べた。そのときにネットワークの全体像を正確に知っている必要がないのかという疑問がでることを予期した。このような疑問に応えるため、以下に述べるような簡単なシミュレーションを行った。図14を参照して説明する。このシミュレーションでは、スタートノードからターゲットノードに向かって情報を伝達する際に、本実施例の経路選択装置10を用いてターゲットノードまでの最短経路を選択し、その最短経路に沿って情報を発信する。経路選択装置10が最短経路を求める際に用いたネットワークの全体像情報では、ネットワーク内に障害が発生中のノードは存在していないものとする。このときに求めた最短経路を構成するリンクを図中に太線で示す。スタートノードから情報を発信した後、その情報がターゲットノードに到着するまでの間に、図14に示すように、先に求めた最短経路上のノードにおいて障害が発生したものとする。このときネットワークの各ノードに対して障害ノード情報を通知する必要が生じる。障害ノード情報は最終的には全ノードに対して伝達され、全てのノードにおけるネットワークの全体像情報が更新されるが、障害発生後すぐの間は、障害が発生したノード周辺のノードでのみ、ネットワークの全体像情報を更新することができ、遠方のノードでは更新することができない。一方、スタートノードが発信した情報を受け取ったノードは、そのノードがその時点までに取得しているネットワークの全体像情報に基づいて求めた、ターゲットノードへの最短経路に従って情報を転送していく。このような前提でシミュレーションを行ったところ、スタートノードが情報を発信した時点では知られていなかった障害を、途中のノードが検知し、以後は更新されたネットワークの全体像情報に基づいて求めた最短経路に従って情報を伝達することが確認出来た。このように、本発明によれば、ネットワークが自律的に障害を回避して、ターゲットノードに情報を伝達することができる。このような効果を得ることができるのは、本発明の経路選択方法によれば非常に高速に経路選択を実行可能であり、このため、従来は頻繁に実行できなかった経路選択計算を各ノードが情報を送信する度毎に実行することができるからである。また、このように本発明の経路選択方法に基づく経路選択処理を、ネットワークの各ノードにおいて、そのノードが情報を送信する度毎に行なうこととすれば、自律的に障害を回避するネットワークを提供することが可能となる。このように、本発明の経路選択方法は、単に動作が高速であることに留まらない非常に有効な経路選択方法である。以上、本発明を実施の形態及び実施例に即して説明したが、本発明はこれに限定されるものではなく、本発明の技術的な意味を変えることなく様々な変形が可能である。例えば、上述の実施形態及び実施例は、ネットワークとして情報ネットワークを挙げて説明したが、本発明はネットワークの種類を問うものではなく、例えば交通ネットワークにおける経路選択に対しても本発明を適用することができる。また、上述の説明では最短経路検索を行なう場合を中心にして説明したが、2番目に短い経路、3番目に短い経路というように様々な経路選択の要請を満足する処理も可能である。また、上述の説明では、ネットワークの結合状態情報としてリンクの距離を用いたため、選択される経路は最短経路となったが、ネットワークの結合状態情報はリンクの距離に限定されるものではない。例えば、リンクを流れる情報流の速度、容量などを結合状態情報として、所望の経路選択を実施することが可能である。また、第1の実施の形態及び実施例1では、ノードを非線形振動子と見立てたが、例えばソリトンを使用して位相を伝達することとしてもよい。この出願は、2010年11月19日に出願された日本出願特願第2010-259345号を基礎とする優先権を主張し、その開示のすべてをここに取り込むものである。