特許 6985669 通信品質劣化推定装置、通信品質劣化推定方法、及びプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6985669

(24)【登録日】2021年11月30日

(45)【発行日】2021年12月22日

(54)【発明の名称】通信品質劣化推定装置、通信品質劣化推定方法、及びプログラム

(51)【国際特許分類】

   H04L 12/70 20130101AFI20211213BHJP

【ＦＩ】

   H04L12/70 100A

【請求項の数】5

【全頁数】12

(21)【出願番号】特願2018-117413(P2018-117413)

(22)【出願日】2018年6月20日

(65)【公開番号】特開2019-220864(P2019-220864A)

(43)【公開日】2019年12月26日

【審査請求日】2020年10月26日

(73)【特許権者】

【識別番号】000004226

【氏名又は名称】日本電信電話株式会社

(73)【特許権者】

【識別番号】504176911

【氏名又は名称】国立大学法人大阪大学

(74)【代理人】

【識別番号】100107766

【弁理士】

【氏名又は名称】伊東 忠重

(74)【代理人】

【識別番号】100070150

【弁理士】

【氏名又は名称】伊東 忠彦

(74)【代理人】

【識別番号】100124844

【弁理士】

【氏名又は名称】石原 隆治

(72)【発明者】

【氏名】田行 里衣

(72)【発明者】

【氏名】池上 大介

(72)【発明者】

【氏名】松田 崇弘

(72)【発明者】

【氏名】滝根 哲哉

【審査官】 中川 幸洋

(56)【参考文献】

【文献】 特開２０１７−１９２０４１（ＪＰ，Ａ）

【文献】 米国特許第０８１４５７４５（ＵＳ，Ｂ１）

【文献】 松田 崇弘 Takahiro MATSUDA，信号処理技術によるネットワークトモグラフィの解法 Signal Processing Approaches for Network Tomography，電子情報通信学会技術研究報告 Ｖｏｌ．１１５ Ｎｏ．３２１ IEICE Technical Report，日本，一般社団法人電子情報通信学会 The Institute of Electronics,Information and Communication Engineers，2015年11月17日，第115巻

【文献】 松田 崇弘 Takahiro MATSUDA，ｌ<SB>１</SB>／ｌ<SB>２</SB>最適化手法を用いたネットワークトモグラフィ Network Tomography with l<SB>1</SB>/l<SB>2</SB> Optimization，電子情報通信学会技術研究報告 Ｖｏｌ．１１０ Ｎｏ．４４９ IEICE Technical Report，日本，社団法人電子情報通信学会 The Institute of Electronics,Information and Communication Engineers，2011年02月24日，第110巻

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｈ０４Ｌ １２／７０

(57)【特許請求の範囲】

【請求項1】

通信網を構成するノードの品質劣化を推定する通信品質劣化推定装置であって、
２点のノード間の品質観測値と、各観測における経路情報と、ノード同士の劣化の類似性を表すグラフの構造を表す行列とを入力する入力手段と、
品質観測値から構成されるベクトルが、各観測の経路情報から構成される行列と各ノードの劣化量の値から構成されるベクトルとの積に近いという関係と、前記グラフの構造を表す行列に基づいた類似ノード間の劣化量の差が小さいという関係とから導出される推定式を解くことにより、各ノードの劣化量の推定値を算出する演算手段と、
前記各ノードの劣化量の推定値、又は、当該推定値に基づいて得られる情報を出力する出力手段と、を備え、
前記演算手段は、各ノードの劣化量の値から構成される縦ベクトルとグラフラプラシアン行列と各ノードの劣化量の値から構成される横ベクトルの積を最小化するような解を選択する推定式から前記推定値を算出する
ことを特徴とする通信品質劣化推定装置。

【請求項2】

前記グラフの構造を表す行列として、接続行列と次数行列の差で表されるグラフラプラシアン行列を使用する
ことを特徴とする請求項１に記載の通信品質劣化推定装置。

【請求項3】

前記演算手段は、グラフ構造を表す行列を用いた制約以外に、劣化するノードが少数であるという仮定に基づいた制約を含む推定式を解くことにより前記推定値を算出する
ことを特徴とする請求項１又は２に記載の通信品質劣化推定装置。

