Matlab小波变换对奇异点的检测.docxVIP

因此可用做小波母函数。 因此可用做小波母函数。 Matlab小波变换对于奇异点的检测 i信号的突变性 突变信号又称奇异信号，突变信号的突变点经常携带比较重要的信息 ，是信号的重要特 征之一。在数字信号处理和数字图像处理中具有非常重要的作用和地位， 信号的突变性检测 是先对原信号在不同尺度上进行“磨光”，再对磨光后信号的一阶或二阶倒数检测其极值点 或过零点。对信号进行磨光处理，主要是为了消除噪声而不是边缘。 传统的信号突变检测方 法是基于傅立叶变换的，由某一函数的傅立叶变换趋近于零的快慢来推断该函数是否具有突 变性，但它只能反映信号的整体突变性， 而对信号的局部突变则无法描述。 这样我们就引入 小波变换算法。 2 ?信号的突变点的检测原理 设h(t)是函数f(t)和g(t)的卷积，即： h(t) f(t) g(t) 则根据傅立叶变换的性质有： F[h(t)] j F[f (t) g(t)] j f( )g() =[j f( )]g( )=f( )[j g()] = F[f(t)] F[g(t)] F[f(t)] F[g(t)] 所以得到：h(t) f(t) g(t) f(t) g(t) 若将函数f(t)看作是信号，g(t)看作是滤波器，那么信号的导数与滤波器的卷积结果可以 看作是滤波器的导数与信号的卷积。例如，如果选 g(t)为高斯函数，则利用其导数可以构造 f(t)的突变点和极值点Morlet小波和Maar f(t)的突变点和极值点 设(t)是一个起平滑作用的低通平稳函数，且满足条件 (t)dt1, lim (t)dt 1, lim (t) 0 /2通常取(t)为高斯函数，即 /2 (t) 假设(t)是二次可导的，并且定义 (1)(t)- (1)(t)- (t)dt te t2/2 1 2 t2/2 .2 0 t)e 则函数⑴ 则函数⑴(t)、 ⑵(t)满足小波的容许条件: ⑴⑴dt 0， (2)(t)dt 0 若记s 1 1，则s(t)表示(t)在尺度因子s下的伸缩。由于小波变换就是将原信 s s 号f (t)同伸缩小波卷积得到的，为此以 ⑴(t),⑵(t)为小波函数定义的卷积型小波变换 为： (1). Ws f (t) f * s0) (t) d s f * s dt (t) s?(f * dt s)(t) (2) ws f(t) f * s2) (t)- f* 2d2s f * s 2 dt (t) 2d2 sdt2(f * s)(t) 由此可见，小波变化 wsi) f (t), ws2) f (t)分别是函数f(t)在尺度s下由(t)平滑后再 取一阶、二阶导数。当 s较小时，用 s(t)对f(t)平滑的结果对 f(t)的突变位置影响不大； 当s较大时，则此平滑过程会将 f (t)的一些细小的突变削去，而只剩下大尺寸的突变。由 此我们可知，当小波函数可看作某一平滑函数的一阶导数时， 信号小波变换模的局部极值点 对应信号的突变点(或边缘)。当小波函数可看作某一平滑函数的二阶导数时，信号小波变 换模的过零点，也对应信号的突变点(或边缘) 。 这就是采用检测小波变换系数模的过零点和局部极值点可检测信号突变点 (或边缘)的 原理。 Matlab小波变换检测奇异点 原始信号是含有奇异点的信号，为确定该奇异点的时间，采用 haar小波进行连续小波变 换后，在对系数进行分析处理。 仿真程序如下： figure(1) plot(cuspamax) xlabel( 时间’);ylabel( 幅值); title(频率突变信号); figure(2) [c,l]=wavedec(cuspamax,5, db6); cfd=zeros(5,1024); for k=1:5 d=detcoef(c,l,k); d=d(o nes(1,2Ak),:); cfd(k,:)=wkeep(d(:),1024) end cfd=cfd(:); I=fin d(abs(cfd)sqrt(eps)); cfd(I)=zeros(size(I)); cfd=reshape(cfd,5,1024); colormap(p in k(64)); img=image(flipud(wcodemat(cfd,64, row ))); set(get(img, parent ), YtickLabel ,[]); title(离散小波变换后系数的绝对值 ) ylabel(层数); figure(3) ccfs=cwt(cuspamax,1:32, haar , plot ); title(连续小波变换系数的绝对值 ) colormap(p in k(64)); ylabel(尺度) xlabel(时间(或者空间)’) 程序的运行结果如下图所示： 图1原始信号的示意图