小波去噪MATLAB代码「建议收藏」

小波去噪MATLAB代码「建议收藏」本文基于北京交通大学陈后金教授的课件。我加以整理,若有冒犯还请谅解1利用MATLAB产生分解与重建滤波器组计算滤波器组的函数[Ld,Hd,Lr,Hr]=wfilters(‘wname’)Ld:分解低通滤波器h0[-n];Hd:分解高通滤波器h1[-n];Lr:分解低通滤波器h0[-n];Hr:分解高通滤波器h1[-n];wfname:小波名eg1:计算db2小波的四…

大家好,又见面了,我是你们的朋友全栈君。

本文基于北京交通大学陈后金教授的课件。我加以整理,若有冒犯还请谅解

1利用MATLAB产生分解与重建滤波器组

计算滤波器组的函数[Ld,Hd,Lr,Hr]=wfilters(‘wname’)

Ld:分解低通滤波器h0[-n];

Hd:分解高通滤波器h1[-n];

Lr:分解低通滤波器h0[-n];

Hr:分解高通滤波器h1[-n];

wfname:小波名

eg1:计算db2小波的四个滤波器,并画出其时域波形。

MATLAB程序如下:
wn=’db2′;
[Ld,Hd,Lr,Hr] = wfilters(wn);
k=0:3;
subplot(221);stem(k,Ld);

title(‘低通分解滤波器Ld’);
subplot(222);stem(k,Lr);
title(‘低通重建滤波器Lr’);
subplot(223);stem(k,Hd);

title(‘高通分解滤波器Hd’);
subplot(224);stem(k,Hd);
title(‘高通重建滤波器Hr’);

2利用MATLAB计算小波函数

[phi,psi,t]=wavefun(‘wname’,lter)

wname:小波名

lter:计算过程的迭代次数

phi:尺度函数ϕ(t)

psi:小波函数ψ(t)

t:尺度函数与小波函数的抽样点

eg2:利用MATLAB计算db2小波的尺度函数与小波函数。

MATLAB程序如下:

lter=20;

wname=’db2′;

[s,w,t]=wavefun(wname,lter);

subolot(211);plot(t,s);

title(‘尺度函数’)

subplot(212);plot(t,w);

title(‘小波函数’)

3 利用MATLAB计算一维DWT和IDWT

计算多级DWT和IDWT的函数为wavedec和waverec,其调用格式为

[C,L]=wavedec(x,N,’wname’)

x=waverec(C,L,’wname’)

其中

wname: 小波名; 
x: 时域信号;
N: 小波变换的级数;
C = [cAN cDN cDN-1 ••• cD1];
L(1)= cAN的长度, L(2)= cDN的长度,
L(N+1)= cD1的长度, L(N+2)= x的长度

3利用部分小波系数重建信号

x=wrcoef(‘type’, C, L, ‘wname’, N)
type=’a’ 由第N级近似分量重建信号
type=‘d’ 由第N级细节分量重建信号
wname: 小波名

若 C = [cA3 cD3 cD2 cD1]
x=wrcoef(‘a’,C,L, ‘wname’,3)
x=IDWT{[cA3 0 0 0]}
x=wrcoef(‘a’,C,L, , ‘wname’,2)
x=IDWT{[cA3 cD3 0 0] }=IDWT{[cA2 0 0]}

eg3 已知某信号的波形如图所示,试计算其5级小波变换系数,
由第5、3、1级小波近似系数重建的信号。

小波去噪MATLAB代码「建议收藏」

MATLAB程序如下:
wname=’db1′;
dwtmode(‘per’)
t=linspace(0,1,1024);
x=20*t.^2.*(1-t).^4.*cos(12*pi*t);
subplot(511);plot(t,x);
axis([0 1 -0.5 0.5]);title(‘Signal’);

[C,L] = wavedec(x,5,wname);
subplot(512);plot(t,C); axis([0 1 -3 3]);

A5=wrcoef(‘a’,C,L,wname,5);
subplot(513);plot(t,A5);
axis([0 1 -0.5 0.5]);
A3=wrcoef(‘a’,C,L,wname,3);
subplot(514);plot(t,A3);
axis([0 1 -0.5 0.5]);

A1=wrcoef(‘a’,C,L,wname,1);
subplot(515);plot(t,A1);
axis([0 1 -0.5 0.5]);

4基于小波的信号去噪

XD = wden(X, TPTR, SORH, SCAL, N, ‘wname’)
其中:
XD: 对噪声信号X去噪后得到的信号;
X: 含噪声信号;
TPTR: 阈值规则,主要有’rigrsure’, ‘heursure’, ‘sqtwolog’, ‘minimaxi’;
SORH: 阈值方法, ‘s’ (soft阈值), ‘h’ (hard阈值);
SCAL: 阈值尺度的调整方法,主要有’one’, ‘ sln’, ‘ mln’ ;
N: 离散小波变换的级数
wname: 小波名

eg4 试利用函数wden对含有噪声的blocks信号进行去噪。

MATLAB程序如下:
snr=5; % 噪声方差
[x, xn] = wnoise(‘blocks’,11,snr);
k=0:length(x)-1;
subplot(311);plot(k,x);
title(‘原信号’);
subplot(312);plot(k,xn);
title(‘含噪信号’);
lev=5;wn=’db1′;
% 利用soft SURE阈值规则去噪
xd1= wden(xn, ‘heursure’, ‘s’, ‘one’, lev, wn);
subplot(313);plot(k,xd1);
title(‘去噪后的结果’);

5基于小波的信号压缩

NC= wthcoef(‘d’, C, L, N)
其中:
‘d’: 表示对DWT系数C中细节(detail)分量进行压缩;
C,L: 由wavedec得到的DWT系数;
N: 若N=[1 2 3]表示将C中1、2和3级细节分量置零。
NC: 由系数C经过压缩后得到的新系数;

eg5试利用函数wthcoef对leleccum信号进行压缩。

MATLAB程序如下:
load leleccum; 
x = leleccum(1001:2024)*0.95/100;
k=0:length(x)-1;
[C,L] = wavedec(x,5,’db3′);
NC1 = wthcoef(‘d’,C,L,[1]);
x1 = waverec(NC1, L, ‘db3’);
NC5 = wthcoef(‘d’,C,L,[1 2 3 4 5]) ; 
x2 = waverec(NC5, L, ‘db3’);
subplot(311);plot(k,x);title(‘原信号’);
subplot(312);plot(k,x1);title(‘2:1压缩’);
subplot(313);plot(k,x2);title(’32:1压缩’);

 

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/148179.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • nginx自动重启脚本_nginx进程数和连接数

    nginx自动重启脚本_nginx进程数和连接数工作中经常要杀掉nginx对应的进程,然后手动命令重启和删除nginx日志,每次都要输入3个或3个以上的命令,效率低。写个脚本万事大吉。#!/bin/shNAME=”nginx-p/home/proj/gamebin/-cconf/depnginx.conf”ID=`ps-ef|grep”$NAME”|grep-v”$0″|grep-v”…

  • java etl 简单实例_东方通ETL开发实例[通俗易懂]

    java etl 简单实例_东方通ETL开发实例[通俗易懂]东方通ETL开发实例下面通过一个简单的“OracleCDC增量抽取”实例,带大家感受一下TIETL的开发思路和强大的数据处理能力。功能需求:同一数据库(localhost_etl)下,现有cdc_source和cdc_target两张表如下:表cdc_source表cdc_target现表cdc_source的增量(增、删、改)需同步到cdc_target中TIETL环境下开发思路:类似于jav…

  • K3 官改新手小白配置阿里DDNS 超级详细「建议收藏」

    K3 官改新手小白配置阿里DDNS 超级详细「建议收藏」K3官改新手小白配置阿里DDNS超级详细写的比较仓促,不对之处请指正,这个是写给小白看的,大神勿喷首先介绍一下什么是DDNSDDNS(DynamicDomainNameServer)是动态域名服务的缩写。DDNS是将用户的动态IP地址映射到一个固定的域名解析服务上,用户每次连接网络的时候客户端程序就会通过信息传递把该主机的动态IP地址传送给位于服务商主机上的服务器程序,服务器…

  • 基于LM331的频率电压转换电路「建议收藏」

    基于LM331的频率电压转换电路「建议收藏」常用的模拟信号的传输方式有电压传输、电流传输和频率传输,其中电压传输的方式最为简单方便,成本最低,但是电压信号在传输的过程中最容易受到干扰,并且传输过程会有损耗,因此不适合远距离传输。将电压转换为电流或者频率后,可以进行远距离传输,且抗干扰能力强,其中电流传输的抗干扰能力最好,传输距离最远,但是成本较高,而频率在距离超过100m时波形会失真,抗干扰能力介于电压和电流之间。本文主要介绍基于LM331的频率-电压转换电路,关于LM331的简介和电压-频率转换电路可以参考《基于LM331的电压频率转换电路》。

  • 最全静态网页模板网站下载_网页首页内容模板

    最全静态网页模板网站下载_网页首页内容模板前言最近好多人私信我:我学了不久前端,想要用一些项目实战来练练手。那么,我以后将定期更新静态网页实战,专用于提供初学者的练习。那么这次博客为大家提供一些模板网站,希望初学者借鉴和学习优质资源:以下网站对于初学者的作用:提供免费的静态网页下载提高自己的css+html5设计水平http://sc.chinaz.com/moban/“站长素材”网站中包括大量的免费模板,可任意下载。…

  • 大芒果魔兽单机版GM命令大全_魔兽335单机大芒果全套文件

    大芒果魔兽单机版GM命令大全_魔兽335单机大芒果全套文件命令:.modifyhp99999999999999增加被选择人物的血量.modifygold99999999增加被选择人物的金币数量.modifyspeed30加速.chea

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号