Идентификация обьекта автоматизации

@Langoliere · Регистрация: 23.01.2013

Студворк — интернет-сервис помощи студентам

Доброго времени суток, вопрос в следующем - имеются экспериментальные данные изменения параметра в ответ на движение регулирующего клапана. Для идентификации объекта запускаю System Identification. Импортирую данные, далее в Estimate выбираю пункт Process Models, получаю передаточную функцию с задержкой и без интегральной части с соответствием почти 95%.
Нет понимания правильно ли я все сделал.Можно ли получить передаточную ф-ю лучше подходящую под данные? Нужно ли это? Возможно во вкладке Estimate нужно было выбирать другие пункты?
Просто график полученной функции выглядит как прямая линия и не совсем соответствует изменениям параметра. Нормально ли это?
В дальнейшем нужно будет выбрать и рассчитать регулятор(ПИ или ПИД) и оценить качество регулирования. Подходит ли эта функция для дальнейшей работы?
Мануалы читал, видео с официального сайта смотрел, понимания все равно не хватает.

Matlab M
1
2
3
4
5
y=[1925,1925,1930.8,1942.8,1957.1,1983.4,2004.5,2021.4,2035.8,2047.1,2056.7,2074.2,2095.3,2110.9,2126.8,2135.9,2141.5,2155.6,2175.3,2190.1,2200]
x=[10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10,10]
t=[0:20]
plot(t,y)
grid on

В архиве файл сессии из System Identification App.

@ida2 · 26.05.2020, 13:08

Сообщение от Langoliere

данные изменения параметра в ответ на движение регулирующего клапана

Судя по тому, что в приведенном вами примере, все элементы вектора x имеют одно и то же значение (а именно - 10), можно предположить, что "движение регулирующего клапана" не происходило. Может быть причина именно в этом?

@Langoliere · 26.05.2020, 13:25 **[ТС]**

Все не совсем так. Вектор Y- изменения параметра со временем. А X это изменение задания регулятора. То есть на 10% он изменил положение с 70 до 80.
То есть клапан открылся на 10% и расход увеличился с 1925 до 2200
Написал это и понял что исходные данные не те по Y. Должно быть

Matlab M
1
y=[0,0,5.8,17.8,32.1,58.4,79.5,96.4,110.8,122.1,131.7,149.2,170.3,185.9,201.8,210.9,216.5,230.6,250.3,265.1,275]

@ida2 · 26.05.2020, 14:04

Сообщение от Langoliere

А X это изменение задания регулятора.

Так вот именно самого изменения Х система и не видит в ваших данных. Изменение задания тоже происходит во времени. Это вы знаете, что 10 (в данном случае) это изменение именно с 70 до 80. А ей бедной, что делать? А что было до этого изменения, вы считаете, что ей знать не нужно? А в случае изменения с 50 до 60 вы тоже напишете 10? Продемонстрируйте это изменение в векторе X, в абсолютных величинах. Например: x=[70, 80, 80, ...], но и У сформируйте соответственно. Или, если в приращениях, то тогда и X и Y (как в вашем предыдущем посте).

@Langoliere · 26.05.2020, 16:06 **[ТС]**

Насколько я понял, нужно указать именно воздействие на объект, разве нет?
Переделал с изменениями регулятора от 70 до 80. Кажется соответствует примерно реальному движению регулятора. Но тут могу ошибаться, может лучше сделать просто линейную зависимость?
Пока выглядит вот так.

Matlab M
1
2
3
4
5
6
7
y=[0,0,5.8,17.8,32.1,58.4,79.5,96.4,110.8,122.1,131.7,149.2,170.3,185.9,201.8,210.9,216.5,230.6,250.3,265.1,275]
x=[70,70,70.5,71,71.8,72.5,73.4,74.4,75.3,75.8,76.2,76.5,77,77.5,78,78.4,78.7,79.2,79.7,79.9,80]
t=[0:20]
plot(t,y)
grid on
plot(t,x)
grid on

