设为首页收藏本站
切换到宽版

矿石收音机论坛

 找回密码
 加入会员

QQ登录

只需一步，快速开始

搜索
矿石收音机论坛»论坛 › ∷技术区∷ › 〓电脑与单片机〓 › cortex-m0的32点FFT程序
返回列表 发新帖
查看: 6806|回复: 14

cortex-m0的32点FFT程序

[复制链接]
yjmwxwx
     
发表于 2018-1-4 15:02:52 | 显示全部楼层 |阅读模式
经过好几天的努力终于写出来了，上次写的太垃圾了。

Radix-2的算法，输入32位数，从0到31。
程序所有文件 fft11.zip (9.58 KB, 下载次数: 261)


所有小数的旋转因子先乘32768,计算完后右移15位，目的是绕开小数计算，

上次很傻居然乘10。写完才知道原来FFT的程序这么简单，都不如驱动数码管的程序复杂。


计算旋转因子的程序，比如FFT的第三步需要4个旋转因子，N就填8

  1. #include <stdio.h>
  2. #include <math.h>
  3. #include <stdlib.h>
  4. void main() {
  5.   float pi,sr,si;

  7.   int rr,ii, n,i;
  8.   n=8;
  9.   pi=3.1415927;   
  10.         for(i=0;i<n/2;i++){
  11.            sr=cos(2*pi*i/n);
  12.            si=-sin(2*pi*i/n);
  13.         rr=sr*32768;
  14.         ii=si*32768;
  15.        printf("0x%04X,0x%04X,",rr,ii);
  16. }
  17.         printf("\n");
  18. }
复制代码



计算器计算的
截图_2018-01-04_14-40-37.png

单片机计算的
截图_2018-01-04_14-38-54.png


从0x20000000开始到0X20000080是实部
从0X20000080开始到0X20000100是虚部

全部顺序排列

和计算器计算的大体对照下感觉差不多，单片机的输出的是十六进制，负数都是0XFFFFXXXX这样的。


