dsp fft 64 256 1024_stm32f103

nashui_sx · 发表于 2019-10-22 10:56:51

本帖最后由 nashui_sx 于 2021-1-17 21:56 编辑

设置点数，只能选这三个#define NPT 1024 //采样点数64 256 1024

dsp_fft_stm32f103.rar (324.35 KB, 下载次数: 404)

2021-1-17 19:58 上传

点击文件名下载附件

STM32F10x_DSP_Lib_V2.0.0_setup.rar (3.11 MB, 下载次数: 1454)

2019-10-22 10:56 上传

点击文件名下载附件
1.原始代码是1024 原始数据是随机的，后面以我测试时候随机数据说明

2.创建测试数据
记住固定采样频率ad采集后数据存在测试数据lBufInArray 采样频率保持与宏Fs一致

由于计算后的频率都是整数，想保留小数可以扩大Fs值，比如例子中采样频率是44800，设置成448000，计算的结果频率就会扩大10倍，得到一位小数，依次类推
void InitBufInArray(void)
{
      u16 i;
      float fx;
      for(i=0; i<NPT; i++)
      {
            fx = 500+                                        //0HZ    交流分量0 直流分量 500
                  1500 * sin(PI2 * i * 350.0 / Fs)  +  1500+ //350HZ 交流分量1500 直流分量1500
                  2700 * sin(PI2 * i * 8400.0 / Fs) +  2700+ //8400HZ  交流分量2700 直流分量2700
                  4000 * sin(PI2 * i * 18725.0 / Fs)+  4000; //18725HZ 交流分量4000 直流分量4000
      //数据特征：交流分量350HZ 1500 8400HZ 2700  18725HZ  4000  直流分量500+1500+2700+4000=8700
      //为了模拟adc都是正数值交流也加上直流分量变为正数
            lBufInArray = ((u16)fx) << 16;
      }
}

3.串口数据分析
原始数据+fft后数据

截取重要4行与原始数据高度吻合
...
  序号频率幅值实部虚部
   0,    0,  8696,    0,  8696
...
   2, 350,  1492, 746, -3
...
48,  8400,  2698,  1349, -1
...
107, 18725,  3996,  1998,    1
...

USART1 OK...
*********************************************************
  8700
13239
  7325
11885
  2764
10171
10848
  9901
12599
  2857
11113
  8201
12625
12247
  4529
11182
  6065
14766
11181
  7431
10339
  5012
15698
  9937
10866
  8872
  5262
15091
  9043
13975
  7311
  6636
13028
  8837
15982
  6253
  8632
  9998
  9358
16404
  6162
10590
  6766
10335
15184
  7206
11907
  4163
11291
12691
  9177
12217
  2852
11708
  9606
11541
11496
  3137
11221
  6718
13578
10051
  4881
  9759
  4700
14601
  8405
  7544
  7592
  3922
14168
  7116
10348
  5269
  4373
12232
  6585
12496
  3456
  5687
  9170
  6929
13406
  2729
  7290
  5696
  7938
12867
  3391
  8591
  2669
  9154
11089
  5352
  9178
852
10028
  8629
  8144
  8941
701
10114
  6209
11036
  8094
  2236
  9234
  4505
13244
  7088
  5044
  7555
  3943
14157
  6444
  8394
  5558
  4592
13518
  6563
11451
  3901
  6159
11508
  7579
13511
  3212
  8103
  8700
  9296
14187
  3888
  9820
  5891
11240
13498
  5948
10836
  3881
12807
11841
  9005
10955
  3242
13456
  9844
12355
10311
  4155
12894
  8165
15163
  9305
  6363
11190
  7285
16698
  8458
  9255
  8770
  7371
16547
  8221
12047
  6310
  8245
14730
  8808
14008
  4532
  9461
11703
10109
14670
  3993
10470
  8229
11712
13943
  4903
10814
  5167
13026
12130
  7051
10283
  3231
13477
  9807
  9855
  8994
  2798
12700
  7640
12518
  7348
  3821
10681
  6178
14262
  5903
  5858
  7793
  5691
14547
  5182
  8222
  4708
  6108
13236
  5492
10193
  2215
  7064
10633
  6809
11237
995
  8041
  7401
  8767
11146
  1417
  8562
  4371
10763
10088
  3424
  8356
  2308
12137
  8527
  6533
  7462
  1701
12387
  7060
  9968
  6218
  2633
11334
  6217
12870
  5152
  4774
  9198
  6286
14542
  4800
  7498
  6551
  7228
14635
  5514
10074
  4160
*********************************************************
  序号频率幅值实部虚部
   0,    0,  8696,    0,  8696
   1, 175, 14, -6, -4
   2, 350,  1492, 746, -3
   3, 525, 13, -6, -3
   4, 700, 10, -4, -3
   5, 875,    8, -4, -1
   6,  1050,    6, -3,    0
   7,  1225,    4, -2,    0
   8,  1400,    6, -3, -1
   9,  1575,    8, -4,    0
