JP2911838B2 - 圧力式波高計 - Google Patents
圧力式波高計Info
- Publication number
- JP2911838B2 JP2911838B2 JP31302896A JP31302896A JP2911838B2 JP 2911838 B2 JP2911838 B2 JP 2911838B2 JP 31302896 A JP31302896 A JP 31302896A JP 31302896 A JP31302896 A JP 31302896A JP 2911838 B2 JP2911838 B2 JP 2911838B2
- Authority
- JP
- Japan
- Prior art keywords
- water pressure
- frequency
- wave height
- water
- equation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Measurement Of Levels Of Liquids Or Fluent Solid Materials (AREA)
Description
【０００１】
【発明の属する技術分野】本発明は圧力式波高計の技術
分野に属する。
分野に属する。
【０００２】
【従来の技術】圧力式波高計は、水中に設置された水圧
計で計測した圧力変化から水表面の高さの変化即ち表面
波を計測するものであり、その原理は水面が高くなれば
水圧も高くなり、水面が低くなれば水圧も低くなるとい
う関係があることを利用したものである。
計で計測した圧力変化から水表面の高さの変化即ち表面
波を計測するものであり、その原理は水面が高くなれば
水圧も高くなり、水面が低くなれば水圧も低くなるとい
う関係があることを利用したものである。
【０００３】ところで水面高の上下の変動と水中圧力の
高低変化との関係は、水面上下動の周期が長い（即ち変
化がゆっくりしている）場合は、水中における単位面積
当たりの圧力は、水面が高くなった場合単位面積当たり
で高くなった水の重量分だけ増加し、水面が低くなった
場合は単位面積当たりで低くなった水の重量分だけ減少
する。従って、水中圧力変化を水柱高の変化で表した値
に一致する。即ち、水柱高の変化に１を掛けた値という
ことになる。
高低変化との関係は、水面上下動の周期が長い（即ち変
化がゆっくりしている）場合は、水中における単位面積
当たりの圧力は、水面が高くなった場合単位面積当たり
で高くなった水の重量分だけ増加し、水面が低くなった
場合は単位面積当たりで低くなった水の重量分だけ減少
する。従って、水中圧力変化を水柱高の変化で表した値
に一致する。即ち、水柱高の変化に１を掛けた値という
ことになる。
【０００４】しかしながら、水面高の変化の周期が短く
（即ち周波数が高く）なってくると、水面高の変化に対
する水圧変化が小さくなってくる。この小さくなってく
る割合即ち減衰率は理論的に定まっており次のようにな
る。今、周波数をｆ（＝１／Ｔ、但しＴは周期）、圧力
センサ直下の水深をＤ、圧力センサの深度をｚ、重力加
速度をｇとすれば、数式３からＬを求め、このＬを数式
４に代入して、減衰率ηが周波数ｆの関数として求めら
れる（詳しくは朝倉書店発行 井島武士著「海岸工学」
３４頁参照）
（即ち周波数が高く）なってくると、水面高の変化に対
する水圧変化が小さくなってくる。この小さくなってく
る割合即ち減衰率は理論的に定まっており次のようにな
る。今、周波数をｆ（＝１／Ｔ、但しＴは周期）、圧力
センサ直下の水深をＤ、圧力センサの深度をｚ、重力加
速度をｇとすれば、数式３からＬを求め、このＬを数式
４に代入して、減衰率ηが周波数ｆの関数として求めら
れる（詳しくは朝倉書店発行 井島武士著「海岸工学」
３４頁参照）
【０００５】
【数３】
【０００６】
【数４】
【０００７】従って、水圧変化の周波数がｆであるとき
はその水圧変化に対応する水柱高の変化の振幅に減衰率
η（ｆ）の逆数Ｈ（ｆ）を補正倍数として乗ずれば正し
い水面高の変化の振幅（即ち波高）が得られることにな
る。そこで、従来は、時間を横軸として水圧センサから
得られた水圧波形からゼロアップクロス法で１波１波分
離してその周期Ｔを計測し、その周期Ｔを用いて数式３
および数式４から補正倍数Ｈ（ｆ）を求め、水圧波形の
振幅に乗じて波高を求めていた。
はその水圧変化に対応する水柱高の変化の振幅に減衰率
η（ｆ）の逆数Ｈ（ｆ）を補正倍数として乗ずれば正し
い水面高の変化の振幅（即ち波高）が得られることにな
る。そこで、従来は、時間を横軸として水圧センサから
得られた水圧波形からゼロアップクロス法で１波１波分
離してその周期Ｔを計測し、その周期Ｔを用いて数式３
および数式４から補正倍数Ｈ（ｆ）を求め、水圧波形の
振幅に乗じて波高を求めていた。
【０００８】
【発明が解決しようとする課題】しかしながら、水圧波
形から周期を正確に読み取ることがむずかしく、補正倍
数が２．５以上では、一般に正しい水圧−波高換算がで
きない。このときの波の周期は約７秒である。従って、
従来の水圧式波高計では７秒以下の周期の波は計測でき
ないという問題があった。
形から周期を正確に読み取ることがむずかしく、補正倍
数が２．５以上では、一般に正しい水圧−波高換算がで
きない。このときの波の周期は約７秒である。従って、
従来の水圧式波高計では７秒以下の周期の波は計測でき
ないという問題があった。
【０００９】本発明の目的は、上記従来の問題点に鑑み
て、ゼロアップクロス法で周期を求めるという方法によ
らずに、水圧波形の含む各周波数成分に対して、それぞ
れの周波数に対応した補正倍数を乗ずる周波数伝達関数
を有するディジタルフィルターを通過させることにより
補正された波高値が得られる圧力式波高計を提供するこ
とにある。
て、ゼロアップクロス法で周期を求めるという方法によ
らずに、水圧波形の含む各周波数成分に対して、それぞ
れの周波数に対応した補正倍数を乗ずる周波数伝達関数
を有するディジタルフィルターを通過させることにより
補正された波高値が得られる圧力式波高計を提供するこ
とにある。
【００１０】
【課題を解決するための手段】本発明は、上記の目的を
達成するために次の手段構成を有する。本発明の圧力式
波高計の構成は、次の各手段を具備するものである。 （イ）水面の波高の変化によって変化する水圧を検知し
アナログ信号で出力する水圧センサ （ロ）水圧センサからの水圧アナログ信号をディジタル
信号に変換するＡ／Ｄ変換器 （ハ）Ａ／Ｄ変換器からの水圧信号の、予め定めた周波
数範囲の周波数成分に対し、数式５の周波数伝達関数
（補正倍率）Ｈ（ｆ）を乗じて出力し、入力信号周波数
ｆに依存して生起する減衰誤差を補正した波高情報のリ
アルタイム出力を確保するディジタルフィルター但し、
数式５のＬは、入力信号周波数をｆとしたとき次の数式
６から得られる値を代入する。
達成するために次の手段構成を有する。本発明の圧力式
波高計の構成は、次の各手段を具備するものである。 （イ）水面の波高の変化によって変化する水圧を検知し
アナログ信号で出力する水圧センサ （ロ）水圧センサからの水圧アナログ信号をディジタル
信号に変換するＡ／Ｄ変換器 （ハ）Ａ／Ｄ変換器からの水圧信号の、予め定めた周波
数範囲の周波数成分に対し、数式５の周波数伝達関数
（補正倍率）Ｈ（ｆ）を乗じて出力し、入力信号周波数
ｆに依存して生起する減衰誤差を補正した波高情報のリ
アルタイム出力を確保するディジタルフィルター但し、
数式５のＬは、入力信号周波数をｆとしたとき次の数式
６から得られる値を代入する。
【００１１】
【数５】
【００１２】
【数６】
【００１３】
【発明の実施の形態】本発明で用いられるディジタルフ
ィルターはその周波数伝達関数が、周波数の関数である
補正倍率になるように設定されているので、ディジタル
信号化された水圧信号を入力するとそれに含まれている
周波数成分毎にその周波数に対応する補正倍率が乗ぜら
れそれらが合成されたと同じ形で自動的に出力され結果
として補正された波高値が得られることになる。
ィルターはその周波数伝達関数が、周波数の関数である
補正倍率になるように設定されているので、ディジタル
信号化された水圧信号を入力するとそれに含まれている
周波数成分毎にその周波数に対応する補正倍率が乗ぜら
れそれらが合成されたと同じ形で自動的に出力され結果
として補正された波高値が得られることになる。
【００１４】ディジタルフィルターの特性は、インパル
ス応答ｈ（ｔ）の時系列データで表される。インパルス
は、ある時刻で値が１で、その前後の値がゼロである時
系列を言う。インパルス応答ｈ（ｔ）は、このインパル
スを入力したときの応答である。もしフィルターのイン
パルス応答ｈ（ｔ）が得られると、次のような出力が計
算される。即ち、入力の時系列をサンプル値の大きさを
持つインパルス列と見なすことができ、サンプル時間間
隔をΔｔとすると、Δｔ時間ずらして、このインパルス
応答を入力することになる。従って、このインパルス列
の各インパルス応答ｈ（ｔ）を時間をずらして加えれ
ば、フィルター出力が得られる。この操作を「畳み込
み」と呼ぶ。結論的に、「ディジタルフィルターを作る
ことは、インパルス応答ｈ（ｔ）を求めること」と言う
ことができる。
ス応答ｈ（ｔ）の時系列データで表される。インパルス
は、ある時刻で値が１で、その前後の値がゼロである時
系列を言う。インパルス応答ｈ（ｔ）は、このインパル
スを入力したときの応答である。もしフィルターのイン
パルス応答ｈ（ｔ）が得られると、次のような出力が計
算される。即ち、入力の時系列をサンプル値の大きさを
持つインパルス列と見なすことができ、サンプル時間間
隔をΔｔとすると、Δｔ時間ずらして、このインパルス
応答を入力することになる。従って、このインパルス列
の各インパルス応答ｈ（ｔ）を時間をずらして加えれ
ば、フィルター出力が得られる。この操作を「畳み込
み」と呼ぶ。結論的に、「ディジタルフィルターを作る
ことは、インパルス応答ｈ（ｔ）を求めること」と言う
ことができる。
【００１５】このインパルス応答ｈ（ｔ）をフーリエ変
換とすると、周波数伝達関数Ｈ（ｆ）が得られる。ま
た、逆に周波数伝達関数Ｈ（ｆ）を与えて、逆フーリエ
変換すると、インパルス応答ｈ（ｔ）が得られる。この
周波数伝達関数Ｈ（ｆ）は、振幅が１である各周波数を
持つ正弦波を、フィルターに入力したときの出力の振幅
の大きさ（倍率）を表すものである。上述のＨ（ｆ）
は、周波数伝達関数Ｈ（ｆ）を表しているので、このＨ
（ｆ）を逆フーリエ変換してインパルス応答ｈ（ｔ）を
求めることにより、ディジタルフィルターを作ることが
できる。
換とすると、周波数伝達関数Ｈ（ｆ）が得られる。ま
た、逆に周波数伝達関数Ｈ（ｆ）を与えて、逆フーリエ
変換すると、インパルス応答ｈ（ｔ）が得られる。この
周波数伝達関数Ｈ（ｆ）は、振幅が１である各周波数を
持つ正弦波を、フィルターに入力したときの出力の振幅
の大きさ（倍率）を表すものである。上述のＨ（ｆ）
は、周波数伝達関数Ｈ（ｆ）を表しているので、このＨ
（ｆ）を逆フーリエ変換してインパルス応答ｈ（ｔ）を
求めることにより、ディジタルフィルターを作ることが
できる。
【００１６】原理的には以上の手法でディジタルフィル
ターを求めるが、実際の計算は次の手法によっている。
従来から、周波数伝達関数Ｈ（ｆ）を与えて、インパル
ス応答ｈ（ｔ）を単純に、逆フーリエ変換で求めるとう
まく行かないことが知られていた。周波数伝達関数Ｈ
（ｆ）を与えて、インパルス応答ｈ（ｔ）を単純に計算
するとインパルス応答ｈ（ｔ）が無限の長さとなる。こ
れを適当な長さで切ると、希望のフィルター特性が得ら
れない。これはＦＦＴ（高速フーリエ変換）を用いてフ
ーリエ変換をし、必要な部分を残して逆変換し、バンド
パスフィルターを通したと考える無謀なことと同じで、
正しくフィルターリングされない。
ターを求めるが、実際の計算は次の手法によっている。
従来から、周波数伝達関数Ｈ（ｆ）を与えて、インパル
ス応答ｈ（ｔ）を単純に、逆フーリエ変換で求めるとう
まく行かないことが知られていた。周波数伝達関数Ｈ
（ｆ）を与えて、インパルス応答ｈ（ｔ）を単純に計算
するとインパルス応答ｈ（ｔ）が無限の長さとなる。こ
れを適当な長さで切ると、希望のフィルター特性が得ら
れない。これはＦＦＴ（高速フーリエ変換）を用いてフ
ーリエ変換をし、必要な部分を残して逆変換し、バンド
パスフィルターを通したと考える無謀なことと同じで、
正しくフィルターリングされない。
【００１７】しかし最近は、周波数伝達関数Ｈ（ｆ）か
らインパルス応答ｈ（ｔ）を計算するとき、窓関数を乗
じる手法「窓関数法」が普及し、もっぱらこの手法が用
いられている。無限長になるインパルス応答ｈ（ｔ）を
窓関数を用いて、インパルス応答ｈ（ｔ）の両端を滑ら
かにゼロに収束させる手法である。このようにすると、
正しくフィルターリングされ、位相ずれのないＦＩＲ型
フィルターが得られる。
らインパルス応答ｈ（ｔ）を計算するとき、窓関数を乗
じる手法「窓関数法」が普及し、もっぱらこの手法が用
いられている。無限長になるインパルス応答ｈ（ｔ）を
窓関数を用いて、インパルス応答ｈ（ｔ）の両端を滑ら
かにゼロに収束させる手法である。このようにすると、
正しくフィルターリングされ、位相ずれのないＦＩＲ型
フィルターが得られる。
【００１８】しかし、フィルター長を有限長にしたた
め、希望の周波数伝達関数Ｈ（ｆ）とは完全に一致しな
い場合がある。一般にこのフィルター長の長さを長くす
ればするほど、希望の周波数伝達関数Ｈ（ｆ）に近づ
く。しかし長くすると、計算量が多くなるため、試行錯
誤しながら最適な長さを決定する。
め、希望の周波数伝達関数Ｈ（ｆ）とは完全に一致しな
い場合がある。一般にこのフィルター長の長さを長くす
ればするほど、希望の周波数伝達関数Ｈ（ｆ）に近づ
く。しかし長くすると、計算量が多くなるため、試行錯
誤しながら最適な長さを決定する。
【００１９】
【実施例】以下、本発明の圧力式波高計の実施例を図面
を参照して説明する。図１は、本発明の第１の実施例を
示すブロック図であり、水圧センサ１、Ａ／Ｄ変換器
２、ディジタルフィルター３および波高データメモリ４
がまとまって水中に設置される。水圧センサ１で検出さ
れた水圧のアナログデータはＡ／Ｄ変換器２でディジタ
ル信号に変換された後、ディジタルフィルター３へ入力
される。
を参照して説明する。図１は、本発明の第１の実施例を
示すブロック図であり、水圧センサ１、Ａ／Ｄ変換器
２、ディジタルフィルター３および波高データメモリ４
がまとまって水中に設置される。水圧センサ１で検出さ
れた水圧のアナログデータはＡ／Ｄ変換器２でディジタ
ル信号に変換された後、ディジタルフィルター３へ入力
される。
【００２０】図５は、ディジタルフィルター３の周波数
特性を示す図である。横軸は周波数、縦軸は補正倍数で
ある。この図は水深１２ｍの水底に水圧センサが置かれ
た場合の特性で、周波数ゼロから約０．２９Ｈｚ（周期
で言うならば約３．５秒）付近の所までの曲線が数式５
で表される曲線である。本実施例では波高計が水底に置
かれているのでＤ＝ｚとなり、その結果数式５の分母は
１となる。数式５に従えばこの曲線は急峻に上昇し大き
な補正倍数となる。しかし大きな補正倍数を乗ずること
は不要なノイズをも増大させることになりかえってＳ／
Ｎを低下させることが考えられるので３０倍で飽和する
特性にした。
特性を示す図である。横軸は周波数、縦軸は補正倍数で
ある。この図は水深１２ｍの水底に水圧センサが置かれ
た場合の特性で、周波数ゼロから約０．２９Ｈｚ（周期
で言うならば約３．５秒）付近の所までの曲線が数式５
で表される曲線である。本実施例では波高計が水底に置
かれているのでＤ＝ｚとなり、その結果数式５の分母は
１となる。数式５に従えばこの曲線は急峻に上昇し大き
な補正倍数となる。しかし大きな補正倍数を乗ずること
は不要なノイズをも増大させることになりかえってＳ／
Ｎを低下させることが考えられるので３０倍で飽和する
特性にした。
【００２１】このディジタルフィルターは、入力波が連
続的に入力されていると、同時に出力波も連続的に出力
されるリアルタイム型フィルターである。従って、入力
波が無限長でもかまわない。このようなディジタルフィ
ルター３を通すことによりその出力は水面の波高データ
となり波高データメモリ４へ格納される。こうして或る
時間あるいは或る期間の波高データが蓄積されると引き
上げられて波高データメモリ４から波高データが読み出
される。
続的に入力されていると、同時に出力波も連続的に出力
されるリアルタイム型フィルターである。従って、入力
波が無限長でもかまわない。このようなディジタルフィ
ルター３を通すことによりその出力は水面の波高データ
となり波高データメモリ４へ格納される。こうして或る
時間あるいは或る期間の波高データが蓄積されると引き
上げられて波高データメモリ４から波高データが読み出
される。
【００２２】図２は、第２の実施例を示すブロック図で
あり、ディジタルフィルター３は水中設置機器の中には
含まれていない。水中設置の水圧センサ１や、Ａ／Ｄ変
換器２の機能は第１の構成の場合と同様であるが、ディ
ジタルデータに変換された水圧データは水圧データメモ
リ５に格納される。こうして或る時間あるいは或る期間
の水圧データが蓄積されると引き上げて水圧データメモ
リ５から水圧データを読み出してディジタルフィルター
３にかけて波高データを得ることになる。
あり、ディジタルフィルター３は水中設置機器の中には
含まれていない。水中設置の水圧センサ１や、Ａ／Ｄ変
換器２の機能は第１の構成の場合と同様であるが、ディ
ジタルデータに変換された水圧データは水圧データメモ
リ５に格納される。こうして或る時間あるいは或る期間
の水圧データが蓄積されると引き上げて水圧データメモ
リ５から水圧データを読み出してディジタルフィルター
３にかけて波高データを得ることになる。
【００２３】図３は、第３の実施例を示すブロック図で
ある。水中設置のディジタルフィルター３からの出力信
号は、水中通信ケーブル６を経由して、陸上に設置され
た波高データメモリ４に転送される。図４は、第４の実
施例を示すブロック図である。水中設置のＡ／Ｄ変換器
２からの出力信号は、水中通信ケーブル６を経由して陸
上に設置されたディジタルフィルター３に入力され、そ
の出力は、波高データとして波高データメモリ４に記憶
される。このように、水中通信ケーブルによるデータ転
送を行うと、リアルタイムに波高のデータが入手でき、
港湾工事や海水浴場の安全管理に利用できる。また、水
中通信ケーブル６の代わりに、ブイと無線通信機を用い
て構成することもできる。
ある。水中設置のディジタルフィルター３からの出力信
号は、水中通信ケーブル６を経由して、陸上に設置され
た波高データメモリ４に転送される。図４は、第４の実
施例を示すブロック図である。水中設置のＡ／Ｄ変換器
２からの出力信号は、水中通信ケーブル６を経由して陸
上に設置されたディジタルフィルター３に入力され、そ
の出力は、波高データとして波高データメモリ４に記憶
される。このように、水中通信ケーブルによるデータ転
送を行うと、リアルタイムに波高のデータが入手でき、
港湾工事や海水浴場の安全管理に利用できる。また、水
中通信ケーブル６の代わりに、ブイと無線通信機を用い
て構成することもできる。
【００２４】図６は、以上の実施例において得られた波
形を示す図である。（ａ）は水圧センサ１から得られる
水圧変動波形である。（ｂ）は（ａ）の波形信号をディ
ジタルフィルター３に通してその出力から得られた波形
である。（ｃ）は正しい波高データが得られると考えら
れている超音波式波高計で測定された同じ場所の水面に
おける波高データである。ここで（ｂ）の波形を（ｃ）
の波形と比較すると、ほぼ正しく復元されていることが
分かる。
形を示す図である。（ａ）は水圧センサ１から得られる
水圧変動波形である。（ｂ）は（ａ）の波形信号をディ
ジタルフィルター３に通してその出力から得られた波形
である。（ｃ）は正しい波高データが得られると考えら
れている超音波式波高計で測定された同じ場所の水面に
おける波高データである。ここで（ｂ）の波形を（ｃ）
の波形と比較すると、ほぼ正しく復元されていることが
分かる。
【００２５】以上の結果により、ディジタルフィルター
を用いることによって、補正倍数が３０程度まで、周期
でいうならば周期３．５秒程度の周期の短い波まで計測
できることになり、従来のゼロアップクロス法で補正倍
数は精々２．５までで７秒以下の周期の波が計測できな
かったのに較べてより正確な復元ができることになる。
また、水圧データ波形から１波１波分離して周期を求め
るという作業も不要となる。
を用いることによって、補正倍数が３０程度まで、周期
でいうならば周期３．５秒程度の周期の短い波まで計測
できることになり、従来のゼロアップクロス法で補正倍
数は精々２．５までで７秒以下の周期の波が計測できな
かったのに較べてより正確な復元ができることになる。
また、水圧データ波形から１波１波分離して周期を求め
るという作業も不要となる。
【００２６】
【発明の効果】以上説明したように、本発明の圧力式波
高計は水圧データから波高データを復元する手段とし
て、波の周波数の関数として理論的に定まっている補正
倍数を周波数伝達関数とするディジタルフィルターを用
い水圧波形データをこのディジタルフィルターを通過さ
せるようにしたので、従来のゼロアップクロス法のよう
に、水圧波形データから１波１波分離して波の周期を割
り出すという作業は必要なくなり、周期が正確に読み取
りにくいという問題もなくなり、また、補正倍数が精々
２．５までで７秒以下の周期の波が計測できなかったの
に較べ、少なくとも補正倍数が３０位まで、周期３．５
秒位迄の短い波まで補正の対象とすることができ、従来
よりも一層正確な波高の復元ができるという利点があ
る。
高計は水圧データから波高データを復元する手段とし
て、波の周波数の関数として理論的に定まっている補正
倍数を周波数伝達関数とするディジタルフィルターを用
い水圧波形データをこのディジタルフィルターを通過さ
せるようにしたので、従来のゼロアップクロス法のよう
に、水圧波形データから１波１波分離して波の周期を割
り出すという作業は必要なくなり、周期が正確に読み取
りにくいという問題もなくなり、また、補正倍数が精々
２．５までで７秒以下の周期の波が計測できなかったの
に較べ、少なくとも補正倍数が３０位まで、周期３．５
秒位迄の短い波まで補正の対象とすることができ、従来
よりも一層正確な波高の復元ができるという利点があ
る。
【図１】本発明の圧力式波高計の第１の実施例のブロッ
ク図である。
ク図である。
【図２】本発明の圧力式波高計の第２の実施例のブロッ
ク図である。
ク図である。
【図３】本発明の圧力式波高計の第３の実施例のブロッ
ク図である。
ク図である。
【図４】本発明の圧力式波高計の第４の実施例のブロッ
ク図である。
ク図である。
【図５】本発明で用いたディジタルフィルターの周波数
特性を示す図である。
特性を示す図である。
【図６】本発明の実施例における、水圧波形図、復元さ
れた波高の波形図および正しい波高データとしての超音
波式波高計で計測された波高データの波形図の比較図で
ある。
れた波高の波形図および正しい波高データとしての超音
波式波高計で計測された波高データの波形図の比較図で
ある。
１ 水圧センサ ２ Ａ／Ｄ変換器 ３ ディジタルフィルター ４ 波高データメモリ ５ 水圧データメモリ ６ 水中通信ケーブル
Claims (1)
- 【請求項１】 次の各手段を具備することを特徴とする
圧力式波高計。 （イ）水面の波高の変化によって変化する水圧を検知し
アナログ信号で出力する水圧センサ （ロ）水圧センサからの水圧アナログ信号をディジタル
信号に変換するＡ／Ｄ変換器 （ハ）Ａ／Ｄ変換器からの水圧信号の、予め定めた周波
数範囲の周波数成分に対し、下記数式１の周波数伝達関
数（補正倍率）Ｈ（ｆ）を乗じて出力し、入力信号周波
数ｆに依存して生起する減衰誤差を補正した波高情報の
リアルタイム出力を確保するディジタルフィルター但
し、数式１のＬは、入力信号周波数をｆとしたとき次の
数式２から得られる値を代入する。 【数１】 【数２】
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP31302896A JP2911838B2 (ja) | 1996-11-08 | 1996-11-08 | 圧力式波高計 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP31302896A JP2911838B2 (ja) | 1996-11-08 | 1996-11-08 | 圧力式波高計 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH10141953A JPH10141953A (ja) | 1998-05-29 |
| JP2911838B2 true JP2911838B2 (ja) | 1999-06-23 |
Family
ID=18036354
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP31302896A Expired - Fee Related JP2911838B2 (ja) | 1996-11-08 | 1996-11-08 | 圧力式波高計 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2911838B2 (ja) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4725883B2 (ja) * | 2005-01-31 | 2011-07-13 | 社団法人海洋調査協会 | 水圧波から表面波への変換手法 |
| CN104457901B (zh) * | 2014-11-28 | 2018-01-05 | 南京信息工程大学 | 一种确定水深的方法及系统 |
| CN106969755B (zh) * | 2017-03-15 | 2019-08-13 | 天津大学 | 一种基于水压力反演海浪波高与周期的方法 |
| CN113155408B (zh) * | 2020-10-24 | 2023-04-18 | 大连理工大学 | 一种基于点压力反馈浪高显示的造波机实验装置及其方法 |
| CN115435757B (zh) * | 2022-11-09 | 2023-03-14 | 海南浙江大学研究院 | 一种利用水下动力压力反演表面波浪波高的方法 |
-
1996
- 1996-11-08 JP JP31302896A patent/JP2911838B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH10141953A (ja) | 1998-05-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110307994B (zh) | 异声检测装置以及异声检测方法 | |
| CA2277119C (en) | Seismic data acquisition and processing using non-linear distortion in a groundforce signal | |
| US5619616A (en) | Vehicle classification system using a passive audio input to a neural network | |
| WILSON | Seiches | |
| CN101877237B (zh) | 降噪装置和噪声确定方法 | |
| US7594423B2 (en) | Knock signal detection in automotive systems | |
| JP3704336B2 (ja) | 波形検出装置とそれを使用した状態監視システム | |
| WO2001033547A1 (en) | Methods and apparatuses for signal analysis | |
| JP2911838B2 (ja) | 圧力式波高計 | |
| GB2297465A (en) | Detecting periodicity or harmonic structure in a waveform | |
| Tang | On the time–frequency analysis of signals that decay exponentially with time | |
| JP4951060B2 (ja) | 渦流量計からの測定信号を処理するための方法 | |
| JP2908358B2 (ja) | 津波計 | |
| US4982433A (en) | Speech analysis method | |
| US5134590A (en) | Method and apparatus for correcting the response of seismic sensors to an emission signal that does not have the same shape as a reference signal | |
| CN106446320A (zh) | 基于imf瞬时能量曲率差值的井架钢结构损伤识别方法 | |
| EP0355994A1 (en) | Cross-correlation apparatus and methods | |
| JP3897398B2 (ja) | 湾水振動を考慮した波浪強制力の測定方法 | |
| Newland | Application of harmonic wavelets to time-frequency mapping | |
| Jeon et al. | Pseudo Wigner-Ville distribution, computer program and its applications to time-frequency domain problems | |
| Teisseire et al. | Measurement of the instantaneous characteristics of natural response modes of a spar platform subjected to irregular wave loading | |
| JP7000981B2 (ja) | 物体検知装置 | |
| JPH06273198A (ja) | 波浪特性抽出装置 | |
| JP3220205B2 (ja) | 雑音抑制装置 | |
| JPH0117599B2 (ja) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| R360 | Written notification for declining of transfer of rights |
Free format text: JAPANESE INTERMEDIATE CODE: R360 |
|
| R371 | Transfer withdrawn |
Free format text: JAPANESE INTERMEDIATE CODE: R371 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090409 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100409 Year of fee payment: 11 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110409 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110409 Year of fee payment: 12 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120409 Year of fee payment: 13 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120409 Year of fee payment: 13 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130409 Year of fee payment: 14 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130409 Year of fee payment: 14 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140409 Year of fee payment: 15 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |