Настройка ПИД регулятора

@Sosnovskyi · Регистрация: 28.06.2013

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

Есть схема в Simulink, с ПИД регулятором. Схема - работает. Выдает кривую.
Однако, этого мало. Нужно посчитать, так называемую, погрешность работы регулятора. Это заключается в том, что нужно высчитать площадь пространства, которое я отметил красным маркером.
Посоветовали делать через отсылку в workspace результатов симуляции, а далее через нахождения площади, через интеграл.
Не могу понять, как это осуществить.
Отсылка в workspace есть, остальное отсутствует.
В приложение архив с симулинк файлом, оригинальное изображение кривой, и что нужно посчитать.

@N_J · 14.11.2013, 15:41

Предлагаю уменьшить шаг дискретизации для начала.
Можно дальше поступить так:
1) подключаешь блок ToWorkspace (допустим, назовешь переменную res), причем в настройках указываешь, что формат не Array и не Structure, а Structure With Time.
2) определяешь max от res.signals.values (пусть будет resmax) и соответствующий этому значению индекс ind
3) далее используешь функцию trapz(x,y), где x = res.time(1:ind), y = res.signals.values(1:ind)

@R2D2 · 14.11.2013, 15:59

Sosnovskyi, блок интегратора поставить что мешает?

@Sosnovskyi · 14.11.2013, 16:26 **[ТС]**

N_J, а можно, примерный код, пожалуйста?..)
R2D2, он же посчитает не конкретну ту, которая нужна, а выдаст значения по всему графику?

@R2D2 · 14.11.2013, 16:57

Ну по идее можно заморочится и сделать так, что бы интегратор интегрировал только до первого максимума. Но если эта площадь нужна как промежуточный результат, то и определение "на глаз" подойдет.

@Hofstadter · 27.11.2013, 08:07

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

N_J, а можно, примерный код, пожалуйста?..)

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

trapz(x,y), где x = res.time(1:ind), y = res.signals.values(1:ind)

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

причем в настройках указываешь, что формат не Array и не Structure, а Structure With Time.

это важно!

мое предложение следующее:
1. если нужно вычислять заштрихованную площадь, то нужно в Matlab загрузить заштрихованную площадь.
2. если нужно вычислить площадь полученной функции в пределах t[0,a], nj vj;yj воспользоваться To Workspace или, что я лично нахожу удобнее воспользоваться интегратором с external reset.
на самом же деле ошибку регулятора никогда не вычисляют по отклику (исходя из моих весьма неглубоких познаний в системах управления). ошибка - это разность. а ошибка регулирования - это разность задания и текущего значения.
3. возможно, ты регулируец не то, что нужно, раз уж Вас интересует площадь?

@Hofstadter · 28.11.2013, 04:11

nj vj;yj = то можно
зачепятался

@N_J 262 / 61 / 5 Регистрация: 29.06.2012 Сообщений: 109
	14.11.2013, 15:41	2
	Предлагаю уменьшить шаг дискретизации для начала. Можно дальше поступить так: 1) подключаешь блок ToWorkspace (допустим, назовешь переменную res), причем в настройках указываешь, что формат не Array и не Structure, а Structure With Time. 2) определяешь max от res.signals.values (пусть будет resmax) и соответствующий этому значению индекс ind 3) далее используешь функцию trapz(x,y), где x = res.time(1:ind), y = res.signals.values(1:ind) 0

@R2D2 939 / 839 / 121 Регистрация: 23.11.2012 Сообщений: 2,489
	14.11.2013, 15:59	3
	Sosnovskyi, блок интегратора поставить что мешает? 0

@Sosnovskyi 0 / 0 / 0 Регистрация: 28.06.2013 Сообщений: 20
	14.11.2013, 16:26 [ТС]	4
	N_J, а можно, примерный код, пожалуйста?..) R2D2, он же посчитает не конкретну ту, которая нужна, а выдаст значения по всему графику? 0

@R2D2 939 / 839 / 121 Регистрация: 23.11.2012 Сообщений: 2,489
	14.11.2013, 16:57	5
	Ну по идее можно заморочится и сделать так, что бы интегратор интегрировал только до первого максимума. Но если эта площадь нужна как промежуточный результат, то и определение "на глаз" подойдет. 0

@Hofstadter 42 / 42 / 0 Регистрация: 27.10.2012 Сообщений: 228
	28.11.2013, 04:11	7
	nj vj;yj = то можно зачепятался 0

Опции темы