48 lines
1.4 KiB
Mathematica
48 lines
1.4 KiB
Mathematica
|
% 从文件读取数据
|
|||
|
|
data = load('E:\Documents\工作资料\01项目资料\03.风电机组三合一采集器项目开发\01.机舱采集器\03.测试文件\ven\Data.txt'); % 假设数据存储在名为data.txt的文件中
|
|||
|
|
|
|||
|
|
% 提取采样率和持续时间
|
|||
|
|
fs = 60000; % 采样率为60KHz
|
|||
|
|
duration = 10; % 持续时间为10秒
|
|||
|
|
|
|||
|
|
% 创建时间向量
|
|||
|
|
t = linspace(0, duration, length(data));
|
|||
|
|
|
|||
|
|
% 绘制时域波形图
|
|||
|
|
figure;
|
|||
|
|
subplot(2,1,1);
|
|||
|
|
plot(t, data);
|
|||
|
|
title('时域波形');
|
|||
|
|
xlabel('时间（秒）');
|
|||
|
|
ylabel('幅值');
|
|||
|
|
ylim([min(data)-1, max(data)+1]); % 空出一定空间保证图形美观
|
|||
|
|
|
|||
|
|
% 计算并绘制频域谱
|
|||
|
|
N = length(data);
|
|||
|
|
f = (-N/2:N/2-1)*fs/N;
|
|||
|
|
f_data = fftshift(abs(fft(data))/N);
|
|||
|
|
|
|||
|
|
subplot(2,1,2);
|
|||
|
|
plot(f, f_data);
|
|||
|
|
title('频域谱');
|
|||
|
|
xlabel('频率（Hz）');
|
|||
|
|
ylabel('幅度');
|
|||
|
|
|
|||
|
|
% 设置频率图的x轴范围
|
|||
|
|
xlim([0, fs/2]);
|
|||
|
|
|
|||
|
|
% 计算指标
|
|||
|
|
rms_value = sqrt(mean(data.^2));
|
|||
|
|
peak_value = max(data);
|
|||
|
|
peak_to_peak = max(data) - min(data);
|
|||
|
|
kurtosis_value = kurtosis(data);
|
|||
|
|
crest_factor = max(data) / sqrt(mean(data.^2));
|
|||
|
|
skewness_value = skewness(data);
|
|||
|
|
|
|||
|
|
% 显示指标计算结果
|
|||
|
|
disp(['有效值（RMS）：', num2str(rms_value)]);
|
|||
|
|
disp(['峰值（Peak）：', num2str(peak_value)]);
|
|||
|
|
disp(['峰峰值（Peak-to-Peak）：', num2str(peak_to_peak)]);
|
|||
|
|
disp(['峭度（Kurtosis）：', num2str(kurtosis_value)]);
|
|||
|
|
disp(['峰值因数（Crest Factor）：', num2str(crest_factor)]);
|
|||
|
|
disp(['歪度（Skewness）：', num2str(skewness_value)]);
|