clc;clear;
filename='TestAcceleration.xlsx'; % Excel数据文件名称
sheet=1; % 数据读写工作表worksheet编号
ReadRange='A4:B40003'; % 读取时域信号数据范围
WriteRange='E4'; % 写入幅频特性数据起始位置
Num=xlsread(filename,sheet,ReadRange);
[m n]=size(Num); % 获取数据维数,主要是获取数据行数
time=Num(:,1);
Acceleration=Num(:,2); % 要进行FFT运算的数据所在列数
plot(time(1:6000),Acceleration(1:6000))
xlabel('时间 [s]');
ylabel('加速度 [m/s2]');
title('加速度时域波形');
Fs=4e4; % 采样频率
NFFT=4096; % 每个分块数据个数
N_Frame=floor(2*m/NFFT-1); % 根据数据长度计算的分块数目,分块间按50%重叠率
temp_Amplitude_Acceleration=zeros(N_Frame,NFFT); % 每一行用来存储每个分块FFT变换后的幅值
w_hann=hann(NFFT,'periodic'); % 汉宁窗,是长度为NFFT的列向量
for i=1:N_Frame
temp=Accele