OpenEdv-开源电子网

 找回密码
 立即注册
正点原子全套STM32/Linux/FPGA开发资料,上千讲STM32视频教程免费下载...
查看: 4087|回复: 10

stm32f103使用dsp库进行FFT计算数据出现很多错误,请大神们给点建议

[复制链接]

8

主题

42

帖子

0

精华

初级会员

Rank: 2

积分
161
金钱
161
注册时间
2014-1-7
在线时间
24 小时
发表于 2019-12-17 17:19:14 | 显示全部楼层 |阅读模式
1金钱
如题:stm32f103使用dsp库进行FFT计算数据出现很多错误,请大神们给点建议:
我用函数产生256个数据组成一组正弦波,运行FFT后发现很多数据出错,请问这种情况如何解决?
下面是我的代码:
u16 fs=1000;
float f=51.1234;
float a=1.2345;
float ph=30.123;
float dc=0.123;
float z1,z2,Z1,Z2,ka,r,fo,ao,pho;
float_t Y1;
uint32_t k;
u16 t=0;
float x[N],xfft[N],y[N];


for(t=0;t<N;t++)
{
  x[t]=a*sin(2*PI*f*t/fs+ph*PI/180)+dc; //波形产生
}  

  cr4_fft_256_stm32(&xfft, &x, N);// FFT计算

  arm_abs_f32(xfft,y,N);
  arm_max_f32(y,N,&Y1,&k);
  k =k+1;
  Z1=xfft[k]-(xfft[k-1]+xfft[k+1])/2;
  Z2=xfft[k+1]-(xfft[k]+xfft[k+2])/2;
  arm_abs_f32(&Z1,&z1,1);
  arm_abs_f32(&Z2,&z2,1);
  ka=z1/z2;
  r = (2-ka)/(1+ka);
  fo=(k-1+r)*fs/N;
  ao=(2*PI*r*(1-r*r)*z1)/(N*sin(r*PI));
波形产生的数据:

波形数据.jpg


FFT计算后的数据:(红框的都是一场数据)

FFT计算数据

FFT计算数据

请问这个有人遇到过吗?这个要怎么处理才能避免呢?



正点原子逻辑分析仪DL16劲爆上市
回复

使用道具 举报

8

主题

42

帖子

0

精华

初级会员

Rank: 2

积分
161
金钱
161
注册时间
2014-1-7
在线时间
24 小时
 楼主| 发表于 2019-12-17 17:20:37 | 显示全部楼层
请各位大神们不吝赐教,多多指点,谢谢
回复

使用道具 举报

6

主题

1127

帖子

0

精华

金牌会员

Rank: 6Rank: 6

积分
1656
金钱
1656
注册时间
2019-8-15
在线时间
102 小时
发表于 2019-12-17 18:56:16 | 显示全部楼层
帮顶                                    
成功没有捷径
回复

使用道具 举报

8

主题

42

帖子

0

精华

初级会员

Rank: 2

积分
161
金钱
161
注册时间
2014-1-7
在线时间
24 小时
 楼主| 发表于 2019-12-18 08:07:07 | 显示全部楼层
自己顶一下,别沉帖了,万一高手来了看不到呢?
回复

使用道具 举报

8

主题

42

帖子

0

精华

初级会员

Rank: 2

积分
161
金钱
161
注册时间
2014-1-7
在线时间
24 小时
 楼主| 发表于 2019-12-18 08:08:32 | 显示全部楼层
自己顶一下,别沉帖了,万一高手来了看不到呢?
回复

使用道具 举报

8

主题

42

帖子

0

精华

初级会员

Rank: 2

积分
161
金钱
161
注册时间
2014-1-7
在线时间
24 小时
 楼主| 发表于 2019-12-18 08:08:59 | 显示全部楼层

自己顶一下,别沉帖了,万一高手来了看不到呢?
回复

使用道具 举报

10

主题

161

帖子

0

精华

高级会员

Rank: 4

积分
750
金钱
750
注册时间
2019-9-5
在线时间
134 小时
发表于 2019-12-18 09:13:16 | 显示全部楼层
float换成double,
回复

使用道具 举报

8

主题

42

帖子

0

精华

初级会员

Rank: 2

积分
161
金钱
161
注册时间
2014-1-7
在线时间
24 小时
 楼主| 发表于 2019-12-18 11:04:07 | 显示全部楼层

好的,我改一下试试!谢谢
回复

使用道具 举报

0

主题

2

帖子

0

精华

新手入门

积分
4
金钱
4
注册时间
2020-8-16
在线时间
0 小时
发表于 2020-8-16 19:06:14 | 显示全部楼层
楼主可以分享下源代码吗???感谢!!
回复

使用道具 举报

1

主题

9

帖子

0

精华

初级会员

Rank: 2

积分
50
金钱
50
注册时间
2020-7-7
在线时间
15 小时
发表于 2020-10-10 11:43:43 | 显示全部楼层
能分享一下F1的DSP库吗
回复

使用道具 举报

0

主题

2

帖子

0

精华

新手入门

积分
5
金钱
5
注册时间
2020-10-10
在线时间
2 小时
发表于 2020-10-13 08:18:25 | 显示全部楼层
想要dsp库
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则



关闭

原子哥极力推荐上一条 /2 下一条

正点原子公众号

QQ|手机版|OpenEdv-开源电子网 ( 粤ICP备12000418号-1 )

GMT+8, 2025-6-20 07:42

Powered by OpenEdv-开源电子网

© 2001-2030 OpenEdv-开源电子网

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