【請求項4】

通信網を構成するノードの品質劣化を推定する通信品質劣化推定装置が実行する通信品質劣化推定方法であって、
２点のノード間の品質観測値と、各観測における経路情報と、ノード同士の劣化の類似性を表すグラフの構造を表す行列とを入力する入力ステップと、
品質観測値から構成されるベクトルが、各観測の経路情報から構成される行列と各ノードの劣化量の値から構成されるベクトルとの積に近いという関係と、前記グラフの構造を表す行列に基づいた類似ノード間の劣化量の差が小さいという関係とから導出される推定式を解くことにより、各ノードの劣化量の推定値を算出する演算ステップと、
前記各ノードの劣化量の推定値、又は、当該推定値に基づいて得られる情報を出力する出力ステップと、を備え、
前記演算ステップにおいて、各ノードの劣化量の値から構成される縦ベクトルとグラフラプラシアン行列と各ノードの劣化量の値から構成される横ベクトルの積を最小化するような解を選択する推定式から前記推定値を算出する
ことを特徴とする通信品質劣化推定方法。

【請求項5】

コンピュータを、請求項１ないし３のうちいずれか１項に記載の通信品質劣化推定装置における各手段として機能させるためのプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、通信ネットワークにおける通信品質状態を管理する技術に関するものである。

【背景技術】

【0002】

モバイル通信ではユーザ環境の多様化、利用方法の多岐化により、様々な箇所で品質劣化が発生する状況である。モバイル通信サービス品質向上のためには通信品質を測定して現在の通信品質状態を把握し、品質劣化検出時にはその要因となる箇所を特定することが重要である。エンドエンド遅延時間やパケット損失率、スループットといった品質を把握する直接的な方法としては、エンドエンド間で試験パケットを送受信することにより品質を能動的に測定する方法がある。

【0003】

非特許文献１では、発着ノードペアでの試験パケットでの品質測定（パケット損失率や遅延）から、ネットワーク内の各リンクでの品質を推定する方法について紹介しており、具体的には、複数の発着ノードペア間でマルチキャスト通信により試験パケットを送出して各リンクのパケット損失率を推定する方法を紹介している。ユニキャスト通信により試験パケット測定を行うと、各リンクのパケット損失率が一意に定まらない場合があるが、パケットペアを送出することでその問題を解決する方法を紹介している。

【0004】

非特許文献２〜４では、複数本のリンクのうち、少数のリンクが品質劣化していると仮定し、リンクの品質推定・品質劣化リンクの特定に関する検討を行っている。非特許文献２では特に品質劣化リンク本数が高々１本のとき、品質劣化リンクを特定可能なネットワークのルーチングに対する条件を導出しており、非特許文献３では必要な測定パス数を導出している。

【0005】

また、近年の端末の高機能化・多様化やネットワーク利用形態の多様化に伴い、特にモバイル端末を用いた通信品質は、そのときのアプリケーションや無線通信環境、端末の位置、アクセス先サーバ、端末のＯＳの種類やバージョン、利用するアクセスネットワーク（ＷｉＦｉ（登録商標）か３Ｇかなど）など要素に大きく依存すると考えられる。そのような状況も考慮した方法が、特許文献１において提案されている。

【0006】

一方で、従来は各ノードに対して同様のモデルや前提を仮定していたのに対して、モバイル通信で影響を与える構成要素(上記で挙げたアプリケーションや無線通信環境、端末の位置、アクセス先サーバ等)は同質なモデルとして扱うことは難しい。それに対して、非特許文献５は、異なる性質を仮定した場合の各要素の推定を提案している。

【先行技術文献】

【非特許文献】

【0007】

【非特許文献1】M. Coates et al., ``Internet tomography,'' IEEE Signal Processing Magazine, May 2002.

【非特許文献2】M. H. Firooz et al., ``Network tomography via compressed sensing,'' Globecom 2010.

【非特許文献3】W. Xu et al., ``Compressive sensing over graphs: how many measurements are needed?'', Forty-Eighth Annual Allerton Conference, 2010.

