为什么我的emd分解要那么长的时间，一个小时多，我的声音信号才4.5s

星海

我的程序如下：2.wav才录了4.5s的声音，请高手指点下？
clear
[Y,FS,NBITS]=wavread('2.wav');
[m,n]=size(Y);
plot(linspace(0,m/FS,m),Y);
title('声音信号')
xlabel('t')
ylabel('Sound')
y= linspace(1,2,500);%添加趋势项
imf=emd(Y);
emd_visu(Y,1:length(Y),imf)
[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray))
grid on

zhangnan3509

采样频率是多少？是不是很高啊？

murhythm

应该是采样频率高吧，数据量大。。。

星海

原帖由 zhangnan3509 于 2008-3-15 21:26 发表
采样频率是多少？是不是很高啊？

我没有设置采样频率，请问在我的程序中，怎样把采样频率降低点呢？初学请见谅！

harechen

fs 是采样频率
可能是 44000
4.5X44000 约20万点DATA

星海

44000是默认的采样频率？

harechen

fs = 44000 是因為人聽到的聲音在22KHz以內，
由於採樣原理，所以兩倍（大概沒想錯吧）

但是因為採樣目的不同，頻率範圍也許可以調整小一些（可能要看錄音器）。

星海

能帮忙用resample函数帮我修改下采样频率吗，我没有找到相关的使用列子，采样频率设置低点，谢谢

星海

谢谢了，问题解决了，非常感谢！

zhangnan3509

请楼主说一下具体解决办法

星海

clear
[Y,FS,NBITS]=wavread('2.wav');
Z=resample(Y,1,100);
[m,n]=size(Z);
figure(1);
plot(linspace(0,100*m/FS,m),Z);
title('声音信号')
xlabel('t')
ylabel('Sound')
y= linspace(1,2,500);%添加趋势项
imf=emd(Z);
emd_visu(Z,1:length(Z),imf)
[A,f,tt]=hhspectrum(imf);
[im,tt]=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray))
grid on


这样就把采样频率降为原来的100分之一了吧

星海

我往语音信号添加高斯噪声后，进行emd分解，怎么样处理imf分量，可以把噪音除去，再重构回原来的语音信号呢？
高手给点意见！

asdf229955

请问这句：y= linspace(1,2,500);%添加趋势项
是干什么用的呀？是在哪里被处理的？

jinnian

用EMDC函数吧，那个速度快很多的。