使用小波分析的各类函数对信号图像进行消噪(含数据以及图像)珍贵版摘要.docVIP

用sym6和db6小波对信号y1进行5层分解。其中信号y1是由正余弦信号y和白噪声信号s构成。选用sure阈值模式。Wden函数是对一维信号的小波进行消噪处理。 xd=wden(x,tptr,sorh,scal,n,’wname’) X即为将要去噪的信号。tptr为所选用的sure阈值模式。sorh为函数选择阈值使用方式，其中s表示软阈值，h为硬阈值。输入参数scal规定阈值处理随噪声水平变化。scal=one,不随噪声水平变化。scal=sln,根据第一层小波分解的噪声水平估计进行调整。scal=mln根据每一层小波分解的噪声水平估计进行调整。 对于函数heursure为启发式阈值，rigrsure为stein无偏估计，sqtwolog 为固定式阈值 minimaxi为极大值极小值阈值。 应用小波分析对信号去噪的程序如下： t=0:0.001:1 f1=5; f2=20; y=3*sin(2*pi*f1*t)+5*cos(5*pi*f2*t) s1=randn(1,length(y)) s=y+s1 subplot(211); plot(t,s); grid on; lev=5; xdH=wden(s,heursure,s,sln,lev,sym6); xdR=wden(s,rigrsure,s,sln,lev,sym6); xdS=wden(s,sqtwolog,s,sln,lev,sym6); xdM=wden(s,minimaxi,s,sln,lev,sym6); subplot(5,2,1); plot(y);title(原始信号) axis([1,2048,-10,10]); subplot(5,2,2); plot(s);title(有噪信号) axis([1,2048,-10,10]); subplot(5,2,3); plot(xdH);xlabel(heursure阈值消噪处理后的信号) axis([1,2048,-10,10]); subplot(5,2,4); plot(xdR);xlabel(rigrsure阈值消噪处理后的信号) axis([1,2048,-10,10]); subplot(5,2,5); plot(xdS);xlabel(sqtwolog阈值消噪处理后的信号) axis([1,2048,-10,10]); subplot(5,2,6); plot(xdM);xlabel(minimaxi阈值消噪处理后的信号) Subplot(5,2,7); plot(y-xdH);xlabel(heursure阈值消噪后与原信号比较) subplot(5,2,8); plot(y-xdR);xlabel(rigrsure阈值消噪后与原信号比较) subplot(5,2,9); plot(y-xdS);xlabel(sqtwolog阈值消噪后与原信号比较) subplot(5,2,10); plot(y-xdM);xlabel(minimaxi阈值消噪后与原信号比较) A=var(y-xdH) B=var(y-xdR) C=var(y-xdS) D=var(y-xdM) A = 0.2415 B = 0.2421 C = 2.2178 D = 1.0703 各信号图形如下： 当采用db6作为小波时： t=0:0.001:1 f1=5; f2=20; y=3*sin(2*pi*f1*t)+5*cos(5*pi*f2*t) s1=randn(1,length(y)) s=y+s1 subplot(211); plot(t,s); grid on; lev=3; xdH=wden(s,heursure,s,sln,lev,db6); xdR=wden(s,rigrsure,s,sln,lev,db6); xdS=wden(s,sqtwolog,s,sln,lev,db6); xdM=wden(s,minimaxi,s,sln,lev,db6); subplot(5,2,1); plot(y);title(原始信号) axis([1,2048,-10,10]); subplot(5,2,2); plot(s);title(有噪信号) axis([1,2048,-10,10]); subplot(5,2,3); plot(xdH);xlabel(heursure阈值消噪处理后的信号) axis([1,2048,-10,10]); subplot(5,2,4); plot(xdR);xlabel(rigrsure阈值消噪处理后的信号) axis([1,2048,-10,10]); subplot(5,2,5); plot(xdS);xlabel(sqtwolog阈值