resid
Вычисление и тестирование невязок
Синтаксис
Описание
resid( вычисляет 1-ступенчатые ошибки предсказания ( невязки) для идентифицированной модели, Data,sys)sys, и строит динамику невязки как одно из следующего, в зависимости от данных в Data:
Для данных временной области,
residстроит графики автокорреляции невязок и перекрестной корреляции невязок с входными сигналами. Корреляции сгенерированы для лагов -25 до 25. Чтобы задать другое значение максимальной задержки, используйтеresidOptions. 99% -доверия область, маркирующая статистически незначительные корреляции, отображается как заштрихованная область вокруг оси X.Для данных частотного диапазона,
residстроит график диаграммы Боде частотной характеристики от входных сигналов до невязок. 99% -доверия область, маркирующая статистически незначительный ответ, показана как область вокруг оси X.
Чтобы изменить параметры отображения, щелкните правой кнопкой мыши график, чтобы получить доступ к контекстному меню. Для получения дополнительной информации о меню смотрите Советы.
resid( вычисляет и строит графики невязок нескольких идентифицированных моделей Data,sys1,...,sysn)sys1..., sysn.
Примеры
График остаточных корреляций модели
Загрузка данных во временной области.
load iddata1
data = z1;Оцените модель ARX.
sys = arx(data,[1 1 0]);
Постройте график автокорреляции невязок и перекрестной корреляции между невязками и входами.
resid(data,sys)
Корреляции вычисляются до максимальной задержки по умолчанию, 25. 99% -доверия область, маркирующая статистически незначительные корреляции, отображается как заштрихованная область вокруг оси X.
Преобразуйте данные в частотный диапазон.
data2 = fft(data);
Вычислите невязки для идентифицированной модели, sys, и данные частотному диапазону. Постройте график остаточной реакции с помощью красных крестов.
resid(data2,sys,'rx')
Для данных частотных диапазонов, resid строит графики Диаграммы Боде, показывающие частотную характеристику от входа до остатков.
Сравнение невязок для нескольких идентифицированных моделей
Загрузка данных во временной области.
load iddata1Оцените модель ARX.
sys1 = arx(z1,[1 1 0]);
Оцените модель передаточной функции.
sys2 = tfest(z1,2);
Постройте график корреляций невязок.
resid(z1,sys1,'b',sys2,'r')
Перекрестная корреляция между невязками sys2 и входы лежат в 99% доверительной полосе для всех лагов.
Задайте максимальную задержку для вычислений остаточной импульсной характеристики
Загрузка данных во временной области.
load iddata1Оцените модель ARX.
sys = arx(z1,[1 1 0]);
Задайте максимальную задержку для вычисления остаточной корреляции.
opt = residOptions('MaxLag',35);Постройте график импульсной характеристики от входа до невязок.
resid(z1,sys,opt,'ir')
Вычисление невязок для системы MIMO
Загрузка данных во временной области.
load iddata7Данные являются набором данных с двумя входами и одним выходом.
Оцените модель ARX.
sys = tfest(z7,2);
Вычислите невязки и их автокорреляции и перекрестные корреляции с входами.
[E,R] = resid(z7,sys);
R матрица корреляций 26 на 3 на 3. Для примера,
R(:,1,1)автокорреляция невязок до запаздывания 25.R(:,1,2)- перекрестная корреляция невязок с первым входом до запаздывания 25.
E является объектом iddata с невязками в качестве выходных данных и входами данных валидации (z7) в качестве входных данных. Можно использовать E для идентификации моделей ошибок и анализа динамики ошибок.
Постройте график ошибки.
plot(E)
Оцените импульсную характеристику между входами и невязками. Постройте график с 3 стандартной доверительной областью отклонения.
I = impulseest(E); showConfidence(impulseplot(I,20),3)
Входные параметры
Выходные аргументы
Совет
При щелчке правой кнопкой мыши по графику открывается контекстное меню, в котором можно получить доступ к следующим опциям:
Systems - выберите системы, чтобы просмотреть остаточную корреляцию или графики отклика. По умолчанию все системы построены.
Show Confidence Region - Просмотр 99% доверительной области, маркирующей статистически незначительные корреляции. Применимо только для корреляционных графиков.
Data Experiment - только для данных нескольких экспериментов. Переключение между данными из различных экспериментов.
Characteristics - Просмотр характеристик данных. Не применяется для корреляционных графиков.
Peak Response - Просмотр максимальной чувствительности данных.
Confidence Region - Просмотр 99% доверительной области, маркирующей статистически незначительный ответ.
Show - Применяется только для графиков частотной характеристики.
Magnitude - Просмотрите величину частотной характеристики.
Phase - Просмотр фазы частотной характеристики.
I/O Grouping - для наборов данных, содержащих более одного входного или выходного канала. Выберите группировку входа и выходных каналов на графике. Не применяется для корреляционных графиков.
None - Стройте графики каналов ввода-вывода в собственных отдельных осях.
All - Группируйте все входные каналы вместе и все выходные каналы вместе.
I/O Selector - для наборов данных, содержащих более одного входного или выходного канала. Выберите подмножество входного и выходного каналов для построения графика. По умолчанию все выходные каналы построены.
Grid - добавить сетки на график.
Normalize - Нормализуйте шкалу Y всех данных на графике. Не применяется для данных частотной характеристики.
Full View - Вернуться к полному представлению. По умолчанию график масштабируется до полного вида.
Initial Condition - Определение обработки начальных условий.
Задайте как одно из следующих:
Estimate - Обрабатывайте начальные условия как параметры оценки.
Zero - Установите все начальные условия на нуль.
Absorb delays and estimate - Поглощает ненулевые задержки в коэффициенты модели и обрабатывает начальные условия как параметры оценки. Используйте эту опцию только для моделей в дискретном времени.
Properties - Откройте диалоговое окно «Property Editor», чтобы настроить атрибуты графика.
Ссылки
[1] Ljung, L. System Identification: Теория для пользователя. Upper Saddle River, NJ: Prentice-Hall PTR, 1999, Section 16.6.