Skip to content

Экспертиза алгоритма

  1. Разбор интервалов частот (Гц):
[0..1]
[1..2]
[2..4]
[4..8]
[8..16]
[16..32]
[32..64]
[64..128]
[128..256]
[256..512]
[512..1024]
[1024..2048]
[2048..4096]
[4096..8192]
[8192..16384]
[16384..32768]

  1. Из этих 16 интервалов выделяем 2 мастер-интервала: 0..127 (8 интервалов) и 128...32768 (также 8 интервалов);

  2. Мастер-интервал 0..127: это зона частот волновых ритмов мозга;

  3. Мастер-интервал 128...32768: это зона частот музыкально-акустического диапазона (согласуется +/- 200 Гц с частотам октавной системы до пятой октавы);

  4. Для отсчета (интервал 2 секунды исходного файла) получаем 8 значений:

sum(20*log10([0..1]))/1     + sum(20*log10([128..256]))/128
sum(20*log10([1..2]))/1     + sum(20*log10([256..512]))/256
sum(20*log10([2..4]))/2     + sum(20*log10([512..1024]))/512
sum(20*log10([4..8]))/4     + sum(20*log10([1024..2048]))/1024
sum(20*log10([8..16]))/8    + sum(20*log10([2048..4096]))/2048
sum(20*log10([16..32]))/16  + sum(20*log10([4096..8192]))/4096
sum(20*log10([32..64]))/32  + sum(20*log10([8192..16384]))/8192
sum(20*log10([64..128]))/64 + sum(20*log10([16384..32768]))/16384

Замечания

Следует изменить интервалы

[0..3]          - дельта-ритм
[4..5]          - нижний тета-ритм
[6..7]          - верхний тета-ритм
[8..13]         - альфа-ритм
[14..16]        - бета-ритм, 1 power
[17..20]        - бета-ритм, 2 power
[21..28]        - бета-ритм, 3 power
[29..63]        - гамма-ритм
[64..127]       - большая октава
[128..255]      - малая октава
[256..511]      - первая октава
[512..1023]     - вторая октава
[1024..2047]    - третья октава
[2048..4095]    - четвертая октава
[4096..8191]    - пятая октава   
[8192..16384]   - спецэффекты

Матчасть

  1. В настоящей реализации алгоритм вычисляет 1 среднее значение для интервала. При восстановлении в амлитудно-временном представленнии (обратное БПФ) с частотой не более 64 Гц (верхний порог частот волновых ритмов мозга) будет получен сигнал с частотой дискретизации не менее 128 Гц (теорема Котельникова).

  2. Исходя и п.2 для восстановления сигнала длиной 2 секунды подтребуется вектор длиной 256 элементов (требуется проверка на модели MATLAB).

Предложения

  1. Проверить достоверность п.2 из раздела Матчасти.

  2. Заполнить вектор соответствующими данными для 8 интервалов - требуют проверки на модели MATLAB следующие варианты: заполняются равные 8 интервалов вектора, каждый интервал заполняется одним и тем же конкретным значением [+ ИЛИ+ ] заполняется начало вектора в соответсвии с частотными интервалами мозговых ритмов (так же одним и тем же конкретным значением).

  3. Восстановить сигнал мозговых ритмов через обратное БПФ.

Edited by Konstantin Narkhov