JP5766896B1 - ピーク周波数検出装置、方法およびプログラム - Google Patents
ピーク周波数検出装置、方法およびプログラム Download PDFInfo
- Publication number
- JP5766896B1 JP5766896B1 JP2015509250A JP2015509250A JP5766896B1 JP 5766896 B1 JP5766896 B1 JP 5766896B1 JP 2015509250 A JP2015509250 A JP 2015509250A JP 2015509250 A JP2015509250 A JP 2015509250A JP 5766896 B1 JP5766896 B1 JP 5766896B1
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- digital data
- peak frequency
- power
- digital
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 title description 19
- 238000005070 sampling Methods 0.000 claims abstract description 84
- 238000001228 spectrum Methods 0.000 claims abstract description 28
- 239000000284 extract Substances 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 abstract description 16
- 238000006243 chemical reaction Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 19
- 238000004364 calculation method Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 6
- 230000006866 deterioration Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000004590 computer program Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 230000010355 oscillation Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 2
- 239000013535 sea water Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/04—Measuring peak values or amplitude or envelope of ac or of pulses
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
- G01R23/165—Spectrum analysis; Fourier analysis using filters
- G01R23/167—Spectrum analysis; Fourier analysis using filters with digital filters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
- G01R23/177—Analysis of very low frequencies
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/50—Systems of measurement, based on relative movement of the target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/356—Receivers involving particularities of FFT processing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/523—Details of pulse systems
- G01S7/526—Receivers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10G—REPRESENTATION OF MUSIC; RECORDING MUSIC IN NOTATION FORM; ACCESSORIES FOR MUSIC OR MUSICAL INSTRUMENTS NOT OTHERWISE PROVIDED FOR, e.g. SUPPORTS
- G10G7/00—Other auxiliary devices or accessories, e.g. conductors' batons or separate holders for resin or strings
- G10G7/02—Tuning forks or like devices
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
- G10L25/90—Pitch determination of speech signals
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- Multimedia (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computational Linguistics (AREA)
- Human Computer Interaction (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Health & Medical Sciences (AREA)
- Signal Processing (AREA)
- Discrete Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
Description
n≧1/(ftg×Ttg)
fs/(n×ftg)≦N≦fs×Ttg
fs>2×n×fch
を満たす。なお本明細書において、fcl〜fchはfcl以上fch以下を意味するものとする。
・ピーク周波数の周波数分解能fn≦ftg
・FFTの時間窓長Tn≦Ttg
fs/(n×ftg)≦N≦fs×Ttg(Nは2のべき乗)
となる。
n=m1×m2×…×mk
fcl(j)≒(m1×m2×…×mj)×fcl
fch(j)≒(m1×m2×…×mj)×fch
であってもよい。
n≧1/(ftg×Ttg)
fs/(n×ftg)≦N≦fs×Ttg
fs>2×n×fch
1.第一実施形態
第一実施形態では、サンプリング周波数fsでサンプリングされた受信デジタルデータ列の、下限値fclと上限値fchとで定められた周波数帯域fcl〜fchにおけるピーク周波数fを、所望の周波数分解能ftgと所望の時間窓長Ttgを満たして検出するピーク周波数検出装置とこれを用いたドップラー計測器について説明する。
(ステップ1.乗数nを定める)
まず、n≧1/(ftg×Ttg)の条件を満たすn(nは、2以上の整数)を定める。例えば、nを、n≧1/(ftg×Ttg)を満たす最小の整数とし、後述するステップ3でNが存在しなければ、nを1増加させて再計算をしても良いし、最初から大きめにしておいても良い。
(ステップ2.FFTのサンプリング周波数fsを選択する)
fs>2×n×fch
となるように選択する。なお、本実施形態では、ピーク周波数検出装置1に入力されるデジタルデータ列のサンプリング周波数がFFTのサンプリング周波数fsとなる。fsが回路上の制約等で上式を満たせない場合は、他の実施形態を使用する。
(ステップ3.FFTのサンプル数Nを定める)
次に、fs/(n×ftg)≦N≦fs×Ttg(Nは2のべき乗)を満たすNを選択する。もし、Nが存在しない場合は、Nが存在するまでnを大きくするか、fsを大きくするか、もしくはその両方を行う。ただし、FFTのサンプリング定理から、
fs>2×n×fch
を満たしていなければならない。
乗数n、fs、Nが設定されたピーク周波数検出装置1に対象とするデジタルデータ列を入力すると、以下のように、所望の周波数分解能ftgと所望の時間窓長Ttgを満たしてピーク周波数が検出される。
送信周波数ftx:120kHz
局部発信周波数floc:137kHz
検出最大速度(水平方向)V:15m/s
検出速度精度(水平方向)V0:0.15m/s
検出対象物の位置の精度(位置分解能)D0:7.5m
Δfp=2×1200Hz=2400Hz
となる。
fmid±(△fp/2)=17000±1200Hz
の信号を通し、次段のA/D変換器28でエイリアシングが発生しないような特性とする。
n≧1/(ftg×Ttg)=1/(12×0.01)=8.3
となる。ここでは、n=12とする。
fs>2×n×fch=2×12×18200=436800Hz
となるため、fs=510kHzとする。
N=4096
となる。
図8は、本発明の第二実施形態としてのピーク周波数検出装置3の構成を示すブロック図である。ピーク周波数検出装置3は、ピーク周波数検出装置1の第一デジタルBPF11の前に、間引き部16が追加された構成となっている。これは、A/D変換後のデジタルデータ列を間引いて減らすことによりサンプリング周波数を下げるためである。ピーク周波数検出装置3に入力されるデジタルデータ列のサンプリング周波数fsが高く、乗数nと第一デジタルBPF11のカットオフ周波数fchがfs>4×n×fchを満たす場合、第二実施形態を適用することが好ましい。
fs>2×n×fch
となるようにする。fsは間引き後のサンプリング周波数である。
P(1):p0,p1,p2,…とし、間引いた後のデジタルデータ列を、Q(1):q0,q1,q2,…としたとき、qi=p(r×i)(i=0,1,2,3,4,・・・)を行う。r=2のときは、q0=p0、q1=p2、q2=p4、q3=p6となる。
間引く方法はこれ以外であっても良い。例えば、r=2のときに、q0=(p0+p1)/2、q1=(p2+p3)/2、q2=(p4+p5)/2、…などとしても良い。
fs/(n×ftg)≦N≦fs×Ttg (Nは2のべき乗)
を満たすNが存在せず、nを大きくすることもできない場合、本実施例は適用できない。この場合は第一実施形態を適用する。
fs=fis/r=10.2MHz/20=510kHz>436.8kHz
となる。
図9は、本発明の第三実施形態を示すピーク周波数検出装置4の構成を示すブロック図である。ピーク周波数検出装置4は、第一実施形態の第二デジタルBPF13の後に間引き部17が追加された構成となっている。第一デジタルBPF部11、n乗部12、第二デジタルBPF13は、入力されるデジタルデータのサンプリング周波数が、第一実施形態ではfsだったものが、fisとなっただけで、動作は同じである。受信デジタルデータ列のサンプリング周波数fsが高く、乗数nと第一デジタルBPF11のカットオフ周波数fchに対して、fs>4×n×fchの場合、計算量を削減するため、一般的には第二実施形態を適用するのが良い。しかし、本実施形態であっても、所望の時間窓長ftg、所望の時間窓長Ttgを満たした処理が可能である。
fs>2×n×fch
となるようにする。間引いた後のデジタルデータ列は、FFT部14に入力され、その後の処置は、第一実施形態と同一である。
fs=fis/r=10.2MHz/20=510kHz>436.8kHz
となる。FFTのサンプリング周波数fsは、fs=510kHzとなる。
図10は、本発明の第四実施形態としてのピーク周波数検出装置5の構成を示すブロック図である。ピーク周波数検出装置5は、第一実施形態のデジタルBPF11の前に、補間部18が追加された構成となっている。補間部18を追加することにより、A/D変換後のデジタルデータ列を補間して増やすことによりサンプリング周波数を上げることができる。ピーク周波数検出装置5に入力されるデジタルデータ列のサンプリング周波数がfsが低く、乗数nと第一デジタルBPF11の高いほうのカットオフ周波数fchに対して、2×fch<fs<2×n×fchの場合、本実施形態を適用する。
vi=u0 (i=0,1,2,3,・・・,(g−1))
vi=u1 (i=g,g+1,g+2,g+3,・・・,(2g−1))
vi=u2 (i=2g,2g+1,2g+2,2g+3,・・・,(3g−1))
・・・
となるように補間する。g=2のときは、v0=u0、v1=u0、v2=u1、v3=u1、v4=u2、v5=u2、…となる。補間する方法はこれ以外であっても良い。例えば、g=2のときに、
v0=u0
v1=(u0+u1)/2
v2=u1
v3=(u1+u2)/2
v4=u2
v5=(u2+u3)/2
・・・
となるように補間しても良い。
2×fch=2×18.2kHz=36.4kHz
2×n×fch=2×12×18.2kHz=436.8kHz
2×fch<fis<2×n×fch
であるから、本実施形態が適用される。
fs=fis×g=42.5kHz×12=510kHz>436.8kHz
となる。
v0〜v11=u0
v12〜v23=u1
v24〜v35=u2
v36〜v47=u3
・・・
となるように補間を行う。
図12Aは、A/D変換器28への入力信号で、
y=sin(2π17000t)+2
とする。つまり、アナログフィルタで取りきれなかった直流分2が残っている、振幅1、周波数17kHzのsin波形とする。この入力信号を、サンプリング周波数が17kHzの2.5倍である42.5kHzでA/D変換したものが、ピーク周波数検出装置1に入力されるデジタルデータ列に相当する。
図13は、本発明の第五実施形態としてのピーク周波数検出装置6の構成を示すブロック図である。ピーク周波数検出装置6は、第一実施形態のn乗部12と第二デジタルBPF13を多重べき乗部19に置き換えたものである。
n=m1×m2×・・・×mk
が成立するように選択される。
また、
fcl(j)≒(m1×m2×…×mj)×fcl
fch(j)≒(m1×m2×…×mj)×fch
と設定する。
とする。すなわち、bi= (ai)4 (i=0,1,2,3,4,…)となる。
ここまで説明したピーク周波数検出装置1、3−6に、パラメータ設定部を追加してもよい。図15は、第一実施形態のピーク周波数検出装置1にパラメータ設定部20を追加した構成を示すブロック図である。パラメータ設定部20は、プロセッサ、メモリ、入出力機構を備えるコンピュータであって、図示しないキーボードやマウスやタッチパネルディスプレイ等の図示しない操作部を用いたユーザの入力に応じてピーク周波数検出装置1、3−6にパラメータの値を設定する。パラメータは、これまで説明したとおり、次のようなものである。
受信デジタルデータ列のサンプリング周波数fis
デジタルデータ列のサンプリング周波数fs
所望の周波数分解能ftg
所望の時間窓長Ttg
n乗部の乗数n(nは2以上の整数)
第一デジタルBPFの周波数帯域 ほぼfcl〜fch(fcl<fch)
第二デジタルBPFの周波数帯域 ほぼn×fcl〜n×fch
FFTのサンプリング周波数fs
FFTのサンプリング数N
以上説明した本発明の実施形態によると、FFTによる周波数解析の際に制約となる周波数分解能(f0)と時間窓長(T0)との間にあるf0=1/T0と言う相反関係の問題を回避し、前記所望の周波数分解能ftgと前記所望の時間窓長Ttgでもって、信号波のピーク周波数の検出を可能にすることができる。受信デジタルデータ列のサンプリング周波数は、2×fch より大きければピーク周波数計算が可能となる。そして、近似計算やカーブフィッティング、平均化などの処理を必要としないため、ピーク周波数の計算精度の劣化は無い。
尚、本発明の技術的範囲は、上述した実施例に限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変更を加え得ることは勿論である。
Claims (11)
- 所定の周波数帯域(fcl〜fch)においてパワースペクトルが最大となるピーク周波数を検出するピーク周波数検出装置であって、
デジタルデータ列の各要素をn乗(nは2以上の整数)するn乗部と、
サンプリング周波数fsと周波数分解能ftgと時間窓長Ttgに応じて決まるN(Nは2のべき乗の整数)個のn乗されたサンプリング周波数fsのデジタルデータ列に対して高速フーリエ変換を行って得られるパワースペクトルの最大値に対応する周波数を仮想ピーク周波数として導出するFFT部と、
前記仮想ピーク周波数を1/n倍した値をデジタルデータ列のピーク周波数として出力する1/n倍部と、
を備え、
n≧1/(ftg×Ttg)
fs/(n×ftg)≦N≦fs×Ttg
fs>2×n×fch
を満たすピーク周波数検出装置。 - 前記所定の周波数帯域に含まれる周波数のデジタルデータ列を抽出する第一デジタルバンドパスフィルタを備え、
前記n乗部には、前記第一デジタルバンドパスフィルタの出力が入力される、
請求項1に記載のピーク周波数検出装置。 - サンプリング周波数fisのデジタルデータ列を1/r(rは2以上の整数)に間引いてサンプリング周波数をfsにする間引き部を備え、
前記第一デジタルバンドパスフィルタには前記間引き部の出力が入力される、
ことを特徴とする請求項2に記載のピーク周波数検出装置。 - 前記デジタルデータ列をg倍(gは2以上の整数)に補間してサンプリング周波数をfsにする補間部を備え、
前記第一デジタルバンドパスフィルタには前記補間部の出力が入力される、
ことを特徴とする請求項2記載のピーク周波数検出装置。 - n乗されたN個の前記デジタルデータ列から第二の周波数帯域に含まれるデジタルデータ列を抽出する第二デジタルバンドパスフィルタを備え、
前記FFT部には、前記第二デジタルバンドパスフィルタによって抽出されたデジタルデータ列が入力され、
前記第二の周波数帯域は、ほぼn×fcl〜n×fchである、
請求項1に記載のピーク周波数検出装置。 - 前記第二デジタルバンドパスフィルタで抽出されたデジタルデータ列を1/r(rは2以上の整数)に間引いてサンプリング周波数をfsとする間引き部を備え、
前記FFT部には、前記間引き部の出力が入力される、
請求項5に記載のピーク周波数検出装置。 - 前記特定の周波数帯域に含まれる周波数のデジタルデータ列を抽出する第一デジタルバンドパスフィルタと、
前記n乗部の出力から第二の周波数帯域に含まれるデジタルデータ列を抽出する第二デジタルバンドパスフィルタとを備え、
前記n乗部には、前記第一デジタルバンドパスフィルタの出力が入力され、
前記FFT部には、前記第二デジタルバンドパスフィルタの出力が入力され、
前記第二の周波数帯域は、ほぼn×fcl〜n×fchである、
請求項1記載のピーク周波数検出装置。 - 前記n乗部に代えて、入力されるデジタルデータ列をmj乗(mjは2以上の整数)するべき乗部(j)(j=1,2,…,k)と前記べき乗部(j)の出力から特定の周波数帯域fcl(j)〜fch(j)の信号を抽出するデジタルバンドパスフィルタ(j)とを備えるべき乗ブロック(j)をk(kは2以上の整数)段備える多重べき乗部を備え、
n=m1×m2×…×mk
fcl(j)≒(m1×m2×…×mj)×fcl
fch(j)≒(m1×m2×…×mj)×fch
である請求項1〜4のいずれか一項に記載のピーク周波数検出装置。 - ユーザの指示を受け付ける操作部と、
前記指示に応じたn、fs、Nの少なくともいずれかを設定するパラメータ設定部と、
を備える請求項1から8のいずれか一項に記載のピーク周波数検出装置。 - 所定の周波数帯域(fcl〜fch)においてパワースペクトルが最大となるピーク周波数を検出するピーク周波数検出方法であって、
デジタルデータ列の各要素をn乗(nは2以上の整数)し、
サンプリング周波数fsと周波数分解能ftgと時間窓長Ttgに応じて決まるN(Nは2のべき乗の整数)個のn乗されたサンプリング周波数fsのデジタルデータ列に対して高速フーリエ変換を行って得られるパワースペクトルの最大値に対応する周波数を仮想ピーク周波数として導出し、
前記仮想ピーク周波数を1/n倍した値をデジタルデータ列のピーク周波数として出力する、
ことを含み、
n≧1/(ftg×Ttg)
fs/(n×ftg)≦N≦fs×Ttg
fs>2×n×fch
を満たすピーク周波数検出方法。 - 所定の周波数帯域(fcl〜fch)においてパワースペクトルが最大となるピーク周波数を検出するピーク周波数検出プログラムであって、
デジタルデータ列の各要素をn乗(nは2以上の整数)するn乗部と、
サンプリング周波数fsと周波数分解能ftgと時間窓長Ttgに応じて決まるN(Nは2のべき乗の整数)個のn乗されたサンプリング周波数fsのデジタルデータ列に対して高速フーリエ変換を行って得られるパワースペクトルの最大値に対応する周波数を仮想ピーク周波数として導出するFFT部と、
前記仮想ピーク周波数を1/n倍した値をデジタルデータ列のピーク周波数として出力する1/n倍部と、
してコンピュータを機能させ、
n≧1/(ftg×Ttg)
fs/(n×ftg)≦N≦fs×Ttg
fs>2×n×fch
を満たすピーク周波数検出プログラム。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2014/068432 WO2016006079A1 (ja) | 2014-07-10 | 2014-07-10 | ピーク周波数検出装置、方法およびプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP5766896B1 true JP5766896B1 (ja) | 2015-08-19 |
JPWO2016006079A1 JPWO2016006079A1 (ja) | 2017-04-27 |
Family
ID=53888036
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2015509250A Active JP5766896B1 (ja) | 2014-07-10 | 2014-07-10 | ピーク周波数検出装置、方法およびプログラム |
Country Status (6)
Country | Link |
---|---|
US (1) | US9857399B2 (ja) |
EP (1) | EP3168638A4 (ja) |
JP (1) | JP5766896B1 (ja) |
KR (1) | KR101785714B1 (ja) |
CN (1) | CN106415310B (ja) |
WO (1) | WO2016006079A1 (ja) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107402326B (zh) * | 2017-07-20 | 2019-08-23 | 南京理工大学 | 一种改进s变换的有限窗长时频分析方法 |
CN107589424B (zh) * | 2017-08-31 | 2021-12-17 | 努比亚技术有限公司 | 超声波采样方法、装置及计算机可读存储介质 |
CN114371342B (zh) * | 2022-03-21 | 2022-05-27 | 国仪量子(合肥)技术有限公司 | Fpga及基于其的实时信号测频方法以及锁相放大器 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01110351A (ja) * | 1987-10-23 | 1989-04-27 | Aloka Co Ltd | 超音波ドプラ診断装置 |
JP2808954B2 (ja) * | 1991-11-13 | 1998-10-08 | 国際電信電話株式会社 | 無変調信号検出及び周波数引き込み装置 |
US5431169A (en) * | 1992-08-03 | 1995-07-11 | Olympus Optical Co., Ltd. | Ultrasonic diagnosing apparatus |
KR970003966B1 (ko) * | 1994-11-22 | 1997-03-24 | 삼성전자 주식회사 | 윈도우필터를 이용한 직접확산통신시스템의 수신기 |
JP3271504B2 (ja) * | 1996-02-02 | 2002-04-02 | 三菱電機株式会社 | 周波数推定回路およびそれを用いたafc回路 |
US5799038A (en) * | 1996-04-30 | 1998-08-25 | Advantest Corporation | Method for measuring modulation parameters of digital quadrature-modulated signal |
JPH09298572A (ja) * | 1996-04-30 | 1997-11-18 | Advantest Corp | Psk信号の搬送波周波数推定器 |
US7239674B2 (en) * | 2003-06-04 | 2007-07-03 | Honeywell Federal Manufacturing & Technologies, Llc | Method of differential-phase/absolute-amplitude QAM |
KR20070056818A (ko) * | 2005-11-30 | 2007-06-04 | 주식회사 유컴테크놀러지 | 무선주파수인식 시스템 |
US20080052335A1 (en) * | 2006-08-01 | 2008-02-28 | Gee Edward C | Systems and methods for time domain to frequency domain conversion using frequency shifting |
JP4352082B2 (ja) * | 2007-06-18 | 2009-10-28 | 株式会社東芝 | 周波数同期回路、方法、プログラム及びこれらを用いた受信装置 |
JP5089460B2 (ja) * | 2008-01-16 | 2012-12-05 | 三菱電機株式会社 | 伝搬遅延時間測定装置及びレーダ装置 |
US8390508B1 (en) | 2010-04-05 | 2013-03-05 | Raytheon Company | Generating radar cross-section signatures |
JP5094937B2 (ja) | 2010-09-21 | 2012-12-12 | 三菱電機株式会社 | 周波数解析装置 |
JP2012247304A (ja) * | 2011-05-27 | 2012-12-13 | Sonic Corp | 短時間信号のピークパワースペクトルを検出する方法及び装置 |
JP5840868B2 (ja) * | 2011-05-27 | 2016-01-06 | 株式会社ソニック | 周波数検出方法及び装置 |
JP5670836B2 (ja) | 2011-05-27 | 2015-02-18 | 株式会社ソニック | フーリエ変換でのサンプル数を削減した、短時間信号のピークパワースペクトルを検出する方法及び装置 |
-
2014
- 2014-07-10 EP EP14897177.3A patent/EP3168638A4/en not_active Withdrawn
- 2014-07-10 US US15/314,614 patent/US9857399B2/en not_active Expired - Fee Related
- 2014-07-10 KR KR1020167033030A patent/KR101785714B1/ko active IP Right Grant
- 2014-07-10 JP JP2015509250A patent/JP5766896B1/ja active Active
- 2014-07-10 WO PCT/JP2014/068432 patent/WO2016006079A1/ja active Application Filing
- 2014-07-10 CN CN201480079481.8A patent/CN106415310B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
JPWO2016006079A1 (ja) | 2017-04-27 |
EP3168638A4 (en) | 2018-06-13 |
US20170205448A1 (en) | 2017-07-20 |
KR101785714B1 (ko) | 2017-11-06 |
CN106415310A (zh) | 2017-02-15 |
EP3168638A1 (en) | 2017-05-17 |
CN106415310B (zh) | 2018-02-02 |
WO2016006079A1 (ja) | 2016-01-14 |
US9857399B2 (en) | 2018-01-02 |
KR20170009878A (ko) | 2017-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2178082B1 (en) | Cyclic signal processing method, cyclic signal conversion method, cyclic signal processing device, and cyclic signal analysis method | |
JP5670836B2 (ja) | フーリエ変換でのサンプル数を削減した、短時間信号のピークパワースペクトルを検出する方法及び装置 | |
AU2018208683B2 (en) | Flow meter configuration and calibration | |
CN107210046B (zh) | 用于处理和分析信号的方法,以及实现这种方法的装置 | |
RU2011134419A (ru) | Блок модулированных фильтров с малым запаздыванием | |
JP2013541024A (ja) | 信号を解析し、瞬時周波数および短時間フーリエ変換を提供するための方法、ならびに信号を解析するためのデバイス | |
JP5766896B1 (ja) | ピーク周波数検出装置、方法およびプログラム | |
JP5840868B2 (ja) | 周波数検出方法及び装置 | |
JP4734057B2 (ja) | 移動物体検出装置 | |
JP5492606B2 (ja) | 演算装置、及び演算装置を備えた流量計 | |
JP2012247304A (ja) | 短時間信号のピークパワースペクトルを検出する方法及び装置 | |
Huang et al. | Resolution doubled co-prime spectral analyzers for removing spurious peaks | |
JP2000181472A (ja) | 信号分析装置 | |
JP6241131B2 (ja) | 音響用フィルタ装置、音響用フィルタリング方法、およびプログラム | |
JP2002303645A (ja) | 周波数計測装置、周波数計測方法およびレーダ装置 | |
Jarrot et al. | A time-frequency characterization framework for signals issued from underwater dispersive environments | |
White et al. | Signal processing techniques | |
JP7056739B2 (ja) | 波源方向推定装置、波源方向推定方法、およびプログラム | |
WO2021193637A1 (ja) | 基本周波数推定装置、アクティブノイズコントロール装置、基本周波数の推定方法及び基本周波数の推定プログラム | |
KR100987306B1 (ko) | 신호 처리 장치 및 그 방법 | |
JP3964095B2 (ja) | 大気温度測定方法および装置 | |
JPH11160422A (ja) | 電波高度計 | |
JPH07231244A (ja) | 狭帯域ディジタルフィルタ及びスペクトル分析器及びドップラ周波数推定器及び地殻同定器及び地磁気同定器及び地殻推定器及び地磁気推定器 | |
JP5487062B2 (ja) | 雑音除去装置 | |
Nelson | High-resolution correlation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20150519 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20150617 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5766896 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |