cmsis 的fft库,实现中的一些疑问.

alca_bigV 2017-08-28 09:46:36

1)在matlab中用如下命令生成一个正弦波的采样点,512个点.t=0:1/6656:511/6656;%生成采样时间序列
s1=sin(54*2*pi*t);%生成采样点
fftsa1=fft(s1,512);%fft变换
2)得到的图形如下:结果正确

3)命令行输出元素如下,结果正确.

4)把s1导入板上,用fft算法(从arm 的 cmsis库里得来)进行运算.输出一个1024的数据结构,..生成的数据如下:


5)我就不理解 了,第四步输出的数据,应该和第3步对应的数据 一致,即第4个元素的实部和虚部应该是 1 和 0,但这里,所有的512个元素的实部与虚部都完全与matlab中的fft转换后得到的数据不一样,。。
我开始怀疑是我的程序有问题,用另一套fft算法进行运算,得出的结果虽然和4看到的结果不同,但进行模计算,再进行归一化处理后,结果完全 一致(全部归一化后,得到一个512个float型 的数组,全部是1),..
所以肯定是我哪里理解 有问题,...在arm上运算后得到的数据 不应该是和matlab第三步显示的 结果一样?第三步的结果,除了第四个元素为1 ,0,以外,其他 都 为零,完全 能理解 ,因为除了第四点对应的频率 外,其他 频率 都是0向量.
(感觉用了fft运算后得到的结果是:一个圆,分成128份,从0度转到360度,计算出了每一时刻所对应的半径的x和y坐标,这样计算模值后,就全部 为1)..
如何理解,才能和matlab中结果 3 对应上?(5)我就不理解 了,第四步输出的数据,应该和第3步对应的数据 一致,即第4个元素的实部和虚部应该是 1 和 0,但这里,所有的512个元素的实部与虚部都完全与matlab中的fft转换后得到的数据不一样,。。
我开始怀疑是我的程序有问题,用另一套fft算法进行运算,得出的结果虽然和4看到的结果不同,但进行模计算,再进行归一化处理后,结果完全 一致(全部归一化后,得到一个512个float型 的数组,全部是1),..
所以肯定是我哪里理解 有问题,...在arm上运算后得到的数据 不应该是和matlab第三步显示的 结果一样?第三步的结果,除了第四个元素为1 ,0,以外,其他 都 为零,完全 能理解 ,因为除了第四点对应的频率 外,其他 频率 都是0向量.
(感觉用了fft运算后得到的结果是:一个圆,分成128份,从0度转到360度,计算出了每一时刻所对应的半径的x和y坐标,这样计算模值后,就全部 为1)..
如何理解,才能和matlab中结果 3 对应上?(偏置4所对应的实部与虚部,应该有数值,其他都为零)?
...全文
716 6 打赏 收藏 转发到动态 举报
写回复
用AI写文章
6 条回复
切换为时间正序
请发表友善的回复…
发表回复
Simon_Lancaster 2018-07-30
  • 打赏
  • 举报
回复
你碰上一装13的,此人不懂MATLAB与嵌入式系统!
赵4老师 2017-08-30
  • 打赏
  • 举报
回复
理解讨论之前请先学会如何观察! 调试的过程就是观察。
alca_bigV 2017-08-30
  • 打赏
  • 举报
回复
引用 3 楼 zhao4zhong1 的回复:
代码功能归根结底不是别人帮自己看或讲解或注释出来的;而是被自己静下心来花足够长的时间和精力亲自动手单步或设断点或对执行到某步获得的中间结果显示或写到日志文件中一步一步分析出来的。 提醒:再牛×的老师也无法代替学生自己领悟和上厕所! 单步调试和设断点调试是程序员必须掌握的技能之一。
这话说的太对了...看得出,你应该很牛b.. 不过,我有问代码功能 吗?或者,你有看我说的东西了吗? 我好像只是 说,我的理解思路有问题,想请人指点一下哪一个理解的地方出错好不? 上面是有很多代码和数据,但那些对于搞过这些的,一眼扫过都可以,请教人总要把步骤说清楚,别人才能理解 我的思路 哪里出错,不是? 从matlab和cortex里面的代码调试,我好像一个字没提 需要如何调试... 我只是告诉了别人,结果是这样,过程是这样,而已.
赵4老师 2017-08-29
  • 打赏
  • 举报
