How can I calculate Kurtosis of Sound data?

Question

Serhat Sayim 2021-1-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/707713-how-can-i-calculate-kurtosis-of-sound-data

评论： dpb 2021-1-3

I have sound data of people who have asthma.I made separate plots windows thanks to the start and end indexes of wheezing and non-wheezing sounds. How can I calculate the kurtosis value of each wheeze or nonwheeze windows.

wheezing starting indexes

wheezing ending indexes

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Star Strider 2021-1-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/707713-how-can-i-calculate-kurtosis-of-sound-data#answer_589803

在 MATLAB Online 中打开

Not certain what you want.

Try this:

startidx = [2239                                              
            19576
            30338
            46537
            56438
            69869
            97067
            107264
            123725
            134174];
endidx = [3239
          24791
          31829
          49212
          57662
          72755
          99823
          108708
          126971
          136379];
framelen = endidx - startidx;
fl_mean = mean(framelen);
fl_varv = var(framelen);
fl_kurt = kurtosis(framelen);

.

9 个评论
显示 7更早的评论隐藏 7更早的评论

Star Strider 2021-1-3

在 MATLAB Online 中打开

As always, my pleasure!

Try this for both of them:

for k = 1:numel(startidx)
    ktsis_wheeze(k) = kurtosis(signal(startidx(k):endidx(k)));
end
for k = 1:numel(startidx)-1
    idxrng(k,:) = [endidx(k) startidx(k+1)];                                % Information To Show Indices (Delete)
    ktsis_nonwheeze(k) = kurtosis(signal(endidx(k):startidx(k+1)));
end

The first loop is the same as previously, I just re-named the variable.

See if the second loop does what you want.

(Note that ‘ktsis_nonwheeze’ it is 1 element shorter than the ‘ktsis_wheeze’ vector because I do not know how your data are organised, so I do not know whether to include any other indices.)

dpb 2021-1-3

在 MATLAB Online 中打开

Kin=arrayfun(@(2239,3239)kurtosis(X(2239:3239)),iStart,iEnd);

That's not the code I gave you...just plug your variables in where I told you:

Kin=arrayfun(@(i1,i2)kurtosis(X(i1:i2)),iStart,iEnd);

Use whatever is your array/variable name for the data in place of "X". You plotted the data so you have the variable already, use it again.

请先登录，再进行评论。

Answer 2

dpb 2021-1-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/707713-how-can-i-calculate-kurtosis-of-sound-data#answer_589793

在 MATLAB Online 中打开

Kin=arrayfun(@(i1,i2)kurtosis(X(i1:i2)),iStart,iEnd);

where X are your data and the two indexing arrays above for the "in group" sections.

The "out group" indices can be derived from the above in group values by adjustment of point numbers by one and adding the first, last points of the X vector. The same logic/functional form will then work for those sections.

Will leave as "exercise for Student" to consider the details...

2 个评论
显示无隐藏无

Serhat Sayim 2021-1-3

Let me explain clearly. For example I have a wheeze plot window between 2239 and 3239. How can I calculate the kurtosis value between these values and for this window?

dpb 2021-1-3

编辑：dpb 2021-1-3

That's exactly the code I gave you for in-between the two indices -- you have to do a minimal amount of adjustment of those two arrays to get the start-stop for the other segments from those values.

请先登录，再进行评论。