waveletfunc1.m

使用小波分析进行转速计算的程序

%%%采样点1000MS测转速程序
clear all;
close all;
rotatespeed=zeros(6,5);
load 004.dat;
for j=1:6;
    m=j*10000;
x=X004(1:m);

figure(1);
subplot(6,1,1);
plot(x,'b');
title('原始信号数据和各层近似');
xlabel('时间（ms）');
ylabel('激光CCD读数');


[c,l]=wavedec(x,5,'db3');
for i=1:5
    decom=wrcoef('a',c,l,'db3',i);
    subplot(6,1,i+1);
    plot(decom);
    ylabel(['a',num2str(i)]);
    %%%%%%%%%%%%%%%%%%下面进行fft求取该层的分解获得的频率
    
    ffts=fft(decom);
    absffts=abs(ffts);
    
    N = length(ffts);
    ffts(1) = [];
    power = abs(ffts(1:N/2)).^2;
    nyquist = 1/2;
    freq = (1:N/2)/(N/2)*nyquist;
   
    period = 1./freq;
    
    [mp,index] = max(power);
    t=period(index);
    n_zhuansu=60*1000/t
    rotatespeed(j,i)=n_zhuansu;
end
end
%%%%%%%%%%%%%%%%%%%
figure(2);
subplot(6,1,1);
plot(x,'b');
title('原始信号数据和各层近似');
xlabel('时间（ms）');
ylabel('激光CCD读数');

[c,l]=wavedec(x,5,'db3');
for i=1:5
    decom=wrcoef('d',c,l,'db3',i);
    subplot(6,1,i+1);
    plot(decom);
    ylabel(['d',num2str(i)]);
end
    
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%