Exponenta Banner

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) элементами является ковариация между действительными и мнимыми частями ответа.

Примеры

свернуть все

Скрипт Open Live Script

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

sys=[01s+1s-1s+21]

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.

Скрипт Open Live Script

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

sys=[01s+1s-1s+21]

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)−1B

В дискретное время частотная характеристика является передаточной функцией, оцененной в точках на модульном круге, которые соответствуют действительным частотам. freqresp сопоставляет действительные частоты w(1)..., w(N) к точкам на модульном круге использование преобразования z=ejωTs. Ts является шагом расчета. Функция возвращает значения передаточной функции в получившихся значениях z. Для моделей с незаданным шагом расчета, freqresp использование Ts = 1.

Алгоритмы

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

D+C(jω−A)−1B,ω=ω1,…,ωN

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

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

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

Ссылки

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

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

| | | | | | (System Identification Toolbox)

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

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

Поддержка