#基于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);