表面肌实验报告

gend('原始信号','去噪后信号');title('train样本滤波前后信号对比');subplot(4,1,2);plot(sample_train(:,2),'r');holdon;plot(new_train(:,2),'b');subplot(4,1,3);plot(sample_train(:,3),'r');holdon;plot(new_train(:,3),'b');subplot(4,1,4);plot(sample_train(:,4),'r');holdon;plot(new_train(:,4),'b');%%%%%trin样本特征提取%%%%%%Threshold=0.006;WAMP=zeros(Channel,Count);forc=1:Countforn=1:Channelforw=1:Window-1if(abs(new_train((c-1)*WinLap+w,n)-new_train((c-1)*WinLap+w+1,n))>Threshold)WAMP(n,c)=WAMP(n,c)+1;endendendendFeature_train=WAMP;%figure(3)%plot(Feature_train);%%%%%test样本小波去噪%%%%%%fori=1:Channelx=sample_test(:,i);[thr,sorh,keepapp]=ddencmp('den','wv',x);[c,l]=wavedec(x,3,'db1');a3=appcoef(c,l,'db1',3);d3=detcoef(c,l,3);d2=detcoef(c,l,2);d1=detcoef(c,l,1);s4=wdencmp('gbl',c,l,'db1',3,thr,sorh,keepapp);new_test(:,i)=s4;end