Определим спектральную плотность и корреляционную функцию выходного процесса 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) постоянная спектральная плотность
заменяется на спектральную плотность
.
или с использованием интеграла Дюамеля для
будем иметь
(13.6.5)