Matlab 信号小波阈值去噪实例_长笛人倚楼Gloria

http://blog.sina.com.cn/u/2214742602

首页博文目录关于我

个人资料

微博

加好友发纸条

写留言加关注

博客等级：
博客积分：

博客访问：
关注人气：
获赠金笔：0支
赠出金笔：0支
荣誉徽章：

正文字体大小：大中小

Matlab 信号小波阈值去噪实例

(2014-02-19 17:25:25)

标签：

小波

阈值

去噪

ddencmp

wdencmp

分类：信号处理

截取自 http://blog.sina.com.cn/s/blog_6163bdeb0102dw7a.html

本实例中去噪过程使用的函数是ddencmp 和 wdencmp。详细介绍可以参考 matlab help 或者 http://blog.sina.com.cn/s/blog_84024a4a0101fn1a.html

1. 载入信号并进行小波分解和重构

clear all;

close all;

clc;

load leleccum

% 截取信号

s = leleccum(1:3920);

l_s = length(s);

[C,L] = wavedec(s,3,'db1');

cA3 = appcoef(C,L,'db1',3);

cD3 = detcoef(C,L,3);

cD2 = detcoef(C,L,2);

cD1 = detcoef(C,L,1);

A3 = wrcoef('a',C,L,'db1',3);

% 从C中重建1、2、3层细节

D1 = wrcoef('d',C,L,'db1',1);

D2 = wrcoef('d',C,L,'db1',2);

D3 = wrcoef('d',C,L,'db1',3);

figure

subplot(2,2,1); plot(A3);

title( 'Approximation A3' )

subplot(2,2,2); plot(D1);

title( 'Detail D1' )

subplot(2,2,3); plot(D2);

title( 'Detail D2' )

subplot(2,2,4); plot(D3);

title( 'Detail D3' )

http://s4/middle/002pSPISzy6GIc1gHgTd3&690信号小波阈值去噪实例" TITLE="Matlab 信号小波阈值去噪实例" />

2. 使用近似信号进行重构，实现粗糙去噪

figure

subplot(2,1,1);plot(s);title('Original'); axis off

subplot(2,1,2);plot(A3);title('Level 3 Approximation');

axis off

http://s2/middle/002pSPISzy6GIc1iNlnf1&690信号小波阈值去噪实例" TITLE="Matlab 信号小波阈值去噪实例" />

3. 阈值去噪

[thr,sorh,keepapp] = ddencmp('den','wv',s);

clean = wdencmp('gbl',C,L,'db1',3,thr,sorh,keepapp);

figure

subplot(2,1,1); plot(s(2000:3920)); title('Original')

subplot(2,1,2); plot(clean(2000:3920)); title('denoised')

http://s1/middle/002pSPISzy6GIc1lQMU90&690信号小波阈值去噪实例" TITLE="Matlab 信号小波阈值去噪实例" />

来自为知笔记(Wiz)

阅读┊ 收藏 ┊ 喜欢 ▼ ┊打印┊举报/Report

前一篇：[转载]MATLAB中小波阈值去噪

后一篇：Matlab help 小波降噪

新浪BLOG意见反馈留言板　欢迎批评指正