Разработка системы управления предприятием

 

Министерство образования и науки РФ

ФГБО ВПО Сибирский государственный индустриальный университет

кафедра систем информатики и управления







Курсовая работа







Выполнил ст. гр. ИП-09 Пятаев А.М.

Проверил: к.т.н, доцент Огнев С.П.












Новокузнецк 2012

Содержание


Задание

Постановка задачи анализа и синтеза СУ

. Анализ объекта управления

.1. Построение математической модели ОУ

.2. Статические и динамические хар-ки ОУ

.3. Анализ частотных хар-к

. Синтез системы управления

.1. Определение критерия управления, ограничений, внешних факторов, выбор типового закона управления, расчет настроечных параметров

.2. Построение структурной схемы системы управления

. Анализ системы управления

.1. Анализ устойчивости

.2. Анализ качества управления в типовой системе

.3. Анализ чувствительности СУ к изменениям параметров модели объекта исследования

. Программный модуль

4.1. Листинг

Задание


Доход от импорта продукции для предприятия колеблется от 3,5 (летний и зимний период) до 7,8 (весенний и осенний период) млн. руб. в месяц.

Увеличение инвестиций на 25% в производство увеличивает доход на 38% вне зависимости от сезона. Какие инвестиции необходимы для стабилизации среднегодового дохода предприятия на уровне 70 млн. руб.?


Постановка задачи анализа и синтеза СУ


Дано:

1)Объект управления - импортирующее предприятие, которое характеризуется каналом управления (по каналу управления инвестиции - доход).

2)Основное входное воздействие на объект:

: объем инвестиций

, %/мес.

) Основное выходное воздействие объекта:

: доход

, млн.руб./мес.

4)Режимы работы объекта:

а) Номинальный режим характеризуется:


= 100%/мес.


=5,65 млн.р./мес.


б) Исследовательский режим характеризуется:


==125% мес.

==7,8 млн.р./мес.


в) Заданный режим характеризуется:


y* = 5,83 млн.р./мес.


5)Производственные характеристики:

Определяют время перехода объекта в новый режим и оцениваются:

временем распределения инвестиций:


=0,5 мес.


временем изменения дохода от импорта (сезонное запаздывание):


=6 мес.


6)Ограничения:

а) По входному воздействию:


==25%/мес.

=180%/мес.


б) По выходному воздействию:


=2,09 млн.р./мес.

=8,475 млн.р./мес.


в) Ограничения по регламенту работы объекта:

включают технологию производства, нормы качества продукции и сырья.

Примем условия линейности объекта в рассматриваемых диапазонах инвестиций и производительности.

) Внешние факторы:

w, характеризуют эквивалентные изменения номинальных инвестиций, например, за счет аварий, приобретения нового оборудования, брака и т.д.

Примем:


==30%/мес.


8)Линейные методы и алгоритмы анализа и синтеза системы управления.

9)Критерии исследования и управления:

показатели качества управления, характеризующие отклонения от заданной производительности: статическая ошибка, время управления, максимальное динамическое отклонение, перерегулирование, степень затухания колебаний, модульная интегральная ошибка .

Требуется разработать систему, отвечающую всем требованиям и критериям.

1. Анализ объекта управления


Цель: исследование основных динамических характеристик предприятия по заданному каналу управления, результаты которого достаточны для синтеза управляющей системы (СУ).


1.1 Построение математической модели ОУ


На основе режимов работы функционирования объекта с соответствующими параметрами динамики построим экспериментальную характеристику временного регламента.



Для определения средних заданных инвестиций, необходимых для определения заданного дохода y* решим пропорцию:


=103%


По виду экспериментальной характеристики в исследовательском режиме определено:

объект характеризуется начальным запаздыванием вследствие распределения инвестиций =0,5 мес.

доход от импорта увеличивается только в следующем сезоне (через период в 6 мес.), т.е. через =6 мес.

состояние объекта характеризуется установившимся периодическим режимом.

Исходя из таких предпосылок, объект можно описать следующей структурой:



Данную структуру можно описать через типовые звенья:

распределение инвестиций характеризует запаздывание по денежным средствам и описывается звеном запаздывания или транспортной моделью:



где время запаздывания


мес. Отсюда:


рост дохода до нового режима соответствует типовому колебательному звену или транспортной модели:



при (незатухающие колебания), где характеризует соотношение изменения дохода к соответствующему изменению объема инвестиций:



постоянная времени характеризует время роста дохода в исслед режиме:


мес. Отсюда:



Структура объекта последовательная:



Модель в дифференциальной форме:

Из уравнения Лапласа получим:



Модель в дискретной форме:


где - период контроля продаж, который определяется из условия:

Отсюда:

0,1 мес. 3,1 дня


Далее:



Окончательно получаем:



1.2 Статические и динамические характеристики ОУ


Статическая характеристика ОУ

На основании экспериментальной характеристики объект характеризуется следующими выражениями:


РежимыНоминальныйИсследовательскийЗаданный100%/мес.125%/мес.103%/мес.y5,65 млн.р./мес.7,8 млн.р./мес.5,83 млн.р./мес.


В диапазоне входных воздействий характеристика является квазилинейной, что обуславливает применение линейных методов управления.

Переходная характеристика ОУ

Рассчитывается исходя из постоянных инвестиций в предприятия при следующих условиях:

Начальные условия принимаем равными нулевым:


0%/мес.

= 0 млн.р./мес.


За входное воздействие принимаем увеличение инвестиций на 25%:


25%/мес.


Для расчета используем дискретную модель объекта:



