[2018年最新整理]小波变换及其应用2014.ppt

School of mathematics physics 华北电力大学数理学院 给定正交尺度函数 ，记正交小波函数为 ，因而两尺度关系为： 定义： 由 构成的函数集 称为小波包。 School of mathematics physics 华北电力大学数理学院 * * * * School of mathematics physics 华北电力大学数理学院 一维信号降噪 function yanshi10 load noisbump; x=noisbump; [c,l]=wavedec(x,5,sym6); [thr,sorh,keepapp]=ddencmp(den,wv,x); %全局降噪的阀值thr。 [thr1,nkeep]=wdcbm(c,l,2); %使用birge-massart策略确定各层降噪的阀值。 %全局软阈值降噪信号xc,保留一层以上的所有分解系数[cxc,lxc] [xc,cxc,lxc,perf0,perfl2] =wdencmp(gbl,c,l,sym6,5,thr,s,1); %分层软阀值降噪xc; [xc1,cxc1,lxc1,perf01,perfl21] =wdencmp(lvd,c,l,sym6,5,thr1,s); subplot(3,1,1);plot(x);title(原始信号,fontsize,8); subplot(3,1,2);plot(xc);title(‘全局软阈值降噪信号，能量成分[perf01,perfl21],fontsize,8); subplot(3,1,3);plot(xc1);title(分层软阀值降噪信号，能量成分[perf01,perfl21],fontsize,8); School of mathematics physics 华北电力大学数理学院 School of mathematics physics 华北电力大学数理学院 School of mathematics physics 华北电力大学数理学院 二维信号降噪 load woman randn(seed,2055615866); x=X+15*randn(size(X)); %加噪 [thr,sorh,keepapp]=ddencmp(den,wv,x); %全局阀值降噪 [xc,cxc,lxc,perf0,perfl2]=wdencmp(‘gbl’,x,‘sym4’,2,thr,sorh,keepapp); %显示 colormap(pink(255)),sm=size(map,1); subplot(221),image(wcodemat(X,sm));title(X); subplot(222),image(wcodemat(x,sm));title(x); subplot(223),image(wcodemat(xc,sm));title(xc); [c,s]=wavedec2(x,3,sym4);%三层二维小波分解 [th1,nkeep1]=wdcbm2(c,s,1.5,2.7*prod(s(1,:))); [xc1,cxc1,lxc1,perf01,perfl21]=wdencmp(lvd,c,s,sym4,3,th1,s); subplot(224),image(wcodemat(xc1,sm));title(xc1); School of mathematics physics 华北电力大学数理学院 School of mathematics physics 华北电力大学数理学院 应用2：小波变换用于信号奇异点检测 信号的奇异点（突变点）常包含重要的故障信息（对于图像它是边缘）。提取此突变点的位置及奇异性（光滑度）是信号分析与处理的重要工作之一 School of mathematics physics 华北电力大学数理学院 函数奇异性的描述，李氏（Lipschitz）指数α： 连续不可导函数：称函数f(t)在点 的李氏指数为（ ），如果满足： 不连续函数：称函数f(t)在点 的李氏指数α（ ），如果其原函数F(t)的李氏指数为α+1。 n阶连续导数，但n+1阶不可导函数：称函数f(t)在点 的李氏指数为α（ ），如果满足： School of mathematics physics 华北电力大学数理学院 例如： 李氏指数α=1 李氏指数α=0 李氏指数α=-1 School of mathematics physics 华北电力大