【非特許文献4】T. Matsuda et al., ``Link quality classifier with compressed sensing based on l1-l2 optimization,'' IEEE Communication Letters, Vol. 15, No. 10, Oct. 2011.

【非特許文献5】H. Kinsho et al., `` Heterogeneous Delay Tomography Based on Graph Fourier Transform in Mobile Networks,'' CQR 2017.

【特許文献】

【0008】

【特許文献1】特許第６０６１８３８号

【発明の概要】

【発明が解決しようとする課題】

【0009】

上述したような異なる要素（ノード等）からなるネットワークにおいて、ユーザ端末から計測を行ったクラウドソーシング等の手段で収集したエンドエンドの測定データを入力とする場合に、推定したい要素の計測値が得られるとは限らないという課題がある。

【0010】

本発明は上記の点に鑑みてなされたものであり、異なる要素からなるネットワークにおいてスケーラブルに監視可能な情報から、通信品質劣化検出時に劣化箇所を特定し、当該劣化箇所の品質推定値を推定することを可能とする技術を提供することを目的とする。

【課題を解決するための手段】

【0011】

開示の技術によれば、 通信網を構成するノードの品質劣化を推定する通信品質劣化推定装置であって、
２点のノード間の品質観測値と、各観測における経路情報と、ノード同士の劣化の類似性を表すグラフの構造を表す行列とを入力する入力手段と、
品質観測値から構成されるベクトルが、各観測の経路情報から構成される行列と各ノードの劣化量の値から構成されるベクトルとの積に近いという関係と、前記グラフの構造を表す行列に基づいた類似ノード間の劣化量の差が小さいという関係とから導出される推定式を解くことにより、各ノードの劣化量の推定値を算出する演算手段と、
前記各ノードの劣化量の推定値、又は、当該推定値に基づいて得られる情報を出力する出力手段と、を備え、
前記演算手段は、各ノードの劣化量の値から構成される縦ベクトルとグラフラプラシアン行列と各ノードの劣化量の値から構成される横ベクトルの積を最小化するような解を選択する推定式から前記推定値を算出する
ことを特徴とする通信品質劣化推定装置が提供される。

【発明の効果】

【0012】

開示の技術によれば、異なる要素からなるネットワークにおいてスケーラブルに監視可能な情報から、通信品質劣化検出時に劣化箇所を特定し、当該劣化箇所の品質推定値を推定することを可能とする技術が提供される。

【図面の簡単な説明】

【0013】

【図1】実施例１から３の関係を表す図である。

【図2】通信品質劣化推定装置の構成例を示す図である。

【図3】通信品質劣化推定装置のハードウェア構成例を示す図である。

【図4】実施例４におけるグラフ構成を示す図である。

【図5】実施例４におけるグラフ構成を示す図である。

【図6】実施例４における基地局の劣化を示す図である。

【図7】実施例４におけるサーバの劣化を示す図である。

【図8】実施例４における与えた真値と推定結果を示す図である。

【発明を実施するための形態】

【0014】

以下、本発明の実施の形態として、実施例１〜実施例４を説明する。実施例１〜実施例４のうちの実施例１〜３に関し、図１にそれらの関係を示す。

【0015】

本発明の実施の形態では、通信品質劣化推定装置１００が開示される。図１に示すように、実施例１では、入力の準備について説明する。図１に示すように、入力の準備は、エンドエンド観測値、経路情報（観測フローが通過した要素の情報）、要素間の関係性を表す情報（グラフ構造）に基づいて実行される。実施例２では、１種類の構成要素からなる品質推定について説明する。実施例３では、２種類の構成要素からなる品質推定について説明する。

【0016】

（装置構成例）
実施例１〜実施例４を説明するにあたり、実施例１〜実施例４における処理を実行する通信品質劣化推定装置１００について説明する。図２は、通信品質劣化推定装置１００の構成図である。図２に示すように、通信品質劣化推定装置１００は、入力情報生成部１１０、情報管理部１２０、及び品質ベクトル推定部１３０を有する。

【0017】

入力情報生成部１１０は、エンドエンド観測値、経路情報（観測フローが通過した要素の情報）、要素間の関係性を表す情報（グラフ構造）を入力とし、これらに基づいて、ノード同士の劣化の類似性を表すグラフの構造を表す行列等を、情報管理部１２０への入力情報として生成する。情報管理部１２０は、入力情報に基づいて推定式を決定し、当該推定式を品質ベクトル推定部１３０に入力する。品質ベクトル推定部１３０は、入力された推定式から品質ベクトルを推定し、当該品質ベクトル（各ノードの劣化量の推定値）を出力する。品質ベクトル推定部１３０は、当該推定値に基づいて得られる情報（例：グラフィカルな画像、統計値等）を出力してもよい。具体的な処理内容については各実施例で説明する。

【0018】

通信品質劣化推定装置１００は、コンピュータに内蔵されるＣＰＵやメモリ等のハードウェア資源を用いて、当該品質劣化推定装置１００で実施される処理に対応するプログラムを実行することによって実現することが可能である。上記プログラムは、コンピュータが読み取り可能な記録媒体（可搬メモリ等）に記録して、保存したり、配布したりすることが可能である。また、上記プログラムをインターネットや電子メール等、ネットワークを通して提供することも可能である。

【0019】

図３は、本実施の形態における上記コンピュータのハードウェア構成例を示す図である。図３のコンピュータは、それぞれバスＢで相互に接続されているドライブ装置１５０、補助記憶装置１５２、メモリ装置１５３、ＣＰＵ１５４、インタフェース装置１５５、表示装置１５６、及び入力装置１５７等を有する。

【0020】

当該コンピュータでの処理を実現するプログラムは、例えば、ＣＤ−ＲＯＭ又はメモリカード等の記録媒体１５１によって提供される。プログラムを記憶した記録媒体１５１がドライブ装置１５０にセットされると、プログラムが記録媒体１５１からドライブ装置１５０を介して補助記憶装置１５２にインストールされる。但し、プログラムのインストールは必ずしも記録媒体１５１より行う必要はなく、ネットワークを介して他のコンピュータよりダウンロードするようにしてもよい。補助記憶装置１５２は、インストールされたプログラムを格納すると共に、必要なファイルやデータ等を格納する。

【0021】

メモリ装置１５３は、プログラムの起動指示があった場合に、補助記憶装置１５２からプログラムを読み出して格納する。ＣＰＵ１５４は、メモリ装置１５３に格納されたプログラムに従って、当該品質劣化推定装置１００に係る機能を実現する。インタフェース装置１５５は、ネットワークに接続するためのインタフェースとして用いられる。表示装置１５６はプログラムによるＧＵＩ（Ｇｒａｐｈｉｃａｌ Ｕｓｅｒ Ｉｎｔｅｒｆａｃｅ）等を表示する。入力装置１５７はキーボード及びマウス、ボタン、又はタッチパネル等で構成され、様々な操作指示を入力させるために用いられる。

【0022】

以下、実施例１〜実施例４を説明する。

【0023】

（実施例１）
実施例１では、入力の準備に関して説明する。本発明の実施の形態では、品質劣化推定の対象となるネットワークを構成する要素（ノード）は、品質劣化についての周囲の機器との依存性や品質劣化の時間変化について依存性等を有するという特徴を持つ要素群ａ＿１，ａ＿２，…，ａ＿Ｎ_Ａ（例：アクセスネットワーク機器等）と、劣化の発生が全体のうちの少数であるという特徴を持つ要素群ｂ＿１，ｂ＿２，…，ｂ＿Ｎ_Ｂ（例：サーバ・コアネットワークの機器群等）からなる。

【0024】

まず、上記のネットワークから、通信フローｍにおけるエンドエンド通信品質情報に関する値ｚ＿ｍを収集する。収集する値ｚ＿ｍについては、例えば、ネットワークにおけるあるサーバからデータをダウンロード又はアップロードしているとし、該サーバ−該端末間での片側の通信品質でもよいし、あるサーバとの往復の通信品質でもよい。通信品質の指標としては、例えば、遅延、ラウンドトリップタイム、パケット損失率等を用いることができる。

【0025】