Расчет ведем в приращениях, т.е. результирующее выходное значение ОУ получается за счет суммирования и расчетного .


Таблица расчета переходной характеристики

lty(l)-0,99y(l-2)1,98y(l-1)0,00086V(l-5)V(l)y(l) рез00000005,6510,10000255,6520,20000255,6530,30000255,6540,40000255,6550,50000255,6560,60,021288000,0212875255,67128870,70,06343700,042149250,0212875255,71343780,80,125818-0,021074630,1256047650,0212875255,77581890,90,207604-0,062802380,2491189270,0212875255,8576041010,307784-0,124559460,4110560090,0212875255,957784111,10,425172-0,2055280,6094124090,0212875256,075172121,20,558422-0,30470620,8418403710,0212875256,208422131,30,706042-0,420920191,10567490,0212875256,356042141,40,866414-0,552837451,3979635850,0212875256,516414151,51,037805-0,698981791,7154989970,0212875256,687805161,61,218391-0,85774952,0548533150,0212875256,868391171,71,406276-1,027426662,4124148060,0212875257,056276181,81,599506-1,20620742,7844257850,0212875257,249506191,91,796096-1,392212893,1670216460,0212875257,4460962021,994047-1,583510823,5562705820,0212875257,644047212,12,191366-1,778135293,9482135730,0212875257,841366222,22,386085-1,974106794,3389042480,0212875258,036085232,32,576284-2,169452124,7244482230,0212875258,226284242,42,760105-2,362224115,1010415270,0212875258,410105252,52,935774-2,550520765,4650077330,0212875258,585774262,63,101617-2,732503875,8128334490,0212875258,751617272,73,256073-2,906416726,1412018240,0212875258,906073282,83,39771-3,070600916,4470237460,0212875259,04771292,93,525242-3,223511876,7274664620,0212875259,175242280282,527668-2,427604314,9339852560,0212875258,17766828128,12,559078-2,466992635,004783530,0212875258,20907828228,22,585871-2,502391765,0669752360,0212875258,23587128328,32,607824-2,533487625,1200245220,0212875258,25782428428,42,624768-2,560012265,1634923210,0212875258,27476828528,52,636581-2,581746165,1970397680,0212875258,28658128628,62,643198-2,598519885,2204305930,0212875258,29319828728,72,644605-2,61021535,2335324540,0212875258,29460528828,82,640838-2,616766235,2363172220,0212875258,29083828928,92,631989-2,618158615,228860220,0212875258,281989290292,618196-2,614430115,2113384350,0212875258,26819629129,12,599646-2,605669225,1840277340,0212875258,24964629229,22,576573-2,592013875,1472991130,0212875258,22657329329,32,549252-2,573649565,1016140370,0212875258,19925229429,42,517999-2,550807025,0475189220,0212875258,16799929529,52,483167-2,523759464,9856388180,0212875258,13316729629,62,445138-2,492819414,9166703770,0212875258,09513829729,72,404326-2,458335194,8413741670,0212875258,054326



Исходя из переходной характеристики, определили:

) Характеристика соответствует экспериментальной с некоторой погрешностью, но в целом модель адекватная.

) Объект характеризуется транспортным запаздыванием 0,5 мес., в течение которого доход от продаж не изменяется.

) Рост дохода происходит за общее время 6,5 мес.

) Объект является колебательным с периодом Т = 6 мес.

Импульсная характеристика ОУ

Начальные условия принимаем равными нулевым:


0%/мес.

= 0 млн.р./мес.


За входное воздействие принимаем единовременное в течение месяца увеличение инвестиций на 25%:


25%/мес.


В последующие месяцы инвестиции остаются нулевыми.


Для расчета используем дискретную модель объекта:



Расчет ведем в приращениях, т.е. результирующее выходное значение ОУ получается за счет суммирования и расчетного .


Таблица расчета импульсной характеристики

lty(l)-0,99y(l-2)1,98y(l-1)0,00086V(l-5)V(l)y(l) рез00000005,6510,10000255,6520,20000255,6530,30000255,6540,40000255,6550,50000255,6560,60,021288000,0212875255,67128870,70,06343700,042149250,0212875255,71343780,80,125818-0,021074630,1256047650,0212875255,77581890,90,207604-0,062802380,2491189270,0212875255,8576041010,307784-0,124559460,4110560090,0212875255,957784111,10,425172-0,2055280,6094124090,021287506,075172121,20,558422-0,30470620,8418403710,021287506,208422131,30,706042-0,420920191,10567490,021287506,356042141,40,866414-0,552837451,3979635850,021287506,516414151,51,037805-0,698981791,7154989970,021287506,687805161,61,197104-0,85774952,054853315006,847104171,71,342839-1,027426662,370265556006,992839181,81,473688-1,185132782,65882102007,123688191,91,588492-1,329410512,917902719007,2384922021,686263-1,458951363,145214573007,336263212,11,766194-1,572607293,338801164007,416194222,21,827663-1,669400583,497063876007,477663232,31,870241-1,748531943,618773323007,520241242,41,893691-1,809386663,703077942007,543691252,51,89797-1,851538973,749508736007,54797262,61,883226-1,874754373,757980134007,53322627927,90,508025-0,527742021,035766932006,158025280280,488006-0,517883471,005889331006,13800628128,10,463307-0,502944670,966251612006,11330728228,20,434222-0,483125810,917347754006,08422228328,30,401086-0,458673880,859759458006,05108628428,40,36427-0,429879730,79414945006,0142728528,50,324179-0,397074720,721254047005,97417928628,60,281248-0,360627020,641875058005,93124828728,70,235934-0,320937530,556871108005,88593428828,80,188713-0,278435550,467148487005,83871328928,90,140077-0,233574240,373651607005,790077290290,090527-0,18682580,27735318005,74052729129,10,040568-0,138676590,179244205005,69056829229,2-0,0093-0,08962210,080323878005,64070229329,3-0,05857-0,04016194-0,018410485005,59142829429,4-0,106770,00920524-0,115973399005,54323229529,5-0,153410,0579867-0,21140095005,49658629629,6-0,198060,10570047-0,303760216005,4519429729,7-0,240280,15188011-0,392158287005,409722



