• formula.
    了解作者
  • WINDOWS
    开发工具
  • 7KB
    文件大小
  • zip
    文件格式
  • 0
    收藏次数
  • 1 积分
    下载积分
  • 10
    下载次数
  • 2017-11-06 09:02
    上传日期
距离评估算法,用于高维特征集的特征选择,得到更好的分类模型
distance.zip
  • distance
  • distanceevaluation.m
    540B
  • cdet_test.m
    35.4KB
  • indis.m
    155B
内容介绍
clc;clear load('normalz.mat') load('normal30xy.mat') load('in30z.mat') load('In30xy.mat') load('combine30z.mat') load('combine30xy.mat') load('ball30z.mat') load('Ball30xy.mat') load('out30xy.mat') load('out30z.mat') b1=ball30xy(:,3);b2=ball30xy(:,5);b3=ballz(:,3); b5=ball30xy(:,2);b6=ball30xy(:,4);b7=ballz(:,2); c1=comxy(:,3);c2=comxy(:,5);c3=comz(:,3); c5=comxy(:,2);c6=comxy(:,4);c7=comz(:,2); i1=inxy(:,3);i2=inxy(:,5);i3=inz(:,3); i5=inxy(:,2);i6=inxy(:,4);i7=inz(:,2); n1=normalxy(:,3);n2=normalxy(:,5);n3=normalz(:,3); n5=normalxy(:,2);n6=normalxy(:,4);n7=normalz(:,2); o1=outxy(:,3);o2=outxy(:,5);o3=outz(:,3); o5=outxy(:,2);o6=outxy(:,4);o7=outz(:,2); n=3000; tr=200;te=3/4*tr; for i=1:tr num=((i-1)*n+1):(n*i); B1(:,i)=b1(num); B2(:,i)=b2(num); B3(:,i)=b3(num); B5(:,i)=b5(num); B6(:,i)=b6(num); B7(:,i)=b7(num); C1(:,i)=c1(num); C2(:,i)=c2(num); C3(:,i)=c3(num); C5(:,i)=c5(num); C6(:,i)=c6(num); C7(:,i)=c7(num); I1(:,i)=i1(num); I2(:,i)=i2(num); I3(:,i)=i3(num); I5(:,i)=i5(num); I6(:,i)=i6(num); I7(:,i)=i7(num); N1(:,i)=n1(num); N2(:,i)=n2(num); N3(:,i)=n3(num); N5(:,i)=n5(num); N6(:,i)=n6(num); N7(:,i)=n7(num); O1(:,i)=o1(num); O2(:,i)=o2(num); O3(:,i)=o3(num); O5(:,i)=o5(num); O6(:,i)=o6(num); O7(:,i)=o7(num); end %RMS均方根值 B1rms=rms(B1);B2rms=rms(B2);B3rms=rms(B3); B5rms=rms(B5);B6rms=rms(B6);B7rms=rms(B7); C1rms=rms(C1);C2rms=rms(C2);C3rms=rms(C3); C5rms=rms(C5);C6rms=rms(C6);C7rms=rms(C7); I1rms=rms(I1);I2rms=rms(I2);I3rms=rms(I3); I5rms=rms(I5);I6rms=rms(I6);I7rms=rms(I7); N1rms=rms(N1);N2rms=rms(N2);N3rms=rms(N3); N5rms=rms(N5);N6rms=rms(N6);N7rms=rms(N7); O1rms=rms(O1);O2rms=rms(O2);O3rms=rms(O3); O5rms=rms(O5);O6rms=rms(O6);O7rms=rms(O7); %平均值 B1me=mean(B1);C1me=mean(C1);I1me=mean(I1);N1me=mean(N1);O1me=mean(O1); B2me=mean(B2);C2me=mean(C2);I2me=mean(I2);N2me=mean(N2);O2me=mean(O2); B3me=mean(B3);C3me=mean(C3);I3me=mean(I3);N3me=mean(N3);O3me=mean(O3); B6me=mean(B6);C6me=mean(C6);I6me=mean(I6);N6me=mean(N6);O6me=mean(O6); B7me=mean(B7);C7me=mean(C7);I7me=mean(I7);N7me=mean(N7);O7me=mean(O7); B5me=mean(B5);C5me=mean(C5);I5me=mean(I5);N5me=mean(N5);O5me=mean(O5); %标准差 B1s=std(B1).^2;C1s=std(C1).^2;I1s=std(I1).^2;N1s=std(N1).^2;O1s=std(O1).^2; B2s=std(B2).^2;C2s=std(C2).^2;I2s=std(I2).^2;N2s=std(N2).^2;O2s=std(O2).^2; B3s=std(B3).^2;C3s=std(C3).^2;I3s=std(I3).^2;N3s=std(N3).^2;O3s=std(O3).^2; B6s=std(B6).^2;C6s=std(C6).^2;I6s=std(I6).^2;N6s=std(N6).^2;O6s=std(O6).^2; B7s=std(B7).^2;C7s=std(C7).^2;I7s=std(I7).^2;N7s=std(N7).^2;O7s=std(O7).^2; B5s=std(B5).^2;C5s=std(C5).^2;I5s=std(I5).^2;N5s=std(N5).^2;O5s=std(O5).^2; %波形指标 B1shape=(n*B1rms)./sum(abs(B1));C1shape=(n*C1rms)./sum(abs(C1));I1shape=(n*I1rms)./sum(abs(I1));N1shape=(n*N1rms)./sum(abs(N1));O1shape=(n*O1rms)./sum(abs(O1)); B2shape=(n*B2rms)./sum(abs(B2));C2shape=(n*C2rms)./sum(abs(C2));I2shape=(n*I2rms)./sum(abs(I2));N2shape=(n*N2rms)./sum(abs(N2));O2shape=(n*O2rms)./sum(abs(O2)); B3shape=(n*B3rms)./sum(abs(B3));C3shape=(n*C3rms)./sum(abs(C3));I3shape=(n*I3rms)./sum(abs(I3));N3shape=(n*N3rms)./sum(abs(N3));O3shape=(n*O3rms)./sum(abs(O3)); B6shape=(n*B6rms)./sum(abs(B6));C6shape=(n*C6rms)./sum(abs(C6));I6shape=(n*I6rms)./sum(abs(I6));N6shape=(n*N6rms)./sum(abs(N6));O6shape=(n*O6rms)./sum(abs(O6)); B7shape=(n*B7rms)./sum(abs(B7));C7shape=(n*C7rms)./sum(abs(C7));I7shape=(n*I7rms)./sum(abs(I7));N7shape=(n*N7rms)./sum(abs(N7));O7shape=(n*O7rms)./sum(abs(O7)); B5shape=(n*B5rms)./sum(abs(B5));C5shape=(n*C5rms)./sum(abs(C5));I5shape=(n*I5rms)./sum(abs(I5));N5shape=(n*N5rms)./sum(abs(N5));O5shape=(n*O5rms)./sum(abs(O5)); %脉冲指标 B1impluse=max(abs(B1),[],1)./sqrt(sum(abs(B1)));C1impluse=max(abs(C1),[],1)./sqrt(sum(abs(C1)));I1impluse=max(abs(I1),[],1)./sqrt(sum(abs(I1))); N1impluse=max(abs(N1),[],1)./sqrt(sum(abs(N1)));O1impluse=max(abs(O1),[],1)./sqrt(sum(abs(O1))); B5impluse=max(abs(B5),[],1)./sqrt(sum(abs(B5)));C5impluse=max(abs(C5),[],1)./sqrt(sum(abs(C5)));I5impluse=max(abs(I5),[],1)./sqrt(sum(abs(I5))); N5impluse=max(abs(N5),[],1)./sqrt(sum(abs(N5)));O5impluse=max(abs(O5),[],1)./sqrt(sum(abs(O5))); B2impluse=max(abs(B2),[],1)./sqrt(sum(abs(B2)));C2impluse=max(abs(C2),[],1)./sqrt(sum(abs(C2)));I2impluse=max(abs(I2),[],1)./sqrt(sum(abs(I2))); N2impluse=max(abs(N2),[],1)./sqrt(sum(abs(N2)));O2impluse=max(abs(O2),[],1)./sqrt(sum(abs(O2))); B3impluse=max(abs(B3),[],1)./sqrt(sum(abs(B3)));C3impluse=max(abs(C3),[],1)./sqrt(sum(abs(C3)));I3impluse=max(abs(I3),[],1)./sqrt(sum(abs(I3))); N3impluse=max(abs(N3),[],1)./sqrt(sum(abs(N3)));O3impluse=max(abs(O3),[],1)./sqrt(sum(abs(O3))); B6impluse=max(abs(B6),[],1)./sqrt(sum(abs(B6)));C6impluse=max(abs(C6),[],1)./sqrt(sum(abs(C6)));I6impluse=max(abs(I6),[],1)./sqrt(sum(abs(I6))); N6impluse=max(abs(N6),[],1)./sqrt(sum(abs(N6)));O6impluse=max(abs(O6),[],1)./sqrt(sum(abs(O6))); B7impluse=max(abs(B7),[],1)./sqrt(sum(abs(B7)));C7impluse=max(abs(C7),[],1)./sqrt(sum(abs(C7)));I7impluse=max(abs(I7),[],1)./sqrt(sum(abs(I7))); N7impluse=max(abs(N7),[],1)./sqrt(sum(abs(N7)));O7impluse=max(abs(O7),[],1)./sqrt(sum(abs(O7))); %峰值指标 B1crest=max(abs(B1),[],1)./B1rms;C1crest=max(abs(C1),[],1)./C1rms;I1crest=max(abs(I1),[],1)./I1rms;N1crest=max(abs(N1),[],1)./N1rms;O1crest=max(abs(O1),[],1)./O1rms; B5crest=max(abs(B5),[],1)./B5rms;C5crest=max(abs(C5),[],1)./C5rms;I5crest=max(abs(I5),[],1)./I5rms;N5crest=max(abs(N5),[],1)./N5rms;O5crest=max(abs(O5),[],1)./O5rms; B2crest=max(abs(B2),[],1)./B2rms;C2crest=max(abs(C2),[],1)./C2rms;I2crest=max(abs(I2),[],1)./I2rms;N2crest=max(abs(N2),[],1)./N2rms;O2crest=max(abs(O2),[],1)./O2rms; B3crest=max(abs(B3),[],1)./B3rms;C3crest=max(abs(C3),[],1)./C3rms;I3crest=max(abs(I3),[],1)./I3rms;N3crest=max(abs(N3),[],1)./N3rms;O3crest=max(abs(O3),[],1)./O3rms; B6crest=max(abs(B6),[],1)./B6rms;C6crest=max(abs(C6),[],1)./C6rms;I6crest=max(abs(I6),[],1)./I6rms;N6crest=max(abs(N6),[],1)./N6rms;O6crest=max(abs(O6),[],1)./O6rms; B7crest=max(abs(B7),[],1)./B7rms;C7crest=max(abs(C7),[],1)./C7rms;I7crest=max(abs(I7),[],1)./I7rms;N7crest=max(abs(N7),[],1)./N7rms;O7crest=max(abs(O7),[],1)./O7rms; %裕度 B1cle=n^2*max(abs(B1),[],1)./(sum(sqrt(abs(B1))).^2);C1cle=n^2*max(abs(C1),[],1)./(sum(sqrt(abs(C1))).^2);I1cle=n^2*max(abs(I1),[],1)./(sum(sqrt(abs(I1))).^2); N1cle=n^2*max(abs(N1),[],1)./(sum(sqrt(abs(N1))).^2);O1cle=n^2*max(abs(O1),[],1)./(sum(sqrt(abs(O1))).^2); B5cle=n^2*max(abs(B5),[],1)./(sum(sqrt(abs(B5))).^2);C5cle=n^2*max(abs(C5),[],1)./(sum(sqrt(abs(C5))).^2);I5cle=n^2*max(abs(I5),[],1)./(sum(sqrt(abs(I5))).^2); N5cle=n^2*max(abs(N5),[],1)./(sum(sqrt(abs(N5))).^2);O5cle=n^2*max(abs(O5),[],1)./(sum(sqrt(abs(O5))).^2); B2cle=n^2*max(abs(B2),[],1)./(sum(sqrt(abs(B2))).^2);C2cle=n^2*max(abs(C2),[],1)./(sum(sqrt(abs(C2))).^2);I2cle=n^2*max(abs(I2),[],1)./(sum(sqrt(abs(I2))).^2); N2cle=n^2*max(abs(N2),[],1)./(sum(sqrt(abs(N2))).^2);O2cle=n^2*max(abs(O2),[],1)./(sum(sqrt(abs(O2))).^2); B3cle=n^2*max(abs(B3),[],1)./(sum(sqrt(abs(B3))).^2);C3cle=n^2*max(abs(C3),[],1)./(sum(sqrt(abs(C3))).^2);I3cle=n^2*max(abs(I3),[],1)./(sum(sqrt(abs(I3))).^2); N3cle=n^2*max(abs(N3),[],1)./(sum(sqrt(abs(N3))).^2);O3cle=n^2*max(abs(O3),[],1)./(sum(sqrt(abs(O3))).^2); B6cle=n^2*max(abs(B6),[],1)./(sum(sqrt(abs(B6))).^2);C6cle=n^2*max(abs(C6),[],1)./(sum(sqrt(abs(C6))).^2);I6cle=n^2*max(abs(I6),[],1)./(sum(sqrt(abs(I6))).^2); N6cle=n^2*max(abs(N6),[],1)./(sum(sqrt(abs(N6))).^2);O6cle=n^2*max(abs(O6),[],1)./(sum(sqrt(abs(O6))).^2); B7cle=n^2*max(abs(B7),[],1)./(sum(sqrt(abs(B7))).^2);C7cle=n^2*max(abs(C7),[],1)./(sum(sqrt(abs(C7))).^2);I7cle=n^2*max(abs(I7),[],1)./(sum(sqrt(abs(I7))).^2); N7cle=n^2*max(abs(N7),[],1)./(sum(sqrt(abs(N7))).^2);O7cle=n^2*max(abs(O7),[],1)./(sum(sqrt(abs(O7))).^2); %峭度 B1kur=kurtosis(B1);B2kur=kurtosis(B2);B3kur=kurtosis(B3);B5kur=kurtosis(B5);B6kur=kurtosis(B6);B7kur=kurtosis(B7); C1kur=kurtosis(C1);C2kur=kurtosi
评论
    相关推荐