时域特征提取程序

matlabbbbbbbbbb 2019-07-15 08:43:37

y=xlsread('Xaxis.xlsx');
y1=zeros(size(y));
for i=1:73
N=length(y(:,i));
[c,l]=wavedec(y(:,i),4,'db4');
[thr1,nkeep]=wdcbm(c,l,4);
y1(:,i)=wdencmp('lvd',c,l,'db4',4,thr1,'s');
end
figure;
p6=zeros(1,73);
for i=1:73
N=length(y1(:,i));
p6(i)=sum(y1(:,i).^4)/N;
end
plot(p6);ylabel('峭度');
figure;
p1=zeros(1,73);
for i=1:73
N=length(y1(:,i));
p1(i)=mean(y1(:,i));
end
plot(p1);ylabel('均值');
figure;
p2=zeros(1,73);
for i=1:73
N=length(y1(:,i));
p2(i)=sqrt(sum(y1(:,i).^2)/N);
end
plot(p2);ylabel('均方根值');
figure;
p3=zeros(1,73);
for i=1:73
p3(i)=(sum(sqrt(abs(y1(:,i))))/N).^2;
end
plot(p3);ylabel('方根幅值');
figure;
p4=zeros(1,73);
for i=1:73
p4(i)=sum(abs(y1(:,i)))/N;
end
plot(p4);ylabel('绝对平均值');
figure;
p5=zeros(1,73);
for i=1:73
p5(i)=sum(y1(:,i).^3)/N;
end
plot(p5);ylabel('斜度');
figure;
p7=zeros(1,73);
for i=1:73
p7(i)=std(y1(:,i));
p7(i)=p7(i).^2;
end
plot(p7);ylabel('方差');
figure;
p8=zeros(1,73);
for i=1:73
p8(i)=max(y1(:,i));
end
plot(p8);ylabel('最大值');
figure;
p9=zeros(1,73);
for i=1:73
p9(i)=min(y1(:,i));
end
plot(p9);ylabel('最小值');
figure;
p10=p8-p9;
plot(p10);ylabel('峰峰值');
figure;
p11=zeros(1,73);
for i=1:73
p11(i)=max(abs(y1(:,i)));
end
plot(p11);ylabel('峰值');
figure;
f1=p2./abs(p1);
plot(f1);ylabel('波形指标');
figure;
f2=p11./p2;
plot(f2);ylabel('峰值指标');
figure;
f3=p11./abs(p1);
plot(f3);ylabel('脉冲指标');
figure;
f4=p11./p3;
plot(f4);ylabel('裕度指标');
figure;
f5=p6./((p2).^4);
plot(f5);ylabel('峰值指标');

...全文