Вычисление и тестирование невязок
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')
Загрузка данных во временной области.
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.