Исходя из импульсной характеристики, определили:

) Объект характеризуется транспортным запаздыванием 0,5 мес.


1.3 Анализ частотных характеристик


Целью является исследование влияния интервалов нагрузки и воздействия на объект по изменению по изменению дохода предприятия и определение рабочего диапазона объекта.

При известных частотных характеристиках типовых моделей, характеристика объекта в целом принимает вид:


где


- характеристика транспортной модели:



- характеристика сезонной модели:



В итоге:



Построим графики частотных характеристик:


wT, месA(w, T)ф(w, T)tp(w, T), мес0Беск.0,086-3,14159-0,162,831850,085149-3,19159-0,5570380,231,415930,082692-3,24159-0,282882320,320,943950,078899-3,29159-0,19149710,415,707960,074138-3,34159-0,145804490,512,566370,0688-3,39159-0,118388920,610,471980,063235-3,44159-0,100111870,78,9759790,057718-3,49159-0,087056840,87,8539820,052439-3,54159-0,077265570,96,9813170,047514-3,59159-0,0696501316,2831850,043-3,64159-0,063557781,15,7119870,038914-3,69159-0,058573131,25,2359880,035246-3,74159-0,054419261,34,8332190,03197-3,79159-0,050904441,44,487990,029054-3,84159-0,047891741,54,188790,026462-3,89159-0,045280741,63,9269910,024157-3,94159-0,042996111,73,6959910,022108-3,99159-0,040980261,83,4906590,020283-4,04159-0,039188391,93,306940,018655-4,09159-0,0375851423,1415930,0172-4,14159-0,036142212,12,9919930,015896-4,19159-0,034836712,22,8559930,014726-4,24159-0,033649892,32,731820,013672-4,29159-0,032566272,42,6179940,012722-4,34159-0,031572952,52,5132740,011862-4,39159-0,03065912,62,416610,011082-4,44159-0,029815542,72,3271060,010374-4,49159-0,029034472,82,2439950,009729-4,54159-0,0283091948,70,1290183,62E-05-27,4916-0,0098525448,80,1287543,61E-05-27,5416-0,0098502448,90,128493,6E-05-27,5916-0,00984794490,1282283,58E-05-27,6416-0,0098456549,10,1279673,57E-05-27,6916-0,0098433749,20,1277073,55E-05-27,7416-0,009841149,30,1274483,54E-05-27,7916-0,0098388449,40,127193,52E-05-27,8416-0,0098365949,50,1269333,51E-05-27,8916-0,0098343549,60,1266773,49E-05-27,9416-0,0098321149,70,1264223,48E-05-27,9916-0,0098298949,80,1261683,47E-05-28,0416-0,009827671940,0323882,28E-06-100,142-0,00900928




Исходя из АЧХ максимальная эффективность соответствует номинальному режиму и определяется номинальной производительностью:



Изменение интервала инвестиций приводит к снижению эффективности по экспоненциальному закону.

Определим рабочую область, которой соответствует снижение эффективности до уровня рентабельности:



Такой эффективности соответствует нагрузка или интервал инвестирования мес. Это интервал кризисной ситуации.

Исходя из ФЧХ, производство характеризуется застойным режимом на всем диапазоне , т.к. , при этом с ростом нагрузки застой усугубляется и рабочему интервалу инвестирования соответствует застойный режим или время застоя


мес.


И тогда любой технологический регламент увеличивается в 0,009 раз.

Например, обеспеченный по условию задачи рост производительности в течение 6 мес., в случае кризиса будет выходить на ту же производительность в течение 6,054 мес.

Таким образом, изменение регламента весьма слабо влияет на временные характеристики ОУ.

Выводы по 1-му разделу:

) Модель объекта управления характеризуется колебательной моделью с запаздыванием.

) Объект статический, с колебательными процессами, характеризующимися самовыравниванием.

) Изменение режимов работы объекта приводит к снижению эффективности, почти не влияет на временные характеристики ОУ. Кризисная ситуация соответствует периоду инвестирования более 4,8 мес.


2. Синтез системы управления


Целью является в соответствии с проведенным анализом объекта управления сформировать вариант СУ на основе типовых законов управления. управление система математический


2.1 Определение критерия управления, ограничений, внешних факторов, выбор типового закона управления, расчет настроечных параметров


Для определения рационального подхода к управлению производительностью по каналу инвестиции - доход приведем следующие обоснования:

) По условию задачи требуется стабилизировать производительность на заданном уровне:


5,83 млн.р./мес.


при этом внешние факторы в виде дополнительных затрат являются второстепенной задачей, поэтому выбираем принцип управления по отклонению, где рассматривается отклонение фактической производительности от заданной:



Таким образом, общая структура СУ:



2) Модель ОУ характеризуется параметрами:


мес.

мес.


Для выбора метода управления рассчитаем динамические показатели:


, где


