声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1083|回复: 2

[小波] 请教小波变换的问题~急!

[复制链接]
发表于 2009-5-2 21:38 | 显示全部楼层 |阅读模式

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

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

x
我对一个语音信号加噪 然后用小波做降噪,用不同的阈值选取方法来去噪,但是程序出现错误,请高手帮忙看看 错误显示如下:
On line 40  ==> c1=[a3 dd1 dd2 dd3];
程序如下
[y,fs,bits]=wavread('matlab.wav');
t=0:1/fs:(size(y)-1)/fs;    %将所加噪声信号的点数调整到与原始信号相同
x2=y+rand(size(y));
sound(x2,fs)

%显示波形
subplot(231);
plot(t,y);
xlabel('样本序号 n');
ylabel('(原始信号)幅值 A');

subplot(232);
plot(x2);
xlabel('样本序号 n');
ylabel('(含噪信号)幅值 A');

%用db1小波对原始信号进行3层分解并提取系数
[c,l]=wavedec(x2,3,'db1');
a3=appcoef(c,l,'db1',3);
d3=detcoef(c,l,3);
d2=detcoef(c,l,2);
d1=detcoef(c,l,1);

%对信号进行强制性消噪处理并图示结果
dd3=zeros(1,length(d3));
dd2=zeros(1,length(d2));
dd1=zeros(1,length(d1));
c1=[a3 dd1 dd2 dd3];
s1=waverec(c1,l,'db1');

subplot(233);
plot(s1);
title('强制消噪后的信号');
xlabel('样本序号 n');
ylabel('幅值 A');
回复
分享到:

使用道具 举报

发表于 2010-4-7 20:02 | 显示全部楼层
把冒号加上试试
发表于 2010-4-7 21:07 | 显示全部楼层
程序经过调整,是因为维度问题。具体调整如下
将c1=[a3 dd1 dd2 dd3];更改为c1=[a3',dd3,dd2,dd1];即可使用了。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-29 17:20 , Processed in 0.056133 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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