freqresp

Частотная характеристика по сетке

Синтаксис

[H,wout] = freqresp(sys) H = freqresp(sys,w) H = freqresp(sys,w,units) [H,wout,covH] = freqresp(idsys,...)

Описание

[H,wout] = freqresp(sys) возвращает частотную характеристику модели sys динамической системы на частотах wout. freqresp команда автоматически определяет частоты на основе динамики sys.

H = freqresp(sys,w) возвращает частотную характеристику на действительной сетке частоты, заданной векторным w.

H = freqresp(sys,w,units) явным образом задает единицы частоты w с units.

[H,wout,covH] = freqresp(idsys,...) также возвращает ковариацию covH из частотной характеристики идентифицированной модели idsys.

Входные параметры

`sys`	Любая модель динамической системы или массив моделей.
`w`	Вектор из действительных частот, на которых можно оценить частотную характеристику. Задайте частоты в модулях `rad/TimeUnit`, где `TimeUnit` единицы измерения времени, заданные в `TimeUnit` свойство `sys`.
`units`	Единицы частот во входном векторе частоты `w`В виде одного из следующих значений: `'rad/TimeUnit'` — радианы на единицу измерения времени заданы в `TimeUnit` свойство `sys` `'cycles/TimeUnit'` — циклы на единицу измерения времени заданы в `TimeUnit` свойство `sys` `'rad/s'` `'Hz'` `'kHz'` `'MHz'` `'GHz'` `'rpm'` Значение по умолчанию: `'rad/TimeUnit'`
`idsys`	Любая идентифицированная модель.

Выходные аргументы

H

Массив, содержащий значения частотной характеристики.

Если sys отдельная модель динамической системы, имеющая Ny выходные параметры и Nu входные параметры, H 3D массив с размерностями Ny- Nu- Nw, где Nw количество точек частоты. Таким образом, H(:,:,k) ответ на частоте w(k) или wout(k).

Если sys массив моделей размера [Ny Nu S1 ... Sn]H массив с размерностями Ny- Nu- Nw- S1-...-by- SN Массив.

Если sys модель данных частотной характеристики (такой как frd, genfrd, или idfrd), freqresp(sys,w) оценивает к NaN для значений w выходить за пределы интервала частоты задано sys.frequency. freqresp команда может интерполировать между частотами в sys.frequency. Однако freqresp не может экстраполировать вне интервала частоты, заданного sys.frequency.

wout

Вектор из частот, соответствующих значениям частотной характеристики в H. Если вы не используете w от входных параметров до freqresp, команда автоматически определяет частоты wout на основе системной динамики. Если вы задаете w, затем wout = w

covH

Ковариация ответа H. Ковариация 5D массив где covH(i,j,k,:,:) содержит ковариационную матрицу 2 на 2 ответа от iвход th к jth выход на частоте w(k). (1,1) элементом этой матрицы 2 на 2 является отклонение действительной части ответа. (2,2) элементом является отклонение мнимой части. (1,2) и (2,1) элементами является ковариация между действительными и мнимыми частями ответа.

Примеры

свернуть все

Вычислите частотную характеристику системы

Этот пример использует:

Control System Toolbox

Скрипт Open Live Script

Создайте следующую систему с 2 выходами, с 2 входами:

$s y s = [\begin{array}{cccccccccccccccccccc} 0 & \frac{1}{s + 1} \\ \frac{s - 1}{s + 2} & 1 \end{array}]$

sys11 = 0; 
sys22 = 1; 
sys12 = tf(1,[1 1]); 
sys21 = tf([1 -1],[1 2]);
sys = [sys11,sys12;sys21,sys22];

Вычислите частотную характеристику системы.

[H,wout] = freqresp(sys);

H 2 2 45 массивами. Каждая запись H(:,:,k) в H матрица 2 на 2, дающая комплексную частотную характеристику всех пар ввода - вывода sys на соответствующей частоте wout(k). Эти 45 частот в wout автоматически выбраны на основе динамики sys.

Вычислите частотную характеристику на заданной сетке частоты

Этот пример использует:

Control System Toolbox

Скрипт Open Live Script

Создайте следующую систему с 2 выходами, с 2 входами:

$s y s = [\begin{array}{cccccccccccccccccccc} 0 & \frac{1}{s + 1} \\ \frac{s - 1}{s + 2} & 1 \end{array}]$

sys11 = 0; 
sys22 = 1; 
sys12 = tf(1,[1 1]); 
sys21 = tf([1 -1],[1 2]);
sys = [sys11,sys12;sys21,sys22];

Создайте логарифмически распределенную сетку 200 точек частоты между 10 и 100 радианами в секунду.

w = logspace(1,2,200);

Вычислите частотную характеристику системы на заданной сетке частоты.

H = freqresp(sys,w);

H 2 2 200 массивами. Каждая запись H(:,:,k) в H матрица 2 на 2, дающая комплексную частотную характеристику всех пар ввода - вывода sys на соответствующей частоте w(k).

Вычислите частотную характеристику и связанную ковариацию

Этот пример использует:

Скрипт Open Live Script

Вычислите частотную характеристику и сопоставленную ковариацию для идентифицированной модели процесса на ее частоте максимальной чувствительности.

Загрузите данные об оценке z1.

load iddata1 z1

Оцените модель процесса SISO с помощью данных.

model = procest(z1,'P2UZ');

Вычислите частоту, на которой модель достигает пикового усиления частотной характеристики. Чтобы получить более точный результат, задайте значение допуска 1e-6.

[gpeak,fpeak] = getPeakGain(model,1e-6);

Вычислите частотную характеристику и сопоставленную ковариацию для model на его частоте максимальной чувствительности.

[H,wout,covH] = freqresp(model,fpeak);

H значение отклика в fpeak частота и wout совпадает с fpeak.

covH 5-мерный массив, который содержит ковариационную матрицу ответа от входа до выхода на частоте fpeak. Здесь covH(1,1,1,1,1) отклонение действительной части ответа и covH(1,1,1,2,2) отклонение мнимой части. covH(1,1,1,1,2) и covH(1,1,1,2,1) элементами является ковариация между действительными и мнимыми частями ответа.

Больше о

свернуть все

Частотная характеристика

В непрерывное время frequency response на частоте ω является значением передаточной функции в s = jω. Для моделей в пространстве состояний этим значением дают

$H (j ω) = D + C {(j ω I - A)}^{- 1} B$

В дискретное время частотная характеристика является передаточной функцией, оцененной в точках на модульном круге, которые соответствуют действительным частотам. freqresp сопоставляет действительные частоты w(1)..., w(N) к точкам на модульном круге использование преобразования $z = e^{j ω T_{s}}$ . _Ts является шагом расчета. Функция возвращает значения передаточной функции в получившихся значениях z. Для моделей с незаданным шагом расчета, freqresp использование _Ts = 1.

Алгоритмы

Для передаточных функций или моделей нулей, полюсов и усиления, freqresp оценивает числитель (числители) и знаменатель (знаменатели) в заданных точках частоты. Для моделей в пространстве состояний непрерывного времени (A, B, C, D), частотная характеристика

$\begin{matrix} D + C {(j ω - A)}^{- 1} B, & ω = \end{matrix} ω_{1}, \dots, ω_{N}$

Для КПД A уменьшается до верхней формы Хессенберга и линейного уравнения (jω − A) X =, B решен в каждой точке частоты, использовав в своих интересах структуру Hessenberg. Сокращение к форме Хессенберга обеспечивает хороший компромисс между КПД и надежностью. См. [1] для получения дополнительной информации об этом методе.

Альтернативы

Использование evalfr оценивать частотную характеристику на отдельных частотах или небольших числах частот. freqresp оптимизирован для средних и крупных векторов из частот.

Ссылки

[1] Laub, A.J., "Эффективные Многомерные Расчеты Частотной характеристики", IEEE^® Transactions на Автоматическом управлении, AC-26 (1981), стр 407-408.

Смотрите также

Представлено до R2006a

Документация

freqresp

Синтаксис

Описание

Входные параметры

Выходные аргументы

Примеры

Вычислите частотную характеристику системы

Вычислите частотную характеристику на заданной сетке частоты

Вычислите частотную характеристику и связанную ковариацию

Больше о

Частотная характеристика

Алгоритмы

Альтернативы

Ссылки

Смотрите также

Документация System Identification Toolbox

Поддержка