- допустимая статистическая ошибка при стабилизации производительности:


млн.р./мес.


- максимальные внешние факторы в виде дополнительных финансовых затрат.

По условию задачи:


Тогда:


Рассчитаем относительное запаздывание, которое характеризует соотношение транспортного запаздывания к инерционности объекта.



В данном случае наиболее подходящим является непрерывный метод управления на базе типовых законов.


2.2 Построение структурной схемы системы управления


Выбор типового алгоритма управления возможен на основе интуитивного подхода, при котором логически определяется необходимость тех или иных составляющих управления.

Оперативность необходима, т.к. надо быстро принимать решения при изменениях объема инвестирования.

Рост дохода сезонный или вялотекущий процесс, поэтому достаточно времени для анализа рационального управления.

Возможно долгосрочное прогнозирование изменения инвестиций.

Таким образом, выбираем ПИД-закон управления.

На основе формального подхода, при котором алгоритм управления выбирается на основе динамических показателей и требований к качеству управления по справочным номограммам [А.Д. Копелович Теория автоматического управления. Краткий справочник].

В качестве требований к типовому качеству управления определяем переходный процесс с 20-ти процентным перерегулированием, т.к. доход допускает небольшие колебания.

Исходя из номограмм, выбирается ПИД-закон:



Таким образом, окончательная структура СУ имеет вид:



) Рассчитаем настроечные коэффициенты алгоритма управления, исходя из справочных данных.


мес.

мес.


Перейдем от полученных настроек к и :


мес.


Характеризует период анализа производительности и составляет 0,0358 мес.


мес.


Характеризует период прогноза на 5,5 мес.



) Алгоритм работы СУ в целом определен следующим циклом управления:

. Формируется новое целевое назначение или задание

. Расчет отклонения фактического состояния от задания :



. Расчет управленческого решения на основе алгоритма управления :



. Расчет фактического входного воздействия на объект с учетом наложения на управление внешних факторов w:



.Расчет фактического состояния объекта y исходя из модели объекта:



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

.1. Учет ограничений по входному воздействию:

Если , то принимаем ;

Если , то принимаем ;

Если , то:



.2. Учет ограничений по выходному состоянию:

Если , то принимаем ;

Если , то принимаем ;

Если , то:


3. Анализ системы управления СУ


3.1 Анализ устойчивости


Система управления является замкнутой, поэтому применяем для исследования критерий Найквиста.

Передаточная функция разомкнутой системы:



Отсюда частотная характеристика разомкнутой системы:



Рассчитаем частотную характеристику УС, исходя из передаточной функции:


p>jw:


Тогда АЧХ УС:



ФЧХ УС:



Окончательно имеем частотную характеристику разомкнутой системы:



Построим годограф разомкнутой системы:


wxy0--0,1-33,728-1,686220,2-16,3656-1,62690,3-10,3998-1,53540,4-7,32296-1,420810,5-5,43528-1,29280,6-4,16698-1,160050,7-3,26901-1,02940,8-2,61214-0,905570,9-2,12086-0,791411-1,7473-0,688251,1-1,45937-0,596351,2-1,23474-0,51531,3-1,05753-0,444281,4-0,91619-0,382271,5-0,80226-0,3282318,5-0,004480,02548418,6-0,003180,02553418,7-0,001890,02552118,8-0,000610,02544418,90,0006570,025305190,001910,02510519,10,0031460,02484519,20,004360,02452519,30,005550,02414819,40,0067140,02371519,50,007850,02322819,60,0089540,02268719,70,0100240,02209719,80,0110580,021457


Годограф не охватывает критическую точку (-1;j0), поэтому разработанная система является устойчивой и работоспособной.

Оценим надежность системы, исходя из запасов устойчивости, и выработаем рекомендации по сохранению такой надежности.

Определим показатели запаса устойчивости:

По амплитуде:



Поэтому система является заведомо устойчивой и надежной при больших колебаниях входных и выходных параметров объекта в диапазоне от номинальных величин, а именно инвестиции могут колебаться в диапазоне:


%/мес., что соответствует условию задачи.


По частоте воздействия:



СУ обладает ограниченной надежностью по всем временным регламентам работы производства, поэтому необходимо рекомендовать стабилизацию во времени всех параметров ОУ, а именно:

) Время инвестирования не должно отклоняться:


мес.


) Время процесса, характеризующее рост дохода:


мес.


В том числе четко соблюдать временные графики поставок, сбыта продукции и другие временные регламенты, влияющие на доход от импорта.


3.2 Анализ качества управления в типовой системе


Цель: расчет переходных процессов в СУ в разных режимах функционирования и определение соответствия показателей качества заданным критериям.

При рассчитанных выше коэффициентах ПИД-закона на практике мы получаем нестабильную систему:



Поэтому коэффициенты ПИД-закона были определены эмпирически:


мес.

мес.


) Рассчитаем переходный процесс при изменении заданного дохода предприятия, который в номинальном режиме 5,65 млн.р./мес. и требуемый y* = 5,83 млн.р./мес.

При этом примем условия, что внешние факторы в виде дополнительных затрат отсутствуют:


= 0%/мес.


За основу принимаем алгоритм работы СУ из пункта 2.2.


Таблица расчета переходного процесса.

lt, месyид(t)dy(t)uп(t)uи(t)uд(t)u(t)w(t)Vид(t)y*(t)005,650000001005,6510,15,67950,15052,25750,150537,62540,033040,0335,8320,25,737910,092091,381350,24259-14,6025-12,97860-12,97865,8330,35,8243570,0056430,0846480,248233-21,6117-21,27880-21,27885,8340,45,937696-0,1077-1,615430,140538-28,3347-29,80960-29,80965,8350,56,076524-0,24652-3,69786-0,10599-34,7071-38,5110-38,5115,8360,66,187627-0,35763-5,36441-0,46361-27,7758-33,60380-33,60385,8370,76,224582-0,39458-5,91873-0,8582-9,2386-16,01550-16,01555,8380,86,180621-0,35062-5,25931-1,2088210,990184,52205404,5220545,8390,96,049657-0,21966-3,29486-1,4284732,740928,01756028,017565,831015,8263870,0036130,054189-1,4248655,8174854,44681054,446815,83111,15,5181870,3118134,677195-1,1130577,050180,61425080,614255,83121,25,1441130,68588710,2883-0,4271693,51841103,37950103,37955,83131,34,7262281,10377216,556580,67661104,4713121,70450121,70455,83141,44,2893551,54064523,109682,217255109,2183134,54530134,54535,83151,53,8607811,96921929,538294,186474107,1435140,86820140,86825,83161,63,4672132,36278735,44186,54926198,39194140,3830140,3835,83171,73,1318152,69818540,472779,24744683,84945133,56970133,56975,83181,82,8731192,95688144,3532112,2043364,67404121,23160121,23165,8349149,15,833,45E-075,18E-0667,79074,94E-0567,79075067,790755,8349249,25,8299995,05E-077,58E-0667,79074E-0567,79074067,790745,8349349,35,8299996,18E-079,26E-0667,79072,81E-0567,79073067,790735,8349449,45,8299996,77E-071,02E-0567,79071,48E-0567,79072067,790725,8349549,55,8299996,81E-071,02E-0567,79071,15E-0667,79071067,790715,8349649,65,8299996,34E-079,51E-0667,7907-1,2E-0567,7907067,79075,8349749,75,8299995,42E-078,13E-0667,7907-2,3E-0567,79068067,790685,8349849,85,834,14E-076,21E-0667,7907-3,2E-0567,79067067,790675,83



Исходя из вида переходных процессов, определим показатели качества управления:

Статическая ошибка:


млн.р./мес.


Время управления:


10 мес.


Максимальное динамическое отклонение:


= 3 млн.р./мес.


Перерегулирование:



Степень затухания колебаний:



Модульная интегральная ошибка:


74,8


Вывод: за первые 2 месяца управления очень сильно колеблется (от -40%/мес. до 140%/мес.), но затем стабилизируется на уровне 67%/мес., что на 37%/мес. меньше номинального значения и на 40% меньше заданного значения без СУ.

) Рассмотрим режим функционирования при условии, что заданная производительность не меняется, т.е. млн.р./мес., но при этом присутствуют дополнительны затраты, определяемые постановкой задачи:



Таблица расчета переходного процесса.

lt, месyид(t)dy(t)uп(t)uи(t)uд(t)u(t)w(t)Vид(t)y*(t)005,650000001005,6510,15,6795-0,0295-0,4425-0,0295-7,375-7,847-30-37,8475,6520,25,73791-0,08791-1,31865-0,11741-14,6025-16,0386-30-46,03865,6530,35,824357-0,17436-2,61535-0,29177-21,6117-24,5188-30-54,51885,6540,45,937696-0,2877-4,31543-0,57946-28,3347-33,2296-30-63,22965,6550,56,076524-0,42652-6,39786-1,00599-34,7071-42,111-30-72,1115,6560,66,12065-0,47065-7,05976-1,47664-11,0316-19,568-30-49,5685,6570,76,063536-0,41354-6,20304-1,8901714,278616,185401-30-23,81465,6580,85,899471-0,24947-3,74207-2,1396441,0162135,1345-305,1345025,6590,95,6236750,0263250,394878-2,1133268,9490967,23065-3037,230655,651015,2323840,4176166,264237-1,695797,82264102,3912-3072,391175,65111,14,7500540,89994613,49919-0,79576120,5825133,2859-30103,28595,65121,24,2045661,44543421,68150,649676136,3719158,7031-30128,70315,65131,33,6269042,02309630,346452,672773144,4157177,4349-30147,43495,65141,43,0507672,59923338,98855,272006144,0342188,2947-30158,29475,65151,52,512143,1378647,06798,409866134,6567190,1345-30160,13455,65161,62,0426043,60739654,1109412,01726117,384183,5122-30153,51225,65490495,653,76E-075,64E-0695,697676,72E-0595,69774-3065,697745,6549149,15,6499996E-079,01E-0695,697675,62E-0595,69774-3065,697745,6549249,25,6499997,67E-071,15E-0595,697674,16E-0595,69773-3065,697735,6549349,35,6499998,65E-071,3E-0595,697672,47E-0595,69771-3065,697715,6549449,45,6499998,93E-071,34E-0595,697676,88E-0695,69769-3065,697695,6549549,55,6499998,51E-071,28E-0595,69768-1E-0595,69768-3065,697685,6549649,65,6499997,48E-071,12E-0595,69768-2,6E-0595,69766-3065,697665,6549749,75,6499995,94E-078,91E-0695,69768-3,8E-0595,69765-3065,697655,6549849,85,654,04E-076,07E-0695,69768-4,7E-0595,69764-3065,697645,65



Исходя из вида переходных процессов, определим показатели качества управления:

Статическая ошибка:


млн.р./мес.


Время управления:


10 мес.


Максимальное динамическое отклонение:


= 4,5 млн.р./мес.