FFT在CORTEX-M0单片机上的程序

  1.         @ stm32f030f4p6 asm
  2.         @fft 32点
  3.          .thumb                 
  4.          .syntax unified
  5. .section .data       
  6.                .equ STACKINIT,                 0x20001000
  7.         .equ xuanzhuanyinzi,                0x20000120
  8. xinxi:        .ascii "fft-32dian-2018-01-04"
  9.         .align 4
  10. shuru:        .int 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15
  11.         .int 16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31
  12. hudie:        .int 0,64,32,96,16,80,48,112,8,72,40,104,24,88,56,120
  13.         .int 4,68,36,100,20,84,52,116,12,76,44,108,28,92,60,124
  14. fft3xuanzhuanyinzi:
  15.         .int 0x8000,0x0000,0x5A82,0xFFFFA57E,0x0000,0xFFFF8000,0xFFFFA57E,0xFFFFA57E
  16. fft4xuanzhuanyinzi:
  17.         .int 0x8000,0x0000,0x7641,0xFFFFCF05,0x5A82,0xFFFFA57E,0x30FB,0xFFFF89BF,0x0000,0xFFFF8000,0xFFFFCF05,0xFFFF89BF,0xFFFFA57E,0xFFFFA57E,0xFFFF89BF,0xFFFFCF05
  18. fft5xuanzhuanyinzi:
  19.         .int        0x8000,0x0000,0x7D8A,0xFFFFE708,0x7641,0xFFFFCF05,0x6A6D,0xFFFFB8E4,0x5A82,0xFFFFA57E,0x471C,0xFFFF9593,0x30FB,0xFFFF89BF,0x18F8,0xFFFF8276,0x0000,0xFFFF8000,0xFFFFE708,0xFFFF8276,0xFFFFCF05,0xFFFF89BF,0xFFFFB8E4,0xFFFF9593,0xFFFFA57E,0xFFFFA57E,0xFFFF9593,0xFFFFB8E4,0xFFFF89BF,0xFFFFCF05,0xFFFF8276,0xFFFFE708

  21.        
  22.                    
  23.                    
  24.        
  25. .section .text
  26. vectors:        
  27.         .word STACKINIT         
  28.         .word _start + 1        
  29.         .word _nmi_handler + 1  
  30.         .word _hard_fault  + 1  
  31.         .word 0
  32.         .word 0   
  33.         .word 0
  34.         .word 0
  35.         .word 0
  36.         .word 0
  37.         .word 0
  38.         .word _svc_handler +1
  39.         .word 0
  40.         .word 0
  41.         .word _pendsv_handler +1
  42.         .align 2
  43. _start:
  44. _shizhong:                             @时钟设置
  45.         ldr r2, = 0x40022000   @FLASH访问控制
  46.         movs r1, # 0x32
  47.         str r1, [r2]
  48.         ldr r0, = 0x40021000 @ rcc
  49.         @0x34时钟控制寄存器 2 (RCC_CR2)
  50.         movs r1, # 0x01
  51.         str r1, [r0, # 0x34]  @ HSI开14M时钟
  52. _dengdai14mshizhongwending:            @等14M时钟稳定
  53.         ldr r1, [r0, # 0x34]
  54.         lsls r1, r1, # 30     @ 左移30位
  55.         bpl _dengdai14mshizhongwending  @ 等待14M时钟稳定

  57. _neicunqingling:                                 @ 0x1000的内存清零
  58.                                          @ 1K=1024BIT=0X400
  59.         ldr r0, = 0x20000000
  60.         movs r1, # 0
  61.         ldr r3, = 0x1000
  62. _neicunqinglingxunhuan:                           @ 内存清零循环
  63.         subs r3, # 4
  64.         str r1, [r0, r3]
  65.         bne _neicunqinglingxunhuan                @



  69.        
  70. _fftjisuan:
  71. _fft1:       
  72.         ldr r0, = shuru
  73.         ldr r1, = 0x20000000
  74.         movs r2, # 16
  75.         ldr r3, = hudie
  76. _fft1xunhuan:
  77.         ldr r6, [r3]
  78.         ldr r7, [r3, # 0x04]
  79.         ldr r6, [r0, r6]
  80.         ldr r7, [r0, r7]
  81.         adds r4, r6, r7
  82.         subs r5, r6, r7
  83.         str r4, [r1]
  84.         str r5, [r1, # 0x04]
  85.         adds r1, r1, # 0x08
  86.         adds r3, r3, # 0x08
  87.         subs r2, r2, # 1
  88.         bne _fft1xunhuan
  89.        
  90. _fft2:
  91.         ldr r0, = 0x20000000
  92.         mov r4, r0
  93.         adds r4, r4, # 0x80
  94.         mov r8, r4
  95. _fft2xunhuan:
  96.         ldr r4, [r0]
  97.         ldr r5, [r0, # 0x08]
  98.         ldr r6, [r0, # 0x04]
  99.         ldr r7, [r0, # 0x0c]
  100.         adds r1, r4, r5
  101.         subs r2, r4, r5
  102.         mov r3, r6                @0r
  103.         movs r4, # 0       
  104.         subs r4, r4, r7                @0i
  105.         mov r5, r6                @3r
  106.         mov r6, r7                @3i
  107.         mov r7, r0
  108.         adds r7, r7, # 0x80
  109.         str r1, [r0]
  110.         str r2, [r0, # 0x08]
  111.         str r3, [r0, # 0x04]
  112.         str r4, [r7, # 0x04]
  113.         str r5, [r0, # 0x0c]
  114.         str r6, [r7, # 0x0c]
  115.         adds r0, r0, # 0x10
  116.         cmp r0, r8
  117.         bne _fft2xunhuan
  118. _fft3:
  119.         ldr r6, = xuanzhuanyinzi
  120.         ldr r2, = fft3xuanzhuanyinzi
  121.         ldr r0, = 0x20000000
  122.         str r2, [r6]
  123.         movs r1, r0
  124.         adds r1, r1, # 0x80
  125.         mov r8, r1
  126.         mov r9, r0
  127.         mov r10, r2
  128.         movs r5, # 0x10
  129.         mov r11, r5
  130.         movs r5, # 4
  131.         mov r12, r5
  132.         bl _fftg
  133. _fft4:
  134.         ldr r6, = xuanzhuanyinzi
  135.         ldr r2, = fft4xuanzhuanyinzi
  136.         ldr r0, = 0x20000000
  137.         str r2, [r6]
  138.         movs r1, r0
  139.         adds r1, r1, # 0x80
  140.         mov r8, r1
  141.         mov r9, r0
  142.         mov r10, r2
  143.         movs r5, # 0x20
  144.         mov r11, r5
  145.         movs r5, # 8
  146.         mov r12, r5
  147.         bl _fftg
  148. _fft5:
  149.         ldr r6, = xuanzhuanyinzi
  150.         ldr r2, = fft5xuanzhuanyinzi
  151.         ldr r0, = 0x20000000
  152.         str r2, [r6]
  153.         movs r1, r0
  154.         adds r1, r1, # 0x80
  155.         mov r8, r1
  156.         mov r9, r0
  157.         mov r10, r2
  158.         movs r5, # 0x40
  159.         mov r11, r5
  160.         movs r5, # 16
  161.         mov r12, r5
  162.         bl _fftg
  163.        
  164.        
  165. tingting:
  166.         b tingting
  167. _fftg:
  168.         push {lr}
  169.         movs r4, # 0
  170.         mov lr, r4
  171. _fftxunhuan:
  172.         mov r0, r9
  173.         mov r1, r8
  174.         mov r2, r10
  175.         mov r3, r11
  176.         ldr r4, [r0]
  177.         ldr r5, [r1]
  178.         ldr r6, [r0, r3]
  179.         ldr r7, [r1, r3]
  180.         ldr r3, [r2]                @sr
  181.         ldr r2, [r2, # 0x04]        @si
  182.         push {r4,r5}
  183.         mov r4, r6
  184.         muls r4, r4, r3
  185.         mov r5, r7
  186.         muls r5, r5, r2
  187.         subs r4, r4, r5                @r       
  188.         asrs r4, r4, # 15
  189.         mov r5, r6
  190.         muls r5, r5, r2
  191.         muls r7, r7, r3
  192.         adds r7, r7, r5                @i
  193.         asrs r7, r7, # 15
  194.         pop {r2,r5}
  195.         adds r3, r2, r4    @0r
  196.         subs r2, r2, r4    @1r
  197.         adds r4, r5, r7
  198.         subs r5, r5, r7
  199.         mov r6, r11
  200.         str r3, [r0]
  201.         str r4, [r1]
  202.         str r2, [r0, r6]
  203.         str r5, [r1, r6]
  204.         adds r0, r0, # 0x04
  205.         adds r1, r1, # 0x04
  206.         mov r9, r0
  207.         mov r8, r1
  208.         mov r2, r10
  209.         adds r2, r2, # 0x08
  210.         mov r10, r2
  211.         mov r4, lr
  212.         adds r4, r4, # 1
  213.         mov lr, r4
  214.         mov r7, r12
  215.         cmp r4, r7
  216.         bne _fftxunhuan
  217.         movs  r4, # 0
  218.         ldr r6, = xuanzhuanyinzi
  219.         ldr r5, = 0x20000080
  220.         ldr r6, [r6]
  221.         mov r10, r6
  222.         mov lr, r4
  223.         mov r7, r11
  224.         adds r0, r0, r7
  225.         adds r1, r1, r7
  226.         mov r9, r0
  227.         mov r8, r1
  228.         cmp r0, r5
  229.         bne _fftxunhuan
  230.         pop {pc}
  231.                
  232. _nmi_handler:
  233.         bx lr
  234. _hard_fault:
  235.         bx lr
  236. _svc_handler:
  237.         bx lr
  238. _pendsv_handler:
  239.         bx lr
复制代码




回复

使用道具 举报

leo121
     
发表于 2018-1-5 00:19:46 | 显示全部楼层
其实arm都提供这个例程的，不用自己写的
回复 支持 反对

使用道具 举报

yjmwxwx
     
 楼主| 发表于 2018-1-5 00:37:54 | 显示全部楼层
la45088d1 发表于 2018-1-4 22:26
建议：ADC定时采样，FFT-128，LED或者TFT显示频谱。这样好玩些。

慢慢学，以前没接触过ADC采样显示频谱，先把用到的几个函数写出来。
回复 支持 反对

使用道具 举报

yjmwxwx
     
 楼主| 发表于 2018-1-5 01:09:56 | 显示全部楼层
leo121 发表于 2018-1-5 00:19
其实arm都提供这个例程的，不用自己写的

有库，CMSIS DSP库里面有，不过我没用过不会用，等把这个库里的函数提取出来研究研究。

有现成收音机IC为啥还有这么多人做矿石机，一个道理就是玩呗。
回复 支持 反对

使用道具 举报

zystbbs
     
发表于 2018-1-5 09:18:28 | 显示全部楼层
算法很重要啊，这么肯定简单了。
回复 支持 反对

使用道具 举报

leo121
     
发表于 2018-1-5 12:54:43 | 显示全部楼层
yjmwxwx 发表于 2018-1-5 01:09
有库，CMSIS DSP库里面有，不过我没用过不会用，等把这个库里的函数提取出来研究研究。

有现成收音机I ...

矿坛牛人多啊
回复 支持 反对

使用道具 举报

yjmwxwx
     
 楼主| 发表于 2018-1-5 15:07:42 | 显示全部楼层
把用到辅助小程序贴出来。

这个是生成位反转表的小程序，结果乘4因为每个数据32位。
  1. //位反转
  2. #include <stdio.h>
  3. void main() {
  4.         int aa,a2,a3,a4,a5,bb,cc,dd;


  7.         bb=256;                                //长度
  8.         dd=bb;
  9.         for(aa=0;aa<bb;aa++){
  10.                 cc=1;
  11.                 a3=0;
  12.                 a4=a3;
  13.                 while(cc<dd){
  14.                 a2=aa & cc;
  15.                 a5=bb>>1;
  16.                 if(a2!=0){
  17.                 a5=a5>>a4;
  18.                 a3=a3 | a5;
  19. }
  20.                 cc=cc<<1;
  21.                 a4=a4+1;

  23.        
  24. }
  25.         printf("0x%04X,",a3*4);               
  26. }       
  27.         printf("\n");
  28. }

  30.        
复制代码


下面这个是生成旋转因子的函数
  1. #include <stdio.h>
  2. #include <math.h>
  3. #include <stdlib.h>
  4. void main() {
  5.   float pi,sr,si;

  7.   int rr,ii, n,i;
  8.   n=256;
  9.   pi=3.1415927;   
  10.         for(i=0;i<n/2;i++){
  11.            sr=cos(2*pi*i/n);
  12.            si=-sin(2*pi*i/n);
  13.         rr=sr*0x8000;
  14.         ii=si*0x8000;
  15.        printf("0x%04X,0x%04X,",rr,ii);
  16. }
  17.         printf("\n");
  18. }
  19.        
复制代码


上边这俩自己写的，写的不太好。
回复 支持 反对

使用道具 举报

yjmwxwx
     
 楼主| 发表于 2018-1-5 15:19:22 | 显示全部楼层
fft64.zip (10.65 KB, 下载次数: 224)

fft128.zip (312.91 KB, 下载次数: 203)

fft256.zip (521.98 KB, 下载次数: 225)


64、128、256的都写出来了，不过好像误差很大。


下面这个是256的，输入0-255

1.png

截图_2018-01-05_14-51-55.png

比较开头一部分，误差有点大


256点的程序

  1.         @ stm32f030f4p6 asm
  2.         @fft 256点
  3.         @2018-01-05
  4.          .thumb                 
  5.          .syntax unified
  6. .section .data       
  7.                .equ STACKINIT,                 0x20001000
  8.         .equ xuanzhuanyinzi,                0x20000804
  9.         .align 4
  10. shuru:
  11.         .int 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126,127,128,129,130,131,132,133,134,135,136,137,138,139,140,141,142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,163,164,165,166,167,168,169,170,171,172,173,174,175,176,177,178,179,180,181,182,183,184,185,186,187,188,189,190,191,192,193,194,195,196,197,198,199,200,201,202,203,204,205,206,207,208,209,210,211,212,213,214,215,216,217,218,219,220,221,222,223,224,225,226,227,228,229,230,231,232,233,234,235,236,237,238,239,240,241,242,243,244,245,246,247,248,249,250,251,252,253,254,255
  12.        
  13. hudieweifanzhuan:
  14.         .int 0x0000,0x0200,0x0100,0x0300,0x0080,0x0280,0x0180,0x0380,0x0040,0x0240,0x0140,0x0340,0x00C0,0x02C0,0x01C0,0x03C0,0x0020,0x0220,0x0120,0x0320,0x00A0,0x02A0,0x01A0,0x03A0,0x0060,0x0260,0x0160,0x0360,0x00E0,0x02E0,0x01E0,0x03E0,0x0010,0x0210,0x0110,0x0310,0x0090,0x0290,0x0190,0x0390,0x0050,0x0250,0x0150,0x0350,0x00D0,0x02D0,0x01D0,0x03D0,0x0030,0x0230,0x0130,0x0330,0x00B0,0x02B0,0x01B0,0x03B0,0x0070,0x0270,0x0170,0x0370,0x00F0,0x02F0,0x01F0,0x03F0,0x0008,0x0208,0x0108,0x0308,0x0088,0x0288,0x0188,0x0388,0x0048,0x0248,0x0148,0x0348,0x00C8,0x02C8,0x01C8,0x03C8,0x0028,0x0228,0x0128,0x0328,0x00A8,0x02A8,0x01A8,0x03A8,0x0068,0x0268,0x0168,0x0368,0x00E8,0x02E8,0x01E8,0x03E8,0x0018,0x0218,0x0118,0x0318,0x0098,0x0298,0x0198,0x0398,0x0058,0x0258,0x0158,0x0358,0x00D8,0x02D8,0x01D8,0x03D8,0x0038,0x0238,0x0138,0x0338,0x00B8,0x02B8,0x01B8,0x03B8,0x0078,0x0278,0x0178,0x0378,0x00F8,0x02F8,0x01F8,0x03F8,0x0004,0x0204,0x0104,0x0304,0x0084,0x0284,0x0184,0x0384,0x0044,0x0244,0x0144,0x0344,0x00C4,0x02C4,0x01C4,0x03C4,0x0024,0x0224,0x0124,0x0324,0x00A4,0x02A4,0x01A4,0x03A4,0x0064,0x0264,0x0164,0x0364,0x00E4,0x02E4,0x01E4,0x03E4,0x0014,0x0214,0x0114,0x0314,0x0094,0x0294,0x0194,0x0394,0x0054,0x0254,0x0154,0x0354,0x00D4,0x02D4,0x01D4,0x03D4,0x0034,0x0234,0x0134,0x0334,0x00B4,0x02B4,0x01B4,0x03B4,0x0074,0x0274,0x0174,0x0374,0x00F4,0x02F4,0x01F4,0x03F4,0x000C,0x020C,0x010C,0x030C,0x008C,0x028C,0x018C,0x038C,0x004C,0x024C,0x014C,0x034C,0x00CC,0x02CC,0x01CC,0x03CC,0x002C,0x022C,0x012C,0x032C,0x00AC,0x02AC,0x01AC,0x03AC,0x006C,0x026C,0x016C,0x036C,0x00EC,0x02EC,0x01EC,0x03EC,0x001C,0x021C,0x011C,0x031C,0x009C,0x029C,0x019C,0x039C,0x005C,0x025C,0x015C,0x035C,0x00DC,0x02DC,0x01DC,0x03DC,0x003C,0x023C,0x013C,0x033C,0x00BC,0x02BC,0x01BC,0x03BC,0x007C,0x027C,0x017C,0x037C,0x00FC,0x02FC,0x01FC,0x03FC
  15. fft3xuanzhuanyinzi:
  16.         .int 0x8000,0x0000,0x5A82,0xFFFFA57E,0x0000,0xFFFF8000,0xFFFFA57E,0xFFFFA57E
  17. fft4xuanzhuanyinzi:
  18.         .int 0x8000,0x0000,0x7641,0xFFFFCF05,0x5A82,0xFFFFA57E,0x30FB,0xFFFF89BF,0x0000,0xFFFF8000,0xFFFFCF05,0xFFFF89BF,0xFFFFA57E,0xFFFFA57E,0xFFFF89BF,0xFFFFCF05
  19. fft5xuanzhuanyinzi:
  20.         .int 0x8000,0x0000,0x7D8A,0xFFFFE708,0x7641,0xFFFFCF05,0x6A6D,0xFFFFB8E4,0x5A82,0xFFFFA57E,0x471C,0xFFFF9593,0x30FB,0xFFFF89BF,0x18F8,0xFFFF8276,0x0000,0xFFFF8000,0xFFFFE708,0xFFFF8276,0xFFFFCF05,0xFFFF89BF,0xFFFFB8E4,0xFFFF9593,0xFFFFA57E,0xFFFFA57E,0xFFFF9593,0xFFFFB8E4,0xFFFF89BF,0xFFFFCF05,0xFFFF8276,0xFFFFE708
  21. fft6xuanzhuanyinzi:
  22.         .int 0x8000,0x0000,0x7F62,0xFFFFF375,0x7D8A,0xFFFFE708,0x7A7D,0xFFFFDAD8,0x7641,0xFFFFCF05,0x70E2,0xFFFFC3AA,0x6A6D,0xFFFFB8E4,0x62F2,0xFFFFAECD,0x5A82,0xFFFFA57E,0x5133,0xFFFF9D0E,0x471C,0xFFFF9593,0x3C56,0xFFFF8F1E,0x30FB,0xFFFF89BF,0x2528,0xFFFF8583,0x18F8,0xFFFF8276,0x0C8B,0xFFFF809E,0x0000,0xFFFF8000,0xFFFFF375,0xFFFF809E,0xFFFFE708,0xFFFF8276,0xFFFFDAD8,0xFFFF8583,0xFFFFCF05,0xFFFF89BF,0xFFFFC3AA,0xFFFF8F1E,0xFFFFB8E4,0xFFFF9593,0xFFFFAECD,0xFFFF9D0E,0xFFFFA57E,0xFFFFA57E,0xFFFF9D0E,0xFFFFAECD,0xFFFF9593,0xFFFFB8E4,0xFFFF8F1E,0xFFFFC3AA,0xFFFF89BF,0xFFFFCF05,0xFFFF8583,0xFFFFDAD8,0xFFFF8276,0xFFFFE708,0xFFFF809E,0xFFFFF375
  23. fft7xuanzhuanyinzi:
  24.         .int 0x8000,0x0000,0x7FD8,0xFFFFF9B9,0x7F62,0xFFFFF375,0x7E9D,0xFFFFED38,0x7D8A,0xFFFFE708,0x7C29,0xFFFFE0E7,0x7A7D,0xFFFFDAD8,0x7884,0xFFFFD4E1,0x7641,0xFFFFCF05,0x73B5,0xFFFFC946,0x70E2,0xFFFFC3AA,0x6DCA,0xFFFFBE32,0x6A6D,0xFFFFB8E4,0x66CF,0xFFFFB3C1,0x62F2,0xFFFFAECD,0x5ED7,0xFFFFAA0B,0x5A82,0xFFFFA57E,0x55F5,0xFFFFA129,0x5133,0xFFFF9D0E,0x4C3F,0xFFFF9931,0x471C,0xFFFF9593,0x41CE,0xFFFF9236,0x3C56,0xFFFF8F1E,0x36BA,0xFFFF8C4B,0x30FB,0xFFFF89BF,0x2B1F,0xFFFF877C,0x2528,0xFFFF8583,0x1F19,0xFFFF83D7,0x18F8,0xFFFF8276,0x12C8,0xFFFF8163,0x0C8B,0xFFFF809E,0x0647,0xFFFF8028,0x0000,0xFFFF8000,0xFFFFF9B9,0xFFFF8028,0xFFFFF375,0xFFFF809E,0xFFFFED38,0xFFFF8163,0xFFFFE708,0xFFFF8276,0xFFFFE0E7,0xFFFF83D7,0xFFFFDAD8,0xFFFF8583,0xFFFFD4E1,0xFFFF877C,0xFFFFCF05,0xFFFF89BF,0xFFFFC946,0xFFFF8C4B,0xFFFFC3AA,0xFFFF8F1E,0xFFFFBE32,0xFFFF9236,0xFFFFB8E4,0xFFFF9593,0xFFFFB3C1,0xFFFF9931,0xFFFFAECD,0xFFFF9D0E,0xFFFFAA0B,0xFFFFA129,0xFFFFA57E,0xFFFFA57E,0xFFFFA129,0xFFFFAA0B,0xFFFF9D0E,0xFFFFAECD,0xFFFF9931,0xFFFFB3C1,0xFFFF9593,0xFFFFB8E4,0xFFFF9236,0xFFFFBE32,0xFFFF8F1E,0xFFFFC3AA,0xFFFF8C4B,0xFFFFC946,0xFFFF89BF,0xFFFFCF05,0xFFFF877C,0xFFFFD4E1,0xFFFF8583,0xFFFFDAD8,0xFFFF83D7,0xFFFFE0E7,0xFFFF8276,0xFFFFE708,0xFFFF8163,0xFFFFED38,0xFFFF809E,0xFFFFF375,0xFFFF8028,0xFFFFF9B9

  26. fft8xuanzhuanyinzi:
  27.         .int 0x8000,0x0000,0x7FF6,0xFFFFFCDC,0x7FD8,0xFFFFF9B9,0x7FA7,0xFFFFF696,0x7F62,0xFFFFF375,0x7F09,0xFFFFF055,0x7E9D,0xFFFFED38,0x7E1D,0xFFFFEA1E,0x7D8A,0xFFFFE708,0x7CE3,0xFFFFE3F5,0x7C29,0xFFFFE0E7,0x7B5D,0xFFFFDDDD,0x7A7D,0xFFFFDAD8,0x798A,0xFFFFD7DA,0x7884,0xFFFFD4E1,0x776C,0xFFFFD1EF,0x7641,0xFFFFCF05,0x7504,0xFFFFCC22,0x73B5,0xFFFFC946,0x7255,0xFFFFC674,0x70E2,0xFFFFC3AA,0x6F5F,0xFFFFC0E9,0x6DCA,0xFFFFBE32,0x6C24,0xFFFFBB86,0x6A6D,0xFFFFB8E4,0x68A6,0xFFFFB64C,0x66CF,0xFFFFB3C1,0x64E8,0xFFFFB141,0x62F2,0xFFFFAECD,0x60EC,0xFFFFAC65,0x5ED7,0xFFFFAA0B,0x5CB4,0xFFFFA7BE,0x5A82,0xFFFFA57E,0x5842,0xFFFFA34C,0x55F5,0xFFFFA129,0x539B,0xFFFF9F14,0x5133,0xFFFF9D0E,0x4EBF,0xFFFF9B18,0x4C3F,0xFFFF9931,0x49B4,0xFFFF975A,0x471C,0xFFFF9593,0x447A,0xFFFF93DC,0x41CE,0xFFFF9236,0x3F17,0xFFFF90A1,0x3C56,0xFFFF8F1E,0x398C,0xFFFF8DAB,0x36BA,0xFFFF8C4B,0x33DE,0xFFFF8AFC,0x30FB,0xFFFF89BF,0x2E11,0xFFFF8894,0x2B1F,0xFFFF877C,0x2826,0xFFFF8676,0x2528,0xFFFF8583,0x2223,0xFFFF84A3,0x1F19,0xFFFF83D7,0x1C0B,0xFFFF831D,0x18F8,0xFFFF8276,0x15E2,0xFFFF81E3,0x12C8,0xFFFF8163,0x0FAB,0xFFFF80F7,0x0C8B,0xFFFF809E,0x096A,0xFFFF8059,0x0647,0xFFFF8028,0x0324,0xFFFF800A,0x0000,0xFFFF8000,0xFFFFFCDC,0xFFFF800A,0xFFFFF9B9,0xFFFF8028,0xFFFFF696,0xFFFF8059,0xFFFFF375,0xFFFF809E,0xFFFFF055,0xFFFF80F7,0xFFFFED38,0xFFFF8163,0xFFFFEA1E,0xFFFF81E3,0xFFFFE708,0xFFFF8276,0xFFFFE3F5,0xFFFF831D,0xFFFFE0E7,0xFFFF83D7,0xFFFFDDDD,0xFFFF84A3,0xFFFFDAD8,0xFFFF8583,0xFFFFD7DA,0xFFFF8676,0xFFFFD4E1,0xFFFF877C,0xFFFFD1EF,0xFFFF8894,0xFFFFCF05,0xFFFF89BF,0xFFFFCC22,0xFFFF8AFC,0xFFFFC946,0xFFFF8C4B,0xFFFFC674,0xFFFF8DAB,0xFFFFC3AA,0xFFFF8F1E,0xFFFFC0E9,0xFFFF90A1,0xFFFFBE32,0xFFFF9236,0xFFFFBB86,0xFFFF93DC,0xFFFFB8E4,0xFFFF9593,0xFFFFB64C,0xFFFF975A,0xFFFFB3C1,0xFFFF9931,0xFFFFB141,0xFFFF9B18,0xFFFFAECD,0xFFFF9D0E,0xFFFFAC65,0xFFFF9F14,0xFFFFAA0B,0xFFFFA129,0xFFFFA7BE,0xFFFFA34C,0xFFFFA57E,0xFFFFA57E,0xFFFFA34C,0xFFFFA7BE,0xFFFFA129,0xFFFFAA0B,0xFFFF9F14,0xFFFFAC65,0xFFFF9D0E,0xFFFFAECD,0xFFFF9B18,0xFFFFB141,0xFFFF9931,0xFFFFB3C1,0xFFFF975A,0xFFFFB64C,0xFFFF9593,0xFFFFB8E4,0xFFFF93DC,0xFFFFBB86,0xFFFF9236,0xFFFFBE32,0xFFFF90A1,0xFFFFC0E9,0xFFFF8F1E,0xFFFFC3AA,0xFFFF8DAB,0xFFFFC674,0xFFFF8C4B,0xFFFFC946,0xFFFF8AFC,0xFFFFCC22,0xFFFF89BF,0xFFFFCF05,0xFFFF8894,0xFFFFD1EF,0xFFFF877C,0xFFFFD4E1,0xFFFF8676,0xFFFFD7DA,0xFFFF8583,0xFFFFDAD8,0xFFFF84A3,0xFFFFDDDD,0xFFFF83D7,0xFFFFE0E7,0xFFFF831D,0xFFFFE3F5,0xFFFF8276,0xFFFFE708,0xFFFF81E3,0xFFFFEA1E,0xFFFF8163,0xFFFFED38,0xFFFF80F7,0xFFFFF055,0xFFFF809E,0xFFFFF375,0xFFFF8059,0xFFFFF696,0xFFFF8028,0xFFFFF9B9,0xFFFF800A,0xFFFFFCDC
  28.                           
  29. .section .text
  30. vectors:        
  31.         .word STACKINIT         
  32.         .word _start + 1        
  33.         .word _nmi_handler + 1  
  34.         .word _hard_fault  + 1  
  35.         .word 0
  36.         .word 0   
  37.         .word 0
  38.         .word 0
  39.         .word 0
  40.         .word 0
  41.         .word 0
  42.         .word _svc_handler +1
  43.         .word 0
  44.         .word 0
  45.         .word _pendsv_handler +1
  46.         .align 2
  47. _start:
  48. _shizhong:                             @时钟设置
  49.         ldr r2, = 0x40022000   @FLASH访问控制
  50.         movs r1, # 0x32
  51.         str r1, [r2]
  52.         ldr r0, = 0x40021000 @ rcc
  53.         @0x34时钟控制寄存器 2 (RCC_CR2)
  54.         movs r1, # 0x01
  55.         str r1, [r0, # 0x34]  @ HSI开14M时钟
  56. _dengdai14mshizhongwending:            @等14M时钟稳定
  57.         ldr r1, [r0, # 0x34]
  58.         lsls r1, r1, # 30     @ 左移30位
  59.         bpl _dengdai14mshizhongwending  @ 等待14M时钟稳定

  61. _neicunqingling:                                 @ 0x1000的内存清零
  62.                                          @ 1K=1024BIT=0X400
  63.         ldr r0, = 0x20000000
  64.         movs r1, # 0
  65.         ldr r3, = 0x1000
  66. _neicunqinglingxunhuan:                           @ 内存清零循环
  67.         subs r3, # 4
  68.         str r1, [r0, r3]
  69.         bne _neicunqinglingxunhuan                @



  73.        
  74. _fftjisuan:
  75. _fft1:       
  76.         ldr r0, = shuru
  77.         ldr r1, = 0x20000000
  78.         movs r2, # 128
  79.         ldr r3, = hudieweifanzhuan
  80. _fft1xunhuan:
  81.         ldr r6, [r3]
  82.         ldr r7, [r3, # 0x04]
  83.         ldr r6, [r0, r6]
  84.         ldr r7, [r0, r7]
  85.         adds r4, r6, r7
  86.         subs r5, r6, r7
  87.         str r4, [r1]
  88.         str r5, [r1, # 0x04]
  89.         adds r1, r1, # 0x08
  90.         adds r3, r3, # 0x08
  91.         subs r2, r2, # 1
  92.         bne _fft1xunhuan
  93.        
  94. _fft2:
  95.         ldr r0, = 0x20000000
  96.         movs r4, # 1
  97.         lsls r4, r4, # 10
  98.         adds r4, r4, r0
  99.         mov r8, r4
  100. _fft2xunhuan:
  101.         ldr r4, [r0]
  102.         ldr r5, [r0, # 0x08]
  103.         ldr r6, [r0, # 0x04]
  104.         ldr r7, [r0, # 0x0c]
  105.         adds r1, r4, r5
  106.         subs r2, r4, r5
  107.         mov r3, r6                @0r
  108.         movs r4, # 0       
  109.         subs r4, r4, r7                @0i
  110.         mov r5, r6                @3r
  111.         mov r6, r7                @3i
  112.         movs r7, # 1
  113.         lsls r7, r7, # 10
  114.         adds r7, r7, r0
  115.         str r1, [r0]
  116.         str r2, [r0, # 0x08]
  117.         str r3, [r0, # 0x04]
  118.         str r4, [r7, # 0x04]
  119.         str r5, [r0, # 0x0c]
  120.         str r6, [r7, # 0x0c]
  121.         adds r0, r0, # 0x10
  122.         cmp r0, r8
  123.         bne _fft2xunhuan
  124. _fft3:
  125.         ldr r6, = xuanzhuanyinzi
  126.         ldr r2, = fft3xuanzhuanyinzi
  127.         ldr r0, = 0x20000000
  128.         str r2, [r6]
  129.         movs r1, # 1
  130.         lsls r1, r1, # 10
  131.         adds r1, r1, r0
  132.         mov r8, r1
  133.         mov r9, r0
  134.         mov r10, r2
  135.         movs r5, # 0x10
  136.         mov r11, r5
  137.         movs r5, # 4
  138.         mov r12, r5
  139.         bl _fftg
  140. _fft4:
  141.         ldr r6, = xuanzhuanyinzi
  142.         ldr r2, = fft4xuanzhuanyinzi
  143.         ldr r0, = 0x20000000
  144.         str r2, [r6]
  145.         movs r1, # 1
  146.         lsls r1, r1, # 10
  147.         adds r1, r1, r0
  148.         mov r8, r1
  149.         mov r9, r0
  150.         mov r10, r2
  151.         movs r5, # 0x20
  152.         mov r11, r5
  153.         movs r5, # 8
  154.         mov r12, r5
  155.         bl _fftg
  156. _fft5:
  157.         ldr r6, = xuanzhuanyinzi
  158.         ldr r2, = fft5xuanzhuanyinzi
  159.         ldr r0, = 0x20000000
  160.         str r2, [r6]
  161.         movs r1, # 1
  162.         lsls r1, r1, # 10
  163.         adds r1, r1, r0
  164.         mov r8, r1
  165.         mov r9, r0
  166.         mov r10, r2
  167.         movs r5, # 0x40
  168.         mov r11, r5
  169.         movs r5, # 16
  170.         mov r12, r5
  171.         bl _fftg

  173. _fft6:
  174.         ldr r6, = xuanzhuanyinzi
  175.         ldr r2, = fft6xuanzhuanyinzi
  176.         ldr r0, = 0x20000000
  177.         str r2, [r6]
  178.         movs r1, # 1
  179.         lsls r1, r1, # 10
  180.         adds r1, r1, r0
  181.         mov r8, r1
  182.         mov r9, r0
  183.         mov r10, r2
  184.         movs r5, # 0x80
  185.         mov r11, r5
  186.         movs r5, # 32
  187.         mov r12, r5
  188.         bl _fftg
  189.        
  190. _fft7:
  191.         ldr r6, = xuanzhuanyinzi
  192.         ldr r2, = fft7xuanzhuanyinzi
  193.         ldr r0, = 0x20000000
  194.         str r2, [r6]
  195.         movs r1, # 1
  196.         lsls r1, r1, # 10
  197.         adds r1, r1, r0
  198.         mov r8, r1
  199.         mov r9, r0
  200.         mov r10, r2
  201.         movs r5, # 1
  202.         lsls r5, r5, # 8
  203.         mov r11, r5
  204.         movs r5, # 64
  205.         mov r12, r5
  206.         bl _fftg

  208. _fft8:
  209.         ldr r6, = xuanzhuanyinzi
  210.         ldr r2, = fft8xuanzhuanyinzi
  211.         ldr r0, = 0x20000000
  212.         str r2, [r6]
  213.         movs r1, # 1
  214.         lsls r1, r1, # 10
  215.         adds r1, r1, r0
  216.         mov r8, r1
  217.         mov r9, r0
  218.         mov r10, r2
  219.         movs r5, # 1
  220.         lsls r5, r5, # 9
  221.         mov r11, r5
  222.         movs r5, # 128
  223.         mov r12, r5
  224.         bl _fftg
  225.        
  226.        
  227. tingting:
  228.         b tingting
  229. _fftg:
  230.         push {lr}
  231.         movs r4, # 0
  232.         mov lr, r4
  233. _fftxunhuan:
  234.         mov r0, r9
  235.         mov r1, r8
  236.         mov r2, r10
  237.         mov r3, r11
  238.         ldr r4, [r0]
  239.         ldr r5, [r1]
  240.         ldr r6, [r0, r3]
  241.         ldr r7, [r1, r3]
  242.         ldr r3, [r2]                @sr
  243.         ldr r2, [r2, # 0x04]        @si
  244.         push {r4,r5}
  245.         mov r4, r6
  246.         muls r4, r4, r3
  247.         mov r5, r7
  248.         muls r5, r5, r2
  249.         subs r4, r4, r5                @r
  250.         asrs r4, r4, # 15
  251.         mov r5, r6
  252.         muls r5, r5, r2
  253.         muls r7, r7, r3
  254.         adds r7, r7, r5    @i
  255.         asrs r7, r7, # 15
  256.         pop {r2,r5}
  257.         adds r3, r2, r4    @0r
  258.         subs r2, r2, r4    @1r
  259.         adds r4, r5, r7
  260.         subs r5, r5, r7
  261.         mov r6, r11
  262.         str r3, [r0]
  263.         str r4, [r1]
  264.         str r2, [r0, r6]
  265.         str r5, [r1, r6]
  266.         adds r0, r0, # 0x04
  267.         adds r1, r1, # 0x04
  268.         mov r9, r0
  269.         mov r8, r1
  270.         mov r2, r10
  271.         adds r2, r2, # 0x08
  272.         mov r10, r2
  273.         mov r4, lr
  274.         adds r4, r4, # 1
  275.         mov lr, r4
  276.         mov r7, r12
  277.         cmp r4, r7
  278.         bne _fftxunhuan
  279.         movs  r4, # 0
  280.         ldr r6, = xuanzhuanyinzi
  281.         ldr r5, = 0x20000400
  282.         ldr r6, [r6]
  283.         mov r10, r6
  284.         mov lr, r4
  285.         mov r7, r11
  286.         adds r0, r0, r7
  287.         adds r1, r1, r7
  288.         mov r9, r0
  289.         mov r8, r1
  290.         cmp r0, r5
  291.         bne _fftxunhuan
  292.         pop {pc}
  293.         .ltorg
  294. _nmi_handler:
  295.         bx lr
  296. _hard_fault:
  297.         bx lr
  298. _svc_handler:
  299.         bx lr
  300. _pendsv_handler:
  301.         bx lr
复制代码



回复 支持 反对

使用道具 举报

小徒弟下山
     
发表于 2018-1-5 15:32:32 来自手机 | 显示全部楼层
一大堆汇编，看的头晕
回复 支持 反对

使用道具 举报

jlfzx
     
发表于 2023-8-1 14:07:52 | 显示全部楼层
我也用的DSP库搞的fft。
不过您这个真厉害！
回复 支持 反对

使用道具 举报

iffi123
     
发表于 2023-8-1 14:36:50 | 显示全部楼层
本帖最后由 iffi123 于 2023-8-1 14:49 编辑

M3内核有ST官方汇编库(64,256,1024)，但fft最好用在有硬件浮点的芯片，速度快很多，84M主频，1024点耗时不到1ms,  现在F3,F4也很便宜
回复 支持 反对

使用道具 举报

bis
     
发表于 2023-8-1 18:05:56 | 显示全部楼层
提示: 作者被禁止或删除 内容自动屏蔽
回复 支持 反对

使用道具 举报

yjmwxwx
     
 楼主| 发表于 2023-8-1 18:10:30 | 显示全部楼层
iffi123 发表于 2023-8-1 14:36
M3内核有ST官方汇编库(64,256,1024)，但fft最好用在有硬件浮点的芯片，速度快很多，84M主频，1024点耗时不 ...

这多年下来我其实还没用过FFT做过什么东西。。
回复 支持 反对

使用道具 举报

yjmwxwx
     
 楼主| 发表于 2023-8-1 18:15:46 | 显示全部楼层
bis 发表于 2023-8-1 18:05
1）CMSIS里面有FFT。

2）M0没有FPU，浮点运算非常慢，CMSIS自带定点运算库。

多年下来感觉这数学计算类的程序写起来好麻烦，自己写的效果不好，速度慢、算的不准、动不动就溢出。。官方库当年也看过，看不懂具体怎么算的，人家写库的经验丰富写出来的东西确实好用。
回复 支持 反对

使用道具 举报

iffi123
     
发表于 2023-8-1 19:45:39 | 显示全部楼层
本帖最后由 iffi123 于 2023-8-1 19:50 编辑
yjmwxwx 发表于 2023-8-1 18:10
这多年下来我其实还没用过FFT做过什么东西。。


最多的就是频谱显示了，比如sdr, 音频播放这些，我用过官方库，不过包含相关dsp文件后，代码增加很多，主要那个系数表奇大无比，增加约100K
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 加入会员

本版积分规则

小黑屋|手机版|矿石收音机 ( 蒙ICP备05000029号-1 )

蒙公网安备 15040402000005号

GMT+8, 2025-4-27 21:42

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表