(57)【要約】【課題】多数チャンネル毎に時間変化に対する振幅確率分布のグラフ表示を基本として、その表示範囲やグラフ形態等の表示条件を容易に変更して観察可能にする。【解決手段】周波数分析部１００で被測定信号の周波数成分を分析し、レベル検出部２００で各周波数成分（各チャンネル）の振幅をそれぞれ検出し、ＡＰＤ部３００が各周波数成分について単位時間当たりの振幅確率であって所定確率範囲に入る振幅確率を時間経過毎に求める。表示制御部５００は、各周波数成分毎に一方の軸を振幅の大きさ、他方の軸を時間とするグラフ上に、所定確率範囲に入る振幅確率の分布を、ＡＰＤ部３００の出力に応答して表示部６００に表示させる構成とした。【選択図】      図１
2009133683.tif 000002
【特許請求の範囲】
【請求項１】
被測定信号の周波数成分を分析する周波数分析部（１００）と、分析された各周波数成分の振幅をそれぞれ検出するレベル検出部（２００）と、検出された各周波数成分について単位時間当たりの振幅確率を時間経過毎に求めるＡＰＤ部（３００）と、該ＡＰＤ部が出力する各周波数成分の振幅確率を時間経過に対応して記憶する記憶部（４００）と、表示部（６００）と、操作部（７００）と、該操作部からの指示を基に、該各周波数成分毎に一方の軸を振幅の大きさ、他方の軸を時間とするグラフとして、該記憶部に記憶された前記各周波数成分の振幅確率を並列に識別可能にして該表示部に表示させる表示制御部（５００）とを備えたことを特徴とするＡＰＤ測定装置。
【請求項２】
前記表示制御部は、更新の都度、表示される時間軸の右端を確率を更新表示するときの現在時間とし、更新の都度、前回更新時の時間との差だけ時間軸をスクロールしたグラフ表示を行うことを特徴とする請求項１に記載のＡＰＤ測定装置。
【請求項３】
前記表示制御部は、前記各周波数成分の少なくとも１つの周波数成分の前記グラフ上で前記操作部により指定された位置に指定された時間範囲のゾーンマーカを付すとともに、拡大指示があったときは該ゾーンマーカにより該指定された時間範囲の目盛りを拡大した拡大グラフを前記表示部に表示させるゾーン拡大表示制御手段（５１２）を備えたことを特徴とする請求項１又は２に記載のＡＰＤ測定装置。
【請求項４】
前記表示制御部は、前記操作部により、前記各周波数成分の少なくとも１つの周波数成分の前記グラフを指定して縮小指示を受けることにより、該グラフの時間軸の目盛りを縮小したグラフを前記表示部に表示させるフルスパン表示制御手段（５１１）を備えたことを特徴とする請求項１〜３のいずれか一つに記載のＡＰＤ測定装置。
【請求項５】
前記ＡＰＤ部は前記各所定確率範囲の振幅確率を求めて複数の異なる前記所定確率範囲に分類し、前記表示制御部は、前記グラフ又は／及び前記拡大グラフ上で該各所定確率範囲毎に識別して視認可能に前記表示部へ表示させることを特徴とする請求項１〜４のいずれか一つに記載のＡＰＤ測定装置。
【請求項６】
前記表示制御部は、前記ＡＰＤ部が求めた各周波数成分毎であって、かつ各所定確率範囲毎に最新の時間における、又は指定した時間における振幅確率を数値表示させることを特徴とする請求項１〜５のいずれか一つに記載のＡＰＤ測定装置。
【請求項７】
前記ＡＰＤ部は、前記操作部により周波数成分の指定された所定確率範囲内で最大の振幅値、その振幅確率、及びその時間を求めるピークサーチ手段（５４０）を有し、さらに、前記表示制御部は、前記表示部に対して、表示されている前記グラフ又は／及び前記拡大グラフ上であって該最大の振幅値とその時間に相当する位置に視認可能なピークマーカを表示させるピークマーカ表示制御部（５３０）を備えたことを特徴とする請求項１〜６のいずれか一つに記載のＡＰＤ測定装置。
【請求項８】
前記ピークマーカ表示制御部は、前記表示部に対して前記ピークマーカを表示させるとともに、前記最大の振幅値とそのときの確率値を数字表示させることを特徴とする請求項７に記載のＡＰＤ測定装置。
【請求項９】
前記表示制御部は、前記操作手段によりグラフ変更指示を受け、前記グラフ又は／及び前記拡大グラフの指定された時間位置の前記振幅確率もしくは該グラフ変更指示時点の前記振幅確率を基に、一方の軸を振幅、他方の軸を確率とする振幅／確率グラフを生成して前記表示部に表示させる振幅／確率グラフ表示制御手段（５６０）を備えたことを特徴とする請求項１〜８のいずれか一つに記載のＡＰＤ測定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、受けた信号の周波数成分を分析し、各周波数成分の大きさが所定の閾値を越える確率（振幅確率分布、或いは、単に時間率と言われる。以下「ＡＰＤ」と言う。）を測定するＡＰＤ測定装置に関する。特に、ＡＰＤの変化を観察しやすくしたＡＰＤ測定装置である。
【背景技術】
【０００２】
ＡＰＤ測定技術は従来からあり、例えば、ＡＰＤを高い振幅分解能と時間分解能で測定する特許文献１の技術があり、それをスペクトラムアナライザ等へ適用し、種々の表示を試みた例として特許文献２の技術がある。
【０００３】
一方、最近、通信方式により、所定複数周波数成分のＡＰＤを並列に同時に測定する必要性が要求されてきている。例えば、非特許文献１に記載のように、地上デジタル放送の規格は電波産業会（ＡＲＩＢ）によってＡＲＩＢ標準規格として、ＯＦＤＭ方式（Ｏｒｔｈｏｇｏｎａｌ  Ｆｒｅｑｕｅｎｃｙ  Ｄｉｖｉｓｉｏｎ  Ｍｕｌｔｉｐｌｅｘ  Ｏｐｅｒａｔｉｏｎ  Ｍｏｄｅ）による変調方式が定められている。そして、このＯＦＤＭ方式では、伝送データを所定の帯域に数千本の低速データに分けて、デジタル変調を行う。例えば、その所定の帯域が５．６ＭＨｚであり、伝送データを１ＫＨｚの搬送波（キャリア）単位で分けて、計５、６００本に分けてデジタル変調している。
【０００４】
したがって、ＯＦＤＭ方式におけるＡＰＤ測定は、信号成分を１ＫＨｚ毎に周波数分析して、各周波数成分のＡＰＤを並列に測定することが望まれる。また、特許文献１，２に記載のＡＰＤで並列に測定可能である。
【０００５】
ＡＰＤの表示については、特許文献２では確率分布を確率の範囲毎に帯表示させることにより、分布を視認しやすいように識別して表示させている。
【０００６】
【特許文献１】特許第３１５６１５２号公報
【特許文献２】特許第３３７４１５４号公報
【非特許文献１】東芝レビューＶＯＬ．５８  Ｎｏ．１２（２００３年）、ｐ２〜６
【発明の開示】
【発明が解決しようとする課題】
【０００７】
ところで、上記の地上デジタル放送でのＡＰＤ測定のように同時に（並列に）ＡＰＤを測定すべき対象とする周波数の数が多いこと、さらに到来信号（或いは雑音）の統計的傾向を調査或いは監視をするために、測定に要する時間が長時間（数時間から２４時間を超える場合がある。）であること等から、表示にあたっても、同時測定している周波数成分の測定結果をできるだけ多く同時表示させること、そして、測定時間全体について結果が表示できることが望ましい。
【０００８】
しかし、表示面積は有限であるため、例えば、時間軸の全範囲を２４時間とすれば、時間軸のデータ分解能の広くなり、細かい観察ができない。さらには時間に対する確率の変化だけでなく、レベル変化に対する確率変化を見たい等の要求がある。細かく観察したい、或いは座標軸を変更して観察したい等の要求は、測定中に生じることがある。例えば、測定時間２４時間経過後にそれらの要求を満たしていたのでは、効率が悪い。
【０００９】
本発明の目的は、多数チャンネル毎に時間変化に対する振幅確率分布のグラフ表示を基本として、その表示範囲やグラフ形態等の表示条件を容易に変更して観察可能にする技術を提供する。
【課題を解決するための手段】
【００１０】
本発明の目的を達成するため、請求項1の発明は、測定信号の周波数成分を分析する周波数分析部（１００）と、分析された各周波数成分の振幅をそれぞれ検出するレベル検出部（２００）と、検出された各周波数成分について単位時間当たりの振幅確率を時間経過毎に求めるＡＰＤ部（３００）と、該ＡＰＤ部が出力する各周波数成分の振幅確率を時間経過に対応して記憶する記憶部（４００）と、表示部（６００）と、操作部（７００）と、該操作部からの指示を基に、該各周波数成分毎に一方の軸を振幅の大きさ、他方の軸を時間とするグラフとして該記憶部に記憶された前記各周波数成分の振幅確率を並列に識別可能にして該表示部に表示させる表示制御部（５００）とを備えた。
【００１１】
請求項２に記載の発明は、請求項1に記載の発明において、前記表示制御部は、更新の都度、表示される時間軸の右端を確率を更新表示するときの現在時間とし、更新の都度、前回更新時の時間との差だけ時間軸をスクロールしたグラフ表示を行う構成とした。
【００１２】
請求項３に記載の発明は、請求項1又は２に記載の発明において、前記表示制御部は、前記各周波数成分の少なくとも１つの周波数成分の前記ＡＰＤ部は複数の異なる前記所定確率範囲に分けて該各所定確率範囲の振幅確率を求め、前記表示制御部は、前記グラフ又は／及び前記拡大グラフ上で該各所定確率範囲毎に識別して視認可能に前記表示部へ表示させる前記グラフ上で前記操作部により指定された位置に指定された時間範囲のゾーンマーカを付すとともに、拡大指示があったときは該ゾーンマー２により該指定された時間範囲の目盛りを拡大した拡大グラフを前記表示部に表示させるゾーン拡大表示制御手段（５１１）を備えた。
【００１３】
請求項４の発明は、請求項１〜３のいずれかに記載の発明において、前記表示制御部は、前記操作部により、前記各周波数成分の少なくとも１つの周波数成分の前記グラフを指定して縮小指示を受けることにより、該グラフの時間軸の目盛りを縮小したグラフを前記表示部に表示させるフルスパン表示制御手段（５１１）を備えた。
【００１４】
請求項５に記載の発明は、請求項１〜４のいずれかに記載の発明において、前記ＡＰＤ部は該各所定確率範囲の振幅確率を求めて複数の異なる前記所定確率範囲に分類し、前記表示制御部は、前記グラフ又は／及び前記拡大グラフ上で該各所定確率範囲毎に識別して視認可能に前記表示部へ表示させる構成とした。
【００１５】
請求項６に記載の発明は、請求項１〜５のいずれかに記載の発明において、前記表示制御部は、前記ＡＰＤ部が求めた各周波数成分毎であって、かつ各所定確率範囲毎に最新の時間における、又は指定した時間における振幅確率を数値表示させる構成とした。
【００１６】
請求項７に記載の発明は、請求項１〜６の何れかに記載の発明において、前記ＡＰＤ部は、前記操作部により指定された周波数成分の指定された所定確率範囲内で最大の振幅値、その振幅確率、及びその時間を求めるピークサーチ手段（５４０）を有し、さらに、前記表示制御部は、前記表示部に対して、表示されている前記グラフ又は／及び前記拡大グラフ上であって該最大の振幅値とその時間に相当する位置に視認可能なピークマーカを表示させるピークマーカ表示制御部（５３０）を備えた。
【００１７】
請求項８に記載の発明は、請求項７に記載の発明において、前記ピークマーカ表示制御部は、前記表示部に対して前記ピークマーカを表示させるとともに、前記最大の振幅値とそのときの確率値を数字表示させる構成とした。
【００１８】
請求項９に記載の発明は、請求項１〜８の何れかに記載の発明において、前記表示制御部は、前記操作手段によりグラフ変更指示を受け、前記グラフ又は／及び前記拡大グラフの指定された時間位置の前記振幅確率もしくは該グラフ変更指示時点の前記振幅確率を基に、一方の軸を振幅、他方の軸を確率とする振幅／確率グラフを生成して前記表示部に表示させる振幅／確率グラフ表示制御手段（５６０）を備えた
【発明の効果】
【００１９】
請求項１に記載の発明によれば、複数の周波数成分（チャンネル）の時間の変化に対する振幅確率の変化を並列に、更新表示できるので、各チャンネルの並列観察が可能である。各チャンネルの時間軸を横軸とし振幅を縦軸とした座標上に複数の振幅確率範囲に分類して識別可能に表示しているので、観察しやすい。識別するために、例えば、複数の振幅確率範囲を色分け表示すれば、混み在った表示画面の中であっても、遠方より傾向を視認することができる。逆に言えば、表示画面を有効に利用可能である。
請求項２に記載の発明によれば、グラフを更新するときに時間軸の右端を更新時刻に直して（スクロールして）表示データを更新するので、全横軸範囲を有効に使用する。
請求項３に記載の発明によれば、所望の範囲を拡大表示できるので精細な観察が可能となる。
請求項４に記載の発明によれば、横軸を縮小して表示できるので、全体の傾向、測定の進行状況が把握しやすい。
請求項５に記載の発明によれば、確率範囲毎に帯状、確率分布を視認可能に表示できる。
請求項６，７，８に記載の発明によれば、数値表示できることから、グラフで読みににくい所望位置の値を定量的に読みとれる。
請求項７に記載の発明によれば、指定した確率範囲内で最大となるレベル及びその時間の位置を容易に知ることができる。
請求項９に記載の発明によれば、指定された時間位置における振幅対確率のグラフに変換して表示させることができる。
【発明を実施するための最良の形態】
【００２０】
本発明に係る実施形態について、図面を参照しながら説明する。図１は、ＡＰＤ測定装置に係る実施形態の全体の機能構成を示す図である。図２は、縦軸が振幅（レベルと言うことがある。）、横軸が時間の座標上での確率分布の表示例を示す図である。図３は、フルスパン表示機能、並びにゾーンマーカ機能及びゾーン拡大表示機能を説明するための図である。図４は、ピークマーカ機能を説明するための図である。図５は、指定マーカ機能を説明するための図である。図６は、振幅確率値を指定時の振幅／確率の表示例である。図７は、振幅値を指定時の振幅／確率の表示例である。図８は、函数尺でなる座標上での振幅／確率の表示例である。
【００２１】
図１に基づいて、動作順に構成・動作を説明する。
図１において、周波数分析部１００は、測定しようとする所望の周波数帯（広帯域）の入力信号を選択し、さらに狭い所望の周波数帯域（狭帯域）に分けて選択するものであって、帯域選択部１１０、Ａ／Ｄ変換手段１２０、及びフィルタバンク１３０で構成される。帯域選択部１１０としては、例えば、高周波の通信周波数帯域におけるノイズを測定しようとするなら、いわば特許文献２に示されるようは、ヘテロダイン方式を採用することができる。つまり、入力信号（ノイズ）を局部発振器でミキシングして低周波数の信号に変換して、ここで所望の周波数帯域（δＦ）のバンドパスフィルタを用いて必要な周波数帯の周波数成分に制限する。なお、測定対象の信号の周波数帯域が低く、Ａ／Ｄ変換器が利用できる周波数帯であれば、帯域選択部１１０は必ずしも必要ではない。
【００２２】
Ａ／Ｄ変換手段１２０は、帯域選択部１１０で帯域制限された低周波数の周波数成分（信号）を受け、クロック発生器（不図示）からの周期のクロックで振幅（レベル）幅ΔＬの細かさ（この細かさに応じた閾値がある。）で標本化（サンプリング）し、デジタルデータに変換する。後に求められる振幅確率は、この標本化時のレベル幅ΔＬの細かさ毎に求められる。
【００２３】
フィルタバンク１３０は、δＦの周波数帯域に亘る低周波信号を受けて、さらにｎ個（以下、ｎ（数値）チャンネル、或いはｎ（数値）ＣＨということがある。）のバンドパスフィルタ群により更に細かい周波数帯域（δＦ／ｎ）に切り分ける。例えば、通信回線のチャンネルが１ｋＨｚ毎であれば中心周波数が１ｋＨｚづつズレ、かつ周波数帯域がほぼ１ｋＨｚのバンドパスフィルタをＢＰＦ１〜ＢＰＦｎのｎ個備える。フィルタバンク１３０はＣＰＵを用い、ソフト的にＦＦＴ処理して行える。
【００２４】
レベル検出部２００は、フィルタバンク１３０のバンドパスフィルタのそれぞれに対応した振幅、つまり各チャンネルの振幅を検出するデテクター検波器ＤＥＴ１〜ＤＥＴｎを有し、それぞれチャンネルの周波数成分（バンドパスフィルタの中心周波数に該当）の振幅（レベル）を検出する。
【００２５】
ＡＰＤ部３００は、検波器ＤＥＴ１〜ＤＥＴｎに対応して振幅確率を測定するＡＰＤ検出器であるＡＰＤ１〜ＡＰＤｎ、クロック発生部３１０、及び範囲分類手段３２０を備える。ＡＰＤ１〜ＡＰＤｎのそれぞれの構成、動作は、同じであって、検波器ＤＥＴ１〜ＤＥＴｎの出力をＬＯＧ変換器で対数変換して振幅の単位を「ｄＢ」変換して振幅確率を測定している。ＡＰＤは、ここでは、従来のＡＰＤを採用できる。例えば、特許文献２に記載のＡＰＤ、或いは、特願２００６−２５３８８９号公報に記載のＡＰＤが使用できる。
【００２６】
ＡＰＤ１〜ＡＰＤｎは、振幅確率（分布）を測定する。振幅確率は、同じ大きさの入力信号が所定の時間範囲（期間：Δｔ）内に受ける回数（発生する回数）を表す。入力信号の大きさを区別する細かさは、この例では、Ａ／Ｄ変換手段１２０の細かさΔＬと同じである。確率を測定する時間範囲Δｔは、クロック発生部３１０からのコントロールに従う。クロック発生部３１０は、パルス発生器とタイマーを用いて、例えば時間範囲Δｔとして、１ｓｅｃ間隔、１０秒間隔、或いは１分間隔（いずれも、これらの時間周期のパルスでも良い。）、等のタイミング信号を発生してＡＰＤ１〜ＡＰＤｎに対し、それらの時間範囲内で受ける入力信号の大きさの頻度を測定させることにより、振幅確率を測定させる。操作部７００からの指示で制御部８００を通してタイマーのカウント時間を変更して時間範囲Δｔを変更できる構成にしても良い。
【００２７】
範囲分類手段３２０は、ＡＰＤ１〜ＡＰＤｎがそれぞれ演算して求めた振幅確率の値に応じて、ある範囲に分類する。これは、例えば、１％の分解能で調査しても良いが、その調査目的に応じた分解能で、調査に適切な範囲に分類するために設けた機能である。この実施形態では、操作部７００からの指示で振幅確率が０から１％、１％から１０％、１０％から５０％、５０％から１００％の４段階に分類している。分類の仕方としては、例えば、ＡＰＤ１〜ＡＰＤｎが求めた振幅確率の値と、分類した範囲に相当する複数の閾値とを比較して、比較結果が属する範囲の分類識別札を付す。操作部７００からの指示で制御部８００を通して所望の閾値に変更することにより、分類範囲を変更できる構成にしても良い。なお、ＡＰＤ１〜ＡＰＤｎからは、早ければ、時間範囲Δｔ毎に、更新された振幅確率値が出力される。
【００２８】
記憶部４００は、ＡＰＤ１〜ＡＰＤｎから出力される振幅値とその振幅確率値とを時間経過ともに記憶する。時間経過は、少なくとも、振幅確率を測定するときの時間範囲Δｔ毎に記憶する。さらに、記憶部４００は、範囲分類手段３２０が分類したときの分類識別札を、該当する振幅確率値に対応して記憶する。なお、上記の範囲分類手段３２０は、記憶部４００の入力側（ＡＰＤ部３００側）には無く、出力側（表示制御部５００側）に在って、記憶されている振幅確率値が出力されるときに、その振幅確率値がどの分類に属するか判定し、その判定された分類識別札を添付して出力する構成にしても良い。
【００２９】
表示制御部５００は、表示を実行するために次の（Ｉ）〜（IV）の要素を備える。つまり、各要素として、記憶部４００からの振幅確率値を受けて表示部６００に、（Ｉ）グラフを表示するための手段として、グラフ表示制御手段５１０，振幅／確率グラフ表示制御手段５６０，及び函数尺表示制御手段５７０を備え、（II）数値表示を行うための数値表示制御手段５２０を備え、（III）マーカを表示するために、ピークマーカ表示制御手段５３０，ピークサーチ手段５４０、及び指定マーカ表示制御手段５５０を備え、並びに（IV）上記、各表示を行うための共通の座標情報記憶部５８０及び表示フォーマット記憶部５９０を備える。表示制御部５００は、上記（Ｉ）〜（IV）の各機能を表すプログラムと記憶手段とプログラムを実行させて各機能を実現させるＣＰＵで構成される。なお、座標情報記憶部５８０が記憶している座標情報は、表示されるグラフ、目盛り、各表示欄等の表示画面の位置を特定するのに必要な情報で、表示制御部５００の各要部が、表示させようとしたときに、表示する位置を特定ために用いられる。以下の、各要素での説明では、この座標についての説明は省略する。
【００３０】
上記（Ｉ）におけるグラフ表示制御手段５１０，振幅／確率グラフ表示制御手段５６０，及び函数尺表示制御手段５７０の機能動作とそれらが出力する各グラフは、スイッチＳで切り替えられる。スイッチＳは、操作部７００からの指示で、時間／レベル確率グラフ表示機能、振幅／確率グラフ表示機能、もしくは函数尺表示機能のいずれかが選択設定されたときに、選択設定に対応して、グラフ表示制御手段５１０、振幅／確率グラフ表示制御手段５６０，もしくは函数尺表示制御手段５７０を機能させる。
【００３１】
グラフ表示制御手段５１０は、図２に示すような本発明における基本的なデータを表示する時間／レベル確率表示手段５１３を有する。図２で、縦方向に下から上に向けた４段に、２３．６１ＭＨｚ、２４．６１ＭＨｚ、２５．６１ＭＨｚ、２６．６１ＭＨｚ、２７．６１ＭＨｚの５チャンネルのデータを表示する例である。表示方法は、各チャンネルとも同じであるので、代表例で、一番上段の２３．６１ＭＨｚのチャンネルについて説明する。時間／レベル確率表示手段５１３は、表示フォーマット記憶部５９０から図２のような時間／レベル確率表示フォーマットを取得し、各チャンネルの振幅確率値及び振幅値等のデータを記憶部４００から取得して表示する。
【００３２】
先ず、時間／レベル確率表示手段５１３が表示させる前には、周波数分析部１００，レベル検出部２００，ＡＰＤ部３００，及び記憶部４００のそれぞれが、操作部７００で指定された測定条件、例えば、上記５チャンネルの周波数、全測定時間（データの希望取得時間、或いは希望監視時間）、全表示時間（全測定時間の一部であって表示したい時間、図２の例では、２５０ｓｅｃ間）、確率測定の時間範囲（例えば、Δｔ＝１ｓｅｃ）、０％から１％、１％から１０％、１０％から５０％、５０％から１００％の４段階の分類範囲、と言う条件で、振幅確率測定を開始し、測定した振幅確率値、振幅値（レベル）、分類範囲札を対応させて記憶を開始している。
【００３３】
そして、時間／レベル確率表示手段５１３は、それらの測定条件、及び測定した振幅確率値、振幅値、分類識別札を記憶部４００から受けて、表示フォーマット記憶部５９０から読み出した時間／レベル確率表示フォーマットに合わせて表示部６００へ表示する。そのとき、図２のＣＨ周波数表示欄１０１０には、チャンネルの周波数を数値表示する。横軸表示欄１０２０には、全表示時間を目盛りを振って表示する。これは、５チャンネルとも共通の横軸目盛りである。縦軸表示欄１０３０には、レベル目盛りを１０ｄＢ単位（これは、画面上における振幅値の表示範囲、及び目盛りは、予め初期条件として設定されている。）で表示する。時間／レベル確率表示手段５１３は、これらの横軸―縦軸で示される座標空間１０４０に、記憶部４００に取得時間に対応して記憶される振幅値とその確率値と該当する分類範囲をプロットし、かつ０％から１％、１％から１０％、１０％から５０％、５０％から１００％の４段階に識別して視認可能に分類して表示する。いわば、確率分布の時間変化を表示する。４段階の識別可能な表示の例として、図２の最下段に図に示すように表示の模様別に識別しているが、色分け等の他の方法でも良い。図２では、帯状に確率分布の分類表示している。その最右端のカレント時間位置１０６０が、現在の時間位置の最新のデータを示している。
【００３４】
以下において、横軸が時間軸である場合の説明には、「目盛り」は物理的な寸法を表し（例えば、１ｃｍ間隔の目盛り）、「時間目盛り」は「目盛り」に値付けされた「時間」（例えば、物理的寸法が１ｃｍであるが、それに１ｓｅｃ、或いは１時間、或いは１６：３８：００というように値付けする時間）を言うものとする。又、「横軸表示欄１０２０」を「横軸」と言うことがあり、横軸が時間のときは、「時間軸」ということがある。「横軸の縮小、拡大」とは、物理的な目盛りの縮小、拡大ではなく、値付けされる時間目盛りの縮小、拡大を示す。言い換えれば、「物理的な単位目盛り当たりの時間が小さい、大きい」を言う。
【００３５】
図２では、数値表示欄１０５０には、［振幅（レベル）値、振幅確率値］を［数値ｄＢ、数値％］で表す表示欄に、つまり、左から、０％から１％、１％から１０％、１０％〜５０％、５０％〜１００％の４段階分の表示欄に並べて表示している。これは、数値表示制御手段５２０が、カレント時間位置１０６０における、現在の振幅値と振幅確率値の数値表示を表示させる。カレント時間位置１０６０は、クロック発生部３１０からのクロックで早くとも時間範囲Δｔ経過毎に更新されて、進行する。スイッチＳでグラフ機能が切り替えられてときは、数値表示制御手段５２０がカレント時間位置１０６０ではなく、後記するように図に示すピークマーカ位置、図５に示す指定マーカ位置で測定された振幅値及び確率値を数値表示欄１０５０に表示することがきる。以下、数値表示制御手段５２０の動作、数値表示欄１０５０の表示内容については、各要素の機能説明なお中で説明する。
【００３６】
図２では、５チャンネル分の表示をしているが、実際の測定が、これ以外の他チャンネルについても同時期に測定し、そのデータが記憶部４００に記憶されていれば、グラフ表示制御手段５１０は、上下方向のスクロール機能、或いは上下方向の頁をめくる機能を有することにより、他チャンネルの確率分布も観測できる。
【００３７】
また、横軸表示欄１０２０には、図２のように、０をスタートとして予定されている全表示時間を表示しているが、これはいわばスタートに対する相対時間で時間目盛りを振ってあるが、取得時刻の絶対時間で表示することもできる。記憶部４００が相対時間、及び時計の時刻の双方に対応して振幅確率値等を記憶する構成しておいて、操作部７００からの指示でいずれかを選択して表示できる構成にしておくことにより、実施できる。
【００３８】
時間／レベル確率表示手段５１３が表示する例として、図２を先に説明したが、図２は、後記するフルスパン表示制御手段５１１からの指示で縮小表示したときの例であるが、横軸を有効に使用するため、この実施形態では、時間／レベル確率表示手段５１３はデフォルト表示としては、図３（ａ）の表示フォーマットで表示を行う。図３（ａ）では、全表示時間は操作部７００からは、［現在時間位置（時間範囲Δｔ置きの時間）―５分間］を全表示時間として指定されており、時間／レベル確率表示手段５１３は、タイマーを有し、その指定を受けて、横軸表示欄１０２０の時間目盛りを計算して求め、その結果１６：３５：００〜現在時間（カレント時間位置１０６０）１６：４０：００を全表示時間としている。つまり、横軸表示欄１０２０の右端の時間メモリ１６：４０：００は現在時間であり、時間／レベル確率表示手段５１３は、時間範囲Δｔ毎に横軸全体における、時間目盛りを都度計算して更新している。そしてこの現在時間の変化とともに振幅確率分布を更新している。したがって、時間／レベル確率表示手段５１３は、横軸の目盛り全体と振幅確率分布とを、時間経過とともに、更新している。図３（ａ）の横軸表示欄１０２０における表示では、目盛りの全領域（フルスケールの領域）が、現在までの確率分布表示に有効に使用可能である。これは長時間測定のときに便利である。図２の場合は、全表示時間と、現在のカレント時間位置１０６０における時間とに差があり、その差の時間だけ表示上、空きスペースが生じる。
【００３９】
グラフ表示制御手段５１０は、基本となるグラフは時間／レベル確率表示手段５１３で表示させているが、そのグラフは時間／レベル確率表示手段５１３と連携して動作する、フルスパン表示制御手段５１１及びゾーン拡大表示制御手段５１２を有する。フルスパン（Ｆｕｌｌ  Ｓｐａｎ）表示とは、横軸表示欄１０２０（横軸）のフルスケール（全目盛り）をスタート時間から測定終了予定の全測定時間として表示することを言う。このような表示にあっては、表示されているデータ波形は、横軸の時間目盛りがもっとも縮小された形態で表示される。ただし、ここでのフルスパン表示制御手段５１１は、後記するように縮小率を最大限にすると前記のように全測定時間を横軸の時間目盛りとして表示（フルスパン表示）するが、その縮小率を最大限となる前の途中段階の縮小率も設定可能にされ、その途中断段階に縮小した横軸の目盛りで表示することができる。操作部７００には、グラフ表示制御手段５１０に時間幅の縮小率（ズーム量）を可変指示するためのフルスパン設定キーと、そのフルスパン設定キーが押されている間の時間をカウントし縮小率を示す電気信号として出力するエンコーダとを有する。また、操作部７００には、ゾーン範囲（観測希望の時間幅）を可変指示するためのゾーン拡大表示設定キー（不図示）を有する。ゾーン設定キーとしては、ゾーンの中心時間Ｔｃを設定するキーとゾーン幅となる時間幅Ｔｗと設定するキー（不図示）を有する。或いはゾーンのスタート時間Ｔｓｔとゾーンのストップ時間Ｔｓｐとを設定するキーであっても良い。  
【００４０】
フルスパン表示制御手段５１１は、例えば、図３（ａ）の表示がされているときにフルスパン設定キーで縮小の指示があったとき、グラフの横軸の全目盛りに縮小指示された全表示時間を割り振り、割り振られた時間に対する振幅確率値等のデータを記憶部４００から読み出して、時間／レベル確率表示手段５１３に対して図３（ｂ）のように横軸の時間目盛りを縮小して表示させる。結果として、全表示時間は広がる。図３（ｂ）では、測定開始の時間０〜２０時間分まで横軸の時間範囲を広げている。その分、図３（ｂ）では図３（ａ）のより全表示時間は広がるが、表示される確率分布データの波形は、図３（ａ）の［カレント時間位置―５分間］で表示されている確率分布データの波形より、横軸が縮小して表示される。横軸の時間範囲は、全測定時間（フルスパン測定：Ｆｕｌｌ  Ｓｐａｎ）まで広げることができる。又、それらの逆動作も可能である。つまり図３（ｂ）の表示中にゾーン拡大表示設定キーで拡大指示することにより、図３（ｃ）とすることができる。
【００４１】
フルスパン表示制御手段５１１は、左側のスケール端の測定開始時間０を基準に、設定される縮小率に応じて縮小する。
【００４２】
ゾーン拡大表示制御手段５１２は、操作部７００のゾーン拡大表示定キーからの指示である、ゾーンの中心時間Ｔｃ、ゾーン時間幅Ｔｗを受けて、それらを表すマーカを生成し、既に表示しているグラフの横軸目盛りの位置にゾーンを示すマーカを付す。図３（ｂ）にゾーンマーカを２本の線で表す表示例を示す。２本の線の間隔がゾーン時間幅Ｔｗで、その中心が中心時間Ｔｃである。
【００４３】
図３（ｂ）の表示状態で、操作部７００のゾーン拡大設定キー（確定キー）を設定して確定させると、ゾーン拡大表示制御手段５１２は、時間／レベル確率表示手段５１３に対して指定されたゾーンマーカの範囲［Ｔｃ−Ｔｗ／２］〜「Ｔｃ＋Ｔｗ／２」を新たな横軸の全表示範囲として目盛りを振り直し、その範囲の確率を拡大したグラフ（拡大グラフ）を表示させる（図３（ｃ）を参照）。
【００４４】
また、図３（ｃ）の状態で、フルスパン表示設定キーで縮小の指示を出すことにより、ゾーン拡大表示制御手段５１２による制御を解除し、フルスパン表示制御手段５１１によって横軸を縮小することにより、表示状態を図３（ｂ）の表示に戻させることができる。
【００４５】
以上のことから、グラフ表示制御手段５１０は、つまり、時間／レベル確率表示制御手段５１３、フルスパン表示制御手段５１１、及びゾーン拡大表示制御手段５１２は、次の（ｉ）〜（iii）の動作を行う。（ｉ）デフォルトにおける表示を行う（図３（ａ）参照）。デフォルト表示設定キーで設定される。（ii）縮小表示を行う（図３（ｂ））。フルスパン（全測定時間）表示を行うことができる。フルスパン設定キーで設定される。（iii）ゾーン拡大表示を行う（図３（ｃ））ゾーン拡大設定キーで設定される。そして、グラフ表示制御手段５１０は、操作部７００で設定される各設定キーの設定を受けて、次の（iv）〜（vi）の切り替え動作をスムースに行う。（iv）上記（ｉ）から（ii）への切り替え動作及びその逆の切り替え動作、（ｖ）上記（ii）から（iii）への切り替え動作及びその逆の切り替え動作、（vi）上記（iii）から（ｉ）への切り替え動作及びその逆の切り替え動作。
【００４６】
数値表示欄１０５０における各確率分類における振幅と振幅確率の数値表示は、数値表示制御手段５２０によって行われるが、上記のように、測定中のデータに基づいてグラフ表示している場合に、フルスパン表示制御、ゾーン拡大表示制御を行うとき、グラフ表示の縮小、ゾーン拡大表示と関係なく、上記と同様、現在時間（カレント時間位置１０６０）の測定された振幅値、振幅確率値を数値表示する。測定しておらず、過去データに基づいてグラフ表示しているときに、グラフ表示の縮小、ゾーン拡大表示を行う場合は、フルスパン表示制御手段５１１、ゾーン拡大表示制御手段５１２は測定中でないことを制御部８００から知って、数値表示を行わないことを数値表示制御手段５２０へ指示することができる。ただ、後記するようにピークマーカ、又は指定マーカが設定されているときは、その指定マーカの位置における振幅値及び振幅確率値の数値表示を行う。尚、測定中でも数値表示を中止にすることは可能である。数値表示制御手段５２０は、表示制御部５００内の各要部と連携して動作する。一部については、各要部を説明するときに合わせて数値表示について説明する。
【００４７】
ピークサーチ手段５４０は、図２の表示がされているときに、操作部７００から、振幅確率値の分類範囲の指定とピークマーカ機能をオンする旨の指示を受けて、表示されている時間範囲内における、指定された分類範囲内の確率振幅値と振幅値とを記憶部４００から読み出して、現在表示されている横軸の全表示時間内、かつ指定された分類範囲で、最大の振幅値を示すときの、対応する振幅確率値、及び時間位置を求める。そして、ピークマーカ表示制御手段５３０は、表示されているグラフにピークサーチ手段５４０が求めた最大の振幅値とそれに対応する時間位置に、視認可能な形状のマーカ（表示フォーマット記憶部５９０に記憶）を付す。又、合わせて、最大の振幅値に対応する振幅確率値を求めておく。
【００４８】
ピークマーカ機能がオンしているとき、数値表示制御手段５２０は、指定された振幅確率値が属する分類範囲を担当する数値表示欄１０５０に、ピークサーチ手段５４０が求めた、最大の振幅値とそれに対応する振幅確率値を、ピーク機能であることを認識可能に数値表示する。そのとき、数値表示制御手段５２０は、指定された分類範囲以外の数値表示欄には、ピークマーカのある時間位置（指定された分類範囲で最大の振幅値を示す時間位置）における振幅値と振幅確率値を表示する。これらのピークマーカの表示例を図４に示す。図４は、指定された分類範囲が１０％＜ＡＰＤ≦５０％の場合で、ＣＨ周波数が２７．６１ＭＨｚの場合は、最大振幅値＝―４５．０ｄＢ、振幅確率値＝１０．２％と、白抜きで数値表示されている。５チャンネルあれば、一度に同じピークマーカ機能を実行することができる。
【００４９】
したがって、ピークマーカは次の（イ）〜（ニ）に示す位置を指すように付される。（イ）指定された分類範囲が０％＜ＡＰＤ≦１％の場合は、この範囲で最大となる振幅値と、そのときの時間位置とで定まる座標を指す。（ロ）指定された分類範囲が１％＜ＡＰＤ≦１０％の場合は、この範囲で最大となる振幅値と、そのときの時間位置とで定まる座標を指す。（ハ）指定された分類範囲が１０％＜ＡＰＤ≦５０％の場合は、この範囲で最大となる振幅値と、そのときの時間位置とで定まる座標を指す。（ニ）指定された分類範囲が５０％＜ＡＰＤ≦１００％の場合は、この範囲で最大となる振幅値と、そのときの時間位置とで定まる座標を指す。
【００５０】
指定マーカ表示制御手段５５０は、操作部７００から現在表示中のグラフの横軸である全表示時間内で時間位置の指定を受けて、表示されているグラフ上の指定された時間位置に、視認可能な形状のマーカ（表示フォーマット記憶部５９０に記憶）を付す。そして記憶部４００からその指定された時間位置における、各分類範囲の振幅値、振幅確率値を読み出して、数値表示制御手段５２０に数値表示させる。表示例を図５に示す。図ででは、一つの時間位置の指定で、各チャンネル同じ時間位置に指定マーカを設定した例である。なお、指定マーカは、操作部７００からの指示で移動させて設定できる。
【００５１】
測定中に、指定マーカ、ピークマーカ等の指示がなければ、指定マーカ表示制御手段５５０は、現在時間位置（カレント位置）にマーカ位置をセットする。この場合は、グラフ上では必ずしも視認可能な形態のマーカを付す必要はない。図２のようにデータの右端がカレント時間位置１０６０と認識できるからである。このとき数値表示は、カレント時間位置１０６０における）（現在時間における）測定データに基づく表示が行われる。
【００５２】
上記の指定マーカ機能及びピークマーカ機能のそれぞれは、グラフ表示制御手段５１０による上記（ｉ）〜（iii）のいずれの表示状態であっても実施できる。つまり、上記のゾーン拡大表示機能やズーム表示機能のそれぞれにおいても実施できる。指定マーカ機能及びピークマーカ機能は、ゾーン拡大表示機能及びズーム表示機能とは取り扱うデータは共有するが動作は独立しているからである。指定マーカ機能とピークマーカ機能の同時併用については、数値表示欄１０５０が限られているためこの実施形態では併用していないが、切り替えて使用すれば可能である。
【００５３】
指定マーカ機能及びピークマーカ機能は、表示が５チャンネル分表示しているとき、それぞれのチャンネルで独立して動作可能である（図５は、異なった時間位置に指定マーカを付した例である。）。ただ、指定マーカ機能及びピークマーカ機能は、観察する設定条件、例えば、指定マーカ機能における時間位置、或いはピークマーカ機能における分類範囲は各チャンネルとも同じく設定することが好ましい。なぜなら例えば、観測の時間位置を設定するとき、各チャンネルで別々な時間で観察する理由がなく、むしろ傾向を知るには同一条件が好ましいからである。
【００５４】
振幅／確率グラフ表示制御手段５６０は、例えば、横軸をレベル（振幅）、縦軸を振幅確率とした座標空間に測定データを表示する（横軸を振幅確率、縦軸をレベルとしても良い。）。具体的には、図５のように所望の時間位置に指定マーカを付して、操作部７００における「振幅／確率表示」キーを設定されることにより、振幅／確率グラフ表示制御手段５６０は、図６のような横軸をレベル、縦軸を振幅確率とした座標空間の表示フォーマットを表示フォーマット記憶部５９０から読み出して表示するとともに、その表示空間に図５における指定マーカで指定された所望の時間位置における振幅確率値、レベル（振幅値）を記憶部４００から受けてプロットすることにより、図６、図７のような振幅／確率のグラフ表示を行う。図６は、確率測定の時間範囲Δｔ＝１ｓｅｃ、で測定経過時間（指定マーカの時間位置）が約６０秒後（図５の指定マーカの時間位置）の振幅／確率のグラフで、同一時間位置における５チャンネル分のデータを識別可能に表示している。
【００５５】
なお、上記例では、振幅／確率グラフ表示制御手段５６０は、図５の表示状態から指定マーカで指定された時間における各チャンネルの振幅／確率のグラフを表示したが、ズーム表示機能で縮小された表示、ゾーン拡大表示機能で拡大された表示の各表示状態で、指定マーカにより指定された時間位置における振幅／確率のグラフを表示することは、上記と同様に行うことができる。又、図５等で指定マーカは指定されていないが、カレント時間位置１０６０で数値表示機能が動作しているとき（表示部６００に指定マーカが表示されてはいないが、実質的に指定マーカがカレント時間位置１０６０に設定されているとき）に、「振幅／確率表示」キーを設定されることにより、カレント時間位置１０６０における振幅／確率のグラフを表示することができる。
【００５６】
さらに、振幅／確率グラフ表示制御手段５６０が振幅／確率のグラフを表示しているときに、振幅確率値又は振幅（レベル）値を指定した指定マーカの設定を操作部７００から受けたときは、指定マーカ表示制御手段５５０は指定された振幅確率値又は振幅（レベル）値の位置にマーカを付すとともに、指定された振幅確率値に該当する振幅値、又は指定された振幅値に該当する振幅確率値を数値表示制御手段５２０により数値表示欄１０７０に表示させる。図６で、振幅確率値（ＡＰＤ）として０．１％を指定して指定マーカを設定したときの例である。数値表示欄１０７０に、各チャンネルでの振幅確率＝１％のときの振幅値を表示している。図７は、振幅（レベル）としてー４０ｄＢを指定して指定マーカを設定したときの例である。数値表示欄１０７０に、各チャンネルでの振幅値＝―４０ｄＢのときの振幅確率値を表示している例である。なお、図６は、図５の時間位置が６０ｓｅｃの時点（カレント時間位置１０６０でない時間）での測定結果に基づく表示であり、確率値、振幅値の更新はない。図５で、カレント時間位置１０６０に指定マーカで指定されて切り替えられたときは、図６の振幅／確率グラフ表示制御手段５６０が、時間範囲Δｔ＝１ｓｅｃ間隔でデータを記憶部４００から振幅値を受けて図６のデータを更新する。図７におけるデータの更新の態様は、図６と同様である。
【００５７】
上記のことから、指定マーカは、次の機能を有する。
・数値表示の時間位置を視覚的に示す（図５，図３等）。
カレント時間位置は、指定マーカで指定しても良いが、
指定しない場合の数値表示は、カレント時間位置を示す。
・時間／振幅（レベル）グラフから表示振幅／確率グラフへの切り替え時の
橋渡しの基準となる時間位置を示す（図５から図６への表示切り替え）。
・振幅値又は振幅確率値を指定し、振幅／確率グラフ上の位置を視覚的に示す。
【００５８】
函数尺表示制御手段５７０は、測定対象とする信号について予測される振幅確率分布に関する函数尺に基づく座標空間を表示フォーマットとして、記憶部４００から受けた振幅確率値をプロットして表示する。このような予測される分布の函数尺を適用した座標に実際の測定した振幅確率値をプロットすれば、測定対象の信号が予測された分布に適合するかどうか判断が容易にできる。一般には、予測される分布としては、測定する信号が雑音であるとすると、レイリー分布、正規分布、指数分布、χ２分布等がある。函数尺を求めることについては、本出願人に係る技術であって特許第３３７４１５４号公報に記載された技術がある。測定対象が雑音として予測される分布がレイリー分布の場合における振幅確率／振幅（レベル）のグラフ表示についての表示例を図８に示す。
函数尺がこのレイリー分布の場合の函数尺表示について、図６と図８との違いで言えば、横軸のスケール、つまり振幅確率のスケールが変換されていて異なる。スケール変換は次のように行われる。
【００５９】
変換前のスケール範囲［ａ−ｂ］にある振幅確率値ｘ（１より小さく０より大）をレイリー分布スケール［ｍ−ｎ］間のＸの位置に表示するとしたとき、位置Ｘは次の式で示される。
Ｘ＝｛（ｎ−ｍ）／（ｆ（ｂ）―ｆ（ａ））｝｛ｆ（ｘ）―ｆ（ａ）｝＋ｍ
ただし、ｆ（ｘ）＝ｌｏｇ１０ｌｏｇ１０（１／ｘ）
したがって、このような計算をして作成したレイリー分布の函数尺を適用した表示フォーマットを表示フォーマット記憶部５９０に予め記憶しておく。
【００６０】
関数尺表示は、一例は、振幅／確率グラフ表示制御手段５６０が、例えば図６を表示している場合に、操作部７００からの指示で函数尺表示機能が設定されたとき函数尺表示制御手段５７０が図８の表示を行う。又、他の例は、グラフ表示制御手段５１０が、例えば、図５の表示をしているときに操作部７００からの指示で函数尺表示機能が設定されたときに（このときは、振幅／確率グラフ表示機能も選択可能である）函数尺表示制御手段５７０が図８の表示を行う。
【００６１】
具体的な動作の一例を説明する。図６のような表示が行われている状態で、操作部７００において「レイリー函数尺表示」機能を設定されることにより、函数尺表示制御手段５７０は、横軸を振幅、縦軸を振幅確率の目盛りを図８のようなレイリー函数尺の表示フォーマットを表示フォーマット記憶部５９０から読み出して表示するとともに、その表示空間に図５における指定マーカで指定された所望の時間位置における振幅確率値、振幅値（レベル）を記憶部４００から受けて、プロットすることにより、図８のようなレイリー分布を函数尺とした振幅／確率のグラフ表示を行う。
【００６２】
また、函数尺表示制御手段５７０が機能中であっても、振幅値、又は振幅確率値を指定した指定マーカ機能を動作させることができる。図８では、振幅確率値が６０％を指定した場合であり、横軸の６０％の位置にマーカが付され、その時の、各チャンネルの振幅値が数値表示されている。これは、数値表示制御手段５２０がマーカで指定されたときに函数尺表示制御手段５７０によってが記憶部４００からサーチして求められた、振幅確率値６０％における各チャンネルの振幅値を、数値表示制御手段５２０が数値表示欄１０５０に表示することによって行われる。
【００６３】
本発明は、上記のように時間／振幅（レベル）座標系における振幅確率のグラフ表示が行え、さらにその表示にあたっては、縮小表示機能、ゾーンマーカ機能又は拡大表示機能、それらにピークマーカ機能、指定マーカ機能が実行できる。また、時間／振幅（レベル）座標系における振幅確率のグラフ表示から、振幅／振幅確率座表系でのグラフ表示、又は函数尺グラフ表示等のグラフ表示が実行でき、さらにそれらの表示において、指定マーカ機能を動作せることができる。いずれも多チャンネルに実施できる。
【図面の簡単な説明】
【００６４】
【図１】ＡＰＤ測定装置に係る実施形態の全体の機能構成を示す図である。
【図２】縦軸が振幅、横軸が時間の座標上での確率分布の表示例を示す図である。
【図３】フルスパン表示機能、並びにゾーンマーカ機能及びゾーン拡大表示機能を説明するための図である。
【図４】ピークマーカ機能を説明するための図である。
【図５】指定マーカ機能を説明するための図である。
【図６】振幅確率値を指定時の振幅／確率の表示例である。
【図７】振幅値を指定時の振幅／確率の表示例である。
【図８】函数尺でなる座標上での振幅／確率の表示例である。
【符号の説明】
【００６５】
１００  周波数分析部、  １１０  帯域選択手段、  １２０  Ａ／Ｄ変換手段、
１３０  フィルタバンク、  ２００  レベル検出部、  ３００  ＡＰＤ部、
３１０  クロック発生部、  ３２０  範囲分類手段、  ４００  記憶部、
５００  表示制御部、  ５１０  グラフ表示制御手段、  ５１１  フルスパン表示制御手段、５１２  ゾーン拡大表示制御手段、  ５１３  時間／レベル確率表示手段、  
５２０  数値表示制御手段、  ５３０  ピークマーカ表示制御手段、
５４０  ピークサーチ手段、  ５５０  指定マーカ表示制御手段、
５６０  振幅／確率グラフ表示制御手段、  ５７０  函数尺表示制御手段、
５８０  座表情報記憶手段、  ５９０  表示フォーマット記憶部、  ６００  表示部、７００  操作部、
１０１０  ＣＨ周波数表示欄、  １０２０  横軸表示欄、
１０３０  縦軸表示欄、  １０４０  座標空間、  １０５０  数値表示欄、
１０６０  カレント時間位置、  １０７０  数値表示欄
Ｓ：スイッチ
【図１】
2009133683.tif 000003
【図２】
2009133683.tif 000004
【図３】
2009133683.tif 000005
【図４】
2009133683.tif 000006
【図５】
2009133683.tif 000007
【図６】
2009133683.tif 000008
【図７】
2009133683.tif 000009
【図８】
2009133683.tif 000010