Перерегулирование:



Степень затухания колебаний:



Модульная интегральная ошибка:


104,9


Вывод: графики переходных характеристик аналогичны предыдущим, где велся расчет по заданному значению. Значения показателей качества управления отличаются незначительно, кроме модульной интегральной ошибки, которая увеличилась на 30.

) Рассчитаем переходный процесс при изменении заданного дохода предприятия, который в номинальном режиме 5,65 млн.р./мес. и требуемый y* = 5,83 млн.р./мес.

При этом примем условия, что внешние факторы в виде дополнительных затрат отсутствуют:


= 0%/мес.


Также введем ограничения на и y из пунктов 6.1. и 6.2. алгоритма работы СУ.


Таблица расчета переходного процесса.

lt, месyфак(t)dy(t)uп(t)uи(t)uд(t)u(t)w(t)Vфак(t)y*(t)005,650000001005,6510,15,67950,15052,25750,150537,62540,033040,0335,8320,25,737910,092091,381350,24259-14,6025-12,97860255,8330,35,8243570,0056430,0846480,248233-21,6117-21,27880255,8340,45,937696-0,1077-1,615430,140538-28,3347-29,80960255,8350,56,076524-0,24652-3,69786-0,10599-34,7071-38,5110255,8360,66,187627-0,35763-5,36441-0,46361-27,7758-33,60380255,8370,76,257243-0,42724-6,40865-0,89086-17,404-24,70350255,8380,86,285091-0,45509-6,82636-1,34595-6,96184-15,13420255,8390,96,271309-0,44131-6,61963-1,787263,4455-4,961390255,831016,216451-0,38645-5,79677-2,1737113,714325,7438390255,83111,16,121478-0,29148-4,37217-2,4651923,743316,905940255,83121,25,98774-0,15774-2,3661-2,6229333,4345628,44554028,445545,83131,35,8169620,0130380,195576-2,6098942,6945740,28026040,280265,83141,45,6112220,2187783,281677-2,3911151,4350352,3256052,32565,83151,55,3729270,4570736,856101-1,9340459,5737364,4958064,49585,83161,65,1047850,72521510,87822-1,2088267,0353176,70471076,704715,83171,74,8127411,01725915,25889-0,1915673,0111388,07845088,078455,83181,84,510131,3198719,798041,12830875,6526196,57897096,578975,83191,94,2104451,61955524,293332,74786474,92141101,96260101,96265,832023,9271181,90288228,543234,65074670,83171104,02570104,02575,83290295,830097-9,7E-05-0,0014567,790160,01071467,79942067,799425,8329129,15,830049-4,9E-05-0,0007467,790110,01189567,80126067,801265,8329229,25,830001-9E-07-1,3E-0567,790110,01210967,8022067,80225,8329329,35,8299554,47E-050,0006767,790150,01139867,80222067,802225,8329429,45,8299168,42E-050,00126367,790240,00986867,80137067,801375,8329529,55,8298850,0001150,00172467,790350,00768467,79976067,799765,8329629,65,8298650,0001350,00202767,790490,00504967,79756067,797565,8329729,75,8298560,0001440,00215867,790630,00218767,79497067,794975,8329829,85,8298590,0001410,00211767,79077-0,0006767,79222067,792225,83



Исходя из вида переходных процессов, определим показатели качества управления:

Статическая ошибка:


млн.р./мес.


Время управления:


8 мес.


Максимальное динамическое отклонение:


= 2,8 млн.р./мес.


Перерегулирование:



Степень затухания колебаний:



Модульная интегральная ошибка:


74,1


Вывод: переходные характеристики при управлении по заданному с ограничениями практически не отличаются от таковых без ограничений. Единственное существенное отличие в сокращении времени управления до 8 мес.

) Рассмотрим режим функционирования при условии, что заданная производительность не меняется, т.е. млн.р./мес., но при этом присутствуют дополнительны затраты, определяемые постановкой задачи:



Также введем ограничения на и y из пунктов 6.1. и 6.2. алгоритма работы СУ.


Таблица расчета переходного процесса.

lt, месyфак(t)dy(t)uп(t)uи(t)uд(t)u(t)w(t)Vфак(t)y*(t)005,650000001005,6510,15,6795-0,0295-0,4425-0,0295-7,375-7,847-30255,6520,25,73791-0,08791-1,31865-0,11741-14,6025-16,0386-30255,6530,35,824357-0,17436-2,61535-0,29177-21,6117-24,5188-30255,6540,45,937696-0,2877-4,31543-0,57946-28,3347-33,2296-30255,6550,56,076524-0,42652-6,39786-1,00599-34,7071-42,111-30255,6560,66,174699-0,5247-7,87048-1,53069-24,5437-33,9449-30255,6570,76,231645-0,58165-8,72468-2,11233-14,2365-25,0735-30255,6580,86,247205-0,59721-8,95808-2,70954-3,89006-15,5577-30255,6590,96,221638-0,57164-8,57457-3,281176,391852-5,46389-30255,651016,15561-0,50561-7,58415-3,7867816,507035,136099-30255,65111,16,050186-0,40019-6,00279-4,1869726,3559816,16623-30255,65121,25,906814-0,25681-3,85221-4,4437835,8428927,54689-30255,65131,35,727308-0,07731-1,15962-4,5210944,8764939,19578-30255,65141,45,5138240,1361762,042635-4,3849253,37151,02872-30255,65151,55,2688370,3811635,717446-4,0037561,2468562,96054-3032,960545,65161,64,9951110,6548899,823335-3,3488668,4314774,90595-3044,905955,65171,74,6956710,95432914,31493-2,3945474,8599486,78033-3056,780335,65181,84,3737691,27623119,14346-1,118380,4755298,50067-3068,500675,65191,94,0328481,61715224,257270,49884785,23018109,9863-3079,98635,652023,6833551,96664529,499682,46549387,37349119,3387-3089,338665,65290295,650164-0,00016-0,0024695,69802-0,0071495,68842-3065,688425,6529129,15,650179-0,00018-0,0026895,69784-0,0036695,6915-3065,69155,6529229,25,650179-0,00018-0,0026995,69767-0,0001195,69487-3065,694875,6529329,35,650166-0,00017-0,0024995,69750,00324895,69825-3065,698255,6529429,45,650142-0,00014-0,0021295,697360,00616495,7014-3065,70145,6529529,55,650108-0,00011-0,0016295,697250,00844795,70408-3065,704085,6529629,65,650068-6,8E-05-0,0010295,697180,00995995,70612-3065,706125,6529729,75,650025-2,5E-05-0,0003895,697160,01063395,70741-3065,707415,6529829,85,6499841,65E-050,00024795,697170,01046695,70789-3065,707895,65




