采样的信号分成多段然后分别FFT求平均值

gumengyinyin · 发表于 2014-7-12 15:43

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

Fs=100;    %采样频率
N=1024;
Y=reshape(y,1024,100);
Yf=0;
for i=1:100
Yf=Yf+fft(Y(:,i));
end
Yf=Yf/100;

f=(0:N-1)*Fs/N;
mag=2*abs(Yf)/N; %幅值
% Pyy=mag.^2;  %功率
subplot(2,1,2);
plot(f(1:N/2),mag(1:N/2))  %频率-幅值图
% plot(f(1:N/2),Pyy(1:N/2))  %功率谱图
xlim([0 30]);
xlabel('f/Hz');
ylabel('幅值/m/s^2');
时域采样了102400个点
现在对其进行频谱分析，把这些点分成100组，对每组求傅里叶变换，求和，发现幅值的大小变量，是数量级上的改变，请问是什么原因呢

westrongmc · 发表于 2014-7-12 18:22

本帖最后由 westrongmc 于 2014-7-12 22:30 编辑

1. help fft
2. 非同步触信号的线性谱不可以做平均处理
3. 对自功率谱进行平均

补充内容 (2014-7-13 10:19):
2. 非同步触发信号的线性谱不可以做平均处理。

lbtv · 发表于 2014-7-13 10:02

复数相加：会不会影响呢！关键是你最好取其实部进行平均
下面是用你的程序对某组数据处理结果数据被分成20段
第一种情况（用你原来的程序）

登录/注册后可看大图

第二种情况（取其模进行平均操作）

登录/注册后可看大图

这可能就是你的结果形成原因！

lbtv · 发表于 2014-7-13 10:03

但是我不知道为什么最高幅值会稍稍变大！

lbtv · 发表于 2014-7-13 10:07

对于我使用的数据相关参数，可以查看
http://forum.vibunion.com/forum. ... amp;authorid=225230

gumengyinyin · 发表于 2014-7-14 16:23

lbtv 发表于 2014-7-13 10:02
复数相加：会不会影响呢！关键是你最好取其实部进行平均
下面是用你的程序对某组数据处理结果数据被分成20 ...

恩恩，我后来是对模进行平均计算的，幅值倒是变化不大，不过对数据分的段数不一样，频率会稍微偏移

账号		自动登录	找回密码
密码			我要加入

[FFT] 采样的信号分成多段然后分别FFT求平均值

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。