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
Application number
JP60144353A
Other languages
English (en)
Inventor
Tadao Imai
今井 忠男
Yukinao Hashizume
橋爪 幸直
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Oki Electric Industry Co Ltd
Original Assignee
Oki Electric Industry Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Oki Electric Industry Co Ltd filed Critical Oki Electric Industry Co Ltd
Priority to JP60144353A priority Critical patent/JPS625477A/ja
Publication of JPS625477A publication Critical patent/JPS625477A/ja
Pending legal-status Critical Current

Links

Landscapes

  • Complex Calculations (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 (産業上の利用分野) 本発明はディジタル信号処理、特にディスクリート・フ
ーリエ変換方式に関するものである。
(従来の技術) 従来、このような分野の技術としては、宮用他著、電子
通信学会編「デジタル信号処理」(P18〜36.コロ
ナ社)に記載されたものがあった。前記文献によれば、
ディスクリート・フーリエ変換手段(以下、DFTと称
す。)においては、窓関数、窓時間及び減衰量に一連の
関係があった。
例えば、方形波窓(窓時間をTとする)では、主ローブ
以外は10数dBの減衰がとれ、極小点が1/T (H
2)に現われた。一方、ハミング窓(あるいはハニング
窓)を用いると、40数dB(あるいは30数dB)の
減衰がとれるが、極小点が2/T (H2)に現われ、
周波数分解能が悪くなり、極小点を同一にするには2T
の窓時間が必要となった。
(発明が解決しようどする問題点) 一般に、一定の周波数分解能で、減衰aを多くとるため
には窓時間を長くとる必要があるという問題点があり、
さらに、同−窓では窓時間を2倍にしても、それ以上の
減衰はほとんどとれず、必要な減衰量によって、方形波
窓(窓時間T)、ハミング又はパニング窓(2T> 、
ブラックマン窓(3T)を選択して使用してしなければ
ならないという問題点があった。
本発明の目的は、DPI−の演算時間が余り問題となら
ない場合、同じ窓時間でダイナミックレンジ(減衰量)
を大きくとることにある。
(問題点を解決するための手段) 本発明では前記問題点を解決するため、入力ディジタル
信号列にCOS核及びSIN核の係数列をそれぞれ掛け
合わせてCO5核乗算値及びSIN核乗算値を算出する
核乗算値算出手段と、該COS核乗算値及びSIN核乗
算値をそれぞれ窓面聞分積分する積分手段と、咳各積分
値を二乗し加算しその平方根をとることにより絶対値を
算出り゛る絶対@算出手段とからなるディスクリート・
フーリエ変換手段を複数段備え、2段目以降の入力ディ
ジタル信号列としては前段の絶対値出力を用い、2段目
以降のCOS核及びSIN核の係数列としては前段のC
O3核乗算値及びSIN核乗算値をそれぞれ用いた。
(作用) 前記構成によれば、2段目以降のDFTにおいては前段
の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は周知のルート(平方根)回路である
。
係数メモリ3及び4は検出しようとする周波数foのC
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を構成
する。
また、同様に乗算器9及び10は核乗算値算出手段を構
成し、加算器16及び17は積分手段を構成し、乗n器
11,12、加算器18及びルート回路20は絶対値算
出手段を構成し、これらは2段目のDFTを構成する。
次に前記回路の動作について詳細に動作を説明する。こ
こで、入力信号列を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毎
に出力される。
一方、前記演算結果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段目の出力である。
ここで、検出しようとする周波数f。が入力した場合の
DFT1段目の演算結果をY。、それ以外の周波数fが
入力した場合のDFT一段目の演算結果をY とし、こ
れらの比−をR1で表わす「 ものとする。
Y、/Yo=Rf         ・・−・(3)こ
の時、減衰量としては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。
=Rf         ・・・・・・(1)となる。
従って、DFT2段目の減衰比はDFT1段目の減衰比
の二乗(対数値でいえば2倍)になっていることがねか
るく但し、2段目の演算結果を得るまでには2丁の時間
を必要とりる。)。
以下同様にDFTを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図に示すように減衰量
が大幅に増えていることがわかる。
なお、前記実施例はハードウェアで構成した例を示した
が、マイクロプロセッサにおけるソフトウェアとして実
現することもできる。
(発明の効果) 以上説明したように本発明によれば、DFTを複数段備
え、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図

Claims (1)

  1. 【特許請求の範囲】 入力ディジタル信号列にCOS核及びSIN核の係数列
    をそれぞれ掛け合わせてCOS核乗算値及びSIN核乗
    算値を算出する核乗算値算出手段と、前記COS核乗算
    値及びSIN核乗算値をそれぞれ窓時間分積分する積分
    手段と、 前記各積分値を二乗し加算しその平方根をとることによ
    り絶対値を算出する絶対値算出手段とからなるディスク
    リート・フーリエ変換手段を複数段備え、 2段目以降の入力ディジタル信号列としては前段の絶対
    値出力を用い、 2段目以降のCOS核及びSIN核の係数列としては前
    段のCOS核乗算値及びSIN核乗算値をそれぞれ用い
    たことを特徴とする ディスクリート・フーリエ変換方式。
JP60144353A 1985-07-01 1985-07-01 デイスクリ−ト・フ−リエ変換方式 Pending JPS625477A (ja)

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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7422905B1 (ja) * 2022-08-02 2024-01-26 三菱電機株式会社 信号処理装置および信号処理方法

Cited By (2)

* Cited by examiner, † Cited by third party
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) ヒヨウホンソクドテイゲンカイロ