|
|||||||||||
| 技術(shù)交流 | 電路欣賞 | 工控天地 | 數(shù)字廣電 | 通信技術(shù) | 電源技術(shù) | 測(cè)控之家 | EMC技術(shù) | ARM技術(shù) | EDA技術(shù) | PCB技術(shù) | 嵌入式系統(tǒng) 驅(qū)動(dòng)編程 | 集成電路 | 器件替換 | 模擬技術(shù) | 新手園地 | 單 片 機(jī) | DSP技術(shù) | MCU技術(shù) | IC 設(shè)計(jì) | IC 產(chǎn)業(yè) | CAN-bus/DeviceNe |
FFT快速付立葉級(jí)數(shù) |
| 作者：firstman 欄目：單片機(jī) |
最近用到快速付立葉級(jí)數(shù)來分析諧波,但找過所有網(wǎng)站竟沒有介紹,哪位高人能介紹一下,用匯編怎樣實(shí)現(xiàn)付立葉級(jí)數(shù)展開! |
| 2樓： | >>參與討論 |
| 作者： firstman 于 2003/9/26 21:15:00 發(fā)布：
51FFT清單 ;12點(diǎn)付立葉分析,UI、UR相當(dāng)于付立葉系數(shù)a1,b1 ;UI=(2/N)*(U0-U6+SQRT(3)/2(U1-U5-U7+U11)+(U2-U4-U8+U10)/2) ;UR=(2/N)*(U3-U9+SQRT(3)/2(U2+U4-U8-U10)+(U1+U5-U7-U11)/2) DATAB EQU 00H ;(12點(diǎn)采樣值浮點(diǎn)數(shù)存放首址) UIREG EQU 30H URREG EQU 33H ;sqrt(3)/2=0.866(00 dd b2) ;============================= FFT:NOP UI:MOV DPTR,#DATAB+3 ;U1 MOV R0,#DAT0 MOV R2,#3 LCALL MOVFFT0 MOV DPTR,#DATAB+15 ;U5 LCALL MOVFFT MOV R0,#DAT0 MOV R1,#DAT4 LCALL FSUB ;U1-U5 MOV DPTR,#DATAB+21 ;U7 LCALL MOVFFT MOV R0,#DAT0 MOV R1,#DAT4 LCALL FSUB ;U1-U5-U7 MOV DPTR,#DATAB+33 ;U11 LCALL MOVFFT MOV R0,#DAT0 MOV R1,#DAT4 LCALL FADD ;U1-U5-U7+U11 MOV DAT4,#00 MOV DAT5,#0DDH MOV DAT6,#0B2H MOV R0,#DAT0 MOV R1,#DAT4 LCALL FMUL ;*0.866 MOV DPTR,#DATAB+6 ;U2 LCALL MOVFFT LCALL FDEC MOV R0,#DAT0 MOV R1,#DAT4 LCALL FADD ;+U2/2 MOV DPTR,#DATAB+12 ;U4 LCALL MOVFFT LCALL FDEC MOV R0,#DAT0 MOV R1,#DAT4 LCALL FSUB ;-U4/2 MOV DPTR,#DATAB+24 ;U8 LCALL MOVFFT LCALL FDEC MOV R0,#DAT0 MOV R1,#DAT4 LCALL FSUB ;-U8/2 MOV DPTR,#DATAB+30 ;U10 LCALL MOVFFT LCALL FDEC MOV R0,#DAT0 MOV R1,#DAT4 LCALL FADD ;+U10/2 MOV DPTR,#DATAB ;U0 LCALL MOVFFT MOV R0,#DAT0 MOV R1,#DAT4 LCALL FADD ;+U0 MOV DPTR,#DATAB+18 ;U6 LCALL MOVFFT MOV R0,#DAT0 MOV R1,#DAT4 LCALL FSUB ;-U6 MOV DAT4,#03H MOV DAT5,#0C0H MOV DAT6,#00H MOV R0,#DAT0 MOV R1,#DAT4 LCALL FDIV ;/6 ;UI MOV R1,#UIREG MOV @R1,DAT0 INC R1 MOV @R1,DAT1 INC R1 MOV @R1,DAT2 ;UI->[UIREG] UR:MOV DPTR,#DATAB+6 ;U2 MOV R0,#DAT0 MOV R2,#3 LCALL MOVFFT0 MOV DPTR,#DATAB+12 ;U4 LCALL MOVFFT MOV R0,#DAT0 MOV R1,#DAT4 LCALL FADD ;U2+U4 MOV DPTR,#DATAB+24 ;U8 LCALL MOVFFT MOV R0,#DAT0 MOV R1,#DAT4 LCALL FSUB ;-U8 MOV DPTR,#DATAB+30 ;U10 LCALL MOVFFT MOV R0,#DAT0 MOV R1,#DAT4 LCALL FSUB ;-U10 MOV DAT4,#00 MOV DAT5,#0DDH MOV DAT6,#0B2H MOV R0,#DAT0 MOV R1,#DAT4 LCALL FMUL ;*0.866 MOV DPTR,#DATAB+3 ;U1 LCALL MOVFFT LCALL FDEC MOV R0,#DAT0 MOV R1,#DAT4 LCALL FADD ;+U1/2 MOV DPTR,#DATAB+15 ;U5 LCALL MOVFFT LCALL FDEC MOV R0,#DAT0 MOV R1,#DAT4 LCALL FADD ;+U5/2 MOV DPTR,#DATAB+21 ;U7 LCALL MOVFFT LCALL FDEC MOV R0,#DAT0 MOV R1,#DAT4 LCALL FSUB ;-U7/2 MOV DPTR,#DATAB+33 ;U11 LCALL MOVFFT LCALL FDEC MOV R0,#DAT0 MOV R1,#DAT4 LCALL FSUB ;-U11/2 MOV DPTR,#DATAB+9 ;U3 LCALL MOVFFT MOV R0,#DAT0 MOV R1,#DAT4 LCALL FADD ;+U3 MOV DPTR,#DATAB+27 ;U9 LCALL MOVFFT MOV R0,#DAT0 MOV R1,#DAT4 LCALL FSUB ;-U9 MOV DAT4,#03H MOV DAT5,#0C0H MOV DAT6,#00H MOV R0,#DAT0 MOV &nbs |
|
|
|
| 免費(fèi)注冊(cè)為維庫(kù)電子開發(fā)網(wǎng)會(huì)員，參與電子工程師社區(qū)討論，點(diǎn)此進(jìn)入 |
Copyright © 1998-2006 m.58mhw.cn 浙ICP證030469號(hào) |