声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1413|回复: 4

[小波] 求教小波降噪程序

[复制链接]
发表于 2009-4-15 11:11 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
各位大侠,求小波降噪程序,最好是用heursure 规则降噪,非常感激!!

[ 本帖最后由 咕噜噜 于 2009-4-15 11:30 编辑 ]

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2009-4-15 20:38 | 显示全部楼层
论坛里有,请参考一下

clear all
clc;
load('c:\temp\1.txt');
s1=X1;

fs=500000;  %采样频率
t=1:max(size(s1));
t=1000*t./fs;
t=t';
plot(s1),title('plot(s1)'),grid on;
%axis([0,51.2,-10000,10000]);
xlabel('时间(t)   单位:ms');
ylabel('电压(V)');
title('原始信号');

%小波分解
[c,l]=wavedec(s1,4,'db4');
[cd1,cd2,cd3,cd4]=detcoef(c,l,[1 2 3 4]);
ca4=appcoef(c,l,'db4',4);
a4=wrcoef('a',c,l,'db4',4);
d4=wrcoef('d',c,l,'db4',4);
d3=wrcoef('d',c,l,'db4',3);
d2=wrcoef('d',c,l,'db4',2);
d1=wrcoef('d',c,l,'db4',1);
%不需要的系数置零
caa4=zeros(size(ca4));

%消噪后信号重构,并画图
cl=[caa4‘,cd4’,cd3‘,cd2’,cd1‘]
s2=waverec(cl,l,'db4');
figure(2);
plot(t,s2);title('消噪后重构信号');grid on;
发表于 2010-4-7 20:12 | 显示全部楼层
fs=500000;  %采样频率
t=1:max(size(s1));
t=1000*t./fs;
t=t';


上面的是什么意思,程序里好像没有用到
发表于 2010-4-8 02:40 | 显示全部楼层
1000点FFT 用的吧,可能原代码里有画频谱图。
发表于 2010-4-11 10:30 | 显示全部楼层
s2=wave(cl,l,'db4')这行出错,是不是上面的那行程序出错啊,cl=['caa4','cd4','cd3','cd2','cd1']。因为这些cd1与cd2的大小就不同。忘各位指点。应该是什么样子的
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-5-13 08:50 , Processed in 0.058188 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表