Исходя из вида переходных процессов, определим показатели качества управления:

Статическая ошибка:


млн.р./мес.


Время управления:


8 мес.


Максимальное динамическое отклонение:


= 3,8 млн.р./мес.


Перерегулирование:



Степень затухания колебаний:



Модульная интегральная ошибка:


104,4


Вывод: переходные характеристики при управлении по отклонению с ограничениями практически не отличаются от таковых без ограничений. Единственное существенное отличие в сокращении времени управления до 8 мес.

3.3 Анализ чувствительности СУ к изменениям параметров модели объекта исследования


Для нестационарных условий производственных систем последним свойственно со временем изменять структуру или параметры, т.е. характеристики объекта плавают во времени, что приводит к изменению надежности и качества. Целью анализа чувствительности является показать влияние такого изменения на общие характеристики СУ.

Так для объекта с моделью:



все 3 характеристики , , могут изменяться во времени:


(, , )


По условию задачи свойства объекта могут изменяться в пределах , поэтому и параметры объекта могут плавать в диапазоне:


мес.

мес.


Порядок расчета чувствительности:

. Расчет влияния коэффициента передачи на надежность. Построим таблицу по полученным данным о запасе устойчивости по и .


-50%-25%+0%+25%+50%90%82%75%70%63%6,90%6,90%6,90%6,90%6,90%

2. Аналогично для .


-50%-25%0%+25%+50%27%66%75%77%80%2,00%4,00%6,90%7,00%7,10%

3. Аналогично для .


-50%-25%0%+25%+50%93%87%75%56%20%11,00%9,00%6,90%4,00%0,50%

4. По сводным таблицам построим графики чувствительности и к параметрам объекта.



. Выводы по анализу чувствительности:

Параметр имеет обратную зависимость с и не влияет на .

Параметр имеет прямую зависимость и с и с .

Параметр имеет обратную зависимость и с и с .

4. Программный модуль


.1 Листинг


unit Unit1;

interface, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,, StdCtrls, ExtCtrls, TeeProcs, TeEngine, Chart, Grids, Series;= class(TForm): TStringGrid;: TChart;: TButton;: TButton;: TButton;: TEdit;: TLabel;: TEdit;: TLabel;: TEdit;: TLabel;: TLineSeries;: TEdit;: TLabel;: TEdit;: TLabel;Button1Click(Sender: TObject);FormActivate(Sender: TObject);Button2Click(Sender: TObject);Button3Click(Sender: TObject);

{ Private declarations }

{ Public declarations };: TForm1;, cs: integer;,ynom,time: real;

{$R *.dfm}TForm1.Button1Click(Sender: TObject);;;TForm1.FormActivate(Sender: TObject);.StringGrid1.Cells[0,0]:='l';.StringGrid1.Cells[1,0]:='t';.StringGrid1.Cells[2,0]:='y(l)';.StringGrid1.Cells[3,0]:='-0,99y(l-2)';.StringGrid1.Cells[4,0]:='1,98y(l-1)';.StringGrid1.Cells[5,0]:='0,00086V(l-5)';.StringGrid1.Cells[6,0]:='V(l)';.StringGrid1.Cells[7,0]:='y(l)результирующее';:=0.1;;TForm1.Button2Click(Sender: TObject);

//Переходная хар-ка.Series1.Clear;.StringGrid1.RowCount:=3;:=1;:=StrToInt(Form1.Edit4.Text);.StringGrid1.Cells[0,cs]:='0';.StringGrid1.Cells[1,cs]:=FloatToStr(dT*StrToFloat(Form1.StringGrid1.Cells[0,cs]));.StringGrid1.Cells[2,cs]:=Form1.Edit1.Text;.StringGrid1.Cells[3,cs]:='0';.StringGrid1.Cells[4,cs]:='0';.StringGrid1.Cells[5,cs]:='0';.StringGrid1.Cells[6,cs]:='0';.StringGrid1.Cells[7,cs]:=Form1.Edit2.Text;:=StrToFloat(Form1.Edit2.Text);.Series1.AddXY(0,ynom);:=cs+1; (cs<=steps) do

begin

//Увеличение кол-ва строк

Form1.StringGrid1.RowCount:=cs+1;

//Увеличение l.StringGrid1.Cells[0,cs]:=IntToStr(StrToInt(Form1.StringGrid1.Cells[0,cs-1])+1);

