JPS625477A - デイスクリ−ト・フ−リエ変換方式 - Google Patents
デイスクリ−ト・フ−リエ変換方式Info
- Publication number
- JPS625477A JPS625477A JP60144353A JP14435385A JPS625477A JP S625477 A JPS625477 A JP S625477A JP 60144353 A JP60144353 A JP 60144353A JP 14435385 A JP14435385 A JP 14435385A JP S625477 A JPS625477 A JP S625477A
- Authority
- JP
- Japan
- Prior art keywords
- kernel
- multiplication value
- sin
- cos
- stage
- 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.)
- Pending
Links
- 108010076504 Protein Sorting Signals Proteins 0.000 claims abstract description 3
- 238000000034 method Methods 0.000 claims description 4
- 238000003775 Density Functional Theory Methods 0.000 abstract 3
- 230000006870 function Effects 0.000 description 5
- 230000015654 memory Effects 0.000 description 5
- 238000013016 damping Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004091 panning Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Landscapes
- Complex Calculations (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
(産業上の利用分野)
本発明はディジタル信号処理、特にディスクリート・フ
ーリエ変換方式に関するものである。
ーリエ変換方式に関するものである。
(従来の技術)
従来、このような分野の技術としては、宮用他著、電子
通信学会編「デジタル信号処理」(P18〜36.コロ
ナ社)に記載されたものがあった。前記文献によれば、
ディスクリート・フーリエ変換手段(以下、DFTと称
す。)においては、窓関数、窓時間及び減衰量に一連の
関係があった。
通信学会編「デジタル信号処理」(P18〜36.コロ
ナ社)に記載されたものがあった。前記文献によれば、
ディスクリート・フーリエ変換手段(以下、DFTと称
す。)においては、窓関数、窓時間及び減衰量に一連の
関係があった。
例えば、方形波窓(窓時間をTとする)では、主ローブ
以外は10数dBの減衰がとれ、極小点が1/T (H
2)に現われた。一方、ハミング窓(あるいはハニング
窓)を用いると、40数dB(あるいは30数dB)の
減衰がとれるが、極小点が2/T (H2)に現われ、
周波数分解能が悪くなり、極小点を同一にするには2T
の窓時間が必要となった。
以外は10数dBの減衰がとれ、極小点が1/T (H
2)に現われた。一方、ハミング窓(あるいはハニング
窓)を用いると、40数dB(あるいは30数dB)の
減衰がとれるが、極小点が2/T (H2)に現われ、
周波数分解能が悪くなり、極小点を同一にするには2T
の窓時間が必要となった。
(発明が解決しようどする問題点)
一般に、一定の周波数分解能で、減衰aを多くとるため
には窓時間を長くとる必要があるという問題点があり、
さらに、同−窓では窓時間を2倍にしても、それ以上の
減衰はほとんどとれず、必要な減衰量によって、方形波
窓(窓時間T)、ハミング又はパニング窓(2T> 、
ブラックマン窓(3T)を選択して使用してしなければ
ならないという問題点があった。
には窓時間を長くとる必要があるという問題点があり、
さらに、同−窓では窓時間を2倍にしても、それ以上の
減衰はほとんどとれず、必要な減衰量によって、方形波
窓(窓時間T)、ハミング又はパニング窓(2T> 、
ブラックマン窓(3T)を選択して使用してしなければ
ならないという問題点があった。
本発明の目的は、DPI−の演算時間が余り問題となら
ない場合、同じ窓時間でダイナミックレンジ(減衰量)
を大きくとることにある。
ない場合、同じ窓時間でダイナミックレンジ(減衰量)
を大きくとることにある。
(問題点を解決するための手段)
本発明では前記問題点を解決するため、入力ディジタル
信号列にCOS核及びSIN核の係数列をそれぞれ掛け
合わせてCO5核乗算値及びSIN核乗算値を算出する
核乗算値算出手段と、該COS核乗算値及びSIN核乗
算値をそれぞれ窓面聞分積分する積分手段と、咳各積分
値を二乗し加算しその平方根をとることにより絶対値を
算出り゛る絶対@算出手段とからなるディスクリート・
フーリエ変換手段を複数段備え、2段目以降の入力ディ
ジタル信号列としては前段の絶対値出力を用い、2段目
以降のCOS核及びSIN核の係数列としては前段のC
O3核乗算値及びSIN核乗算値をそれぞれ用いた。
信号列にCOS核及びSIN核の係数列をそれぞれ掛け
合わせてCO5核乗算値及びSIN核乗算値を算出する
核乗算値算出手段と、該COS核乗算値及びSIN核乗
算値をそれぞれ窓面聞分積分する積分手段と、咳各積分
値を二乗し加算しその平方根をとることにより絶対値を
算出り゛る絶対@算出手段とからなるディスクリート・
フーリエ変換手段を複数段備え、2段目以降の入力ディ
ジタル信号列としては前段の絶対値出力を用い、2段目
以降のCOS核及びSIN核の係数列としては前段のC
O3核乗算値及びSIN核乗算値をそれぞれ用いた。
(作用)
前記構成によれば、2段目以降のDFTにおいては前段
のDFTの演篩結果に前段の核乗輝値を掛け合わせてい
るため、1段のDFTによる減衰比が段数分掛け合わさ
れることになり、対数値でいえば、(DFT1段の減衰
量)×(段数)[dB]の減衰量が得られることになり
、より大きなダイナミックレンジが得られる。
のDFTの演篩結果に前段の核乗輝値を掛け合わせてい
るため、1段のDFTによる減衰比が段数分掛け合わさ
れることになり、対数値でいえば、(DFT1段の減衰
量)×(段数)[dB]の減衰量が得られることになり
、より大きなダイナミックレンジが得られる。
(実施例)
第1図は本発明の一実施例を示すもので、ここではDF
Tを2段構成とした例を示す。同図において、1は入力
端子、2は出力端子、3,4は読み出し専用の係数メモ
リ、5,6,7,8.9゜10.11.12は周知の乗
算器、13,14゜15.16,17.18は周知の加
算器、19゜20は周知のルート(平方根)回路である
。
Tを2段構成とした例を示す。同図において、1は入力
端子、2は出力端子、3,4は読み出し専用の係数メモ
リ、5,6,7,8.9゜10.11.12は周知の乗
算器、13,14゜15.16,17.18は周知の加
算器、19゜20は周知のルート(平方根)回路である
。
係数メモリ3及び4は検出しようとする周波数foのC
OS核及びSIN核の係数、即ちω(t)cos(2π
f t)及びω(t)sin(2πfot)(ω(t
):窓関数)の係、数を、それぞれ窓面@Tをサンプリ
ング時間S で割った数N (−T/S、)だけ記憶℃ している。
OS核及びSIN核の係数、即ちω(t)cos(2π
f t)及びω(t)sin(2πfot)(ω(t
):窓関数)の係、数を、それぞれ窓面@Tをサンプリ
ング時間S で割った数N (−T/S、)だけ記憶℃ している。
乗算器5及び6は、それぞれ入力ディジタル信号列にC
OS核及びSIN核の係数列をそれぞれ掛け合わせてC
OS核乗口値及びSIN核乗算値を算出する核乗算値算
出手段を構成する。加算器13及び14は、それぞれC
OS核乗算値及びSIN核乗算値をそれぞれ窓面聞分積
分する積分手段を構成する。乗算器°?、8、加算器1
5及びルート回路19は各積分値を二乗し加算し、その
平方根をとることにより絶対値を算出する絶対値算出手
段を構成する。而して、これらは1段目のDFTを構成
する。
OS核及びSIN核の係数列をそれぞれ掛け合わせてC
OS核乗口値及びSIN核乗算値を算出する核乗算値算
出手段を構成する。加算器13及び14は、それぞれC
OS核乗算値及びSIN核乗算値をそれぞれ窓面聞分積
分する積分手段を構成する。乗算器°?、8、加算器1
5及びルート回路19は各積分値を二乗し加算し、その
平方根をとることにより絶対値を算出する絶対値算出手
段を構成する。而して、これらは1段目のDFTを構成
する。
また、同様に乗算器9及び10は核乗算値算出手段を構
成し、加算器16及び17は積分手段を構成し、乗n器
11,12、加算器18及びルート回路20は絶対値算
出手段を構成し、これらは2段目のDFTを構成する。
成し、加算器16及び17は積分手段を構成し、乗n器
11,12、加算器18及びルート回路20は絶対値算
出手段を構成し、これらは2段目のDFTを構成する。
次に前記回路の動作について詳細に動作を説明する。こ
こで、入力信号列をXi、係数メモリ5内の係数列を0
1、係数メモリ6内の係数列を3iとする。
こで、入力信号列をXi、係数メモリ5内の係数列を0
1、係数メモリ6内の係数列を3iとする。
入力端子1より入力されたディジタル信号列xiは、乗
算器5及び6において、係数列Ci及びSiと順次掛け
合わされ、COS核乗算値及びSIN 核乗算値(X
t c i及ヒx 1s t >が算出される。これら
の乗算値はそれぞれ加算器13及び14によって窓時間
Tの間、即ちN個積分される(ΣX1C1及びΣxis
’>。これらの積L”I
L=1分値は乗算器7及び8により二乗され、加算器
15で加算され、さらにルート回路19によってその平
方根がとられ、絶対値が算出される。従って、ルート回
路19の出力をYとすると、+(谷X H8+ )2
) 1/2 ・・・・・・(1)し++1 となる。これはDFT1段目の出力であり、窓時間T毎
に出力される。
算器5及び6において、係数列Ci及びSiと順次掛け
合わされ、COS核乗算値及びSIN 核乗算値(X
t c i及ヒx 1s t >が算出される。これら
の乗算値はそれぞれ加算器13及び14によって窓時間
Tの間、即ちN個積分される(ΣX1C1及びΣxis
’>。これらの積L”I
L=1分値は乗算器7及び8により二乗され、加算器
15で加算され、さらにルート回路19によってその平
方根がとられ、絶対値が算出される。従って、ルート回
路19の出力をYとすると、+(谷X H8+ )2
) 1/2 ・・・・・・(1)し++1 となる。これはDFT1段目の出力であり、窓時間T毎
に出力される。
一方、前記演算結果Yは乗算器9及び10に送出され、
乗算器5及び6の出力、即ち1段目の核乗算値(X t
c +及びX;S;>と順次、掛け合わされる(YX
iCi、及(FYXiSi)、前記乗算値はそれぞれ加
算器16及び17によって窓時間Tの間、即ちN個積分
される(ΣYXiL嶌1 Ci及びΣYXiSi)。この積分値は乗算器z1 11及び12により二乗され、加算器18で加算され、
さらにルート回路20によってその平方根がとられて絶
対値にされる。従って、ルート回路20の出力をZとす
ると、 となる。これがDFT2段目の出力である。
乗算器5及び6の出力、即ち1段目の核乗算値(X t
c +及びX;S;>と順次、掛け合わされる(YX
iCi、及(FYXiSi)、前記乗算値はそれぞれ加
算器16及び17によって窓時間Tの間、即ちN個積分
される(ΣYXiL嶌1 Ci及びΣYXiSi)。この積分値は乗算器z1 11及び12により二乗され、加算器18で加算され、
さらにルート回路20によってその平方根がとられて絶
対値にされる。従って、ルート回路20の出力をZとす
ると、 となる。これがDFT2段目の出力である。
ここで、検出しようとする周波数f。が入力した場合の
DFT1段目の演算結果をY。、それ以外の周波数fが
入力した場合のDFT一段目の演算結果をY とし、こ
れらの比−をR1で表わす「 ものとする。
DFT1段目の演算結果をY。、それ以外の周波数fが
入力した場合のDFT一段目の演算結果をY とし、こ
れらの比−をR1で表わす「 ものとする。
Y、/Yo=Rf ・・−・(3)こ
の時、減衰量としては201o(1(Ya / Y 1
)[dB]となる。
の時、減衰量としては201o(1(Ya / Y 1
)[dB]となる。
一方、周波数f。が入力した場合のDFT2段目の演算
結果を7゜、それ以外の周波数fが入力した場合のDF
T2段目の演算結果を7.とすると、前記(2)式より
、 Zo=Yo ・・・・・・(4
)Z、=Y、 ・・・・・
・(5)となる。前記(3)式よりY、−Rf−Yoで
あるから、(5)式は、 Z、 = (R、−Yo) 2・・−・−・(6)と表
わされる。従って、(6)式を(4)式で割ってDF7
2段目の出力におけるZ、と7゜どの比をとると、 Z、/Zo=(R,−Yo) /Y。
結果を7゜、それ以外の周波数fが入力した場合のDF
T2段目の演算結果を7.とすると、前記(2)式より
、 Zo=Yo ・・・・・・(4
)Z、=Y、 ・・・・・
・(5)となる。前記(3)式よりY、−Rf−Yoで
あるから、(5)式は、 Z、 = (R、−Yo) 2・・−・−・(6)と表
わされる。従って、(6)式を(4)式で割ってDF7
2段目の出力におけるZ、と7゜どの比をとると、 Z、/Zo=(R,−Yo) /Y。
=Rf ・・・・・・(1)となる。
従って、DFT2段目の減衰比はDFT1段目の減衰比
の二乗(対数値でいえば2倍)になっていることがねか
るく但し、2段目の演算結果を得るまでには2丁の時間
を必要とりる。)。
の二乗(対数値でいえば2倍)になっていることがねか
るく但し、2段目の演算結果を得るまでには2丁の時間
を必要とりる。)。
以下同様にDFTを3段、4段・・・・・・n段と増や
せば、その減衰比は3乗、4乗、・・・・・・0乗(対
数値でいえば3倍、4倍・・・・・・n倍)となる。
せば、その減衰比は3乗、4乗、・・・・・・0乗(対
数値でいえば3倍、4倍・・・・・・n倍)となる。
第2図は前記実施例により得られた減衰量の周波数特性
を示す。ここで、窓関数はハミング窓(ω(t) −0
,54−0,46cos 2πt/T>であり、図中、
実線はDF71段目の出力、一点鎖線はDF72段目の
出力を示す。また、比較のため、同一の窓関数で窓時間
を単に2倍にした場合の減衰量の周波数特性を第3図に
示す。第3図にa3いて、実線は窓時間が王の場合(前
記実施例の1段目の出力と同一)を、また一点鎖線は窓
時間が2王の場合をそれぞれ示しており、この場合はほ
とんど減衰mが増えていない、即ち特性が改善されてい
ないが、前記実施例によれば第2図に示すように減衰量
が大幅に増えていることがわかる。
を示す。ここで、窓関数はハミング窓(ω(t) −0
,54−0,46cos 2πt/T>であり、図中、
実線はDF71段目の出力、一点鎖線はDF72段目の
出力を示す。また、比較のため、同一の窓関数で窓時間
を単に2倍にした場合の減衰量の周波数特性を第3図に
示す。第3図にa3いて、実線は窓時間が王の場合(前
記実施例の1段目の出力と同一)を、また一点鎖線は窓
時間が2王の場合をそれぞれ示しており、この場合はほ
とんど減衰mが増えていない、即ち特性が改善されてい
ないが、前記実施例によれば第2図に示すように減衰量
が大幅に増えていることがわかる。
なお、前記実施例はハードウェアで構成した例を示した
が、マイクロプロセッサにおけるソフトウェアとして実
現することもできる。
が、マイクロプロセッサにおけるソフトウェアとして実
現することもできる。
(発明の効果)
以上説明したように本発明によれば、DFTを複数段備
え、2段目以降のDFTにおいては前段のDFTの演算
結果に前段の核乗算値を掛け合わせているため、1段の
DFTによる減衰比が段数分掛け合わされることになり
、対数値でいえば、(DFTI段の減衰1)×(段数)
の減衰量が得られることになり、より大きなダイナミッ
クレンジを得ることができる等の利点がある。
え、2段目以降のDFTにおいては前段のDFTの演算
結果に前段の核乗算値を掛け合わせているため、1段の
DFTによる減衰比が段数分掛け合わされることになり
、対数値でいえば、(DFTI段の減衰1)×(段数)
の減衰量が得られることになり、より大きなダイナミッ
クレンジを得ることができる等の利点がある。
図面は本発明の説明に供するもので、第1図は本発明の
ディスクリート・フーリエ変換方式を実現する回路の一
実施例を示すブロック図、第2図は第1図の回路におけ
る減衰量の周波数特性図、第3図は従来のディスクリー
ト・フーリエ変換方式による減衰量の周波数特性図であ
る。 3.4・・・・・・係数メモリ、5,6.7.8.9゜
10.11.12・・・・・・乗算器、13,14゜1
5.16.17.18・・・・・・加算器、19゜20
・・・・・・ルート回路。 特許出願人 沖電気工業株式会社 代理人弁理士 古 1) 精 孝 [dB] 周瀘奴 9乏東方式によ5滅奄量の■し皮収シi生回第3図
ディスクリート・フーリエ変換方式を実現する回路の一
実施例を示すブロック図、第2図は第1図の回路におけ
る減衰量の周波数特性図、第3図は従来のディスクリー
ト・フーリエ変換方式による減衰量の周波数特性図であ
る。 3.4・・・・・・係数メモリ、5,6.7.8.9゜
10.11.12・・・・・・乗算器、13,14゜1
5.16.17.18・・・・・・加算器、19゜20
・・・・・・ルート回路。 特許出願人 沖電気工業株式会社 代理人弁理士 古 1) 精 孝 [dB] 周瀘奴 9乏東方式によ5滅奄量の■し皮収シi生回第3図
Claims (1)
- 【特許請求の範囲】 入力ディジタル信号列にCOS核及びSIN核の係数列
をそれぞれ掛け合わせてCOS核乗算値及びSIN核乗
算値を算出する核乗算値算出手段と、前記COS核乗算
値及びSIN核乗算値をそれぞれ窓時間分積分する積分
手段と、 前記各積分値を二乗し加算しその平方根をとることによ
り絶対値を算出する絶対値算出手段とからなるディスク
リート・フーリエ変換手段を複数段備え、 2段目以降の入力ディジタル信号列としては前段の絶対
値出力を用い、 2段目以降のCOS核及びSIN核の係数列としては前
段のCOS核乗算値及びSIN核乗算値をそれぞれ用い
たことを特徴とする ディスクリート・フーリエ変換方式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60144353A JPS625477A (ja) | 1985-07-01 | 1985-07-01 | デイスクリ−ト・フ−リエ変換方式 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60144353A JPS625477A (ja) | 1985-07-01 | 1985-07-01 | デイスクリ−ト・フ−リエ変換方式 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPS625477A true JPS625477A (ja) | 1987-01-12 |
Family
ID=15360130
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP60144353A Pending JPS625477A (ja) | 1985-07-01 | 1985-07-01 | デイスクリ−ト・フ−リエ変換方式 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS625477A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7422905B1 (ja) * | 2022-08-02 | 2024-01-26 | 三菱電機株式会社 | 信号処理装置および信号処理方法 |
-
1985
- 1985-07-01 JP JP60144353A patent/JPS625477A/ja active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7422905B1 (ja) * | 2022-08-02 | 2024-01-26 | 三菱電機株式会社 | 信号処理装置および信号処理方法 |
| WO2024028966A1 (ja) * | 2022-08-02 | 2024-02-08 | 三菱電機株式会社 | 信号処理装置および信号処理方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US3949206A (en) | Filtering device | |
| AU689439B2 (en) | Digital filter having high accuracy and efficiency | |
| US3971927A (en) | Modular discrete cosine transform system | |
| JPS6196817A (ja) | フイルタ− | |
| JP2957183B2 (ja) | 巡回型ディジタルフィルタ | |
| US4062060A (en) | Digital filter | |
| US4012628A (en) | Filter with a reduced number of shift register taps | |
| SE444730B (sv) | Ljudsyntetisator | |
| JP2529229B2 (ja) | コサイン変換装置 | |
| JP2883494B2 (ja) | デジタルフィルター | |
| US6324222B1 (en) | Digital receiver with polyphase structure | |
| US7290022B2 (en) | Method and filter arrangement for digital recursive filtering in the time domain | |
| Anastassopoulos et al. | Method for designing half-band delta-modulation FIR filters | |
| KR950002072B1 (ko) | 디지틀 대역 분할 주파수 변환에서의 필터 계수 배열 방법 | |
| JPS5990419A (ja) | 2次のデイジタル全域通過回路 | |
| JPH01257274A (ja) | 周波数推定方式 | |
| JPH0226408A (ja) | ディジタルフィルタ | |
| JPS6032364B2 (ja) | デイジタル・フイルタ | |
| JPH0641389Y2 (ja) | デジタルフイルタ | |
| JPS592350B2 (ja) | デイジタルレベル検出装置 | |
| JPH0828646B2 (ja) | ディジタルフィルタ | |
| JPS6318368B2 (ja) | ||
| JPH02272818A (ja) | ディジタル演算回路 | |
| JPH05108693A (ja) | Iirデジタルフイルタ装置 | |
| JPS5842657B2 (ja) | ヒヨウホンソクドテイゲンカイロ |