声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3818|回复: 7

[小波] 请高手帮忙看一下这个小波去噪程序

[复制链接]
发表于 2007-8-17 12:07 | 显示全部楼层 |阅读模式

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

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

x
附件中有本人上传的原始数据。
我想把信号中的低频成分滤掉,想保留30K-150K之间的信号,采样频率为500K。打算用小波去噪,进行四层分解后,把第四层低频信号置零,重构信号,所编程序如下,但是老是报错,请高手指教。
%载入信号,画出波形

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;


报错信息为:
  
??? Error using ==> horzcat
All matrices on a row in the bracketed expression must have the
same number of rows.

Error in ==> p1 at 49
cl=[caa4,cd4,cd3,cd2,cd1]

错误就在程序的倒数第4行,程序中已标出。
请高手帮忙,错在哪里了。

补充一下降噪方法:
       进行四层小波分解后,把第四层低频信号置零,然后重构信号,结果重构时出错。

1.txt

128.5 KB, 下载次数: 147

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2007-8-17 14:43 | 显示全部楼层
在下语句中caa4,cd4,cd3,cd2和cd1都是列数据:
cl=[caa4,cd4,cd3,cd2,cd1]
所以造成错误,改为
cl=[caa4‘,cd4’,cd3‘,cd2’,cd1‘]
就可执行了。

评分

1

查看全部评分

 楼主| 发表于 2007-8-17 15:03 | 显示全部楼层

songzy41 ,谢谢!

原来是这个问题,
非常感谢songzy41 !
你也用小波作东西吗?
有没有联系方式,请多多指教。
发表于 2007-8-17 18:39 | 显示全部楼层
我用小波不算很熟练,但在论坛上有不少搞小波的高手,有问题可拿到论坛上一起讨论。
发表于 2009-8-1 10:33 | 显示全部楼层
楼主为什么要进行单支重构?感觉没什么需要啊!
发表于 2009-8-2 17:24 | 显示全部楼层

回复 板凳 songzy41 的帖子

在下语句中caa4,cd4,cd3,cd2和cd1都是列数据:
cl=[caa4,cd4,cd3,cd2,cd1]
所以造成错误,改为
cl=[caa4‘,cd4’,cd3‘,cd2’,cd1‘]
就可执行了。

改了运行还是有问题  ,画不出figure(2)
请指导,谢谢。
发表于 2011-12-29 21:52 | 显示全部楼层
回复 6 # wangzemin 的帖子

可以啊,我这里能运行,贴不上图
发表于 2012-7-18 09:54 | 显示全部楼层
学习一下
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-14 14:46 , Processed in 0.065812 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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