//Расчет t.StringGrid1.Cells[1,cs]:=FloatToStr(dT*StrToFloat(Form1.StringGrid1.Cells[0,cs]));

//Расчет -0,99y(l-2)cs<3 then.StringGrid1.Cells[3,cs]:=FloatToStr(-0.99*StrToFloat(Form1.StringGrid1.Cells[2,1])).StringGrid1.Cells[3,cs]:=FloatToStr(-0.99*StrToFloat(Form1.StringGrid1.Cells[2,cs-2]));

//Расчет 1,98y(l-1).StringGrid1.Cells[4,cs]:=FloatToStr(1.98*StrToFloat(Form1.StringGrid1.Cells[2,cs-1]));

//Расчет 0,00086V(l-5)cs<6 then.StringGrid1.Cells[5,cs]:=FloatToStr(0.00086*StrToFloat(Form1.StringGrid1.Cells[6,1])).StringGrid1.Cells[5,cs]:=FloatToStr(0.00086*StrToFloat(Form1.StringGrid1.Cells[6,cs-5]));

//Расчет V(l).StringGrid1.Cells[6,cs]:=Form1.Edit3.Text;

//Расчет y(l).StringGrid1.Cells[2,cs]:=FloatToStr(StrToFloat(Form1.StringGrid1.Cells[3,cs])+StrToFloat(Form1.StringGrid1.Cells[4,cs])+(Form1.StringGrid1.Cells[5,cs]));

//Расчет y(l)результирующее.StringGrid1.Cells[7,cs]:=FloatToStr(StrToFloat(Form1.StringGrid1.Cells[2,cs])+ynom);

//Внесение данных на график.Series1.AddXY(StrToFloat(Form1.StringGrid1.Cells[1,cs]),StrToFloat(Form1.StringGrid1.Cells[7,cs]));

//Следующий шаг

cs:=cs+1;

end;

end;TForm1.Button3Click(Sender: TObject);

//Импульсная хар-ка:=StrToFloat(Form1.Edit5.Text);.Series1.Clear;.StringGrid1.RowCount:=3;:=1;:=StrToInt(Form1.Edit4.Text);.StringGrid1.Cells[0,cs]:='0';.StringGrid1.Cells[1,cs]:=FloatToStr(dT*StrToFloat(Form1.StringGrid1.Cells[0,cs]));.StringGrid1.Cells[2,cs]:=Form1.Edit1.Text;.StringGrid1.Cells[3,cs]:='0';.StringGrid1.Cells[4,cs]:='0';.StringGrid1.Cells[5,cs]:='0';.StringGrid1.Cells[6,cs]:='0';.StringGrid1.Cells[7,cs]:=Form1.Edit2.Text;:=StrToFloat(Form1.Edit2.Text);.Series1.AddXY(0,ynom);:=cs+1; (cs<=steps) do

begin

//Увеличение кол-ва строк

Form1.StringGrid1.RowCount:=cs+1;

//Увеличение l.StringGrid1.Cells[0,cs]:=IntToStr(StrToInt(Form1.StringGrid1.Cells[0,cs-1])+1);

//Расчет t.StringGrid1.Cells[1,cs]:=FloatToStr(dT*StrToFloat(Form1.StringGrid1.Cells[0,cs]));

//Расчет -0,99y(l-2)cs<3 then.StringGrid1.Cells[3,cs]:=FloatToStr(-0.99*StrToFloat(Form1.StringGrid1.Cells[2,1])).StringGrid1.Cells[3,cs]:=FloatToStr(-0.99*StrToFloat(Form1.StringGrid1.Cells[2,cs-2]));

//Расчет 1,98y(l-1).StringGrid1.Cells[4,cs]:=FloatToStr(1.98*StrToFloat(Form1.StringGrid1.Cells[2,cs-1]));

//Расчет 0,00086V(l-5)cs<6 then.StringGrid1.Cells[5,cs]:=FloatToStr(0.00086*StrToFloat(Form1.StringGrid1.Cells[6,1])).StringGrid1.Cells[5,cs]:=FloatToStr(0.00086*StrToFloat(Form1.StringGrid1.Cells[6,cs-5]));

//Расчет V(l)StrToFloat(Form1.StringGrid1.Cells[1,cs])<=time then.StringGrid1.Cells[6,cs]:=Form1.Edit3.Text else.StringGrid1.Cells[6,cs]:=FloatToStr(0);

//Расчет y(l).StringGrid1.Cells[2,cs]:=FloatToStr(StrToFloat(Form1.StringGrid1.Cells[3,cs])+StrToFloat(Form1.StringGrid1.Cells[4,cs])+(Form1.StringGrid1.Cells[5,cs]));

//Расчет y(l)результирующее.StringGrid1.Cells[7,cs]:=FloatToStr(StrToFloat(Form1.StringGrid1.Cells[2,cs])+ynom);

//Внесение данных на график.Series1.AddXY(StrToFloat(Form1.StringGrid1.Cells[1,cs]),StrToFloat(Form1.StringGrid1.Cells[7,cs]));

//Следующий шаг

cs:=cs+1;

end;

end;.


Министерство образования и науки РФ ФГБО ВПО Сибирский государственный индустриальный университет кафедра систем информатики и управления

Больше работ по теме:

КОНТАКТНЫЙ EMAIL: [email protected]

Скачать реферат © 2017 | Пользовательское соглашение

Скачать      Реферат

ПРОФЕССИОНАЛЬНАЯ ПОМОЩЬ СТУДЕНТАМ