How to find spectrum envelope from wav file

Question

Daemian 2014-12-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/168113-how-to-find-spectrum-envelope-from-wav-file

评论： Daemian 2015-1-9

How do I measure the energy change of the spectrum envelope 20 - 40 hz from from a wav file and plot it to graph? my recording is in 8000 sampling rate 16bits mono.

I would like to get something like this

Thanks for great help, wish you all merry Christmas and happy new year :)

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Image Analyst 2014-12-28

编辑：Image Analyst 2014-12-28

Daemian, perhaps you overlooked this big hint given in Star's answer "I don’t have your signal, so I can’t test the filter I designed for you with it." HINT, HINT. What do you think you should do now?

Also, please read this.

请先登录，再进行评论。

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

Answer 1

Star Strider 2014-12-26

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/168113-how-to-find-spectrum-envelope-from-wav-file#answer_163352

在 MATLAB Online 中打开

If you have the Signal Processing Toolbox, such a bandpass filter is easy to design. First (to make things easier) convert your 16-bit signed integer signal (that I call ‘x’ here) to double:

xd = double(x);                         % Convert to ‘double’
Fs = 8000;                              % Sampling Frequency
Fn = Fs/2;                              % Nyquist Frequency
Fpb = [20 40]/Fn;                       % Passband
Fsb = [15 50]/Fn;                       % Stopband
[n,Wn] = buttord(Fpb, Fsb, 1, 10);      % Filter Order, With Rp = 1, Rs = 10
[b,a]  = butter(n,Wn);                  % Create Filter Transfer Function
[sos,g] = tf2sos(b,a);                  % Second-Order-Section Implementation
figure(1)
freqz(b,a)                              % Transfer Function Plot
figure(2)                               % Second-Order-Section Plot
freqz(sos)
yd = filtfilt(sos, g, xd);              % Filter ‘xd’ To Get ‘y’

where ‘xd’ is your .wav signal and ‘yd’ is the filtered output. The filtfilt function will filter both channels of your signal at the same time (assuming your mono signal could have two channels, with the same information in both channels), so this code will work regardless of its having 1 or 2 channels. You can then change ‘yd’ to 16-bits with the int16 function if you want to.

Have fun with this!

Merry Christmas (belatedly) and happy New Year to you, too!

15 个评论
显示 13更早的评论隐藏 13更早的评论

Star Strider 2014-12-26

编辑：Star Strider 2014-12-28

在 MATLAB Online 中打开

My pleasure!

I don’t have your signal, so I can’t test the filter I designed for you with it.

I wasn’t sure what you meant by plot (I thought you wanted to plot the transfer function). Again taking ‘yd’ as your filtered signal, this should work:

yd_len = length(yd);                    % Length Of ‘wav’ File Data
Ts = 1/Fs;                              % Sampling Time
Tv = linspace(0, yd_len-1, yd_len)*Ts;  % Time Vector
figure(3)
plot(Tv, yd)                            % Plot Filter Output

You will have to define what you mean by ‘spectrum envelope’ if this doesn’t do what you want. It looks like you want to add your original signal to the output of the filter, in which situation figure(4) here may be what you want:

figure(4)
plot(Tv, yd+xd)                         % Plot Filter Output

This is just a guess.

[26 Dec 2014 - 15:51]

— EDIT —

[28 Dec 2014 - 14:39]

Does this do what you want? If not, what would you want it to do? Please be as specific as possible, and as detailed as necessary.

I have no idea what you mean by ‘energy change of the spectrum envelope’ and can find no online references to anything by that name. If you have a reference, please post the PDF of it, describing the process and what you want to do.

Star Strider 2014-12-28

在 MATLAB Online 中打开

My pleasure.

The problem is that paper gives no information on the signal processing methods it uses to get those plots. It also mentions the physiological tremor frequency band of 8-12 Hz, but does not mention how it interacts with the 20-40 Hz band it mentions elsewhere.

The only reference I can find that has what seems to be a method is behind a paywall (on the Scribd site ‘VOice Stress Analysis’). I do not have access to it.

I have no experience with voice stress analysis. If you provide (and attach here) the ‘VOice Stress Analysis’ paper or another reference that has a detailed description of the methods it uses (in PDF format), I will see if it provides an algorithm I can code in MATLAB.

Without a method or algorithm describing their approach to voice stress analysis, I cannot help. Meanwhile, see if the filter design I posted provides the information you want. If you decide to add a bandpass filter for the 8-12 Hz band, I would design that filter similar to the previous design, only with these passband and stopband changes:

Fpb = [ 8 12]/Fn;                       % Passband
Fsb = [ 5 15]/Fn;                       % Stopband

Everything else is unchanged. Be sure to check the stability of the ‘sos’ implementation of this filter with the freqz function.

Star Strider 2015-1-1

My pleasure!

Happy New Year to you, too!

Daemian 2015-1-9

index.pdf

Hey Star Strider, I found this pdf talking about spectral envelopes I not too sure is it helpful, can you help me take a look?

Thanks!

请先登录，再进行评论。

Answer 2

Youssef Khmou 2014-12-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/168113-how-to-find-spectrum-envelope-from-wav-file#answer_163469

编辑：Youssef Khmou 2014-12-28

在 MATLAB Online 中打开

Hilbert Transformation is used to obtain the envelope of signal, here is an example :

 The envelope is decreasing exponential, 
Fs=80;
F=10;
t=0:1/Fs:4-1/Fs;
x=exp(-t).*real(exp(j*2*pi*F*t));
figure; plot(t,x);
Y=abs(hilbert(x));
hold on;
plot(t,Y,'r');

 fx=fftshift(abs(fft(x))); fx=fx(floor(end/2:end));
 fY=fftshift(abs(fft(Y))); fY=fY(floor(end/2:end));
 figure; plot(fx); hold on
 plot(fY,'r')

2 个评论
显示无隐藏无

Image Analyst 2014-12-30

Cool - I didn't know hilbert() did that. Does it always get the envelope no matter what x looks like (like how fast it oscillates)? Why does it get a little squirrely around t=4?

Youssef Khmou 2014-12-30

编辑：Youssef Khmou 2014-12-30

According to theory yes, oscillation in borders can be interpreted as Gibbs effect.

请先登录，再进行评论。

How to find spectrum envelope from wav file

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

回答（2 个）

15 个评论
显示 13更早的评论隐藏 13更早的评论

2 个评论
显示无隐藏无

另请参阅

类别

标签

Community Treasure Hunt

How to find spectrum envelope from wav file

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

回答（2 个）

15 个评论 显示 13更早的评论隐藏 13更早的评论

2 个评论 显示 无隐藏 无

另请参阅

类别

标签

Community Treasure Hunt

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

15 个评论
显示 13更早的评论隐藏 13更早的评论

2 个评论
显示无隐藏无