I need help on a wavelet transform code

2011 2 17
1 个回答
回答已采纳
2 次查看（30 天）

0 个投票

hi all, i am a masters student and working on speaker identification project. i want to use wavelet transform for feature extraction and could find a code but got an error when i run it... the function “en” in line 35 in the sbc_2.m and the error was as following: ??? Undefined command/function 'en'.
here is the code
energy(1)=en(coef(1:s_no(1)));

请先登录，再回答此问题。

 采纳的回答

zohar
zohar 2011-2-17

0 个投票

Hi Shaikha
Give us somthing about s_no ? (example)

2 个评论
显示 无 隐藏 无

Shaikha Hajri
Shaikha Hajri 2011-2-17
hi zohar,
thnx for ur time. here is the whole code
%****************************************************
% Function to calculate sub-band energy(SBC) parameters
% from enframed signal
% modifided on 17 jan 2008
% *********sarikaya paper algorithm***********
% fs samapling rate 8000
% frame size used 192
%s signal
%
%
%***************************************************
function feature= sbc_2(s,fs)
frames=enframe(s,hamming(192),192);
no_offrames=size(frames);
no_of_frames=no_offrames(:,1);
feature=zeros(no_of_frames,12);
energy=zeros(24,1);
fe_frame=0; % This counter is used to avoid frames with no energy is a sub band
f=statusbar('Extracting SBC Feature');% Creates status bar
for i=1:no_of_frames
energy=ones(24);
f=statusbar((i/no_of_frames),f);% updates status bar
[coef,len]=wavedec(frames(i,:),6,'db3');
s_no=cumsum(len);
%% next two gives the last node in wavelet pacet
energy(1)=en(coef(1:s_no(1)));
energy(2)=en(coef((s_no(1)+1):s_no(2)));
% coef((s_no(1)+1):s_no(2))
%refer diagram in sarikaya paper
[coef1,len1]=wavedec(coef(s_no(2):s_no(3)),1,'db3');
len1=cumsum(len1);
energy(3)=en(coef1(1:len1(1)) );
energy(4)=en(coef1((len1(1)+1):len1(2)));
[coef2,len2]=wavedec(coef((s_no(3)+1):s_no(4)),2,'db3');
len2=cumsum(len2);
energy(5)=en(coef2(1:len2(1)));
energy(6)=en(coef2( (len2(1)+1):len2(2) ) );
[coef3,len3]=wavedec(coef2((len2(2)+1):len2(3)),1,'db3');
len3=cumsum(len3);
energy(7)=en(coef3(1:len3(1)));
energy(8)=en(coef3((len3(1)+1):len3(2)));
% all nodes at level six are over
[coef4,len4]=wavedec(coef(s_no(4):s_no(5)),2,'db3');
len4=cumsum(len4);
energy(9)=en(coef4(1:len4(1)));
energy(10)=en(coef4((len4(1)+1):len4(2)));
[coef5,len5]=wavedec(coef4((len4(2)+1):len4(3)),1,'db3');
len5=cumsum(len5);
energy(11)=en(coef5(1:len5(1)));
energy(12)=en(coef5((len5(1)+1):len5(2)));
[coef7,len7]=wavedec(coef(s_no(5):s_no(6)),3,'db3');
len7=cumsum(len7);
energy(13)=en(coef7(1:len7(1)));
energy(14)=en(coef7((len7(1)+1):len7(2)));
[coef8,len8]=wavedec(coef7((len7(2)+1):len7(3)),1,'db3');
len8=cumsum(len8);
energy(15)=en(coef8(1:len8(1)));
energy(16)=en(coef8((len8(1)+1):len8(2)));
[coef9,len9]=wavedec(coef7((len7(3)+1):len7(4)),2,'db3');
len9=cumsum(len9);
energy(17)=en(coef9(1:len9(1)));
energy(18)=en(coef9((len9(1)+1):len9(2)));
energy(19)=en(coef9((len9(2)+1):len9(3)));
[coef10,len10]=wavedec(coef((s_no(6)+1):s_no(7)),3,'db3');
len10=cumsum(len10);
energy(20)=en(coef10(1:len10(1)));
energy(21)=en(coef10((len10(1)+1):len10(2)));
energy(22)=en(coef10((len10(2)+1):len10(3)));
[coef11,len11]=wavedec(coef10((len(3)+1):len10(4)),1,'db3');
len11=cumsum(len11);
energy(23)=en(coef11(1:len11(1)));
energy(24)=en(coef11((len11(1)+1):len11(2)));
% Taking 12 filter bank equalent
%
%f=rdct(feature); function not woring
% 19 jan DCT modified file in dessai cmtr
if all(energy>0)
fe_frame=fe_frame+1;
log_en=log(energy.*1E+06);
for j=1:12
for k=1:24
feature(fe_frame,j)=feature(fe_frame,j)+log_en(k)*cos((j*(k-0.5)*pi)/24);
end
end
% feature=abs(feature); this needs to be checked
end
end
delete(statusbar)
Walter Roberson
Walter Roberson 2011-2-17
To confirm, you have solved this problem? If so then it would help if you could state what the problem turned out to be.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Continuous Wavelet Transforms 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by