% 从文件读取数据
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)]);