回复
代码功能归根结底不是别人帮自己看或讲解或注释出来的;而是被自己静下心来花足够长的时间和精力亲自动手单步或设断点或对执行到某步获得的中间结果显示或写到日志文件中一步一步分析出来的。 提醒:再牛×的老师也无法代替学生自己领悟和上厕所! 单步调试和设断点调试是程序员必须掌握的技能之一。
alca_bigV 2017-08-28
  • 打赏
  • 举报
回复
引用 1 楼 zhao4zhong1 的回复:
不要做A语言代码修改为B语言代码的无用功。 也不要做用A语言代码直接调用B语言代码库这样复杂、这样容易出错的傻事。 只需让A、B语言代码的输入输出重定向到文本文件,或修改A、B语言代码让其通过文本文件输入输出。 即可很方便地让A、B两种语言之间协调工作。 比如: A将请求数据写到文件a.txt,写完后改名为aa.txt B发现aa.txt存在时,读取其内容,调用相应功能,将结果写到文件b.txt,写完后删除aa.txt,再将b.txt改名为bb.txt A发现bb.txt存在时,读取其内容,读完后删除bb.txt 以上A可以替换为任何一种开发语言或开发环境,B可以替换为任何一种与A不同的开发语言或开发环境。 除非A或B不支持判断文件是否存在、文件读写和文件更名。 但是谁又能举出不支持判断文件是否存在、文件读写和文件更名的开发语言或开发环境呢? 可以将临时文件放在RamDisk上提高效率减少磨损磁盘。 数据的结构很复杂的话,文本文件的格式问题可参考json或xml 共享临时文本文件这种进程之间的通讯方法相比其它方法的优点有很多,下面仅列出我现在能想到的: ·进程之间松耦合 ·进程可在同一台机器上,也可跨机,跨操作系统,跨硬件平台,甚至跨国。 ·方便调试和监视,只需让第三方或人工查看该临时文本文件即可。 ·方便在线开关服务,只需删除或创建该临时文本文件即可。 ·方便实现分布式和负载均衡。 ·方便队列化提供服务,而且几乎不可能发生队列满的情况(除非硬盘空间满) ·…… “跨语言、跨机,跨操作系统,跨硬件平台,跨国,跨*.*的”苦海无边, 回头是“使用共享纯文本文件进行信息交流”的岸!
不好意思,不过,我怎么感觉 你有点 答非所问?还是你讲的理论太高深,我没理解?得仔细思考思考.
赵4老师 2017-08-28
  • 打赏
  • 举报
回复
不要做A语言代码修改为B语言代码的无用功。 也不要做用A语言代码直接调用B语言代码库这样复杂、这样容易出错的傻事。 只需让A、B语言代码的输入输出重定向到文本文件,或修改A、B语言代码让其通过文本文件输入输出。 即可很方便地让A、B两种语言之间协调工作。 比如: A将请求数据写到文件a.txt,写完后改名为aa.txt B发现aa.txt存在时,读取其内容,调用相应功能,将结果写到文件b.txt,写完后删除aa.txt,再将b.txt改名为bb.txt A发现bb.txt存在时,读取其内容,读完后删除bb.txt 以上A可以替换为任何一种开发语言或开发环境,B可以替换为任何一种与A不同的开发语言或开发环境。 除非A或B不支持判断文件是否存在、文件读写和文件更名。 但是谁又能举出不支持判断文件是否存在、文件读写和文件更名的开发语言或开发环境呢? 可以将临时文件放在RamDisk上提高效率减少磨损磁盘。 数据的结构很复杂的话,文本文件的格式问题可参考json或xml 共享临时文本文件这种进程之间的通讯方法相比其它方法的优点有很多,下面仅列出我现在能想到的: ·进程之间松耦合 ·进程可在同一台机器上,也可跨机,跨操作系统,跨硬件平台,甚至跨国。 ·方便调试和监视,只需让第三方或人工查看该临时文本文件即可。 ·方便在线开关服务,只需删除或创建该临时文本文件即可。 ·方便实现分布式和负载均衡。 ·方便队列化提供服务,而且几乎不可能发生队列满的情况(除非硬盘空间满) ·…… “跨语言、跨机,跨操作系统,跨硬件平台,跨国,跨*.*的”苦海无边, 回头是“使用共享纯文本文件进行信息交流”的岸!

5,530

社区成员

发帖
与我相关
我的任务
社区描述
C/C++ 模式及实现
社区管理员
  • 模式及实现社区
加入社区
  • 近7日
  • 近30日
  • 至今
社区公告
暂无公告

试试用AI创作助手写篇文章吧