2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Вейвлет анализ в Matlab
Сообщение15.06.2011, 15:04 
Поставлена задача: преобразовать вейвлетом Морле часть ЭЭГ мозга человека.
Проблема возникает при построении спектрограммы преобразованного сигнала (функция cwt рисует его). Отсчет по оси X идет в отсчетах, а необходимо в секундах (5003 отсчета - 10,006 сек). А по оси Y считает в собственных "a", задаваемых в последней строке "1:4:400", а нужно в Герцах:
Используется синтаксис Matlab M
Fr(i) = BFr*SR/a;

, где BFr - базовая частота вейвлета (BFr = centfrq(w);, SR = 500, а "a" - 1:4:400.
Используется синтаксис Matlab M
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

В получившемся графике получается вся необходимая картина, кроме одного НО: Все подписи по осям и отсчеты по этим осям должны быть в тех единицах, что указаны выше (Просто нужно подменить цифры, оставив график не тронутым)
Если кто сможет объснить, как начертить правильно - тоже хорошо.
Заранее спасибо

 
 
 [ 1 сообщение ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group