Читайте также:
|
|
Предположим, что на вход системы (рис 1.17) поступает помеха f(t), а полезный сигнал x(t)=O. Требуется определить флуктуационную ошибку, вызываемую отработкой системы помехи на ее входе. Спектральная плотность помехи равна . Параметры системы, а именно передаточная функция ее также известна.
Поскольку в этом случае весь сигнал на выходе системы представляет собой сигнал ошибки, спектральная плотность ошибки
, (1.77)
Где —квадрат амплитудно-частотной характеристики (АЧХ) замкнутой системы.
Для облегчения вычисления интеграла (1.76) и приведения его к табличному спектральную плотность входного сигнала представляем, как показано было ранее, в виде
. (1.78)
Подставляя (1.78) в (1.77) и обозначая , получим подынтегральное выражение (1.76), т. е.
. (1.79)
Обычно —рациональная дробь, также может быть представлено в виде рациональной дроби.
Учитывая (1.71) и (1.76), получим выражение для среднего квадрата ошибки в виде табличного интеграла:
(1.80)
где многочлены под интегралом
причем п>т.
Таблицы интегралов до n = 6 приведены в приложении I.
Если помеха действует не на входе системы, то вместо берется передаточная функция Kyf(p), соответствующая месту приложения воздействия f.
Пример 1.4. Определить флуктуационную ошибку
системы с передаточной функцией , если на ее входе действует возмущение в виде белого шума со спектральной плотностью
Для решения задачи находим
а также представляем
Тогда средний квадрат ошибки
Имея в виду, что c(p)=akv и d(p) =kv+p, найдем коэффициенты полиномов числителя и знаменателя:
с0 = akv; do = kv; .
С помощью табличного интеграла окончательно получим:
.
Заметим, что подготовка для вычисления интеграла, сводящаяся к нахождению выражения , может быть упрощена. Так как подынтегральное выражение в формуле (1.80) является результатом возведения в квадрат , то для вычисления интеграла (1.80) коэффициенты полиномов числителя и знаменателя определяют исходя из соотношения
. (1.81)
Таким образом, решение задачи сводится к нахождению произведения .
Рассмотрим другой случай, когда на входе системы помеха отсутствует f(t)= 0, а полезный сигнал представляет собой случайный стационарный процесс со спектральной плотностью .
Методика расчета остается прежней, но заменяется на , a на Fx(p).
Ошибка системы
Zx(p) = Kzx(p)X(p). (1.82)
Спектральная плотность ошибки.
(1.83)
где —квадрат АЧХ ошибки.
Подставляя в выражение (1.83) и представляя в виде произведения сопряженных составляющих , т. е.
,
окончательно получим выражение, аналогичное (1.80):
. (1.84)
Пример 1.5. Вычислить значение среднего квадрата сигнала на выходе схемы рис. 1.18, если на ее вход подается сигнал со спектральной плотностью в2*сек. Параметры схемы: R = 500 ом; С = 0,01 ф.
Передаточная функция системы
.
В соответствии с (1.81) находим
.
Следовательно, , , и значение среднего квадрата сигнала на выходе схемы
в2.
Спектральная плотность напряжения на выходе схемы
.
Рассмотрим более общий случай, когда на систему помимо задающего воздействия x(t) действует одновременно помеха f(t).
Суммарная ошибка
. (1.85)
Спектральная плотность ошибки
(1.86)
Где
и
представляют собой взаимные спектральные плотности полезного сигнала и помехи, а и — частотные характеристики ошибки по полезному сигналу и помехе. При отсутствии корреляции между полезным сигналом и помехой
и формула (1.86) упрощается:
(1.87)
Для помехи, приложенной совместно с задающим воздействием, когда , и при отсутствии корреляции между ними получим
. (1.88)
Средний квадрат ошибки
, (1.89)
Где
(1.90)
составляющая дисперсии ошибки, вызываемая задающим воздействием x(t);
(1.91)
составляющая дисперсии ошибки, вызываемая возмущающим воздействием f(t).
Среднеквадратическое значение суммарной ошибки системы
(1.92)
Среднеквадратическую ошибку системы, определяемую формулой (1.92), не следует смешивать со среднеквадратическим отклонением , которое равно положительному квадратному корню из дисперсии .
Как следует из формулы (1.89), среднее значение квадрата ошибки зависит от структуры системы (вида ее передаточной функции и параметров) и от спектральных плотностей входного сигнала и помехи.
Для минимизации соответствующей составляющей ошибки системы необходимо уменьшать площадь под кривой произведения спектральной плотности входного сигнала на квадрат амплитудно-частотной характеристики.
Заменяя в выражении (1.90) передаточную функцию ошибки на передаточную функцию замкнутой системы Ко(р), получим средний квадрат выходной величины .
Если в задающем сигнале x(t) можно выделить воздействие в виде неслучайной составляющей mx(t), представляющей собой медленно меняющуюся функцию времени, и стационарный центрированный случайный процесс, т. е.
,
то точность системы можно оценить средним квадратом ошибки, равным сумме квадратов динамической и случайной ошибок:
или
Здесь — коэффициент, определяющий удельный вес динамической ошибки;
,
Где ,... — коэффициенты ошибки.
Для случая, когда можно предположить, что скорость изменения задающего воздействия постоянна в течение рассматриваемого интервала времени, т. е. , а помеха— белый шум, в соответствии с (1.89) получим для систем с астатизмом 1-го порядка
.
Пример 1.6. Определить средний квадрат суммарной ошибки САУ с передаточной функцией , если на входе системы действует задающее воздействие со спектральной плотностью и помеха со спектральной плотностью . Ошибка системы определяется формулой (1.89). Вторая составляющая ошибки определена в примере 1.4:
.
Вычислим первую составляющую ошибки (от задающего воздействия).
Передаточная функция ошибки
.
Представим спектральную плотность через сопряженные составляющие:
.
Находим
.
Табличный интеграл
.
Окончательно
.
Из данного выражения следует, что для уменьшения составляющей ошибки от полезного сигнала необходимо увеличение , а для уменьшения составляющей ошибки от помех нужно уменьшать.
Основным достоинством аналитического метода является возможность установления связи между величиной СКО и параметрами системы, что позволяет определять значения -параметров системы, при которых СКО оказывается минимальной.
Дата добавления: 2015-07-15; просмотров: 406 | Нарушение авторских прав
<== предыдущая страница | | | следующая страница ==> |
Аналитический метод расчета | | | Графоаналитический метод расчета |