完善资料让更多小伙伴认识你,还能领取20积分哦, 立即完善>
各位高手,现在用STM32F0的ADC采样一组不规则信号,128个点,200HZ采样;采样数据经过MTLAB验证,FFT转换结果跟keil运行输出结果一致,说明FFT程序正确。现在经过FFT转换输出的128个模数(x.real=sqrt(x.real*x.real+x.img*x.img);),取了64个数据转换功率来分析(根据FFT对称性),把64个功率数据由串口输出,从串口数据来看,在不规则信号幅度低的时候,输出数据如下:(35 10 15 18 20 0F 00 0C 19 1F 15 12 10 0F 06 0C 09 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00);
在不规则信号幅度高的时候,输出数据如下:(35 1C 25 19 29 23 19 16 1B 26 2B 28 2D 29 20 28 13 22 24 18 21 10 21 1B 1C 24 1A 19 1F 13 1A 17 17 18 18 18 1A 14 13 15 17 1A 19 1B 14 13 16 15 16 0D 09 10 14 15 0C 1A 17 10 13 0C 09 06 09); 观察数据发现:输入信号幅度低时,在64个点对应频率0-25HZ以后,各个频点对应功率全部为0;而信号强的输出数据,64个频点对应的功率都大于0。 对于这64个数据,是不是输入信号幅度越高,FFT转换后的对于功率就越大?要分析原始波形的幅度,是不是看这64点功率就可以一一对应判断了?对这64个数据还需要进行什么分析才知道原始波形高点的频率和幅度? //5ms采1个点,128个点,采样频率F=200Hz //FFT结果的128个数据对应频率点是:0,1*200Hz/128,2*200Hz/128,3*200Hz/128, 。。。 ,127*200Hz/128 (0->127) //N=128 //频率点表: // 0.00-> 1.5625 -> 3.125- > 4.6875 -> 6.25 -> 7.8125 -> 9.375 -> 10.9375-> //8PCS // 12.5 -> 14.0625-> 15.625-> 17.1875-> 18.75-> 20.3125-> 21.875-> 23.4375-> //8pcs // 25.00-> 25.5626-> 28.125-> 29.6875-> 31.25-> 32.8125-> 34.375-> 35.9375-> //8PCS // 37.5 -> 39.0625-> 40.625-> 42.1875-> 43.75-> 45.3125-> 46.875-> 48.4375-> //8pcs |
|
相关推荐
8个回答
|
|
请问高频信号,经FFT后,得到的模值和频点,怎么处理才知道原始信号的相关参数?
|
|
|
|
% MTLAB 显示出ADC采样128个点对应波形
%--------------------------------------------------------------------------------- clear; close all; %128点ADC采样最大幅度波形 fs=200; %采样频率 N=128; %采样数据点数 n=0:N-1; %采样序列(0->128) t=n/fs; %时间序列(0->128/200) %--------------------------------------------------------------------------------- xn=[729 729 728 732 733 727 729 727 727 728 730 730 76 2 1 1 729 729 730 727 726 726 727 728 727 708 ... 319 3 1 403 730 597 0 1 2 1 18 1 1 0 0 0 454 727 729 733 734 728 729 723 727 732 643 3 2 3 0 0 2 564 ... 729 726 729 727 729 730 731 730 730 730 729 645 2 0 0 0 0 0 1 0 0 4 0 1 1 1 2 4 730 729 730 730 728 ... 728 727 729 724 726 725 728 726 727 726 728 732 729 724 727 728 730 729 723 726 728 729 727 729 728 ... 729 728 727 732 731 281 1 1 1 0]; %--------------------------------------------------------------------------------- subplot(2,2,1),plot(t,xn); %绘出原始信号图 xlabel('时间/sec'); % x轴为时间 title('最大幅度原始信号'); %标题:最大幅度原始信号图 y=fft(xn); %对输入信号做 FFT 变换 mag=abs(y); %求 FFT 转换结果的模值,即振幅 Ayy = mag/(N/2); %模换算成实际的幅度=该点复数的模值除以N/2 f=(0:N-1)*fs/N; %真实频率 subplot(2,2,3),plot(f,Ayy); %绘出FFT转换后实际的幅度图 xlabel('频率/Hz'); % x轴为频率 title('FFT转换后的实际幅度'); %标题:FFT转换后的实际幅度图 %--------------------------------------------------------------------------------- %%%%%%%% 128点ADC采样小型幅度波形 yn=[449 449 443 444 443 444 443 444 444 434 426 420 420 426 435 443 453 473 491 508 516 526 540 558 574 ... 577 574 571 574 582 588 605 613 611 604 587 569 548 527 507 472 429 388 351 312 284 276 283 298 318 332 ... 345 362 383 410 446 471 484 484 486 501 525 547 564 573 576 573 570 576 588 600 608 598 589 565 544 534 ... 537 562 583 595 587 558 533 509 493 482 482 481 487 493 508 522 541 549 547 543 545 550 551 554 554 565 ... 579 597 602 604 599 594 583 568 545 530 514 501 503 512 524 532 541 536 537 538 537 544 545 538 534]; %--------------------------------------------------------------------------------- subplot(2,2,2),plot(t,yn); %绘出原始信号图 xlabel('时间/sec'); % x轴为时间 title('小型幅度原始信号'); %标题:小型幅度原始信号图 z=fft(yn); %对输入信号做 FFT 变换 maf=abs(z); %求 FFT 转换结果的模值,即振幅 Byy = maf/(N/2); %模换算成实际的幅度=该点复数的模值除以N/2 f=(0:N-1)*fs/N; %真实频率 subplot(2,2,4),plot(f,Byy); %绘出FFT转换后实际的幅度图 xlabel('频率/Hz'); % x轴为频率 title('FFT转换后的实际幅度'); %标题:FFT转换后的实际幅度图 |
|
|
|
请各位不吝赐教,在此感谢
|
|
|
|
不知道楼主有没有对采样数据*2呢?
因为,采样频率要大于信号频率的两倍,才能恢复出实际波形。 以下提供圈圈大神的详细解释: http://www.eeworld.com.cn/Test_a ... /article_12114.html |
|
|
|
本帖最后由 无薪税绵 于 2017-10-25 11:42 编辑
不知道楼主有没有对采样数据*2呢? 因为,采样频率要大于信号频率的两倍,才能恢复出实际波形。 以下提供圈圈大神的详细解释: http://www.eeworld.com.cn/Test_a ... /article_12114.html 再提供一份:FFT后的物理意义 http://www.xuebuyuan.com/539160.html 刚发现,这个解释才牛: https://zhuanlan.zhihu.com/wille/19763358 |
|
|
|
请问楼主是用DSP库的吗?STM32F0的DSP库在哪里可以下载?
|
|
|
|
楼主讲得详细
|
|
|
|
|
|
|
|
你正在撰写答案
如果你是对答案或其他答案精选点评或询问,请使用“评论”功能。
STM32进行DMA数据转运,为什么只有A组最后一位数据被转运到B组去?
1148 浏览 2 评论
1497 浏览 3 评论
1129 浏览 0 评论
嵌入式开发 Win11安装ST-Link 提示 设备描述符请求失败 怎么处理?
1266 浏览 0 评论
1553 浏览 1 评论
小黑屋| 手机版| Archiver| 电子发烧友 ( 湘ICP备2023018690号 )
GMT+8, 2024-8-16 17:00 , Processed in 1.287677 second(s), Total 87, Slave 71 queries .
Powered by 电子发烧友网
© 2015 www.ws-dc.com
关注我们的微信
下载发烧友APP
电子发烧友观察
版权所有 © 湖南华秋数字科技有限公司
电子发烧友 (电路图) 湘公网安备 43011202000918 号 电信与信息服务业务经营许可证:合字B2-20210191 工商网监 湘ICP备2023018690号