Купить Matlab  |  Mathematica  |  Mathcad  |  Maple  |  Statistica  |  Другие пакеты Поиск по сайту
Internet-класс  |  Примеры  |  Методики  |  Форум  |  Download
https://hub.exponenta.ru/


10 ноутбуков
В.П.Дьяконова
 
Вычисление амплитудно-частотной и фазо-частотной характеристик линейной системы по ее переходной характеристике

Расчет в аналитическом виде

Переходную характеристику линейной цепи легко снять с помощью генератора прямоугольных импульсов и осциллографа. Однако снятие амлитудно-частотной и фазо-частотных характеристик линейной цепи намного сложнее, ибо приходится измерять коэффициенты передачи и фазовые сдвиги многократно для ряда частот. Поэтому представляет интерес получение амплитудно-частотной A(w) и фазо-частотной j(w) характеристик по известной переходной характеристике a(t).

Вначале рассмотрим аналитический метод решения данной задачи. Зададим переходную характеристику линейной цепи a(t) и вычислим ее производную:

[Graphics:Images/index_gr_1.gif]
[Graphics:Images/index_gr_2.gif]
[Graphics:Images/index_gr_3.gif]
[Graphics:Images/index_gr_4.gif]

Построим график переходной характеристики:

[Graphics:Images/index_gr_5.gif]

[Graphics:Images/index_gr_6.gif]

[Graphics:Images/index_gr_7.gif]

Чтобы избежать больших погрешностей в вычислениях важно, что бы в интервале времен от 0 до t0 переходная характеристика достигала предельного значения (в нашем случае 1). Найдем косинусные и синусные коэффициенты:

[Graphics:Images/index_gr_9.gif]
[Graphics:Images/index_gr_10.gif]

Теперь можно вычислить амплитудно-частотную и фазо-частотную характеристики данной системы:

[Graphics:Images/index_gr_11.gif]
[Graphics:Images/index_gr_12.gif]
[Graphics:Images/index_gr_13.gif]
[Graphics:Images/index_gr_14.gif]
[Graphics:Images/index_gr_15.gif]
[Graphics:Images/index_gr_16.gif]
[Graphics:Images/index_gr_17.gif]

[Graphics:Images/index_gr_18.gif]

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

Теперь можно построить графики амплитудно-частотной и фазо-частотной характеристик:

[Graphics:Images/index_gr_19.gif]
[Graphics:Images/index_gr_20.gif]

[Graphics:Images/index_gr_21.gif]

[Graphics:Images/index_gr_22.gif]

Приближенное вычисление

В ряде случаев переходная характеристика задается рядом дискретных отсчетов - например экспериментальных. Сформируем такую таблицу, пока на основе аналитической зависимости:

[Graphics:Images/index_gr_23.gif]
[Graphics:Images/index_gr_24.gif]
[Graphics:Images/index_gr_25.gif]
[Graphics:Images/index_gr_26.gif]
[Graphics:Images/index_gr_27.gif]
[Graphics:Images/index_gr_28.gif]

[Graphics:Images/index_gr_29.gif]

[Graphics:Images/index_gr_30.gif]

Для численного вычисления амплитудно-частотной и фазо-частотной характеристик теперь можно использовать следующие формулы:

[Graphics:Images/index_gr_31.gif]
[Graphics:Images/index_gr_32.gif]

Ниже даны примеры вычислений и графики этих зависимостей

[Graphics:Images/index_gr_33.gif]
[Graphics:Images/index_gr_34.gif]
[Graphics:Images/index_gr_35.gif]
[Graphics:Images/index_gr_36.gif]
[Graphics:Images/index_gr_37.gif]
[Graphics:Images/index_gr_38.gif]
[Graphics:Images/index_gr_39.gif]

[Graphics:Images/index_gr_40.gif]

[Graphics:Images/index_gr_41.gif]

Поскольку функция Arg определена в интервале от -p до p, то при выходе из этого интервала фазовый угол на графике дает скачок. На самом деле это означает, что часть графика фазо-частотной характеристики после скачка надо просто перенести вниз.

Внимание: возведите выражение для переходной характеристики при аналитических вычислениях в куб и исполните данный ноутбук. В этом случае один и тот же пример будет вычислен двумя методами - аналитическим и приближенным. Сравните полученные характеристики и убедитесь в их идентичности.

В начало страницы

| На первую страницу | Поиск | Купить Matlab

Исправляем ошибки: Нашли опечатку? Выделите ее мышкой и нажмите Ctrl+Enter


Copyright © 1993-2024. Компания Softline. Все права защищены.

Дата последнего обновления информации на сайте: 04.03.17
Сайт начал работу 01.09.00