CC BY
46
Радиофизика
Вестник Нижегородского университета им. Н.И. Лобачевского, 2010, № 4 (1), с. 62-66
УДК 537.86:519.2
ВЛИЯНИЕ ФОРМЫ ПОТЕНЦИАЛА НА СВОЙСТВА СТОХАСТИЧЕСКОГО РЕЗОНАНСА В МОНОСТАБИЛЬНЫХ СИСТЕМАХ
© 2010 г. Н.В. Агудов, А.В. Кричигин
Нижегородский госуниверситет им. Н.И. Лобачевского [email protected]
Поступила в редакцию 02.12.2009
Рассматривается нелинейная инерционная моностабильная система, на вход которой поступает аддитивная смесь гармонического сигнала и белого шума. При определенных параметрах усиление мощности сигнала и отношение сигнал-шум на выходе системы могут быть немонотонными функциями интенсивности входного шума, то есть возможно возникновение стохастического резонанса и антирезонанса. На основании исследования двух систем, описываемых кусочно-линейными моностабильны-ми потенциалами, проводится анализ влияния формы потенциала на поведение усиления мощности сигнала и отношения сигнал-шум.
Ключевые слова: стохастический резонанс, моностабильная система, броуновское движение.
Введение
С точки зрения классической радиотехники принято считать, что наличие шума в системе только ухудшает ее характеристики. Исследования, проведенные в последние годы, показывают, что шум в нелинейных системах может стать причиной новых явлений, наблюдение которых без шума принципиально невозможно (см., напр., [1]). Наличие флуктуаций в нелинейных системах может не только приводить к отклонениям от среднего значения, но и играть конструктивную роль, вызывая увеличение упорядоченности в системах [2], подавление внутренних шумов с помощью внешнего шумового сигнала [3], синхронизацию фаз в системах с несколькими степенями свободы [4]. Кроме того в последние годы растет интерес к экспериментальным и теоретическим исследованиям индуцированных шумом эффектов и явлений в биологических системах (см., напр., [5]).
Внешнее воздействие, в зависимости от вида, величины и способа внесения, а также в зависимости от состояния системы, может приводить к различным явлениям (например, индуцированный шумом хаос [6], стохастический резонанс [7] и когерентный резонанс [8]). Самым простым и наглядным из них является стохастический резонанс. Явление стохастического резонанса представляет собой увеличение отклика нелинейной бистабильной системы на слабый внешний периодический сигнал при добавлении аддитивного шума определенной оптимальной интенсивности. При этом внешнее
периодическое воздействие настолько мало, что само по себе не может переводить бистабильную систему из одного метастабильного состояния в другое. Однако интенсивность шума и период сигнала можно подобрать таким образом, чтобы они «помогали» друг другу при таких переходах. Это приводит к тому, что отношение сигнал-шум (ОСШ) на выходе такой системы имеет отчетливо выраженный максимум при некотором оптимальном уровне шума на входе. Именно наличие максимума в функции ОСШ от входной интенсивности шума и является классическим проявлением стохастического резонанса. При этом анализ множества бистабильных систем показал, что немонотонное поведение ОСШ сопровождается аналогичным немонотонным поведением функции усиления мощности сигнала (УМС) в системе от интенсивности шума на входе [9].
Наряду с классическим стохастическим резонансом в бистабильных инерционных системах в последнее время возрос интерес к данному явлению в различных других системах [10-12], и в частности в моностабильных [13-18]. В работах [13, 14] рассматривались системы, управляемые мультипликативным шумом наряду с аддитивным. В работе [15] рассмотрено не ОСШ, а УМС в зависимости от интенсивности аддитивного шума на входе моностабиль-ной системы и предложено иное определение стохастического резонанса. Стохастическим резонансом назван эффект, при котором УМС имеет максимум. Кроме этого, было высказано предположение о возможности существования
Жх
Ж
ЖФ(х)
Жх
+ 5(0 + ^(0 ,
(1)
Рис. 1. Жестко-мягкий потенциальный профиль
стохастического антирезонанса, то есть наличия минимума такой зависимости. В работах [16-18] рассмотрены нелинейные инерционные моностабильные системы, на вход которых поступает аддитивная смесь гармонического сигнала и шума. На основе теории линейного отклика были впервые получены зависимости ОСШ на выходе системы как функции интенсивности входного шума. Обнаружено, что для некоторых моностабильных систем существуют режимы работы, при которых наблюдается явление стохастического резонанса в классическом понимании. Отмечено, что механизм стохастического резонанса в моноста-бильной системе принципиально отличен от механизма классического стохастического резонанса, так как не существует потенциального барьера, разделяющего метастабильные состояния. Также была аналитически подтверждена возможность существования в моноста-бильных системах максимума или минимума УМС как функции интенсивности входного шума, что тоже отличается от того, что ранее наблюдалось для бистабильных систем. Вместе с тем необходимо отметить, что хотя различие стохастического резонанса в бистабильных и моностабильных системах неоднократно отмечалось, сами свойства этого явления в мо-ностабильных системах исследованы недостаточно.
В данной работе на примере двух видов мо-ностабильных потенциалов исследуются некоторые свойства стохастического резонанса в моностабильных системах. В частности, анализируется влияние формы моностабильного потенциального профиля на возникновение немонотонного поведения функций ОСШ и УМС.
Постановка задачи
В данной работе рассматривается инерционная моностабильная система, возбуждаемая аддитивной смесью сигнала и шума, описываемая уравнением Ланжевена:
где ^(?) - белый гауссовский шум: <^(?)> = 0, <^(?) ^(?+х)> = 2ц5(т), Ф(х) - постоянный во времени потенциальный профиль, характеризующий саму систему, 5^) - входной сигнал, х(?) -выходной сигнал. В качестве входного сигнала рассмотрим гармонический:
5(0 = А СО8(ю0 t) .
В настоящей работе исследуются системы, описываемые моностабильными кусочно-линейными потенциальными профилями, изображенными на рис. 1, 2. Отметим, что основное отличие данных потенциалов заключено в различных формах профилей вблизи минимума потенциалов. Следовательно, при больших значениях интенсивности входного шума характеристики обеих систем должны совпадать.
Потенциальный профиль, изображенный на рис. 1, назовем жестко-мягким, так как он состоит из жесткого и мягкого потенциала:
0, |х| < Ь,
ф(х)=<
к(X - Ь)+ Н, |х| > Ь.
Мягкость или жесткость потенциала означает, что при удалении от состояния равновесия сила, возвращающая систему обратно в минимум потенциала, растет медленнее или быстрее, чем в линейной системе (1), соответствующей параболическому потенциальному профилю Ф(х) = ах2. Потенциальный профиль, изображенный на рис. 2, назовем мягко-мягким, так как он составлен из двух мягких потенциалов с разными наклонами потенциальных стенок:
ф(х) =
Н
—|х|, |х| < ь,
к(X - ь)+ Н, |х| > Ь.
Если входной сигнал мал (А << Ь/ц), то, согласно теории линейного отклика, спектр выходного процесса х(^ будет состоять из шумовой платформы Ох(о) и слагаемого, соответст-
вующего гармоническои компоненте на частоте входного сигнала [9]:
а 2
Бх (ю) = в* (ю) +—[8(ю - ®о)+5(ю + ®о)],
где а - амплитуда сигнала на выходе.
В данноИ работе рассматриваются следующие характеристики:
2
УМС
п = -
2
А
и ОСШ на выходе системы
2
а
Р =
26* (ю)'
(2)
(3)
Методы решения
В работах [16, 17] предложен метод отыскания оценок функций (2) и (3) для систем вида (1), описываемых произвольными потенциальными профилями Ф(х). В частности, этот метод позволяет получить соответствующие выражения и для мо-ностабильных систем, рассматриваемых в данной работе. В случае жестко-мягкой системы (рис. 1):
25 Ь4С4 (д)
1
9(1 + УеР) (5ду2Сі(д)) + (ь2С2(д)юо) ’
(4)
Р =
5 пА 2 С14 (д )
6 дС 2 (д)(і +Iе р)
(5)
С1(д)=у 3ер + 3(2 + 2у + у2),
С2 (д) = 2у5ер + 15у (1 + у )2 е~р +
+ 5(15 + 12у + 3у2 + 2у3 + 2у4); в случае мягко-мягкоИ системы (рис. 2):
П =
д 214р2Сз4 (д)
(в - у (1 - ер| (дСз (д))2 + (ьС (дК }
у 1 - ^ уд^з
р п А 2ву С з2 (д )
2 дС 2 (д - у(1 - е С (д)= 2 + 2У + У2 - 2(1 - ев)+ 2Р + Р
(6)
(7)
С4 (д)=
у
5 + 4у + у2 1 + у
в3
у5 ру2
2 + 2у + 3у2 2(1 + р + у + ру)е_р
Р3у2
Р3 У2
,(1+у)2(1 -е_р)-
РУ 4
4 + 4р + р3 - 5еР + (1 + р)2 е-р р5
где в = Н/д, у = кЬ/д.
Отметим, что наряду с оценками (4) - (7) в работах [16, 18] были вычислены также функции п и р для данных потенциальных профилей на основе точных выражений лапласовских образов функций корреляций соответствующих случайных процессов. Здесь эти выражения не приводятся в силу их громоздкости, но они использованы для верификации результатов анализа оценок (4) - (7).
Анализ характеристик моностабильных систем
Как было показано в работах [16-18], в рассматриваемых системах существуют режимы, при которых функции УМС и ОСШ от интенсивности шума на входе могут быть монотонными либо немонотонными и иметь минимум и/или максимум. В данной работе впервые получены кривые значений параметров, при которых в системах наступает смена режимов поведения этих функций. Данные кривые построены по точкам, при этом каждая точка рассчитывается отдельно. Необходимо сразу отметить, что кривые, построенные данным образом, достаточно близки к прямым линиям, отклонения от прямых линий составляет не более 1%.
На рис. 3 приведена плоскость параметров жестко-мягкой системы: по оси абсцисс отложена глубина потенциальной ямы Н, а по оси ординат наклон к внешних стенок потенциального профиля. Изображены линии смены режимов работы для различных значений ширины потенциальной ямы: кривая 1 соответствует Ь = 1, кривая 2 - Ь = 2. Область выше каждой линии соответствует монотонному поведению ОСШ (5). Область ниже отвечает режиму работы, при котором в системе имеет место явление классического стохастического резонанса, то есть наличие максимума ОСШ (5). Таким образом, кривые на рисунке соответствуют параметрам, при которых ОСШ (5) имеет точку перегиба.
Как видно из приведенного рисунка, в жестко-мягкой системе могут существовать два режима работы: первый соответствует минимуму УМС (стохастический антирезонанс - по определению, данному в [15]) и монотонному спаданию ОСШ, а второй соответствует минимуму УМС и, одновременно, максимуму ОСШ («классический» стохастический резонанс). По графику эмпирически построим выражение «прямой» смены режимов работы:
к * 0.2682Н.
I
2
+
2 4 Л 8 10
Рис. 3. Области режимов работы жестко-мягкой системы
В мягко-мягкой системе в отличие от жестко-мягкой в зависимости от параметров имеет место смена режимов не только по фунции ОСШ, но и по функци УМС. На рис. 4 изображена плоскость параметров, на которой показаны эти режимы работы: по оси абсцисс отложена глубина потенциальной ямы Н, по оси ординат - наклон стенок основного потенциала к, при этом Ь = 1. Пунктирная линия соответствует границам режимов (изменению функции УМС (6) от интенсивности входного шума). Область ниже линии соответствует монотонному, а область выше линии - немонотонному поведению данной функции. Примерное выражение этой линии:
к * 3.268Н.
Ь
Сплошная линия на рис. 4 соответствует смене режимов работы (изменению поведения ОСШ (7) как функции от интенсивности шума на входе). При этом область выше линии соответствует монотонному поведению ОСШ (7), а область ниже - стохастическому резонансу, то есть наличию максимума ОСШ (7). Примерное выражение «прямой» смены режимов:
к * 0.1613Н .
Ь
Как видно из приведенного рисунка, в мягко-мягкой системе существуют три режима работы: первый соответствует немонотонному поведению УМС, но монотонному спаданию ОСШ, при втором режиме обе функции ведут себя монотонно, а третий режим соответствует монотонному росту УМС и, одновременно, наличию максимума ОСШ («классическому» стохастическому резонансу).
Анализируя рис. 3 и 4, можно отметить, что при одинаковом коэффициенте наклона внешних стенок потенциального профиля к явление «классического» стохастического резонанса в обеих системах наступает при близких значениях высоты и ширины потенциальной ямы. Сле-
Рис. 4. Области режимов работы мягко-мягкой системы
довательно, форма потенциальной ямы в данном случае значительной роли не играет, важны только ее характерные размеры. С другой стороны, форма потенциальной ямы существенно влияет на функцию УМС. Так, при жесткой потенциальной яме УМС ведет себя немонотонно при любых значениях параметров системы (рис. 3), а при мягкой потенциальной яме - не при любых (рис. 4).
Заключение
В данной работе на основе методов, предложенных и реализованных в работах [16-18], проанализированы выходные характеристики двух моностабильных систем: УМС и ОСШ как функции интенсивности входного шума. Рассмотренные моностабильные системы отличаются лишь формой потенциальной ямы вблизи минимума потенциала. Для каждой из них на плоскостях параметров построены границы смены режимов работы, то есть значения параметров, при которых указанные функции имеют точку перегиба.
Рассмотренные примеры впервые демонстрируют, что свойства стохастического резонанса в моностабильных системах имеют намного большее разнообразие по сравнению с классическим вариантом этого явления, наблюдаемым в бистабильных системах, в которых два локально устойчивых состояния разделены потенциальным барьером. В классическом явлении стохастического резонанса немонотонное поведение функций УМС и ОСШ происходит одновременно, в одной и той же области параметров системы. Поэтому некоторые авторы (как, например, в [15]) стали отождествлять в принципе разные определения стохастического резонанса: немонотонность
функции УМС и немонотонность функции ОСШ. В данной работе впервые показано, что в моностабильных системах такого однообразия не наблюдается и обе функции могут иметь
максимум (и даже минимум) в различных областях параметров системы.
Показано, что в жестко-жесткой системе могут существовать два принципиально различных режима работы, а в мягко-мягкой системе -три. Сделан вывод о том, что наличие явления стохастического резонанса определяется не формой потенциальной ямы, а ее характерными размерами: глубиной и шириной. Однако форма потенциальной ямы является критической при изменении режимов работы, связанных с немонотонностью поведения УМС.
Работа выполнена при поддержке РФФИ (грант № 08-02-01259-а).
Список литературы
1. Хорстхемке В., Лефевр Р. Индуцированные шумом переходы. М.: Мир, 1987.
2. Гинзбург С.Л., Пустовойт М.А. Индуцированная шумом сверхчувствительность к слабым сигналам // Письма в ЖЭТФ. 1998. Т. 67. № 8. С. 592-596.
3. Vilar J., Rubi J. // Phys. Rev. Lett. 2001. V. 86. P. 950.
4. Горностырев Ю.Н., Ждахин Д.И., Кацнельсон М.И., Трефилов А.В. Стохастический резонанс между предельными циклами. Пружинный маятник в термостате // Письма в ЖЭТФ. 1999. Т. 69. № 8. С. 585-589.
5. Spagnolo B., Spezia S., Curcio L. et al. Noise effects in two different biological systems // Eur. Phys. J. B. 2009. V. 69. P. 133-146.
6. Gao J.B., Tung W., Rao N. // Phys. Rev. Lett. 2002. V. 89. № 25. P. 254101.
7. Анищенко В.С., Нейман А.Б., Мосс Ф., Ши-манский-Гайер Л. // Успехи физических наук. 1999. Т. 169. № 1. С. 7-38.
8. Lai Y.C., Liu Z. // Phys. Rev. Lett. 2001. V. 86. P. 4737.
9. Anishchenko V.S., Astakhov V.V., Neiman A.B., et al. Nonlinear dynamics of chaotic and stochastic systems. Berlin: Springer, 2002.
10. Berdichevsky V., Gitterman M. // Phys. Rev. E. 1999. V. 60. P. 1494.
11. Stocks N.G. // Phys. Rev. Lett. 2000. V. 84. P. 2310.
12. Kosko B., Mitaim S. // Phys. Rev. E. 2004. V. 70. P. 031911.
13. Li J. // Phys. Rev. E. 2002. V. 66. P. 031104.
14. Guo F., Zhou Y.R., Jiang S.Q., Gu T.X. // J. Phys. A: Math. Gen. 2006. V. 39. P. 13861-13868.
15. Evstigneev M., Reimann P., Pankov V., Prince R.H. // Europhys. Lett. 2004. V. 65. № 1. P. 7-12.
16. Агудов Н.В., Кричигин А.В. Стохастический резонанс и антирезонанс в моностабильных системах // Известия вузов. Радиофизика. 2008. Т. 51. № 10. С. 899-913.
17. Агудов Н.В., Кричигин А.В. Преобразование сигнала и белого шума нелинейной инерционной системой // Актуальные проблемы статистической радиофизики (Малаховский сборник). 2006. Т. 5. С. 103-135.
18. Agudov N., Krichigin A. Stochastic resonance
in an overdamped monostable system // arXiv.org, Physics, Condensed Matter. 2008. 0808.2095v1
http://arxiv.org/abs/0808.2095v1.
INFLUENCE OF THE POTENTIAL PROFILE SHAPE ON THE PROPERTIES OF STOCHASTIC RESONANCE IN MONOSTABLE SYSTEMS
N. V. Agudov, A. V. Krichigin
A nonlinear inertial monostable system is considered with the additive mixture of a harmonic signal and a white noise supplied to its input. Under certain parameters the signal power amplification and the signal-to-noise ratio at the system output can be nonmonotonic functions of the input noise intensity, i.e. there is a possibility of stochastic resonance and antiresonance to arise. An analysis of the potential profile shape impact on the signal power amplification and signal-to-noise ratio is carried out on the basis of two systems described by piecewise linear monostable potentials.
Keywords: stochastic resonance, monostable system, Brownian motion.
