How can I calculate Kurtosis of Sound data?
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
Serhat Sayim
am 3 Jan. 2021
Kommentiert: dpb
am 3 Jan. 2021
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
2239
19576
30338
46537
56438
69869
97067
107264
123725
134174
wheezing ending indexes
3239
24791
31829
49212
57662
72755
99823
108708
126971
136379
0 Kommentare
Akzeptierte Antwort
Star Strider
am 3 Jan. 2021
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 Kommentare
Star Strider
am 3 Jan. 2021
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
am 3 Jan. 2021
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.
Weitere Antworten (1)
dpb
am 3 Jan. 2021
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 Kommentare
Siehe auch
Kategorien
Mehr zu Graphics Performance finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!