Matalb读取TDMS文件并提取13个频谱特征

 提取了十三个特征

clc;clear ;close all
nor=convertTDMS(0,'Tr.tdms');%读取文件
data1=nor.Data.MeasuredData(3).Data;    %选择几通道
   
Fs = 2560;%采样频率
N=length(data1);
data1 = abs(data1);
%频谱均值
F1 = mean(data1);
%谱幅值标准差
F2 = std(data1);
%谱幅值偏度
F3 = sum((data1-F1).^3)/((N-1)*F2^3);
%谱幅值峭度
F4 = F3 = sum((data1-F1).^4)/((N-1)*F2^4);
%计算重心频率
f = Fs* [1:N ]/N;
py = 2 * data1 /N;
% FC=sum(f.*py)/sum(data1)
F = sum(data1);
q = sum(f.*py);
F5=sum(q)/F;

%谱频率均方根值
FZ=sum(f.^2.*py);
F6 = sqrt(sum(FZ)/F);

%谱根4/2矩比
Fj=su