Оценка распределения модуля мгновенных значений фонограммы
11.01.2023
Здравствуйте! Получил задание на разработку алгоритма и программы, реализующих оценку распределения модуля мгновенных значений фонограммы.
1) Разработать методику, алгоритм и программу оценки распределния модуля мгновенных значений фонограммы;
2) Разработать алгоритм определения моды распределния.
Может ли кто-нибудь, пожалуйста, помочь с программой в MATLAB, ибо я полный профан в этой области
Лучший ответ
В DSP тулбоксе есть аудиофайлик, его можно использовать для примера.
%%audio = dsp.AudioFileReader('guitartune.wav', ...'SamplesPerFrame',44100);x = audio();%Модули мгновенной амплитуды с вырезанной паузой (нулевыми значениями)x=abs(x(x~=0));%Гистограммаhistogram(x,'Normalization','pdf');hold on%Пусть это будет логнормальное распределениеpd = fitdist(x,'Lognormal');%Можно pdf построить дополнитльно к гистограммеx_vals = 0:0.01:0.5;y = pdf(pd,x_vals);plot(x_vals,y,'LineWidth',2)%Модаmd = mode(x);%Дорисовываем моду на графикеxline(md,'LineWidth',2,'LineStyle','--')hold offОтветы
Your writing is really informative, especially because it's so meaningful and updated. Thanks for sharing this wonderful post! Retro Bowl game
Вы бы сначала алгоритм подробный расписали - какие входные данные, какие формулы, тогда может и подсказать что-то получится.
Хотелось еще бы понять, что имеется в виду под оценкой распределения модулей? Нужно закон распределения получить? Есть ли какие-то сведения о возможном характере распределения?
Фонограмма - это звук какой-то?
Модуль - это на подобие выпрямленного сигнала?
Из того, что написано думается можно гистограмму для начала (ф-ции hist, bar, и т.п.)
А вот алгоритм "реализующий оценку распределения" - это либо что-то очень простое заумно написано, либо реально что-то крутое, но тогда нужна еще информация об этом всем.
В DSP тулбоксе есть аудиофайлик, его можно использовать для примера.
%%audio = dsp.AudioFileReader('guitartune.wav', ...'SamplesPerFrame',44100);x = audio();%Модули мгновенной амплитуды с вырезанной паузой (нулевыми значениями)x=abs(x(x~=0));%Гистограммаhistogram(x,'Normalization','pdf');hold on%Пусть это будет логнормальное распределениеpd = fitdist(x,'Lognormal');%Можно pdf построить дополнитльно к гистограммеx_vals = 0:0.01:0.5;y = pdf(pd,x_vals);plot(x_vals,y,'LineWidth',2)%Модаmd = mode(x);%Дорисовываем моду на графикеxline(md,'LineWidth',2,'LineStyle','--')hold off