matlab用fft还原真实幅值,matlab仿真FFT结果幅值比实际的略有降低?

该博客通过MATLAB处理Excel文件中的加速度时域信号,使用FFT进行频率分析。通过对数据分块、加窗处理并计算平均幅值,得出幅度-频率曲线,最后将结果写回Excel文件。在过程中注意到FFT结果的幅值可能低于实际值。
摘要由CSDN通过智能技术生成

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

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值