基于matlab的声发射信号量化诊断系统

[X,textdata]=xlsread('C:\Users\Administrator\Desktop\shuju2.xls'); %读取信号特征向量矩阵 x=zscore(X); %数据标准化 [idx,C,sumd,D]=kmeans(x,4,'Replicates',5); Cluster1=X(idx==1,:); %第一类信号 Cluster2=X(idx==2,:); %第二类信号 Cluster3=X(idx==3,:); %第三类信号 Cluster4=X(idx==4,:); %第四类信号 figure(1); %分类后幅值与峰频分布图 plot(Cluster1(:,1),Cluster1(:,2), 'rs','Markersize',12); hold on plot(Cluster2(:,1),Cluster2(:,2), 'b^','Markersize',12); plot(Cluster3(:,1),Cluster3(:,2), 'cp','Markersize',12); plot(Cluster4(:,1),Cluster4(:,2), 'bo','Markersize',12); axis([35 90 0 400]); xlabel('Amplitude,dB'); ylabel('Peak frequency,kHz'); figure(2); %分类后持续时间与峰频 plot(Cluster1(:,3),Cluster1(:,2), 'rs', 'Markersize',12); %（红色正方形） hold on plot(Cluster2(:,3),Cluster2(:,2), 'b^', 'Markersize',12); % (蓝色三角形标注) plot(Cluster3(:,3),Cluster3(:,2), 'cp', 'Markersize',12); % (墨绿色五角星标注) plot(Cluster4(:,3),Cluster4(:,2), 'bo', 'Markersize',12); % (蓝色圆标注) axis([0 4000 0 400]); xlabel('Duration, μs'); ylabel('Peak frequency, kHz'); figure(3); %分类后持续时间与幅值 plot(Cluster1(:,3),Cluster1(:,1), 'rs', 'Markersize',12); hold on plot(Cluster2(:,3),Cluster2(:,1), 'b^', 'Markersize',12); plot(Cluster3(:,3),Cluster3(:,1), 'cp', 'Markersize',12); plot(Cluster4(:,3),Cluster4(:,1), 'bo', 'Markersize',12); axis([0 4000 0 100]); xlabel('Duration, μs'); ylabel(' Amplitude, dB ');