エンドエンド観測値である値ｚ＿ｍを収集することに関しては、通信品質劣化推定装置１００の外部で行うこととしてよい。エンドエンド観測値を取得する際に、観測フローが通過した要素の情報である経路情報が得られる。

【0026】

推定を行いたい要素群ａ＿１，ａ＿２，…，ａ＿Ｎ_Ａの劣化（遅延等）をｘ＿１，ｘ＿２，…，ｘ＿Ｎ_Ａとし、要素群ｂ＿１，ｂ＿２，…，ｂ＿Ｎ_Ｂの劣化（遅延等）をｙ＿１，ｙ＿２，…，ｙ＿Ｎ_Ｂとする。これらの要素からなるベクトルをそれぞれベクトルｘ，ベクトルｙとする。

【0027】

ここで、要素間の関係性を表す情報として、ａ＿１，ａ＿２，…，ａ＿Ｎ_Ａの関係性を表すグラフ構造を作成する。このグラフ構造の作成については、通信品質劣化推定装置１００の外部で行うこととしてよい。

【0028】

グラフ構造の作成においては、各要素を頂点として、例えば距離の近い頂点同士を辺で接続する。ここで、時間変化の関係性を扱いたい場合には、各タイムスロットにおける同じ機器を別の頂点として定義し、それらを辺で接続する。この辺は、距離に関連する関数等で重み付けを行うこともできる。

【0029】

上述したエンドエンド観測値、経路情報、及びグラフ構造が、通信品質劣化推定装置１００の入力情報生成部１１０に与えられる。

【0030】

入力情報生成部１１０は、グラフ構造の重み付け接続行列をＷ（Ｎ_Ａ×Ｎ_Ａ行列）として生成する。また、入力情報生成部１１０は、このグラフの次数行列Ｃ（Ｎ_Ａ×Ｎ_Ａ行列）を用いて、グラフラプラシアンＬ（Ｎ_Ａ×Ｎ_Ａ行列）をＬ＝Ｗ−Ｃとして生成する。グラフラプラシアンＬは、ノード同士の劣化の類似性を表すグラフの構造を表す行列である。

【0031】

なお、エンドエンド観測値、経路情報、及びグラフラプラシアンＬが通信品質劣化推定装置１００の外部で作成され、通信品質劣化推定装置１００の入力手段（入力情報生成部１１０等）から通信品質劣化推定装置１００に入力されることとしてもよい。

【0032】

（実施例２）
実施例２では、実施例１により得られた情報に基づいて、１種類の構成要素からなる品質推定を実行する処理について説明する。

【0033】

まず、入力情報生成部１１０は、経路情報に基づいて、ａ＿１，ａ＿２，…，ａ＿Ｎ_Ａの要素群からなるネットワークの観測において、各観測値が通過した要素を表す０，１の要素からなる行列であるＲ_Ａを生成する。観測値、Ｒ_Ａ、及びＬは情報管理部１２０に入力される。

【0034】

ａ＿１，ａ＿２，…，ａ＿Ｎ_Ａの要素群からなるネットワークの観測において、観測値ｚ＿ｍからなるベクトルをｚとすると、Ｒ_Ａ及び誤差項εを用いることで、ベクトルｚとベクトルｘとの関係を下記の式（１）で表わすことができる。

【0035】

ｚ＝Ｒ_Ａｘ＋ε 式（１）
上記の式（１）は、各観測の品質観測値から構成されるベクトルが、各観測の経路情報から構成される行列と各ノード（要素）の劣化量の値から構成されるベクトルとの積に近いという関係に基づく式である。

【0036】

情報管理部１２０は、上記の式（１）とＬに基づいて、推定したｘをｘ'とした推定式として、下記の式（２）を決定する。

【0037】

ｘ'＝ａｒｇｍｉｎ＿｛ｘ｝｛１／２｜｜ｚ-Ｒ_Ａｘ｜｜^２＋ηｘＬｘ^Ｔ｝ 式（２）
上記の式（２）において、｜｜＊｜｜はｌ_２ノルムを表し、＊^Ｔは転置を表す。ηは最適化パラメータである。

【0038】

