当前位置:   article > 正文

关于matlab fft之后幅度矫正(搬运)_matlab信号幅相偏差矫正

matlab信号幅相偏差矫正

已经有网友解答在用matlab过程中使用FFT()函数求解信号频谱过程中,相应的幅值要做修正的问题,即参加FFT运算的点数为N时,最后幅值要除以N乘以2(直流分量时除以N即可),那么关于最后计算的结果为何比我们预期的要小,这个问题有两个方面原因,第一个原因是离散信号的点数,一般我们要求是2的幂指数例如L个,若不满足要求我们会补充零值,可以从能量守恒的角度考虑补充了信息点数使得信号长度增加,但因为是零值所以能量不变,但FFT()作用后我们频点数往往也是这个L,而且我们补充的零值是要参与其中计算的,所以这种能量不变却是频点数增加的做法,最后的结果比理论值小,这个是主要原因。第二个因素不是回答这个问题的,但往往初次使用matlab的FFT()函数会人为的设置参加的运算点数,这个点数设置的过大比如刚好比信号离散点数大了一倍,那么最后的频谱幅值刚好为理论的一半,特别要注意这个问题。
追问: 我知道不是的 是跟采样点数有个比例关系 但是这个比例关系也有误差啊
追答:假设原始信号的峰值为A,采样点数N,那么fft后的每个点(除了第一个点直流分量之外)的模值(幅值)就是A的N/2倍。而第一个点就是直流分量,它的模值就是直流分量的N倍。幅值的偏差要么是你点数与峰值相除小数点后边的精度问题,要么就是你的采样频率不够大,即采样的点数不够多,造成的精度问题,这问题很常见,通常都是加大采样频率,采样的点数越多越精确,千万别相信采样定理,只将采样频率增大一倍,实际中这样根本没法满足需求的精度,搞个5、6倍还差不多,越多越好,不过计算的时间可就增大了,所以你可以试试,满足你精度要求的采样频率到底多少合适。
我用3sin(300pit);采样t=0:0.001:1;%设定信号时间为0到1秒,频率为1000HZ来计算,发现仍有许多t点的值,在峰值仍然达不到幅值3,所以这么多幅值累计的误差很可观,通过N/2倍计算幅值也是2.8左右,都是离散取点不能完全描述连续信号的问题,所以误差是难免的,除了增大点数别无他法。

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/神奇cpp/article/detail/949013
推荐阅读
相关标签
  

闽ICP备14008679号