Matlab的时域数据的FFT计算及频域平均_matlab怎么频域多次平均

#基于Matlab的FFT计算
针对采集到的较长时间时域信号，往往需要对整个时域进行FFT计算，本代码例以随机振动数据FFT处理计算为例，通过加hann窗，计算每个时间块内的数据计算FFT，时间块设置50%重叠，最后对所有块内的FFT数据进行overall平均模式计算，最终获得整个时域内的FFT幅频特性。

x=data; %输入时域数据
n=length(x);
fs=2000; %时域信号的采样频率
L=fs*1.024; %设置FFT块的数据长度
nb=fix(n/(L/2));
l=1:L/2:nb*(L/2)+1;
windows=hann(L); %加hann窗，随机振动常用
f=(0:(L-1)/2)/L*fs; %根据对称性取半，获取各频率点
FFTA=zeros(L/2,length(l)-2);
for i=1:1:length(l)-2  %针对每个时间块内的数据进行FFT计算
    x1=x(l(i):l(i+2)-1).*windows; 
    y=fft(x1,L);
    P2 = abs(y/L);
    P1 = P2(1:L/2);
    P1(2:end-1) = 2*P1(2:end-1);
    P1=P1*2;   %加hann窗的幅值恢复系数为2
    FFTA(:,i)=P1;
end
FOA=rms(FFTA,2); %所有FFT块同一频率下幅值计算有效值
plot(f,FOA);
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20