Связь между спектральными плотностями и корреляционными функциями входного и выходного процессов для одностепенного осциллятора




Определим спектральную плотность и корреляционную функцию выходного процесса y (t) для линейного осциллятора через спектральную плотность и корреляционную функцию входного центрированного случайного процесса х (t). Перепишем для удобства уравнение (13.2.1) для одностепенного осциллятора

(13.6.1)

Корреляционная функция выходного процесса y (t) по определению имеет вид

или с использованием интеграла Дюамеля для будем иметь (13.6.2)

Запишем выражение для спектральной плотности Sy (w) выходного процесса y, как преобразование Фурье корреляционной функции Ky (t) (13.6.2)

(13.6.3)

Введем два сомножителя в (13.6.3), произведение которых равно 1,

(13.6.4)

 

(13.6.5)

Проведем замену переменных t=t – x+h Þ t =t+ x – h и, соответственно, dt=dt.

Тогда (13.6.5) преобразуется к следующему виду

(13.6.6)

В (13.6.6) входят передаточная функция H (w) и сопряженная с ней функция H* (w).

(13.6.7)

Окончательно связь между спектральными плотностями выходного Sy (w) и входного Sх (w) процессов для одностепенного осциллятора запишется так

(13.6.8)

Отметим, что корреляционную функцию Ky (t) можно получить путем применения обратного преобразования Фурье к спектральной плотности Sy (w)

(13.6.9)

Известно, что связана с дисперсией и среднеквадратичным значением простым соотношением

(13.6.10)

Тогда с учетом (13.6.9) и (13.6.10) среднеквадратичное значение выходной величины будет связано со спектральной плотностью входной величины следующим равенством

(13.6.11)

Типовой задачей для осциллятора является следующая задача.

Задан осциллятор с массой m, жесткостью k и логарифмическим декрементом d. На массу осциллятора действует внешняя случайная сила со спектральной плотностью . Определить среднеквадратичное усилие в пружине осциллятора (собственная частота осциллятора )

Правая часть уравнений (13.6.1) будет выглядеть так

(13.6.12)

усилие в пружине осциллятора

(13.6.13)

Если в уравнение (13.6.1) подставить значение из (13.6.13), а в правую часть из (13.6.12), тогда уравнение для усилия в пружине N приобретет вид

(13.6.14)

Структура уравнения для усилия в пружине (13.6.14) очевидно аналогична структуре уравнений (13.6.1) для перемещения . При соотношении между правыми частями уравнений (13.6.1) и (13.6.14)

(13.6.15)

спектральная плотность будет связана со спектральной плотностью равенством

(13.6.16)

Передаточная функция будет иметь вид (13.3.5).

Поэтому далее для определения среднеквадратичного значения усилия в пружине будем пользоваться формулами (13.6.11) и (13.6.16).

(13.6.16а)

Остановимся на оценке дисперсии (13.6.10) в случае постоянной спектральной плотности входного воздействия с использованием теории вычетов.

С учетом представления для передаточных функций в виде простых дробей (13.4.18) выражение для дисперсии выходной величины (13.6.10) запишется так

(13.6.17)

Комплексно сопряженная величина передаточной функции в (13.6.17) используется без разложения на простые дроби, так как у нее нет полюсов в верхней полуплоскости.

Интеграл в (13.6.17) будем также вычислять путем интегрирования по контуру С, расположенному в верхней полуплоскости (см. рисунок 1), с использованием теории вычетов.

При ограниченной спектральной плотности Sx (w) интеграл по верхней полуокружности равен нулю, так как модуль подынтегральной функции z (w) на полуокружности большого радиуса R мажорируется выражением

(13.6.18)

В области, ограниченной контуром С, имеется два простых полюса w 1 и w 2. Поэтому с использованием теории вычетов для интеграла (13.6.17) получим следующее значение

(13.6.19)

На основе формулы для вычетов в простых полюсах (13.4.23) в случае белого шума при постоянной плотности будем иметь

(13.6.20)

После выполнения алгебраических выкладок в скобках формулы (13.6.20) с учетом значений для и из (13.4.19) получим следующее выражение для дисперсии при постоянном значении спектральной плотности входного воздействия

(13.6.21)

С учетом связи между обобщенной спектральной плотностью и спектральной плотностью внешней случайной силы (13.6.16) получаем окончательное выражение для дисперсии усилия в пружине осциллятора в зависимости от постоянной спектральной плотности внешней силы

(13.6.22)

Отметим, что в случае плавно меняющейся спектральной плотности в окрестности резонансной частоты w 0 формула (13.6.22) позволяет получить приближенную оценку дисперсии. При этом в формуле (13.6.22) постоянная спектральная плотность заменяется на спектральную плотность .

 



Поделиться:




Поиск по сайту

©2015-2024 poisk-ru.ru
Все права принадлежать их авторам. Данный сайт не претендует на авторства, а предоставляет бесплатное использование.
Дата создания страницы: 2019-06-03 Нарушение авторских прав и Нарушение персональных данных


Поиск по сайту: