leeting918 发表于 2009-4-15 11:11

求教小波降噪程序

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

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

news 发表于 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();
xlabel('时间（t）   单位：ms');
ylabel('电压（V)');
title('原始信号');

%小波分解
=wavedec(s1,4,'db4');
=detcoef(c,l,);
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=
s2=waverec(cl,l,'db4');
figure(2);
plot(t,s2);title('消噪后重构信号');grid on;

woaizhoulu 发表于 2010-4-7 20:12

fs=500000;%采样频率
t=1:max(size(s1));
t=1000*t./fs;
t=t';


上面的是什么意思，程序里好像没有用到

sparrow 发表于 2010-4-8 02:40

1000点FFT 用的吧，可能原代码里有画频谱图。

youxue 发表于 2010-4-11 10:30

s2=wave(cl,l,'db4')这行出错，是不是上面的那行程序出错啊，cl=['caa4','cd4','cd3','cd2','cd1']。因为这些cd1与cd2的大小就不同。忘各位指点。应该是什么样子的
页: [1]
查看完整版本: 求教小波降噪程序