📄 fft.lis
字号:
.module fft.c
.area text(rom, con, rel)
0000 .dbfile D:\HJN\Design\AVR\fft\fft.c
0000 .dbfunc e four1 _four1 fV
0000 ; theta -> y+52
0000 ; wpi -> y+48
0000 ; wpr -> y+44
0000 ; wtemp -> y+40
0000 ; m -> y+36
0000 ; istep -> y+32
0000 ; n -> y+28
0000 ; mmax -> y+24
0000 ; wr -> y+20
0000 ; wi -> y+16
0000 ; tempi -> y+12
0000 ; tempr -> y+8
0000 ; i -> y+4
0000 ; j -> y+0
0000 ; isign -> R10,R11
0000 ; nn -> y+64
0000 ; data -> R12,R13
.even
0000 _four1::
0000 3A93 st -y,r19
0002 2A93 st -y,r18
0004 0E940000 xcall push_gset4
0008 6801 movw R12,R16
000A E897 sbiw R28,56
000C FE01 movw R30,R28
000E EC5B subi R30,188 ; addi 68
0010 FF4F sbci R31,255
0012 A080 ldd R10,z+0
0014 B180 ldd R11,z+1
0016 .dbline -1
0016 .dbline 9
0016 ; /*data是输入和输出（复数），nn是FFT的点数 isign指示FFT变换方向1为正变换，
0016 ; -1为反变换，程序出自《Numerical Recipes in C》 */
0016 ;
0016 ;
0016 ; #include <math.h>
0016 ; #define SWAP(a,b) tempr=(a);(a)=(b);(b)=tempr
0016 ;
0016 ; void four1(float data[], unsigned long nn, int isign)
0016 ; {
0016 .dbline 14
0016 ; unsigned long n,mmax,m,j,istep,i;
0016 ; double wtemp,wr,wpr,wpi,wi,theta;
0016 ; float tempr,tempi;
0016 ;
0016 ; n=nn << 1; //nn乘以2赋给n,作为取样点数
0016 FE01 movw R30,R28
0018 E05C subi R30,192 ; addi 64
001A FF4F sbci R31,255
001C 2080 ldd R2,z+0
001E 3180 ldd R3,z+1
0020 4280 ldd R4,z+2
0022 5380 ldd R5,z+3
0024 220C lsl R2
0026 331C rol R3
0028 441C rol R4
002A 551C rol R5
002C FE01 movw R30,R28
002E 248E std z+28,R2
0030 358E std z+29,R3
0032 468E std z+30,R4
0034 578E std z+31,R5
0036 .dbline 15
0036 ; j=1;
0036 41E0 ldi R20,1
0038 50E0 ldi R21,0
003A 60E0 ldi R22,0
003C 70E0 ldi R23,0
003E FE01 movw R30,R28
0040 4083 std z+0,R20
0042 5183 std z+1,R21
0044 6283 std z+2,R22
0046 7383 std z+3,R23
0048 .dbline 16
0048 41E0 ldi R20,1
004A 50E0 ldi R21,0
004C 60E0 ldi R22,0
004E 70E0 ldi R23,0
0050 FE01 movw R30,R28
0052 4483 std z+4,R20
0054 5583 std z+5,R21
0056 6683 std z+6,R22
0058 7783 std z+7,R23
005A 79C1 xjmp L5
005C L2:
005C .dbline 16
005C ; for (i=1;i<n;i+=2) {
005C .dbline 17
005C ; if (j > i) {
005C FE01 movw R30,R28
005E 2480 ldd R2,z+4
0060 3580 ldd R3,z+5
0062 4680 ldd R4,z+6
0064 5780 ldd R5,z+7
0066 FE01 movw R30,R28
0068 6080 ldd R6,z+0
006A 7180 ldd R7,z+1
006C 8280 ldd R8,z+2
006E 9380 ldd R9,z+3
0070 2614 cp R2,R6
0072 3704 cpc R3,R7
0074 4804 cpc R4,R8
0076 5904 cpc R5,R9
0078 08F0 brlo X0
007A F6C0 xjmp L6
007C X0:
007C .dbline 17
007C .dbline 18
007C ; SWAP(data[j],data[i]);
007C FE01 movw R30,R28
007E 2080 ldd R2,z+0
0080 3180 ldd R3,z+1
0082 4280 ldd R4,z+2
0084 5380 ldd R5,z+3
0086 44E0 ldi R20,4
0088 50E0 ldi R21,0
008A 60E0 ldi R22,0
008C 70E0 ldi R23,0
008E 5A92 st -y,R5
0090 4A92 st -y,R4
0092 3A92 st -y,R3
0094 2A92 st -y,R2
0096 8A01 movw R16,R20
0098 9B01 movw R18,R22
009A 0E940000 xcall empy32u
009E F801 movw R30,R16
00A0 EC0D add R30,R12
00A2 FD1D adc R31,R13
00A4 2080 ldd R2,z+0
00A6 3180 ldd R3,z+1
00A8 4280 ldd R4,z+2
00AA 5380 ldd R5,z+3
00AC FE01 movw R30,R28
00AE 2086 std z+8,R2
00B0 3186 std z+9,R3
00B2 4286 std z+10,R4
00B4 5386 std z+11,R5
00B6 .dbline 18
00B6 FE01 movw R30,R28
00B8 2480 ldd R2,z+4
00BA 3580 ldd R3,z+5
00BC 4680 ldd R4,z+6
00BE 5780 ldd R5,z+7
00C0 44E0 ldi R20,4
00C2 50E0 ldi R21,0
00C4 60E0 ldi R22,0
00C6 70E0 ldi R23,0
00C8 5A92 st -y,R5
00CA 4A92 st -y,R4
00CC 3A92 st -y,R3
00CE 2A92 st -y,R2
00D0 8A01 movw R16,R20
00D2 9B01 movw R18,R22
00D4 0E940000 xcall empy32u
00D8 F801 movw R30,R16
00DA EC0D add R30,R12
00DC FD1D adc R31,R13
00DE 2080 ldd R2,z+0
00E0 3180 ldd R3,z+1
00E2 4280 ldd R4,z+2
00E4 5380 ldd R5,z+3
00E6 FE01 movw R30,R28
00E8 6080 ldd R6,z+0
00EA 7180 ldd R7,z+1
00EC 8280 ldd R8,z+2
00EE 9380 ldd R9,z+3
00F0 44E0 ldi R20,4
00F2 50E0 ldi R21,0
00F4 60E0 ldi R22,0
00F6 70E0 ldi R23,0
00F8 9A92 st -y,R9
00FA 8A92 st -y,R8
00FC 7A92 st -y,R7
00FE 6A92 st -y,R6
0100 8A01 movw R16,R20
0102 9B01 movw R18,R22
0104 0E940000 xcall empy32u
0108 F801 movw R30,R16
010A EC0D add R30,R12
010C FD1D adc R31,R13
010E 2082 std z+0,R2
0110 3182 std z+1,R3
0112 4282 std z+2,R4
0114 5382 std z+3,R5
0116 .dbline 18
0116 FE01 movw R30,R28
0118 2480 ldd R2,z+4
011A 3580 ldd R3,z+5
011C 4680 ldd R4,z+6
011E 5780 ldd R5,z+7
0120 44E0 ldi R20,4
0122 50E0 ldi R21,0
0124 60E0 ldi R22,0
0126 70E0 ldi R23,0
0128 5A92 st -y,R5
012A 4A92 st -y,R4
012C 3A92 st -y,R3
012E 2A92 st -y,R2
0130 8A01 movw R16,R20
0132 9B01 movw R18,R22
0134 0E940000 xcall empy32u
0138 F801 movw R30,R16
013A EC0D add R30,R12
013C FD1D adc R31,R13
013E DE01 movw R26,R28
0140 1896 adiw R26,8
0142 2D90 ld R2,x+
0144 3D90 ld R3,x+
0146 4D90 ld R4,x+
0148 5C90 ld R5,x
014A 2082 std z+0,R2
014C 3182 std z+1,R3
014E 4282 std z+2,R4
0150 5382 std z+3,R5
0152 .dbline 19
0152 ; SWAP(data[j+1],data[i+1]);
0152 41E0 ldi R20,1
0154 50E0 ldi R21,0
0156 60E0 ldi R22,0
0158 70E0 ldi R23,0
015A FE01 movw R30,R28
015C 2080 ldd R2,z+0
015E 3180 ldd R3,z+1
0160 4280 ldd R4,z+2
0162 5380 ldd R5,z+3
0164 240E add R2,R20
0166 351E adc R3,R21
0168 461E adc R4,R22
016A 571E adc R5,R23
016C 44E0 ldi R20,4
016E 50E0 ldi R21,0
0170 60E0 ldi R22,0
0172 70E0 ldi R23,0
0174 5A92 st -y,R5
0176 4A92 st -y,R4
0178 3A92 st -y,R3
017A 2A92 st -y,R2
017C 8A01 movw R16,R20
017E 9B01 movw R18,R22
0180 0E940000 xcall empy32u
0184 F801 movw R30,R16
0186 EC0D add R30,R12
0188 FD1D adc R31,R13
018A 2080 ldd R2,z+0
018C 3180 ldd R3,z+1
018E 4280 ldd R4,z+2
0190 5380 ldd R5,z+3
0192 FE01 movw R30,R28
0194 2086 std z+8,R2
0196 3186 std z+9,R3
0198 4286 std z+10,R4
019A 5386 std z+11,R5
019C .dbline 19
019C 41E0 ldi R20,1
019E 50E0 ldi R21,0
01A0 60E0 ldi R22,0
01A2 70E0 ldi R23,0
01A4 FE01 movw R30,R28
01A6 2480 ldd R2,z+4
01A8 3580 ldd R3,z+5
01AA 4680 ldd R4,z+6
01AC 5780 ldd R5,z+7
01AE 240E add R2,R20
01B0 351E adc R3,R21
01B2 461E adc R4,R22
01B4 571E adc R5,R23
01B6 44E0 ldi R20,4
01B8 50E0 ldi R21,0
01BA 60E0 ldi R22,0
01BC 70E0 ldi R23,0
01BE 5A92 st -y,R5
01C0 4A92 st -y,R4
01C2 3A92 st -y,R3
01C4 2A92 st -y,R2
01C6 8A01 movw R16,R20
01C8 9B01 movw R18,R22
01CA 0E940000 xcall empy32u
01CE F801 movw R30,R16
01D0 EC0D add R30,R12
01D2 FD1D adc R31,R13
01D4 2080 ldd R2,z+0
01D6 3180 ldd R3,z+1
01D8 4280 ldd R4,z+2
01DA 5380 ldd R5,z+3
01DC 41E0 ldi R20,1
01DE 50E0 ldi R21,0
01E0 60E0 ldi R22,0
01E2 70E0 ldi R23,0
01E4 FE01 movw R30,R28
01E6 6080 ldd R6,z+0
01E8 7180 ldd R7,z+1
01EA 8280 ldd R8,z+2
01EC 9380 ldd R9,z+3
01EE 640E add R6,R20
01F0 751E adc R7,R21
01F2 861E adc R8,R22
01F4 971E adc R9,R23
01F6 44E0 ldi R20,4
01F8 50E0 ldi R21,0
01FA 60E0 ldi R22,0
01FC 70E0 ldi R23,0
01FE 9A92 st -y,R9
0200 8A92 st -y,R8
0202 7A92 st -y,R7
0204 6A92 st -y,R6
0206 8A01 movw R16,R20
0208 9B01 movw R18,R22
020A 0E940000 xcall empy32u
020E F801 movw R30,R16
0210 EC0D add R30,R12
0212 FD1D adc R31,R13
0214 2082 std z+0,R2
0216 3182 std z+1,R3
0218 4282 std z+2,R4
021A 5382 std z+3,R5
021C .dbline 19
021C 41E0 ldi R20,1
021E 50E0 ldi R21,0
0220 60E0 ldi R22,0
0222 70E0 ldi R23,0
0224 FE01 movw R30,R28
0226 2480 ldd R2,z+4
0228 3580 ldd R3,z+5
022A 4680 ldd R4,z+6
022C 5780 ldd R5,z+7
022E 240E add R2,R20
0230 351E adc R3,R21
0232 461E adc R4,R22
0234 571E adc R5,R23
0236 44E0 ldi R20,4
0238 50E0 ldi R21,0
023A 60E0 ldi R22,0
023C 70E0 ldi R23,0
023E 5A92 st -y,R5
0240 4A92 st -y,R4
0242 3A92 st -y,R3
0244 2A92 st -y,R2
0246 8A01 movw R16,R20
0248 9B01 movw R18,R22
024A 0E940000 xcall empy32u
024E F801 movw R30,R16
0250 EC0D add R30,R12
0252 FD1D adc R31,R13
0254 DE01 movw R26,R28
0256 1896 adiw R26,8
0258 2D90 ld R2,x+
025A 3D90 ld R3,x+
025C 4D90 ld R4,x+
025E 5C90 ld R5,x
0260 2082 std z+0,R2
0262 3182 std z+1,R3
0264 4282 std z+2,R4
0266 5382 std z+3,R5
0268 .dbline 20
0268 ; }
0268 L6:
0268 .dbline 21
0268 ; m=n >> 1; //将n除以2的值赋给m
0268 FE01 movw R30,R28
026A 248C ldd R2,z+28
026C 358C ldd R3,z+29
026E 468C ldd R4,z+30
0270 578C ldd R5,z+31
0272 5694 lsr R5
0274 4794 ror R4
0276 3794 ror R3
0278 2794 ror R2
027A FE01 movw R30,R28
027C 24A2 std z+36,R2
027E 35A2 std z+37,R3
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -