删除高频分量matlab,关于FFT之后怎么提出该频率分量

本帖最后由 领养的羚羊 于 2015-5-27 17:10 编辑

就是一个很简单的信号，图2是时域图1是频域。然后在网上找到这样的话：

假设FFT之后某点n用复数a+bi表示，那么这个复数的模就是An=根号a*a+b*b，相位就是Pn=atan2(b,a)。

根据以上的结果，就可以计算出n点(n≠1，且n<=N/2)对应的信号的表达式为：

An/(N/2)*cos(2*pi*Fn*t+Pn)，即2*An/N*cos(2*pi*Fn*t+Pn)。对于n=1点的信号，是直流分量，幅度即为A1/N。

然后按照这段话，那我的信号中两个频域峰值所对应的信号分量应该分别是：

2*1065.8*cos(2*3.1416*311*t-1.1216)/2200;

2*1065.8*cos(2*3.1416*1891*t-1.1216)/2200;

(幅值1065.8，两个点的频率分别是311和1891)，根据公式计算出的相位是1.1216，单位是弧度，2200是点数。然后我发现不管是这两个信号中的单独任意一个还是两个信号之和都和源信号相差甚远。会是什么原因呢？

另外，t是我定义的数列，1、2、3、4………………2200。另外在想是不是和FFT变换时给的采样率fs有关系？

程序：

t=(1:2210);

x3=sin(80*t);

%以下为对信号进行整周期截取%

for i=50:2:length(x3)

j=i*0.5;

c13=x3(