!!!
!!
!!
!!
!!
!
!
"""
""
""
""
""
"
"
###
##
##
##
##
#
#
$$$
$$
$$
$$
$$
$
$
%%%
%%
%%
%%
%%
%
%
! " # $ % &
ＥＬＥＣＴＲＯＮＩＣ ＭＥＡＳＵＲＥＭＥＮＴ ＴＥＣＨＮＯＬＯＧＹ
'３２( '２)
２００９*２+
x©8W犉狅狌狉犻犲狉+W8ªe
Nz{ $|, h}D
（
Y!L'®)%&Bºâ'(
，
Y!6)7ºâ23æ8C9: Y! ２０００７２）
1 2：PZＦｏｕｒｉｅｒ*S÷5Ä?]0²Hæ-g~R，×jH88PZＦｏｕｒｉｅｒ*（ＤＦＲＦＴ）3
H8©+Ä
。
:@È..%Ç»~ZPÅ8ＤＦＲＦＴÄ，"9ZPＤＦＲＦＴ8¨>?WfßPZÚ
8ＤＦＲＦＴ，>FP"ＤＦＴＨｅｒｍｉｔｅ:¦¤?¨
A>1&88ＦｏｕｒｉｅｒP*（ＩＤＦＴ），YSäU<P
ZÚ8!0²æ5Ä%ØＩＤＦＴ，mu5Ä<7H*$，WH[\ßÄ$，Tg~SPZÚ8
$âã¡#?qrñIJ
，
ÝÞ@é¥.TÄ8Hh
。
345
：
PZＦｏｕｒｉｅｒ*；8PZＦｏｕｒｉｅｒ*；ZPÅ
6789:
：ＴＮ９１１．７ ;<=>?：Ａ
犇犻狊犮狉犲狋犲犳狉犪犮狋犻狅狀犪犾犉狅狌狉犻犲狉狋狉犪狀狊犳狅狉犿犪犾犵狅狉犻狋犺犿
狏犻犪犳狉犪犮狋犻狅狀犪犾犱狅犿犪犻狀犱犲犮狅犿狆狅狊犻狋犻狅狀
ＸｉｅＷｅｉｑｕｎ ＬｉｕＺｈｏｎｇｊｉｅ ＭａＳｈｉｗｅｉ
（ＳｃｈｏｏｌｏｆＭｅｃｈａｔｒｏｎｉｃａｌＥｎｇｉｎｅｅｒｉｎｇ＆Ａｕｔｏｍａｔｉｏｎ，ＳｈａｎｇｈａｉＵｎｉｖｅｒｓｉｔｙ，
ＳｈａｎｇｈａｉＫｅｙＬａｂｏｒａｔｏｒｙｏｆＰｏｗｅｒＳｔａｔｉｏｎＡｕｔｏｍａｔｉｏｎＴｅｃｈｎｏｌｏｇｙ，Ｓｈａｎｇｈａｉ２０００７２）
犃犫狊狋狉犪犮狋：ＦｒａｃｔｉｏｎａｌＦｏｕｒｉｅｒＴｒａｎｓｆｏｒｍｃａｎｂｅｗｉｄｅｌｙｕｓｅｄｉｎｓｃｉｅｎｔｉｆｉｃｃａｌｃｕｌａｔｉｏｎａｎｄｅｎｇｉｎｅｅｒｉｎｇ．Ｈｏｗｅｖｅｒｉｔｌａｃｋｓ
ｆａｓｔａｌｇｏｒｉｔｈｍｓｆｏｒｄｉｓｃｒｅｔｅＦｒａｃｔｉｏｎａｌＦｏｕｒｉｅｒＴｒａｎｓｆｏｒｍ（ＤＦＲＦＴ）．Ｉｎｔｈｉｓｐａｐｅｒ，ａＤＦＲＦＴａｌｇｏｒｉｔｈｍｕｓｉｎｇｆｒａｃｔｉｏｎａｌ
ｄｏｍａｉｎｄｅｃｏｍｐｏｓｉｔｉｏｎｉｓｐｒｏｐｏｓｅｄ．Ｗｉｔｈｔｈｉｓｍｅｔｈｏｄ，ｔｈｅＤＦＲＦＴａｔａｎｙｆｒａｃｔｉｏｎａｌｄｏｍａｉｎｏｆａｓｉｇｎａｌｃａｎｂｅｏｂｔａｉｎｅｄｂｙ
ａｗｅｉｇｈｔｅｄｃｏｍｂｉｎａｔｉｏｎｏｆｉｔｓＤＦＲＦＴｓａｔｓｐｅｃｉａｌｆｒａｃｔｉｏｎａｌｄｏｍａｉｎｓ，ｔｈｅｗｅｉｇｈｔｉｎｇｃｏｅｆｆｉｃｉｅｎｔｓａｒｅｏｂｔａｉｎｅｄｆｒｏｍａｎ
ｉｎｖｅｒｓｅＤｉｓｃｒｅｔｅＦｏｕｒｉｅｒＴｒａｎｓｆｏｒｍ（ＩＤＦＴ）ｃｏｍｐｕｔａｔｉｏｎｏｆａｓｅｒｉｅｓｃｏｎｓｉｓｔｅｄｏｆＤＦＴＨｅｒｍｉｔｅｅｉｇｅｎｖａｌｕｅｓａｎｄａｐｐｅｎｄｅｄ
ｚｅｒｏｓ．Ｉｎｓｅａｒｃｈｉｎｇｏｐｔｉｍａｌｆｒａｃｔｉｏｎａｌｄｏｍａｉｎ，ｉｔｄｏｅｓｎｔｔｏｃａｌｃｕｌａｔｅａｌｌｏｆｔｈｅｔｒａｎｓｆｏｒｍｋｅｒｎｅｌｓａｎｄｊｕｓｔｎｅｅｄｔｏｄｏａｎ
ＩＤＦＴｃｏｍｐｕｔａｔｉｏｎｏｎｃｅ．Ｈｅｎｃｅ，ｔｈｅｃｏｍｐｕｔａｔｉｏｎａｌｂｕｒｄｅｎｓｃａｎｂｅｒｅｄｕｃｅｄｅｆｆｉｃｉｅｎｔｌｙ．Ｉｔｃａｎｂｅｕｓｅｄｆｏｒｍｕｌｔｉ
ｃｏｍｐｏｎｅｎｔｓｉｇｎａｌｄｅｔｅｃｔｉｏｎａｎｄｆｉｌｔｅｒｉｎｇｉｎｆｒａｃｔｉｏｎａｌＦｏｕｒｉｅｒｄｏｍａｉｎｓ．Ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｒｅｓｕｌｔｓｍａｎｉｆｅｓｔｅｄｔｈｅ
ｅｆｆｅｃｔｉｖｅｎｅｓｓｏｆｔｈｉｓｍｅｔｈｏｄ．
犓犲狔狑狅狉犱狊：ｆｒａｃｔｉｏｎａｌＦｏｕｒｉｅｒｔｒａｎｓｆｏｒｍ；ｄｉｓｃｒｅｔｅｆｒａｃｔｉｏｎａｌＦｏｕｒｉｅｒｔｒａｎｓｆｏｒｍ；ｆｒａｃｔｉｏｎａｌｄｏｍａｉｎｄｅｃｏｍｐｏｓｉｔｉｏｎ
０ @ A
PZＦｏｕｒｉｅｒ*（ＦＲＦＴ）¤Ｆｏｕｒｉｅｒ*8æ:H
í
，
WÌ¢âãS§;ϸ8.ä
，
Ã~+WXH
âã8§;h
，
SY0Å
、
$"Ô÷
、
T÷âãñ
Ä
、
}âÙÚHæ-8g~R
［１］。
"Ti8g~
ijＦＲＦＴ88Hí<©+Ä，º*Î8P
ZＦｏｕｒｉｅｒ*（ＤＦＲＦＴ）Ä8Ë"`$8I>［２］。
û_.8Ä=BN8H ＤＦＲＦＴ0ºNÄ、
ＤＦＲＦＴp«NÄ、:¦*$NＤＦＲＦＴÄ，0ºN
?p«NÄÞd<?ＦＲＦＴëÁ8@é［３４］，
:¦*$N ＤＦＲＦＴ Ä»~Jâã88
Ｆｏｕｒｉｅｒ*（ＤＦＴ）~:¦¤8PZ]9:，Y~
ＤＦＴＨｅｒｍｉｔｅ:¦*$^¢<:¦*$，WX?
ＦＲＦＴëÁ8@é［５６］。פ，Ã~:¦*$NＤＦＲＦＴÄ
，
ÌØZP§
，
ɲ³u5Ä<Îg8$~
，
Y
5Ä$~Eâã8@
，
5Ä+8
。
:@È.
.%Ç»~ZPÅ8ＤＦＲＦＴÄ，ÝÞiÍÈ..Ã
~TÄÎ`Fâã²há;vw8¡#@é
。
１ 犉犚犉犜+«
~犉φîＦＲＦＴÄ"，y§Úâã狊（狋）8ＦＲＦＴ
9:¢
：
犉φ｛狊（狋）｝＝∫犽φ（狋，狌）狊（狋）ｄ狋 （１）
í²
：φ¤.äNd，$/P犽φX２π¢!)，î¢
·３６·
更多技术文章，论文请登录www.srvee.com
内容版权归作者所有
'３２( ) * 0 1 2 3
犓φ（狋，狌）＝
δ（狌－狋） φ＝０
犲－犼２π狌狋 φ＝π／２
δ（狌＋狋） φ＝π
ｅｊ２πｕｔ φ＝３π／２
１－ｊ犮狅狋槡 φｅｊ
（犮狅狋φ·ｔ
２
＋犮狅狋φ·ｕ
２
－２犮狊犮φ·ｔ·ｕ），
<
烅
烄
烆 |
（２）
Ä"犉φ HXhÐ：
.ä犉２狀π ＝犐；I@*
犉π／２＝犉；.äë¨犉α犉β＝犉α＋β；２π.ä犉２π＝犐，<²犐
¢HÄ"
。
φc§，ＦＲＦＴ®î.âã§Ú
î;Úî8¯c!0
。
?φ∈ ［０，２π］，»φ＝
犪π／２，犪¢iP，犪∈［０，４］，®5ＦＲＦＴ8*$rX４¢
!)
。
P÷V
，ＦＲＦＴ¤ＦＴ8æ:Hí，WÌ¢
hYÄ"^~TJ/PV
，
}!49?ÓTibP
犪8æ:cYÄ"，8犪ZＦＲＦＴr¤ＦＴÄ"8犪
Ø]
。
®5
，犉０¤HÄ"，犉１¤ＤＦＴ*Ä"。
２ BC;¬d+犇犉犚犉犜e
"TＤＦＴ8*$~;H１，犼，－１，－ｊ«４
:¦¤
，
Y<:¦*$d¤k%8¤"４*$ϸA
>
。
"°¾,_..9T ＤＦＴ $~:¦*$8
ＤＦＲＦＴÄ［３］，+9:¢ＤＦＴ$~:¦*$8PZ
]
，
BYg~ＤＦＴＨｅｒｍｉｔｅ:¦*$AàＤＦＲＦＴ8:
¦*$
，
<*Ä"î¢
：
犉犪＝犉２φπ犞Ｔ＝
∑
犖－１
犽＝０
犲－犼犽φ狏犽狏犽Ｔ 犖¢aP
∑
犖－２
犽＝０
犲－犼犽φ狏犽狏犽Ｔ＋犲－犼犖φ狏犖狏犖Ｔ 犖¢¼
烅
烄
烆 P
（３）
í²
：犖î8âã1&8àd；犞¤"ＤＦＴＨｅｒｍｉｔｅ
:¦*$7A>8°~
，狏犽¤ＤＦＴ8'犽 Ｈｅｒｍｉｔｅ
:¦*$
，
犖¢aP§，犞＝［狏０狘狏１狘…狘狏犖－２狘狏犖－１］，
犖¢¼P§，犞＝［狏０狘狏１狘…狘狏犖－２狘狏犖］；犇２φ／π¤ＤＦＲＦＴ
8*$~
，
+¤ÎN¢ＤＦＲＦＴ:¦¤犲－ｊ犽犪π／２
8ÎN~
，
:¦¤8夫yÑ１7î，<²犿
W?f߬P
，
>ß8¤犖¢¼P§UG%:¦
¤HH9A
。
|１ 犇犉犚犉犜PiC;£®£¯
犖 :¦¤
４犿 犲－犼犽犪π／２，犽＝０，１，２，…，（４犿－２），４犿
４犿＋１ 犲－犼犽犪π／２，犽＝０，１，２，…，（４犿－１），４犿
４犿＋２ 犲－犼犽犪π／２，犽＝０，１，２，…，４犿，（４犿＋２）
４犿＋３ 犲－犼犽犪π／２，犽＝０，１，２，…，（４犿＋１），４犿＋２
®5ÎTàd¢ 犖 88âã1&狓，<犪Z
ＤＦＲＦＴ¢：
犡犪＝犉犪｛狓｝＝
∑
犖－１
犽＝０
犲－犼犽犪π／２狏犽狏犽Ｔ狓 犖¢aP
∑
犖－２
犽＝０
犲－犼犽犪π／２狏犽狏犽Ｔ狓＋犲－犼犖犪π／２狏犖狏犖Ｔ狓 犖 ¢¼
烅
烄
烆 P
（４）
Ç9T:¦*$8ＤＦＲＦＴÄS犪?犖 %9§，
<ＤＦＴＨｅｒｍｉｔｅ:¦*$ＤＦＲＦＴ$~ÉWXÍ3
5ÄC
，
פ²³5ÄＤＦＲＦＴ$~Eâã1&8@，
í
（４）(H8~Eæ$8@ßÄú8d¢犗（犖２）。
ì¢ê8¤
，
ZP犪§，m-¤ ＤＦＲＦＴÀ¤
ＦＲＦＴÄÉ
u5Ä<*$，ßCuSäU
<PZÚ8!0²¹``t¨5Ä$
，
ÎTＦＲＦＴS±
Ñâã¡#
、
Uoqr?Mígwg~²Ë9U<P
ZÚ¤dÃ8
。
３ 犇犉犚犉犜+W8ªe
¢.YcＤＦＲＦＴ85Ä，².fßZP8ＤＦＲＦＴ
ÉWX}!%9ZPＤＦＲＦＴ8¨>?î。"１
ßＤＦＲＦＴ$~:¦¤S犖¢¼P§<å¤H9A，S
°&犖¢aP?¼P8DEwH&¦［７］。
Ñé狓¤àd¢犖 88âã1&，Y犖 ¢aP，
犱犽＝犲－ｊ犽犪π／２îＤＦＲＦＴ$~'犽:¦¤（犽＝０，１，
…，犖－１），1&犱＝［犱０，犱１，…，犱犖－１，０，…，０］，<²犱犖－１
GÁ
（犕－犖）
，YＭ≥Ｎ。Õ狓8犪ZＤＦＲＦＴW
X>?Hí
犡犪＝∑
犕－１
犿＝０
犮犿，犪犡犿，犫 （５）
y<¨>FP犮犿，犪W"1&犱8ＩＤＦＴÈ.
犮犿，犪＝ １犕∑
犕－１
犽＝０
犱犽ｅ犼（２π／犕）ｋｍ （６）
Y犫＝４／犕。
Ñé狓¤àd¢犖 88âã1&，Y犖¢¼P，
îＤＦＲＦＴ$~'犽:¦¤（犽＝０，１，…，犖－２），
犱犖－１＝ｅ－犼犖犪π／２î'（犖－１）:¦¤，1&犱＝［犱０，犱１，
…，犱犖－２，０，犱犖－１，０，…，０］，<²犱犖－１GÁ （犕－犖－１）
，
Y犕≥犖＋１。Õ狓8犪ZＤＦＲＦＴWX>
?Hí
犡犪＝∑
犕－１
犿＝０
犮犿，犪犡犿，犫 （７）
y¨>FP犮犿，犪W"1&犱8ＩＤＦＴÈ.
犮犿，犪＝ １犕∑
犕－１
犽＝０
犱犽ｅ犼（２π／犕）ｋｍ （８）
Y犫＝４／犕。
VoÝí¥
，
%âã8fßZＤＦＲＦＴÉWX}
!Tâã8Õv9ZＤＦＲＦＴ8¨>?5Ä，¨
>FPWX"ＤＦＲＦＴ$~8:¦¤（ＤＦＴＨｅｒｍｉｔｅ:
¦¤
）
?¨
7A>1&8ＩＤＦＴà，kÇÄ
¢ＤＦＲＦＴ8ZPÅÄ。S¨>úPd8DE
·４６·
更多技术文章，论文请登录www.srvee.com
内容版权归作者所有
;
Ｆｏｕｒｉｅｒ`#?>¯§?fg '２)
，
35Ä.âãS9ZP8ＤＦＲＦＴBý<，
ZP犪§，æ
5Ä%ØＩＤＦＴ（WX~ＦＦＴ©+
5Ä
［８］），
á~9ZP8ＤＦＲＦＴ，¶@BW
ＤＦＲＦＴ。m³<|ＤＦＲＦＴÄ®«u5Ä<7H
*$
。
¯°
，
S
¹¸U<PZÚ§
，
Ñâã¡#
?qrg~²
，
Ç9TZPÅ8ＤＦＲＦＴÄWH
[\5Ä$
。
４ pq/ß
há;âãæ-g~TY1
、
Uó?³â]ã
#FG
，
×+Î<|âãV¤%ÇHIvwâã
，¯
°S
i8g~²SS
¡#â㲤 )Hhá;âã
，
Bª¥qrñÄ
。
"ＦＲＦＴ89:Wß，ＦＲＦＴ$/P
8á;{Ehá;âã8á;{%µ§
，
há;
âã8ＦＲＦＴý´jò/P¦，$/P8á;{
?ÓT.äNdBPZÚ8ZP犪。7XÃ~±%U<
ZP8ＦＲＦＴWXý§Ú（犪＝０）²8há;âã
>PZÚ8òâã
，
1¡#há;âã8û
8
，
BWXÎHHhá;vw8$â㪥qr
ñÄ
。
（ａ）PZÚ犪＝０．５
（ｂ）PZÚ犪＝１．１５
（ｃ）PZÚ犪＝１．５
v１ âãSdlPZÚ8ＤＦＲＦＴ
（
û¢ZPÅÄ
，
¢:¦*$Ä
）
4Jâã¢Jhá;vwµ¶8`Fâã¢
ｅ－（狋－３０）
２／２０＋０．１ｅ犼（狋
２／１０－２狋），
»~ ＭＡＴＬＡＢç0wij.
ＤＦＲＦＴ8ZPÅÄ?9T:¦*$8Ä，BÎT
â㪥*ñÄ
，
<²â㻫àd犖＝２０５，:¦¤?
Á
1&àd?犕＝２２０。ÝÞ@é¥，}!Øá¬
ＤＦＲＦＴ8ZP犪，WX¸S犪＝１．１５§âã>·
¥Y8
，
ë¢犪＝１．１５r¤U<PZÚ，Ñv１
（ｂ）7î，<²²ST)ô１００～１２０·¸.j¤8âã
>·¢há;âã
，
ST)ô８０～１００·¸¤8â
ã>·¢`Fâã
。
５ L {
9TZPÅ8ＤＦＲＦＴÄ»~9ZPＤＦＲＦＴ
¨>?Hí5ÄfßZP8ＤＦＲＦＴ，¨ >FP"ＤＦＴ
Ｈｅｒｍｉｔｅ:¦¤?¨
7A>1&8ＩＤＦＴ，TÄ
dæWXE9T:¦*$8ＤＦＲＦＴÄëÁ8@
é
，
YSäU<PZÚ8!0²æ5Ä%Ø
ＩＤＦＴ，m³<|ＦＲＦＴ®«u5Ä<7H*
$
，
WXH[\ßÄ$
，
wTg~SPZ
Ú8$âã¡#?qrñIJ
。
; <
［１］ ＮＡＭＩＡＳＶ．ＴｈｅｆｒａｃｔｉｏｎａｌｏｒｄｅｒＦｏｕｒｉｅｒｔｒａｎｓｆｏｒｍ
ａｎｄｉｔｓａｐｐｌｉｃａｔｉｏｎｔｏｑｕａｎｔｕｍｍｅｃｈａｎｉｃｓ［Ｊ］．Ｊｏｕｒｎａｌ
ｏｆＩｎｓｔｉｔｕｔｅｏｆ ＭａｔｈｅｍａｔｉｃＡｐｐｌｉｃａｔｉｏｎ，１９８０，２５：
２４１２６５．
［２］ ＣＡＮＤＡＮＣ，ＯＺＡＫＴＡＳ Ｈ Ｍ．Ｔｈｅ Ｄｉｓｃｒｅｔｅ
ＦｒａｃｔｉｏｎａｌＦｏｕｒｉｅｒＴｒａｎｓｆｏｒｍ［Ｊ］．ＩＥＥＥＴｒａｎｓａｃｔｉｏｎｓ
ｏｎＳｉｇｎａｌＰｒｏｃｅｓｓｉｎｇ，２０００，４８：１３２９１３３７．
［３］ ＯＺＡＫＴＡＳＨＭ，ＡＮＫＡＮ Ｏ，ＫＵＴＡＹ Ａ，ｅｔａｌ．
Ｄｉｇｉｔａｌ Ｃｏｍｐｕｔａｔｉｏｎ ｏｆ ｔｈｅ Ｆｒａｃｔｉｏｎａｌ Ｆｏｕｒｉｅｒ
Ｔｒａｎｓｆｏｒｍ ［Ｊ］．ＩＥＥＥ Ｔｒａｎｓａｃｔｉｏｎ ｏｎ Ｓｉｇｎａｌ
Ｐｒｏｃｅｓｓｉｎｇ，１９９６，４４：２１４１２１５０． （EFG１０４H）
·５６·
更多技术文章，论文请登录www.srvee.com
内容版权归作者所有
浙公网安备 33021202002483