ｘＬｘ^Ｔ（各ノードの劣化量の値から構成される縦ベクトルとグラフラプラシアン行列と各ノードの劣化量の値から構成される横ベクトルの積）を最小化することは、構成したグラフにおける接続された頂点間の差の２乗を最小化することに対応する。つまり、推定式に、ｘＬｘ^Ｔの項を加えることで、グラフ構造を表す行列に基づいた類似ノード間の劣化量の差が小さいという関係を反映できる。

【0039】

式（２）は品質ベクトル推定部１３０に入力され、品質ベクトル推定部１３０が式（２）を解くことによって、ベクトルｘを推定する。品質ベクトル推定部１３０はベクトルｘを出力する。あるいは、品質ベクトル推定部１３０はベクトルｘから得られたグラフィカルな画像、統計値等を出力してもよい。

【0040】

（実施例３）
実施例３では、実施例１により得られた情報に基づいて、２種類の構成要素からなる品質推定を実行する処理について説明する。

【0041】

まず、入力情報生成部１１０は、経路情報に基づいて、ａ＿１，ａ＿２，…，ａ＿Ｎ_Ａとｂ＿１，ｂ＿２，…，ｂ＿Ｎ_Ｂの要素群からなるネットワークの観測において、各観測値が通過した要素を表す０，１の要素からなる行列であるＲ_Ａ及びＲ_Ｂを生成する。観測値、Ｒ_Ａ、Ｒ_Ｂ、及びＬは情報管理部１２０に入力される。

【0042】

ａ＿１，ａ＿２，…，ａ＿Ｎ_Ａとｂ＿１，ｂ＿２，…，ｂ＿Ｎ_Ｂの要素群からなるネットワークの観測において、観測値ｚ＿ｍからなるベクトルをｚとすると、Ｒ_Ａ及びＲ_Ｂ及び誤差項εを用いることで、ベクトルｚとベクトルｘ及びベクトルｙとの関係を下記の式（３）で表わすことができる。

【0043】

ｚ＝Ｒ_Ａｘ＋Ｒ_Ｂｙ＋ε 式（３）
上記の式（３）は、各観測の品質観測値から構成されるベクトルが、各観測の経路情報から構成される行列と各ノード（要素）の劣化量の値から構成されるベクトルとの積に近いという関係に基づく式である。

【0044】

情報管理部１２０は、上記の式（３）とＬに基づいて、推定したｘ、ｙをｘ'、ｙ'とした推定式として、下記の式（４）を決定する。

【0045】

｛ｘ'，ｙ'｝＝ａｒｇｍｉｎ＿｛ｘ，ｙ｝｛１／２｜｜ｚ-Ｒ_Ａｘ−Ｒ_Ｂｙ｜｜^２＋η_１ｘＬｘ^Ｔ＋η_２｜ｙ｜｝ 式（４）
上記の式（４）において、｜｜＊｜｜、｜＊｜はそれぞれｌ_２ノルム、ｌ_１ノルムを表し、＊^Ｔは転置を表す。η_１，η_２は最適化パラメータである。

【0046】

実施例２と同様に、ｘＬｘ^Ｔ（各ノードの劣化量の値から構成される縦ベクトルとグラフラプラシアン行列と各ノードの劣化量の値から構成される横ベクトルの積）を最小化することは、構成したグラフにおける接続された頂点間の差の２乗を最小化することに対応する。つまり、推定式に、ｘＬｘ^Ｔの項を加えることで、グラフ構造を表す行列に基づいた類似ノード間の劣化量の差が小さいという関係を反映できる。また、推定式における｜ｙ｜の項は、グラフ構造を表す行列を用いた制約以外の別の制約の例であり、劣化するノードが少数であるという仮定に基づいた制約を示している。

【0047】

式（４）は品質ベクトル推定部１３０に入力され、品質ベクトル推定部１３０が式（４）を解くことによって、ベクトルｘ、及びベクトルｙを推定する。品質ベクトル推定部１３０はベクトルｘ、ベクトルｙを出力する。あるいは、品質ベクトル推定部１３０はベクトルｘ、ベクトルｙから得られたグラフィカルな画像、統計値等を出力してもよい。

【0048】

