Поставлена задача: преобразовать вейвлетом Морле часть ЭЭГ мозга человека.
Проблема возникает при построении спектрограммы преобразованного сигнала (функция cwt рисует его). Отсчет по оси X идет в отсчетах, а необходимо в секундах (5003 отсчета - 10,006 сек). А по оси Y считает в собственных "a", задаваемых в последней строке "1:4:400", а нужно в Герцах:
, где BFr - базовая частота вейвлета (BFr = centfrq(w);, SR = 500, а "a" - 1:4:400.
load signals.txt;
l = input ('Введите номер датчика от 1 до 21 ');
t = 0;
w = 'Morl';
BFr = centfrq(w);
for i=1:5003
sensor(i) = signals(i,l);
time(i) = t;
t = t+0.002;
end;
Sc = cwt(sensor,1:4:400,w,'plot');
Сигнал:
http://narod.ru/disk/16068707001/signals.txt.htmlВ получившемся графике получается вся необходимая картина, кроме одного НО: Все подписи по осям и отсчеты по этим осям должны быть в тех единицах, что указаны выше
(Просто нужно подменить цифры, оставив график не тронутым)Если кто сможет объснить, как начертить правильно - тоже хорошо.
Заранее спасибо