Диссертация (1143463), страница 10
Текст из файла (страница 10)
2.10:Пример вейвлет-преобразования динамики трансмембранного потенциала при изменении интенсивности света: Верхняя панель- хаотическая динамика потенциала, безразмерный потенциал vs время; модуль вейвлет-преобразования (средняя панель) и соответствующие максимумы модуля как функция медленно меняющейся интенсивности (нижняя панель)ся, так как траектория вышла на устойчивый предельный цикл. Периодглавной гармоники достаточно слабо увеличивается до = 19.5 по сравнению с начальным участком.
В то же время качественно вид колебанийменяется, а именно, они приобретают более релаксационный характер,то есть медленное движение перемежается быстрым. Это хорошо видно на диаграмме полного вейвлет преобразования(Рис. 2.8, средняя панель) где острые спайки соответствуют релаксаторам. Надо отметить,что ненулевая неопределенность вейвлета Морле в частотном пространстве позволяет определить единственный (главный) период осцилляций(сплошная линия на нижней панели рис. 2.8).
Это довольно серьезноепреимущество перед методом Фурье-преобразования, который имеет нулевую неопределенность, что приводит к детектированию более высоких67Рис. 2.11:Двухпериодические колебания трансмембранного потенциала (верхняяпанель), модуль вейвлет-преобразования (средняя панель) и соответствующие максимумы модуля как функция медленно меняющейся интенсивности (нижняя панель)периодических компонент, ограничивая тем самым бифуркационный анализ.
Последняя часть на рис.2.8 показывает возникновение новых паттернов при > 0.7024775, где линия периодов становится неустойчивой.Характерный период для релаксационных колебаний между значениями = 19.5 и более высокими значениями.При дальнейшем уменьшении ключевого параметра до = 0.70155появляются максимумы, соответствующие периодам = 19.5 и =21.25. Такие дублеты можно объяснить релаксационным характером колебаний, однако форма колебаний близка к синусоидальной в интервале = 0.70155 < < 0.70109 и главным периодом остается становится наибольший из двух = 21.25.
Однако при достижении значения < 0.70109 происходит бифуркация удвоения периода: представленныеоба модуль CWT и его максимумы (рис. 2.9, средняя и нижняя панели)68Рис. 2.12: Переход двух-периодического аттрактора к однопериодическому релаксационному. Вехняя панель: динамика потенциала при изменении интенсивности света,модуль вейвлет-преобразования (средняя панель) и соответствующие максимумы модуля как функция медленно меняющейся интенсивности (нижняя панель)позволяют довольно подробно проанализировать данный процесс. Видно, что тут возникают два главных периода = 21.25 и = 42.5 ссопровождающими их пунктирными линиями, что указывает на релаксационные свойства колебаний, как и было в предыдущем случае. Такжеможно увидеть, что второй максимум не возникает сразу, как видно изрис. 2.9, средняя панель, максимум = 21.25 расщепляется в точке < 0.70109.
Расщепление растет до устойчивых значений периода, соответствующего устойчивой величине удвоенного периода. Такое поведениеуказывает на “мягкую” потерю устойчивости в бифуркационной точке.Сосуществование малого и удвоенного периодов сосуществует в оченьузком интервале 0.7008381 < < 0.70041.При интенсивности равной 0.70041 поведение системы приобретаетболее сложный характер: решение состоит из нерегулярных “прыжков”69Рис. 2.13:Финальный этап эволюции колбаний – релаксационные однопериодические колебания (верхняя панель), модуль их вейвлет-преобразования (средняя панель) и и соответствующие максимумы модуля как функция медленно меняющейсяинтенсивности (нижняя панель)между колебаниями с периодами = 21.25 и = 42.5, как показанона рис. 2.10, средняя и нижняя панели.
Такое поведение может соответствовать хаотической перемежаемости. Дальнейшее уменьшение интенсивности приводит к стабилизации колебаний, и при < 0.69904 полныйпериод разделяется на две периодические компоненты. Надо отметить,что характер движения носит релаксационный характер, на что указывает наличие множества точек максимумов внизу (см. рис. 2.10, нижняяпанель). По сравнению с рис. 2.9, средняя панель, они возникают благодаря форме индивидуальных колебаний, характерной для крайне сильныхрелаксаторов, то есть с почти вертикальным графиком быстрой фазыколебания.Такой устойчивый двухпериодический аттрактор существует в довольно широком диапазоне параметра, см.
рис. 2.11, до значений =70Рис. 2.14: Фазовые портреты, соответствующие Рис. 2.7: трансмембранный потенциал vs. производная потенциала по времени (левая панель) и действительная часть vs.мнимой части модуля вейвлета, соответствующего максимальному периоду (праваяпанель).0.69445. Ниже этого значения происходит новая бифуркация: удвоенныйпериод исчезает, как показано на рис.2.12, средняя и нижняя панели.Анализ полного вейвлет-преобразования показывает, что этот переходосуществляется до величины = 0.69388.При дальнейшем понижении интенсивности до = 0, как показанона рис. 2.13-2.12, преобразование периодических паттернов происходитпостепенно: от гармонических до слабо релаксационных.
используя максимумы из рис. 2.12 (нижняя панель), и полную диаграмму из рис. 2.13(средняя панель), движение может быть охарактеризовано как почтигармонические колебания с двумя ведущими Фурье-компонентами.Анализ нелинейных динамических режимов, рассмотренных выше,может быть также проиллюстрирован сравнением фазовых портретовдинамики потенциала и вейвлетных компонент (рис. 2.14–2.20). В част71Рис.
2.15: Фазовые портреты, соответствующие Рис. 2.8: трансмембранный потенциал vs. производная потенциала по времени (левая панель) и действительная часть vs.мнимой части модуля вейвлета, соответствующего максимальному периоду (праваяпанель).ности, на графиках изображены действительная часть вейвлет-преобразования как функция его комплексно-сопряженной части в шкале значений, соответствующих максимумам определенным выше. Например,рис. 2.14 показывает интервал роста гармонических колебаний соответствующему на рис.
2.7. Можно видеть, что фазовый портрет представля˙ ) собой раскручивающующая эллипсоидальную спираль. Тает ((),()же самая картина наблюдается и для фазового портрета вейвлет-компонент. Полуоси обоих предельных эллипсов почти те же самые, различиятолько наблюдаются в центральной части, что объяснимо тем, что вейвлет преобразование “уничтожает” постоянную часть сигнала, то есть72Рис. 2.16: Фазовые портреты, соответствующие Рис. 2.9: трансмембранный потенциал vs. производная потенциала по времени (левая панель) и действительная часть vs.мнимой части модуля вейвлета, соответствующего максимальному периоду (праваяпанель).сдвигает абсциссу к нулю. Мнимая часть соответствует скорости ˙ , которая есть первая производная, и она не сдвигается.Следующий момент, требующий обсуждения – это предельныйцикл, который также качественно воспроизводится см.
рис. 2.15. Здесьесть только количественные различия. Это связано с тем, что фазоваявейвлет-траектория, представленая на рис. 2.15 (правая панель) показывает динамику только одной вейвлет-компоненты с фиксированнымпериодом. Это эквивалентно гармоническим колебаниям с тем же самым периодом, что и у релаксационных колебаний, где быстрые интервалы движений соответствуют прямым линиям динамики потенциала нарис.
2.15 (левая панель). Как обсуждалось уже выше, быстрое движениеприводит к появлению спайков на диаграмме полного вейвлет-преобразования. Однако, анализируя фазовые портреты (рис. 2.15), можно уви73Рис. 2.17:Фазовые портреты, соответствующие Рис. 2.10: трансмембранный потенциал vs. производная потенциала по времени (левая панель) и действительная частьvs.
мнимой части модуля вейвлета, соответствующего максимальному периоду (правая панель).деть, что они только влияют на форму и размер траектории, не меняяпри этом топологии. То же самое демонстрируют и последующие частианализируемого поведения. Изображая фазовые портреты вейвлетов, используя двухпериодические слагаемые с периодами, соответствующимидвум главным максимумам модуля вейвлет-преобразования, как показано на рис. 2.9(нижняя панель), можно видеть, что качественные топологические портреты динамики потенциала (рис.
2.16, левая) и вейвлетгенерируемых фазовых кривых(рис. 2.16, справа) совпадают с высокойстепенью точности. Оба портрета содержат три плотно намотанных цикла, разделенных областями с “редкими” траекториями. Надо отметить,что топология неоднородности хорошо совпадает как на вейвлетном, таки на истинном фазовых портретах.Типичен вид хаотического поведения представленного на соот74Рис. 2.18:Фазовые портреты, соответствующие Рис.
2.11: трансмембранный потенциал vs. производная потенциала по времени (левая панель) и действительная частьvs. мнимой части модуля вейвлета, соответствующего максимальному периоду (правая панель).ветствующих фазовых портретах см. рис. 2.17, который есть аналогрис. 2.10. После стабилизации пережаемости, двухпериодические релаксационные колебания, описанные на рис. 2.11, имеют типичныйвид: двухпетлевые фазовые кривые, представленные на рис.
2.18,˙ ) и (Re ((1 , ) + (2 , )),которые идентичны в обоих ((), ()−Im ((1 , ) + (2 , ))).“Мягкое исчезновение” удвоенного периода (Fig. 2.12) ясно представлено постепенным слипанием фазовых кривых на рис. 2.19. В итоге движение приобретает характер неэллиптического предельного цикла (рис. 2.20, левая панель). Вейвлет-преобразование показывает присутствие двух главных гармоник в колебаниях. Однако, высшая из нихвоздействует только на форму фазового портрета, здесь нет второй петли, что приводит к заключению, что эта часть решения может быть75Рис.