clc
clear
data = importdata('E:文件所在路径\文件名.csv');
timedata = data.data(:,1); % 时间列数据
data1 = data.data(:,2); % 通道数据更改加1
fs=8000;%采样频率
N=320000;%采样点数
dataSize = size(timedata); % 数据大小比如说320000个,如果数据不固定就用datasize
% 原始数据
figure(1)
plot(timedata, data1);
title('原始数据')
% 求FFT
ffty=fft(data1,N);
ffty_abs=abs(ffty);
ffty_abs=ffty_abs*2/N;
ffty_abs(1,1)=ffty_abs(1,1)/2;
n=1:1:N;
Fn=(n-1)*fs/N;
figure(2)
plot(Fn(1:N/2),ffty_abs(1:N/2));
title('FFT')
grid on
% 求RMS
rmsData = zeros(1, (dataSize(1) - 100)/100);%每100个求一次RMS所以是1到3199
j=1;
for i = 1:100:dataSize - 100
rmsData(j) = rms(data1(i:i+99));
j=j+1;
end
figure(3)
plot(rmsData);
title('RMS值')
关于读取EXCEL数据进行FFT+RMS处理
最新推荐文章于 2024-05-20 22:43:36 发布