10,  1750,    6, -3, -1
11,  1925,    6, -3,    0
12,  2100,    6, -3, -1
13,  2275,    5, -2, -2
14,  2450,    7, -3, -2
15,  2625,    8, -3, -3
16,  2800,    2, -1,    0
17,  2975,    6, -3,    0
18,  3150,    6, -3,    0
19,  3325,    6, -3,    0
20,  3500,    6, -3,    0
21,  3675,    6, -3,    0
22,  3850,    4, -2,    0
23,  4025,    4, -2,    0
24,  4200,    6, -3,    0
25,  4375,    6, -3,    0
26,  4550,    6, -3,    0
27,  4725,    4, -2,    0
28,  4900,    4, -2, -1
29,  5075,    2, -1,    0
30,  5250,    4, -2,    0
31,  5425,    0,    0,    0
32,  5600,    4, -2, -1
33,  5775,    6, -3,    0
34,  5950,    2, -1, -1
35,  6125,    6, -3, -1
36,  6300,    2, -1, -1
37,  6475,    4, -2, -1
38,  6650,    2, -1,    0
39,  6825,    4, -2, -1
40,  7000,    4, -2, -1
41,  7175,    6, -3,    0
42,  7350,    4, -2, -1
43,  7525,    4, -2,    0
44,  7700,    4, -2, -1
45,  7875,    2, -1,    0
46,  8050,    4, -2, -1
47,  8225,    2, -1,    0
48,  8400,  2698,  1349, -1
49,  8575,    4, -2,    0
50,  8750,    0,    0,    0
51,  8925,    2, -1, -1
52,  9100,    2, -1, -1
53,  9275,    2, -1, -1
54,  9450,    2, -1, -1
55,  9625,    2, -1, -1
56,  9800,    4, -2, -1
57,  9975,    2, -1, -1
58, 10150,    6, -3, -1
59, 10325,    4, -2,    0
60, 10500,    2, -1,    0
61, 10675,    2, -1,    0
62, 10850,    4, -2, -1
63, 11025,    4, -2, -1
64, 11200,    0,    0,    0
65, 11375,    2, -1,    0
66, 11550,    0,    0,    0
67, 11725,    2,    0, -1
68, 11900,    4, -2, -1
69, 12075,    2, -1,    0
70, 12250,    2, -1,    1
71, 12425,    2, -1,    1
72, 12600,    2, -1,    0
73, 12775,    2, -1,    1
74, 12950,    2, -1,    0
75, 13125,    4, -2,    0
76, 13300,    4, -2,    0
77, 13475,    2, -1,    0
78, 13650,    2, -1, -1
79, 13825,    4, -2, -1
80, 14000,    2, -1,    1
81, 14175,    4, -2,    1
82, 14350,    2, -1,    1
83, 14525,    4, -2,    1
84, 14700,    4, -2,    0
85, 14875,    4, -2,    1
86, 15050,    4, -2,    0
87, 15225,    2,    0,    1
88, 15400,    4, -2,    1
89, 15575,    4, -2,    1
90, 15750,    4, -2,    1
91, 15925,    2, -1,    1
92, 16100,    2, -1,    1
93, 16275,    2,    0, -1
94, 16450,    4, -2,    1
95, 16625,    2, -1,    1
96, 16800,    2, -1, -1
97, 16975,    4, -2,    1
98, 17150,    2, -1,    0
99, 17325,    4, -2,    0
100, 17500,    2, -1,    0
101, 17675,    4, -2,    0
102, 17850,    2, -1,    1
103, 18025,    2, -1,    0
104, 18200,    4, -2,    1
105, 18375,    4, -2,    0
106, 18550,    4, -2,    1
107, 18725,  3996,  1998,    1
108, 18900,    2, -1,    1
109, 19075,    2, -1,    1
110, 19250,    4, -2,    0
111, 19425,    4, -2,    1
112, 19600,    2,    0,    1
113, 19775,    2, -1,    1
114, 19950,    0,    0,    0
115, 20125,    4, -2,    1
116, 20300,    2,    0,    1
117, 20475,    4, -2,    1
118, 20650,    2, -1,    1
119, 20825,    2,    0,    1
120, 21000,    2, -1,    1
121, 21175,    2, -1,    0
122, 21350,    4, -2,    1
123, 21525,    4, -1,    2
124, 21700,    0,    0,    0
125, 21875,    2, -1,    1
126, 22050,    2, -1,    1
127, 22225,    2,    0,    1
*********************************************************

anlog · 发表于 2019-10-27 01:46:55

mark，以备后查

迷途 · 发表于 2020-8-16 19:04:56

感谢楼主分享！！！

nashui_sx · 发表于 2021-1-17 21:57:27

总有人问，更新了些说明备注

帐号		自动登录	找回密码
密码			立即注册

dsp fft 64 256 1024_stm32f103

原子哥极力推荐 /2