（実施例４）
通信品質劣化推定装置１００による具体的な処理の例を実施例４として説明する。実施例４では、１タイムスロットを任意の時間幅で設定して、２０タイムスロットの期間において観測が実施されているとする。ネットワークは２５個の基地局と１５個のサーバから構成されるものとする。基地局は、実施例１で説明したアクセスネットワーク機器等の例である。

【0049】

実施例４では、通信品質劣化推定装置１００が、各タイムスロットにおける各基地局及び各サーバの遅延量をエンドエンドの遅延から推定する。推定する変数は、２０×２５＋２０×１５＝８００個となる。ここで、エンドエンドの遅延が、基地局の遅延とサーバの遅延の和で表されるとする。

【0050】

次に、基地局（２０タイムスロット×２５箇所）の関係性を表すグラフの構成を行う。各基地局を頂点として、基地局間の距離に関して辺を接続したグラフを図４に示す。また、図５で示すように各タイムスロットが隣り合っている物理的位置が同じ基地局についてもグラフで接続する。

【0051】

図６は、基地局に与えられる劣化量（推定の対象とする量）を示す。図６に示すように、ｔ＝１、４、７１０、１３、１６のそれぞれにおいて、図示のように各基地局に劣化量が与えられている。

【0052】

図７は、サーバに与えられる劣化量（推定の対象とする量）を示す。図７に示すように、ｔ＝１、４、７１０、１３、１６のそれぞれにおいて、図示のように各サーバ（横軸）に劣化量（縦軸）が与えられている。

【0053】

図６、７に示すように基地局とサーバの劣化を与えた場合に、通信品質劣化推定装置１００が実施例３の方法で推定を行った結果を図８に示す。ここでは、全ての基地局で観測値を持つわけではなく、図８（中央）で値を持つ箇所でのみ観測値を持つものとする。全基地局で観測値を持たない場合にも、図８（右）に示す通り、良い推定値を得ることができる。

【0054】

（実施の形態のまとめ）
以上、説明したように、本実施の形態によれば、通信網を構成するノードの品質劣化を推定する通信品質劣化推定装置であって、２点のノード間の品質観測値と、各観測における経路情報と、ノード同士の劣化の類似性を表すグラフの構造を表す行列とを入力する入力手段と、品質観測値から構成されるベクトルが、各観測の経路情報から構成される行列と各ノードの劣化量の値から構成されるベクトルとの積に近いという関係と、前記グラフの構造を表す行列に基づいた類似ノード間の劣化量の差が小さいという関係とから導出される推定式を解くことにより、各ノードの劣化量の推定値を算出する演算手段と、前記各ノードの劣化量の推定値、又は、当該推定値に基づいて得られる情報を出力する出力手段とを備えることを特徴とする通信品質劣化推定装置が提供される。

【0055】

入力情報生成部１１０は入力手段の例である。情報管理部１２０及び品質ベクトル推定部１３０は演算手段の例であり、品質ベクトル推定部１３０は出力手段の例である。

【0056】

前記グラフの構造を表す行列として、例えば、接続行列と次数行列の差で表されるグラフラプラシアン行列を使用する。

【0057】

前記演算手段は、各ノードの劣化量の値から構成される縦ベクトルとグラフラプラシアン行列と各ノードの劣化量の値から構成される横ベクトルの積を最小化するような解を選択する推定式から前記推定値を算出することとしてもよい。

【0058】

前記演算手段は、グラフ構造を表す行列を用いた制約以外に、劣化するノードが少数であるという仮定に基づいた制約を含む推定式を解くことにより前記推定値を算出することとしてもよい。

【0059】

以上、本実施の形態について説明したが、本発明はかかる特定の実施形態に限定されるものではなく、特許請求の範囲に記載された本発明の要旨の範囲内において、種々の変形・変更が可能である。

【符号の説明】

【0060】

１００ 通信品質劣化推定装置
１１０ 入力情報生成部
１２０ 情報管理部
１３０ 品質ベクトル推定部
１５０ ドライブ装置
１５１ 記録媒体
１５２ 補助記憶装置
１５３ メモリ装置
１５４ ＣＰＵ
１５５ インターフェース装置
１５６ 表示装置
１５７ 入力装置

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】