[Dmitry_DM]

Здравствуйте. Вот уже несколько дней ломаю голову над поставленной задачей. В общем: есть значения случайной величины - сэмплы из wav файла. Надо подсчитывать кол-во значений случайной величины, которые попадают в j-ый интервал (N[j]).
Nj - абсолютная частота для j-го интервала;
j - номер интервала, j меняется от 1 до m;
m - кол-во интервалов;
Считается m так:

Код:

m:=trunc(3.3*log10(NumSamples)+1);

delX - длина интервала m (кол NumSamples, которые попали в m)
Считается delX так:

Код:

delX:=NumSamples / m;

То есть это все семплы разбились на интервалы, например NumSamples=1079, а m=11, delX=98. Получается 98*11 = 1078.
NumSamples - кол-во значений случайной величины;
С этим все ясно.
Вернемся к задаче.. (выделена красным)
Для этого надо создать цикл по i , которое изменять от 1 до NumSamples. В цикле проверять, в какой интервал попало i -е значение и при этом изменять соответствующий счетчик - абсолютную частоту.
Вот такой у меня цикл для всех операций с сэмплами:

Код:

   for I := 1 to NumSamples do
      begin
        for J := 1 to WavInfo.WaveFormat.nChannels do        // так как бывает 2 канала. 
           begin
////////////////////////Здесь все происходит
           end;
      end;

Вот формула: (обведена красным)

Но пока Yj отложим на потом.
И на этом я застрял. Не знаю как это все осмыслить и как должно считаться. Как найти Nj??? Еще добавлю, что Nj по количеству должно быть столько же, сколько и m. Пока все. Надеюсь на вашу помощь. Заранее благодарен! Для меня это очень важно..

[poli-smen]

Сообщение от Dmitry_DM

В общем: есть значения случайной величины - сэмплы из wav файла. Надо подсчитывать кол-во значений случайной величины, которые попадают в j-ый интервал (N[j]).

Написано очень много, единственно только не написано по какому критерию необходимо разбрасывать сэмплы по интервалам.

Сообщение от Dmitry_DM

То есть это все семплы разбились на интервалы, например NumSamples=1079, а m=11, delX=98. Получается 98*11 = 1078.

А что с 1079-ым сэмплом? Он не учитывается?

[Dmitry_DM]

Сообщение от poli-smen

Написано очень много, единственно только не написано по какому критерию необходимо разбрасывать сэмплы по интервалам.

Да вот это я не понимаю.. Мне кажется просто по-порядку, как идут в файле.
Можно как-то сделать пример и посмотреть что вышло? И еще - Nj должны быть ТОЛЬКО положительного значения. И их кол-во равняться кол m столбиков.

Сообщение от poli-smen

А что с 1079-ым сэмплом? Он не учитывается?

Это я сам не знаю. Всегда в первом m интервале на 1 семпл меньше, чем в остальных.
Ниже мой юнит, весь проект нету смысла выкладывать.

[poli-smen]

Сообщение от Dmitry_DM

Да вот это я не понимаю..

Так как же можно решать задачу, если не известно её условие? (не метод решения, а именно условие)

Сообщение от Dmitry_DM

Мне кажется просто по-порядку, как идут в файле.

Ну тогда к примеру для случая NumSamples=1079 получим, что в нулевой интервал попадают сэмплы с 0 по 97, в первый интервал сэмплы с 98 по 195, во второй интервал сэмплы с 196 по 293 и т.д. И если подсчитывать кол-во значений случайной величины, которые попадают в j-ый интервал (N[j]), то получим что в нулевом интервале будет 98 сэмплов (т.е. N[0]=98) и в первом интервале будет 98 сэмплов (т.е. N[1]=98) и во втором интервале будет 98 сэмплов (т.е. N[2]=98) и т.д.

[Dmitry_DM]

Сообщение от poli-smen

Так как же можно решать задачу, если не известно её условие? (не метод решения, а именно условие)
Ну тогда к примеру для случая NumSamples=1079 получим, что в нулевой интервал попадают сэмплы с 0 по 97, в первый интервал сэмплы с 98 по 195, во второй интервал сэмплы с 196 по 293 и т.д. И если подсчитывать кол-во значений случайной величины, которые попадают в j-ый интервал (N[j]), то получим что в нулевом интервале будет 98 сэмплов (т.е. N[0]=98) и в первом интервале будет 98 сэмплов (т.е. N[1]=98) и во втором интервале будет 98 сэмплов (т.е. N[2]=98) и т.д.

Да в том то и дело, что условие точно так и не стало мне понятно.. Вот я и спросил, есть ли какие-то варианты подсчетов.. А 98(данный пример) это дельта Х(delX). Это длина одного столбца. Как еще можно что-нибудь сделать, при этом использовать цикл, и что бы значения были разными и маленькими, потому что для построения гистограммы нужно откладывать границы (это по X). Левая граница первого интервала - это минимальное значение случайной величины, когда сделать нормировку(то есть разделить мин и макс на макс амплитуду) (оно всегда в районе (-1)) , правая граница первого интервала - это минимальное значение случайной величины плюс длина интервала m, m=3,3log*NumSamples+1, где NumSamples - кол-во значений случайной величины. Таким образом все столбики, (при NumSamples=1079 их 11) должны по Х откладываться в диапазоне от -1 до 1. Из этого выплывает, что в -1 до 1 не могут влезть столбики, шириной 98. Может произвести нормировку всех сэмплов? (поделить их в цикле на макс ампдитуду). Я так подумал, потому что на самом первом графике, на рис., было -128 + 128, а потом зачеркнулось и поставилось -1 +1. И логично предположить, что Y тоже не будут огромными. Гистограмма должна получится, как на прикрепленном мною рисунке, чуть выше формулы в прямоугольнике. Так же с той картинки видно, что 1<=j<=m. Значит j и m это разные значения. И еще видно, что один столбец выше, другой - ниже. Не все одинаковые. Это все что известно.. Вот так я ломаю голову уже долго..

[Dmitry_DM]

Ну что, poli-smen, есть идеи?