小波滤波器构造和消噪程序Word格式.docx
- 文档编号:21735039
- 上传时间:2023-02-01
- 格式:DOCX
- 页数:8
- 大小:15.41KB
小波滤波器构造和消噪程序Word格式.docx
《小波滤波器构造和消噪程序Word格式.docx》由会员分享,可在线阅读,更多相关《小波滤波器构造和消噪程序Word格式.docx(8页珍藏版)》请在冰豆网上搜索。
h1=sin(a)*cos(b);
h2=-sin(a)*sin(b);
h3=cos(a)*sin(b);
low_construct=[h0,h1,h2,h3];
L_fre=4;
%滤波器长度
low_decompose=low_construct(end:
-1:
1);
%确定h0(-n),低通分解滤波器
fori_high=1:
L_fre;
%确定h1(n)=(-1)^n,高通重建滤波器
if(mod(i_high,2)==0);
coefficient=-1;
else
coefficient=1;
end
high_construct(1,i_high)=low_decompose(1,i_high)*coefficient;
end
high_decompose=high_construct(end:
%高通分解滤波器h1(-n)
L_signal=100;
%信号长度
n=1:
L_signal;
%信号赋值
f=10;
t=0.001;
y=10*cos(2*pi*50*n*t).*exp(-20*n*t);
figure
(1);
plot(y);
title('
原信号'
);
check1=sum(high_decompose);
%h0(n)性质校验
check2=sum(low_decompose);
check3=norm(high_decompose);
check4=norm(low_decompose);
l_fre=conv(y,low_decompose);
%卷积
l_fre_down=dyaddown(l_fre);
%抽取,得低频细节
h_fre=conv(y,high_decompose);
h_fre_down=dyaddown(h_fre);
%信号高频细节
figure
(2);
subplot(2,1,1)
plot(l_fre_down);
小波分解的低频系数'
subplot(2,1,2);
plot(h_fre_down);
小波分解的高频系数'
l_fre_pull=dyadup(l_fre_down);
%0差值
h_fre_pull=dyadup(h_fre_down);
l_fre_denoise=conv(low_construct,l_fre_pull);
h_fre_denoise=conv(high_construct,h_fre_pull);
l_fre_keep=wkeep(l_fre_denoise,L_signal);
%取结果的中心部分,消除卷积影响
h_fre_keep=wkeep(h_fre_denoise,L_signal);
sig_denoise=l_fre_keep+h_fre_keep;
%信号重构
compare=sig_denoise-y;
%与原信号比较
figure(3);
subplot(3,1,1)
ylabel('
y'
%原信号
subplot(3,1,2);
plot(sig_denoise);
sig\_denoise'
%重构信号
subplot(3,1,3);
plot(compare);
compare'
%原信号与消噪后信号的比较
2.消噪
此函数用于消噪处理
%角度赋值
%此处赋值使滤波器系数恰为db9
%分解的高频系数采用db9较好,即它的消失矩较大
%分解的有用信号小波高频系数基本趋于零
%对于噪声信号高频分解系数很大,便于阈值消噪处理
[l,h]=wfilters('
db10'
'
d'
low_construct=l;
L_fre=20;
%滤波器长度
coefficient=1;
%原始信号赋值
y=10*cos(2*pi*50*n*t).*exp(-30*n*t);
zero1=zeros(1,60);
%信号加噪声信号产生
zero2=zeros(1,30);
noise=[zero1,3*(randn(1,10)-0.5),zero2];
y_noise=y+noise;
subplot(2,1,1);
plot(y_noise);
受噪声污染的信号'
%h0(n),性质校验
l_fre=conv(y_noise,low_decompose);
h_fre=conv(y_noise,high_decompose);
%消噪处理
fori_decrease=31:
44;
ifabs(h_fre_down(1,i_decrease))>
=0.000001
h_fre_down(1,i_decrease)=(10^-7);
%消噪后信号重构
%平滑处理
forj=1:
2
fori=60:
70;
sig_denoise(i)=sig_denoise(i-2)+sig_denoise(i+2)/2;
end;
%与原信号比较
%消噪后信号
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 滤波器 构造 程序