xiaobojiangzao.m

自己编写的比较全面的故障诊断matlab函数程序

%**************************************************************************************************************************************************************
%   本程序为函数功能，目的是对输入信号进行小波降噪，采用低通滤波的方法，将高于阀值的细节信号去掉
%      
%   中国北方发动机研究所试验测试中心  苏志霄   2006年10月
%**************************************************************************************************************************************************************
%
function R=xiaobojiangzao(data,n,m);           %n为小波分解的尺度,m为希望保留的细节信号阶数，m<n,R为输出的信号
len=max(size(data));
R=zeros(n+1,len);
[c,l]=wavedec(data,n,'db3');
%对分解结构[c,l]中的第n逼近部分进行重构
s(1,:)=wrcoef('a',c,l,'db3',n);
%对分解结构[c,l]中的各层细节部分进行重构
for i=1:n;
    s(i+1,:) = wrcoef('d',c,l,'db3',n+1-i);
end;
R=s(1,:);
for i=1:m-1;
    R=R+R(i+1);
end;