В приложении новая сессия с новыми исходными значениями

@ida2 · 26.05.2020, 22:12

Ну, положим, так уже краше. Но остаются вопросы.
1) Воздействие на объект - положение задвижки (кажется), изменяется с 70 до 80. При этом, одновременно с изменением положения задвижки, изменяется расход. Но вначале, вы писали, что расход изменяется с 1925 до 2200, а тут с 0. Может стоило бы все привести к одним единицам - или все в абсолютных, или все в приращениях?
2) Судя по тому, что весь процесс обрывается при первом достижении задвижкой положения 80, вы считаете, что переходный процесс на этом заканчивается. Если это так - замечательно, но если нет, то полученная модель не в полной мере будет соответствовать реальности. Если есть возможность, добавьте еще несколько точек, соответствующих положению задвижки 80, и убедитесь, что значение расхода перестало изменяться.

@Langoliere · 27.05.2020, 09:54 **[ТС]**

Если сделать все в приращениях то соответствие передаточной функции падает и коэфициенты странно выглядят. Сессия в файле variant3

Matlab M
1
2
3
4
5
6
7
y=[0,0,5.8,17.8,32.1,58.4,79.5,96.4,110.8,122.1,131.7,149.2,170.3,185.9,201.8,210.9,216.5,230.6,250.3,265.1,275]
x=[0,0,0.5,1,1.8,2.5,3.4,4.4,5.3,5.8,6.2,6.5,7,7.5,8,8.4,8.7,9.2,9.7,9.9,10]
t=[0:20]
plot(t,y)
grid on
plot(t,x)
grid on

Если все в абсолютных значениях, то коэфициенты вроде адекватные, но задержки нет, хотя судя по данным должна быть. Соответствие тоже меньше 95. Сессия в файле variant4

Matlab M
1
2
3
4
5
6
7
y=[1925,1925,1930.8,1942.8,1957.1,1983.4,2004.5,2021.4,2035.8,2047.1,2056.7,2074.2,2095.3,2110.9,2126.8,2135.9,2141.5,2155.6,2175.3,2190.1,2200]
x=[70,70,70.5,71,71.8,72.5,73.4,74.4,75.3,75.8,76.2,76.5,77,77.5,78,78.4,78.7,79.2,79.7,79.9,80]
t=[0:20]
plot(t,y)
grid on
%plot(t,x)
%grid on

Вообще для идентификации подобных процессов Process Models используется? Или другие лучше использовать?

@Centurio · 27.05.2020, 12:15

Сообщение от Langoliere

задержки нет, хотя судя по данным должна быть

Она там и есть, небольшая.

@Langoliere · 27.05.2020, 14:58 **[ТС]**

В итоговой передаточной функции обьекта после идентификации задержки нет.
Два файла сессии забыл скинуть, можно посмотреть что получается в итоге

Новые блоги и статьи Все статьи Все блоги /
делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .
Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .

@ida2 164 / 134 / 61 Регистрация: 16.05.2015 Сообщений: 372
	26.05.2020, 22:12
	Ну, положим, так уже краше. Но остаются вопросы. 1) Воздействие на объект - положение задвижки (кажется), изменяется с 70 до 80. При этом, одновременно с изменением положения задвижки, изменяется расход. Но вначале, вы писали, что расход изменяется с 1925 до 2200, а тут с 0. Может стоило бы все привести к одним единицам - или все в абсолютных, или все в приращениях? 2) Судя по тому, что весь процесс обрывается при первом достижении задвижкой положения 80, вы считаете, что переходный процесс на этом заканчивается. Если это так - замечательно, но если нет, то полученная модель не в полной мере будет соответствовать реальности. Если есть возможность, добавьте еще несколько точек, соответствующих положению задвижки 80, и убедитесь, что значение расхода перестало изменяться. 0