% 假设您已经有了时域数据,存储在名为"signal"的向量中
data=load('C:\Users\jiema\Desktop\ys2.txt');
% data=load('C:\Users\jiema\Downloads\运输2.txt');
x=data(:,1);
y=data(:,2);%读取Y轴数据
signal =y
rms1=rms(y)
% 假设您已经有了时域数据,存储在名为"signal"的向量中
% 假设您已有一个频域信号的幅度谱数据
f = fft(signal);
frequency_spectrum = abs(f);
pha = angle(f);
time_domain_signal2 = ifft(frequency_spectrum);
rms2=rms(time_domain_signal2)
%将1hz以下幅值改为0 1Hz/800Fs*点数,傅里叶变换对称,所以尾巴也变成0
frequency_spectrum(1:1021) = 0
frequency_spectrum(815774-1021:815774) = 0
% 将频域信号转换为时域信号
% f2 = frequency_spectrum.*exp(1j*pha); % 根据幅值谱和相位谱重建频域信号
% time_domain_signal = ifft(f2);
time_domain_signal = ifft(frequency_spectrum);
rms3=rms(time_domain_signal)
% 计算时域信号的平方
squared_signal = abs(time_domain_signal).^2;
% 计算信号的平均能量
average_energy = mean(squared_signal);
% 计算RMS值
rms_value = sqrt(average_energy);
disp(['RMS:', num2str(rms_value)]);
disp(['RMS1 原数据:', num2str(rms1)]);
disp(['RMS2 全频:', num2str(rms2)]);
disp(['RMS >1Hz频段值:', num2str(rms3)]);
用matlabl直接滤波,输出大于1Hz频带RMS
最新推荐文章于 2024-03-10 15:22:09 发布
![](https://img-home.csdnimg.cn/images/20240711042549.png)