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

КурсоваяПомощь в написанииУзнать стоимостьмоей работы

Рассмотрим режим функционирования при условии, что заданная производительность не меняется, т. е. млн.р./мес., но при этом присутствуют дополнительны затраты, определяемые постановкой задачи: Министерство образования и науки РФ ФГБО ВПО Сибирский государственный индустриальный университет кафедра систем информатики и управления Курсовая работа Выполнил ст. гр. ИП-09 Пятаев А. М. Изменение режимов… Читать ещё >

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

4.1. Листинг

Задание

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

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

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

Дано:

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

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

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

%/мес.

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

: доход

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

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

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

= 100%/мес.

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

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

==125% мес.

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

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

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

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

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

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

=0,5 мес.

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

=6 мес.

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

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

==25%/мес.

=180%/мес.

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

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

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

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

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

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

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

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

Примем:

==30%/мес.

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

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

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

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

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

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

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

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

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

=103%

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

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

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

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

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

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

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

где время запаздывания мес. Отсюда:

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

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

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

мес. Отсюда:

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

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

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

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

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

Отсюда:

0,1 мес. 3,1 дня Далее:

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

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

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

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


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

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

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

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

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

0%/мес.

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

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

25%/мес.

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

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

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


l	t	y (l)	-0,99y (l-2)	1,98y (l-1)	0,00086V (l-5)	V (l)	y (l) рез
							5,65
	0,1						5,65
	0,2						5,65
	0,3						5,65
	0,4						5,65
	0,5						5,65
	0,6	0,21 288			0,212 875		5,671 288
	0,7	0,63 437		0,4 214 925	0,212 875		5,713 437
	0,8	0,125 818	— 0,2 107 463	0,125 604 765	0,212 875		5,775 818
	0,9	0,207 604	— 0,6 280 238	0,249 118 927	0,212 875		5,857 604
		0,307 784	— 0,12 455 946	0,411 056 009	0,212 875		5,957 784
	1,1	0,425 172	— 0,205 528	0,609 412 409	0,212 875		6,75 172
	1,2	0,558 422	— 0,3 047 062	0,841 840 371	0,212 875		6,208 422
	1,3	0,706 042	— 0,42 092 019	1,1 056 749	0,212 875		6,356 042
	1,4	0,866 414	— 0,55 283 745	1,397 963 585	0,212 875		6,516 414
	1,5	1,37 805	— 0,69 898 179	1,715 498 997	0,212 875		6,687 805
	1,6	1,218 391	— 0,8 577 495	2,54 853 315	0,212 875		6,868 391
	1,7	1,406 276	— 1,2 742 666	2,412 414 806	0,212 875		7,56 276
	1,8	1,599 506	— 1,2 062 074	2,784 425 785	0,212 875		7,249 506
	1,9	1,796 096	— 1,39 221 289	3,167 021 646	0,212 875		7,446 096
		1,994 047	— 1,58 351 082	3,556 270 582	0,212 875		7,644 047
	2,1	2,191 366	— 1,77 813 529	3,948 213 573	0,212 875		7,841 366
	2,2	2,386 085	— 1,97 410 679	4,338 904 248	0,212 875		8,36 085
	2,3	2,576 284	— 2,16 945 212	4,724 448 223	0,212 875		8,226 284
	2,4	2,760 105	— 2,36 222 411	5,101 041 527	0,212 875		8,410 105
	2,5	2,935 774	— 2,55 052 076	5,465 007 733	0,212 875		8,585 774
	2,6	3,101 617	— 2,73 250 387	5,812 833 449	0,212 875		8,751 617
	2,7	3,256 073	— 2,90 641 672	6,141 201 824	0,212 875		8,906 073
	2,8	3,39 771	— 3,7 060 091	6,447 023 746	0,212 875		9,4 771
	2,9	3,525 242	— 3,22 351 187	6,727 466 462	0,212 875		9,175 242
		2,527 668	— 2,42 760 431	4,933 985 256	0,212 875		8,177 668
	28,1	2,559 078	— 2,46 699 263	5,478 353	0,212 875		8,209 078
	28,2	2,585 871	— 2,50 239 176	5,66 975 236	0,212 875		8,235 871
	28,3	2,607 824	— 2,53 348 762	5,120 024 522	0,212 875		8,257 824
	28,4	2,624 768	— 2,56 001 226	5,163 492 321	0,212 875		8,274 768
	28,5	2,636 581	— 2,58 174 616	5,197 039 768	0,212 875		8,286 581
	28,6	2,643 198	— 2,59 851 988	5,220 430 593	0,212 875		8,293 198
	28,7	2,644 605	— 2,6 102 153	5,233 532 454	0,212 875		8,294 605
	28,8	2,640 838	— 2,61 676 623	5,236 317 222	0,212 875		8,290 838
	28,9	2,631 989	— 2,61 815 861	5,22 886 022	0,212 875		8,281 989
		2,618 196	— 2,61 443 011	5,211 338 435	0,212 875		8,268 196
	29,1	2,599 646	— 2,60 566 922	5,184 027 734	0,212 875		8,249 646
	29,2	2,576 573	— 2,59 201 387	5,147 299 113	0,212 875		8,226 573
	29,3	2,549 252	— 2,57 364 956	5,101 614 037	0,212 875		8,199 252
	29,4	2,517 999	— 2,55 080 702	5,47 518 922	0,212 875		8,167 999
	29,5	2,483 167	— 2,52 375 946	4,985 638 818	0,212 875		8,133 167
	29,6	2,445 138	— 2,49 281 941	4,916 670 377	0,212 875		8,95 138
	29,7	2,404 326	— 2,45 833 519	4,841 374 167	0,212 875		8,54 326

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

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

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

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

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

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

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

0%/мес.

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

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

25%/мес.

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

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

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


l	t	y (l)	-0,99y (l-2)	1,98y (l-1)	0,00086V (l-5)	V (l)	y (l) рез
							5,65
	0,1						5,65
	0,2						5,65
	0,3						5,65
	0,4						5,65
	0,5						5,65
	0,6	0,21 288			0,212 875		5,671 288
	0,7	0,63 437		0,4 214 925	0,212 875		5,713 437
	0,8	0,125 818	— 0,2 107 463	0,125 604 765	0,212 875		5,775 818
	0,9	0,207 604	— 0,6 280 238	0,249 118 927	0,212 875		5,857 604
		0,307 784	— 0,12 455 946	0,411 056 009	0,212 875		5,957 784
	1,1	0,425 172	— 0,205 528	0,609 412 409	0,212 875		6,75 172
	1,2	0,558 422	— 0,3 047 062	0,841 840 371	0,212 875		6,208 422
	1,3	0,706 042	— 0,42 092 019	1,1 056 749	0,212 875		6,356 042
	1,4	0,866 414	— 0,55 283 745	1,397 963 585	0,212 875		6,516 414
	1,5	1,37 805	— 0,69 898 179	1,715 498 997	0,212 875		6,687 805
	1,6	1,197 104	— 0,8 577 495	2,54 853 315			6,847 104
	1,7	1,342 839	— 1,2 742 666	2,370 265 556			6,992 839
	1,8	1,473 688	— 1,18 513 278	2,65 882 102			7,123 688
	1,9	1,588 492	— 1,32 941 051	2,917 902 719			7,238 492
		1,686 263	— 1,45 895 136	3,145 214 573			7,336 263
	2,1	1,766 194	— 1,57 260 729	3,338 801 164			7,416 194
	2,2	1,827 663	— 1,66 940 058	3,497 063 876			7,477 663
	2,3	1,870 241	— 1,74 853 194	3,618 773 323			7,520 241
	2,4	1,893 691	— 1,80 938 666	3,703 077 942			7,543 691
	2,5	1,89 797	— 1,85 153 897	3,749 508 736			7,54 797
	2,6	1,883 226	— 1,87 475 437	3,757 980 134			7,533 226
	27,9	0,508 025	— 0,52 774 202	1,35 766 932			6,158 025
		0,488 006	— 0,51 788 347	1,5 889 331			6,138 006
	28,1	0,463 307	— 0,50 294 467	0,966 251 612			6,113 307
	28,2	0,434 222	— 0,48 312 581	0,917 347 754			6,84 222
	28,3	0,401 086	— 0,45 867 388	0,859 759 458			6,51 086
	28,4	0,36 427	— 0,42 987 973	0,79 414 945			6,1 427
	28,5	0,324 179	— 0,39 707 472	0,721 254 047			5,974 179
	28,6	0,281 248	— 0,36 062 702	0,641 875 058			5,931 248
	28,7	0,235 934	— 0,32 093 753	0,556 871 108			5,885 934
	28,8	0,188 713	— 0,27 843 555	0,467 148 487			5,838 713
	28,9	0,140 077	— 0,23 357 424	0,373 651 607			5,790 077
		0,90 527	— 0,1 868 258	0,27 735 318			5,740 527
	29,1	0,40 568	— 0,13 867 659	0,179 244 205			5,690 568
	29,2	— 0,0093	— 0,896 221	0,80 323 878			5,640 702
	29,3	— 0,5 857	— 0,4 016 194	— 0,18 410 485			5,591 428
	29,4	— 0,10 677	0,920 524	— 0,115 973 399			5,543 232
	29,5	— 0,15 341	0,579 867	— 0,21 140 095			5,496 586
	29,6	— 0,19 806	0,10 570 047	— 0,303 760 216			5,45 194
	29,7	— 0,24 028	0,15 188 011	— 0,392 158 287			5,409 722

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

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

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

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

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

где

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

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

В итоге:

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


w	T, мес	A (w, T)	ф (w, T)	tp (w, T), мес
	Беск.	0,086	— 3,14 159	;
0,1	62,83 185	0,85 149	— 3,19 159	— 0,557 038
0,2	31,41 593	0,82 692	— 3,24 159	— 0,28 288 232
0,3	20,94 395	0,78 899	— 3,29 159	— 0,1 914 971
0,4	15,70 796	0,74 138	— 3,34 159	— 0,14 580 449
0,5	12,56 637	0,0688	— 3,39 159	— 0,11 838 892
0,6	10,47 198	0,63 235	— 3,44 159	— 0,10 011 187
0,7	8,975 979	0,57 718	— 3,49 159	— 0,8 705 684
0,8	7,853 982	0,52 439	— 3,54 159	— 0,7 726 557
0,9	6,981 317	0,47 514	— 3,59 159	— 0,6 965 013
	6,283 185	0,043	— 3,64 159	— 0,6 355 778
1,1	5,711 987	0,38 914	— 3,69 159	— 0,5 857 313
1,2	5,235 988	0,35 246	— 3,74 159	— 0,5 441 926
1,3	4,833 219	0,3 197	— 3,79 159	— 0,5 090 444
1,4	4,48 799	0,29 054	— 3,84 159	— 0,4 789 174
1,5	4,18 879	0,26 462	— 3,89 159	— 0,4 528 074
1,6	3,926 991	0,24 157	— 3,94 159	— 0,4 299 611
1,7	3,695 991	0,22 108	— 3,99 159	— 0,4 098 026
1,8	3,490 659	0,20 283	— 4,4 159	— 0,3 918 839
1,9	3,30 694	0,18 655	— 4,9 159	— 0,3 758 514
	3,141 593	0,0172	— 4,14 159	— 0,3 614 221
2,1	2,991 993	0,15 896	— 4,19 159	— 0,3 483 671
2,2	2,855 993	0,14 726	— 4,24 159	— 0,3 364 989
2,3	2,73 182	0,13 672	— 4,29 159	— 0,3 256 627
2,4	2,617 994	0,12 722	— 4,34 159	— 0,3 157 295
2,5	2,513 274	0,11 862	— 4,39 159	— 0,306 591
2,6	2,41 661	0,11 082	— 4,44 159	— 0,2 981 554
2,7	2,327 106	0,10 374	— 4,49 159	— 0,2 903 447
2,8	2,243 995	0,9 729	— 4,54 159	— 0,2 830 919
48,7	0,129 018	3,62E-05	— 27,4916	— 0,985 254
48,8	0,128 754	3,61E-05	— 27,5416	— 0,985 024
48,9	0,12 849	3,6E-05	— 27,5916	— 0,984 794
	0,128 228	3,58E-05	— 27,6416	— 0,984 565
49,1	0,127 967	3,57E-05	— 27,6916	— 0,984 337
49,2	0,127 707	3,55E-05	— 27,7416	— 0,98 411
49,3	0,127 448	3,54E-05	— 27,7916	— 0,983 884
49,4	0,12 719	3,52E-05	— 27,8416	— 0,983 659
49,5	0,126 933	3,51E-05	— 27,8916	— 0,983 435
49,6	0,126 677	3,49E-05	— 27,9416	— 0,983 211
49,7	0,126 422	3,48E-05	— 27,9916	— 0,982 989
49,8	0,126 168	3,47E-05	— 28,0416	— 0,982 767
	0,32 388	2,28E-06	— 100,142	— 0,900 928

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

мес.

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

где

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

млн.р./мес.

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

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

Тогда:

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

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

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

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

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

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

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

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

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

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

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

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

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

мес.

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

мес.

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

мес.

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

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

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

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

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

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

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

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

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

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

Если, то:

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

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

Если, то:

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

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

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

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

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

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

p>jw:

Тогда АЧХ УС:

ФЧХ УС:

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

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


w	x	y
	;	;
0,1	— 33,728	— 1,68 622
0,2	— 16,3656	— 1,6269
0,3	— 10,3998	— 1,5354
0,4	— 7,32 296	— 1,42 081
0,5	— 5,43 528	— 1,2928
0,6	— 4,16 698	— 1,16 005
0,7	— 3,26 901	— 1,0294
0,8	— 2,61 214	— 0,90 557
0,9	— 2,12 086	— 0,79 141
	— 1,7473	— 0,68 825
1,1	— 1,45 937	— 0,59 635
1,2	— 1,23 474	— 0,5153
1,3	— 1,5 753	— 0,44 428
1,4	— 0,91 619	— 0,38 227
1,5	— 0,80 226	— 0,32 823
18,5	— 0,448	0,25 484
18,6	— 0,318	0,25 534
18,7	— 0,189	0,25 521
18,8	— 0,61	0,25 444
18,9	0,657	0,25 305
	0,191	0,25 105
19,1	0,3 146	0,24 845
19,2	0,436	0,24 525
19,3	0,555	0,24 148
19,4	0,6 714	0,23 715
19,5	0,785	0,23 228
19,6	0,8 954	0,22 687
19,7	0,10 024	0,22 097
19,8	0,11 058	0,21 457

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

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

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

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

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

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

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

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

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

мес.

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

мес.

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

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

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

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

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

мес.

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

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

= 0%/мес.

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

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


l	t, мес	yид (t)	dy (t)	uп (t)	uи (t)	uд (t)	u (t)	w (t)	Vид (t)	y*(t)
		5,65								5,65
	0,1	5,6795	0,1505	2,2575	0,1505	37,625	40,033		40,033	5,83
	0,2	5,73 791	0,9 209	1,38 135	0,24 259	— 14,6025	— 12,9786		— 12,9786	5,83
	0,3	5,824 357	0,5 643	0,84 648	0,248 233	— 21,6117	— 21,2788		— 21,2788	5,83
	0,4	5,937 696	— 0,1077	— 1,61 543	0,140 538	— 28,3347	— 29,8096		— 29,8096	5,83
	0,5	6,76 524	— 0,24 652	— 3,69 786	— 0,10 599	— 34,7071	— 38,511		— 38,511	5,83
	0,6	6,187 627	— 0,35 763	— 5,36 441	— 0,46 361	— 27,7758	— 33,6038		— 33,6038	5,83
	0,7	6,224 582	— 0,39 458	— 5,91 873	— 0,8582	— 9,2386	— 16,0155		— 16,0155	5,83
	0,8	6,180 621	— 0,35 062	— 5,25 931	— 1,20 882	10,99 018	4,522 054		4,522 054	5,83
	0,9	6,49 657	— 0,21 966	— 3,29 486	— 1,42 847	32,7409	28,1 756		28,1 756	5,83
		5,826 387	0,3 613	0,54 189	— 1,42 486	55,81 748	54,44 681		54,44 681	5,83
	1,1	5,518 187	0,311 813	4,677 195	— 1,11 305	77,0501	80,61 425		80,61 425	5,83
	1,2	5,144 113	0,685 887	10,2883	— 0,42 716	93,51 841	103,3795		103,3795	5,83
	1,3	4,726 228	1,103 772	16,55 658	0,67 661	104,4713	121,7045		121,7045	5,83
	1,4	4,289 355	1,540 645	23,10 968	2,217 255	109,2183	134,5453		134,5453	5,83
	1,5	3,860 781	1,969 219	29,53 829	4,186 474	107,1435	140,8682		140,8682	5,83
	1,6	3,467 213	2,362 787	35,4418	6,549 261	98,39 194	140,383		140,383	5,83
	1,7	3,131 815	2,698 185	40,47 277	9,247 446	83,84 945	133,5697		133,5697	5,83
	1,8	2,873 119	2,956 881	44,35 321	12,20 433	64,67 404	121,2316		121,2316	5,83
	49,1	5,83	3,45E-07	5,18E-06	67,7907	4,94E-05	67,79 075		67,79 075	5,83
	49,2	5,829 999	5,05E-07	7,58E-06	67,7907	4E-05	67,79 074		67,79 074	5,83
	49,3	5,829 999	6,18E-07	9,26E-06	67,7907	2,81E-05	67,79 073		67,79 073	5,83
	49,4	5,829 999	6,77E-07	1,02E-05	67,7907	1,48E-05	67,79 072		67,79 072	5,83
	49,5	5,829 999	6,81E-07	1,02E-05	67,7907	1,15E-06	67,79 071		67,79 071	5,83
	49,6	5,829 999	6,34E-07	9,51E-06	67,7907	— 1,2E-05	67,7907		67,7907	5,83
	49,7	5,829 999	5,42E-07	8,13E-06	67,7907	— 2,3E-05	67,79 068		67,79 068	5,83
	49,8	5,83	4,14E-07	6,21E-06	67,7907	— 3,2E-05	67,79 067		67,79 067	5,83

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

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

млн.р./мес.

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

10 мес.

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

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

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

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

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

74,8

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

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

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


l	t, мес	yид (t)	dy (t)	uп (t)	uи (t)	uд (t)	u (t)	w (t)	Vид (t)	y*(t)
		5,65								5,65
	0,1	5,6795	— 0,0295	— 0,4425	— 0,0295	— 7,375	— 7,847	— 30	— 37,847	5,65
	0,2	5,73 791	— 0,8 791	— 1,31 865	— 0,11 741	— 14,6025	— 16,0386	— 30	— 46,0386	5,65
	0,3	5,824 357	— 0,17 436	— 2,61 535	— 0,29 177	— 21,6117	— 24,5188	— 30	— 54,5188	5,65
	0,4	5,937 696	— 0,2877	— 4,31 543	— 0,57 946	— 28,3347	— 33,2296	— 30	— 63,2296	5,65
	0,5	6,76 524	— 0,42 652	— 6,39 786	— 1,599	— 34,7071	— 42,111	— 30	— 72,111	5,65
	0,6	6,12 065	— 0,47 065	— 7,5 976	— 1,47 664	— 11,0316	— 19,568	— 30	— 49,568	5,65
	0,7	6,63 536	— 0,41 354	— 6,20 304	— 1,89 017	14,27 861	6,185 401	— 30	— 23,8146	5,65
	0,8	5,899 471	— 0,24 947	— 3,74 207	— 2,13 964	41,1 621	35,1345	— 30	5,134 502	5,65
	0,9	5,623 675	0,26 325	0,394 878	— 2,11 332	68,94 909	67,23 065	— 30	37,23 065	5,65
		5,232 384	0,417 616	6,264 237	— 1,6957	97,82 264	102,3912	— 30	72,39 117	5,65
	1,1	4,750 054	0,899 946	13,49 919	— 0,79 576	120,5825	133,2859	— 30	103,2859	5,65
	1,2	4,204 566	1,445 434	21,6815	0,649 676	136,3719	158,7031	— 30	128,7031	5,65
	1,3	3,626 904	2,23 096	30,34 645	2,672 773	144,4157	177,4349	— 30	147,4349	5,65
	1,4	3,50 767	2,599 233	38,9885	5,272 006	144,0342	188,2947	— 30	158,2947	5,65
	1,5	2,51 214	3,13 786	47,0679	8,409 866	134,6567	190,1345	— 30	160,1345	5,65
	1,6	2,42 604	3,607 396	54,11 094	12,1 726	117,384	183,5122	— 30	153,5122	5,65
		5,65	3,76E-07	5,64E-06	95,69 767	6,72E-05	95,69 774	— 30	65,69 774	5,65
	49,1	5,649 999	6E-07	9,01E-06	95,69 767	5,62E-05	95,69 774	— 30	65,69 774	5,65
	49,2	5,649 999	7,67E-07	1,15E-05	95,69 767	4,16E-05	95,69 773	— 30	65,69 773	5,65
	49,3	5,649 999	8,65E-07	1,3E-05	95,69 767	2,47E-05	95,69 771	— 30	65,69 771	5,65
	49,4	5,649 999	8,93E-07	1,34E-05	95,69 767	6,88E-06	95,69 769	— 30	65,69 769	5,65
	49,5	5,649 999	8,51E-07	1,28E-05	95,69 768	— 1E-05	95,69 768	— 30	65,69 768	5,65
	49,6	5,649 999	7,48E-07	1,12E-05	95,69 768	— 2,6E-05	95,69 766	— 30	65,69 766	5,65
	49,7	5,649 999	5,94E-07	8,91E-06	95,69 768	— 3,8E-05	95,69 765	— 30	65,69 765	5,65
	49,8	5,65	4,04E-07	6,07E-06	95,69 768	— 4,7E-05	95,69 764	— 30	65,69 764	5,65

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

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

млн.р./мес.

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

10 мес.

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

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

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

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

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

104,9

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

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

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

= 0%/мес.

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

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


l	t, мес	yфак (t)	dy (t)	uп (t)	uи (t)	uд (t)	u (t)	w (t)	Vфак (t)	y*(t)
		5,65								5,65
	0,1	5,6795	0,1505	2,2575	0,1505	37,625	40,033		40,033	5,83
	0,2	5,73 791	0,9 209	1,38 135	0,24 259	— 14,6025	— 12,9786			5,83
	0,3	5,824 357	0,5 643	0,84 648	0,248 233	— 21,6117	— 21,2788			5,83
	0,4	5,937 696	— 0,1077	— 1,61 543	0,140 538	— 28,3347	— 29,8096			5,83
	0,5	6,76 524	— 0,24 652	— 3,69 786	— 0,10 599	— 34,7071	— 38,511			5,83
	0,6	6,187 627	— 0,35 763	— 5,36 441	— 0,46 361	— 27,7758	— 33,6038			5,83
	0,7	6,257 243	— 0,42 724	— 6,40 865	— 0,89 086	— 17,404	— 24,7035			5,83
	0,8	6,285 091	— 0,45 509	— 6,82 636	— 1,34 595	— 6,96 184	— 15,1342			5,83
	0,9	6,271 309	— 0,44 131	— 6,61 963	— 1,78 726	3,4455	— 4,96 139			5,83
		6,216 451	— 0,38 645	— 5,79 677	— 2,17 371	13,71 432	5,743 839			5,83
	1,1	6,121 478	— 0,29 148	— 4,37 217	— 2,46 519	23,7433	16,90 594			5,83
	1,2	5,98 774	— 0,15 774	— 2,3661	— 2,62 293	33,43 456	28,44 554		28,44 554	5,83
	1,3	5,816 962	0,13 038	0,195 576	— 2,60 989	42,69 457	40,28 026		40,28 026	5,83
	1,4	5,611 222	0,218 778	3,281 677	— 2,39 111	51,43 503	52,3256		52,3256	5,83
	1,5	5,372 927	0,457 073	6,856 101	— 1,93 404	59,57 373	64,4958		64,4958	5,83
	1,6	5,104 785	0,725 215	10,87 822	— 1,20 882	67,3 531	76,70 471		76,70 471	5,83
	1,7	4,812 741	1,17 259	15,25 889	— 0,19 156	73,1 113	88,7 845		88,7 845	5,83
	1,8	4,51 013	1,31 987	19,79 804	1,128 308	75,65 261	96,57 897		96,57 897	5,83
	1,9	4,210 445	1,619 555	24,29 333	2,747 864	74,92 141	101,9626		101,9626	5,83
		3,927 118	1,902 882	28,54 323	4,650 746	70,83 171	104,0257		104,0257	5,83
		5,830 097	— 9,7E-05	— 0,145	67,79 016	0,10 714	67,79 942		67,79 942	5,83
	29,1	5,830 049	— 4,9E-05	— 0,74	67,79 011	0,11 895	67,80 126		67,80 126	5,83
	29,2	5,830 001	— 9E-07	— 1,3E-05	67,79 011	0,12 109	67,8022		67,8022	5,83
	29,3	5,829 955	4,47E-05	0,67	67,79 015	0,11 398	67,80 222		67,80 222	5,83
	29,4	5,829 916	8,42E-05	0,1 263	67,79 024	0,9 868	67,80 137		67,80 137	5,83
	29,5	5,829 885	0,115	0,1 724	67,79 035	0,7 684	67,79 976		67,79 976	5,83
	29,6	5,829 865	0,135	0,2 027	67,79 049	0,5 049	67,79 756		67,79 756	5,83
	29,7	5,829 856	0,144	0,2 158	67,79 063	0,2 187	67,79 497		67,79 497	5,83
	29,8	5,829 859	0,141	0,2 117	67,79 077	— 0,67	67,79 222		67,79 222	5,83

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

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

млн.р./мес.

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

8 мес.

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

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

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

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

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

74,1

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

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

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

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


l	t, мес	yфак (t)	dy (t)	uп (t)	uи (t)	uд (t)	u (t)	w (t)	Vфак (t)	y*(t)
		5,65								5,65
	0,1	5,6795	— 0,0295	— 0,4425	— 0,0295	— 7,375	— 7,847	— 30		5,65
	0,2	5,73 791	— 0,8 791	— 1,31 865	— 0,11 741	— 14,6025	— 16,0386	— 30		5,65
	0,3	5,824 357	— 0,17 436	— 2,61 535	— 0,29 177	— 21,6117	— 24,5188	— 30		5,65
	0,4	5,937 696	— 0,2877	— 4,31 543	— 0,57 946	— 28,3347	— 33,2296	— 30		5,65
	0,5	6,76 524	— 0,42 652	— 6,39 786	— 1,599	— 34,7071	— 42,111	— 30		5,65
	0,6	6,174 699	— 0,5247	— 7,87 048	— 1,53 069	— 24,5437	— 33,9449	— 30		5,65
	0,7	6,231 645	— 0,58 165	— 8,72 468	— 2,11 233	— 14,2365	— 25,0735	— 30		5,65
	0,8	6,247 205	— 0,59 721	— 8,95 808	— 2,70 954	— 3,89 006	— 15,5577	— 30		5,65
	0,9	6,221 638	— 0,57 164	— 8,57 457	— 3,28 117	6,391 852	— 5,46 389	— 30		5,65
		6,15 561	— 0,50 561	— 7,58 415	— 3,78 678	16,50 703	5,136 099	— 30		5,65
	1,1	6,50 186	— 0,40 019	— 6,279	— 4,18 697	26,35 598	16,16 623	— 30		5,65
	1,2	5,906 814	— 0,25 681	— 3,85 221	— 4,44 378	35,84 289	27,54 689	— 30		5,65
	1,3	5,727 308	— 0,7 731	— 1,15 962	— 4,52 109	44,87 649	39,19 578	— 30		5,65
	1,4	5,513 824	0,136 176	2,42 635	— 4,38 492	53,371	51,2 872	— 30		5,65
	1,5	5,268 837	0,381 163	5,717 446	— 4,375	61,24 685	62,96 054	— 30	32,96 054	5,65
	1,6	4,995 111	0,654 889	9,823 335	— 3,34 886	68,43 147	74,90 595	— 30	44,90 595	5,65
	1,7	4,695 671	0,954 329	14,31 493	— 2,39 454	74,85 994	86,78 033	— 30	56,78 033	5,65
	1,8	4,373 769	1,276 231	19,14 346	— 1,1183	80,47 552	98,50 067	— 30	68,50 067	5,65
	1,9	4,32 848	1,617 152	24,25 727	0,498 847	85,23 018	109,9863	— 30	79,9863	5,65
		3,683 355	1,966 645	29,49 968	2,465 493	87,37 349	119,3387	— 30	89,33 866	5,65
		5,650 164	— 0,16	— 0,246	95,69 802	— 0,714	95,68 842	— 30	65,68 842	5,65
	29,1	5,650 179	— 0,18	— 0,268	95,69 784	— 0,366	95,6915	— 30	65,6915	5,65
	29,2	5,650 179	— 0,18	— 0,269	95,69 767	— 0,11	95,69 487	— 30	65,69 487	5,65
	29,3	5,650 166	— 0,17	— 0,249	95,6975	0,3 248	95,69 825	— 30	65,69 825	5,65
	29,4	5,650 142	— 0,14	— 0,212	95,69 736	0,6 164	95,7014	— 30	65,7014	5,65
	29,5	5,650 108	— 0,11	— 0,162	95,69 725	0,8 447	95,70 408	— 30	65,70 408	5,65
	29,6	5,650 068	— 6,8E-05	— 0,102	95,69 718	0,9 959	95,70 612	— 30	65,70 612	5,65
	29,7	5,650 025	— 2,5E-05	— 0,38	95,69 716	0,10 633	95,70 741	— 30	65,70 741	5,65
	29,8	5,649 984	1,65E-05	0,247	95,69 717	0,10 466	95,70 789	— 30	65,70 789	5,65

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

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

млн.р./мес.

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

8 мес.

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

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

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

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

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

104,4

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

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

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

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

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

(, ,)

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

мес.

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

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


— 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. По сводным таблицам построим графики чувствительности и к параметрам объекта.

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

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

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

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

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

4.1 Листинг

unit Unit1;

interface

uses

Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,

Dialogs, StdCtrls, ExtCtrls, TeeProcs, TeEngine, Chart, Grids, Series;

type

TForm1 = class (TForm)

StringGrid1: TStringGrid;

Chart1: TChart;

Button1: TButton;

Button2: TButton;

Button3: TButton;

Edit1: TEdit;

Label1: TLabel;

Edit2: TEdit;

Label2: TLabel;

Edit3: TEdit;

Label3: TLabel;

Series1: TLineSeries;

Edit4: TEdit;

Label4: TLabel;

Edit5: TEdit;

Label5: TLabel;

procedure Button1Click (Sender: TObject);

procedure FormActivate (Sender: TObject);

procedure Button2Click (Sender: TObject);

procedure Button3Click (Sender: TObject);

private

{ Private declarations }

public

{ Public declarations }

end;

var

Form1: TForm1;

steps, cs: integer;

dT, ynom, time: real;

implementation

{$R *.dfm}

procedure TForm1. Button1Click (Sender: TObject);

begin

Close;

end;

procedure TForm1. FormActivate (Sender: TObject);

begin

Form1.StringGrid1.Cells[0,0]: ='l';

Form1.StringGrid1.Cells[1,0]:='t';

Form1.StringGrid1.Cells[2,0]:='y (l)';

Form1.StringGrid1.Cells[3,0]:='-0,99y (l-2)';

Form1.StringGrid1.Cells[4,0]:='1,98y (l-1)';

Form1.StringGrid1.Cells[5,0]:='0,00086V (l-5)';

Form1.StringGrid1.Cells[6,0]:='V (l)';

Form1.StringGrid1.Cells[7,0]:='y (l)результирующее';

dT:=0.1;

end;

procedure TForm1. Button2Click (Sender: TObject);

begin

//Переходная хар-ка

Form1.Series1.Clear;

Form1.StringGrid1.RowCount:=3;

cs:=1;

steps:=StrToInt (Form1.Edit4.Text);

Form1.StringGrid1.Cells[0,cs]: ='0';

Form1.StringGrid1.Cells[1,cs]:=FloatToStr (dT*StrToFloat (Form1.StringGrid1.Cells[0,cs]));

Form1.StringGrid1.Cells[2,cs]:=Form1.Edit1.Text;

Form1.StringGrid1.Cells[3,cs]:='0';

Form1.StringGrid1.Cells[4,cs]:='0';

Form1.StringGrid1.Cells[5,cs]:='0';

Form1.StringGrid1.Cells[6,cs]:='0';

Form1.StringGrid1.Cells[7,cs]:=Form1.Edit2.Text;

ynom:=StrToFloat (Form1.Edit2.Text);

Form1.Series1.AddXY (0,ynom);

cs:=cs+1;

while (cs<=steps) do

begin

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

Form1.StringGrid1.RowCount:=cs+1;

//Увеличение l

Form1.StringGrid1.Cells[0,cs]: =IntToStr (StrToInt (Form1.StringGrid1.Cells[0,cs-1])+1);

//Расчет t

Form1.StringGrid1.Cells[1,cs]: =FloatToStr (dT*StrToFloat (Form1.StringGrid1.Cells[0,cs]));

//Расчет -0,99y (l-2)

if cs<3 then

Form1.StringGrid1.Cells[3,cs]: =FloatToStr (-0.99*StrToFloat (Form1.StringGrid1.Cells[2,1]))

else

Form1.StringGrid1.Cells[3,cs]:=FloatToStr (-0.99*StrToFloat (Form1.StringGrid1.Cells[2,cs-2]));

//Расчет 1,98y (l-1)

Form1.StringGrid1.Cells[4,cs]: =FloatToStr (1.98*StrToFloat (Form1.StringGrid1.Cells[2,cs-1]));

//Расчет 0,00086V (l-5)

if cs<6 then

Form1.StringGrid1.Cells[5,cs]: =FloatToStr (0.86*StrToFloat (Form1.StringGrid1.Cells[6,1]))

else

Form1.StringGrid1.Cells[5,cs]:=FloatToStr (0.86*StrToFloat (Form1.StringGrid1.Cells[6,cs-5]));

//Расчет V (l)

Form1.StringGrid1.Cells[6,cs]: =Form1.Edit3.Text;

//Расчет y (l)

Form1.StringGrid1.Cells[2,cs]: =FloatToStr (StrToFloat (Form1.StringGrid1.Cells[3,cs])+StrToFloat (Form1.StringGrid1.Cells[4,cs])+

StrToFloat (Form1.StringGrid1.Cells[5,cs]));

//Расчет y (l)результирующее

Form1.StringGrid1.Cells[7,cs]: =FloatToStr (StrToFloat (Form1.StringGrid1.Cells[2,cs])+ynom);

//Внесение данных на график

Form1.Series1.AddXY (StrToFloat (Form1.StringGrid1.Cells[1,cs]), StrToFloat (Form1.StringGrid1.Cells[7,cs]));

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

cs:=cs+1;

end;

procedure TForm1. Button3Click (Sender: TObject);

begin

//Импульсная хар-ка

time:=StrToFloat (Form1.Edit5.Text);

Form1.Series1.Clear;

Form1.StringGrid1.RowCount:=3;

cs:=1;

steps:=StrToInt (Form1.Edit4.Text);

Form1.StringGrid1.Cells[0,cs]: ='0';

Form1.StringGrid1.Cells[1,cs]:=FloatToStr (dT*StrToFloat (Form1.StringGrid1.Cells[0,cs]));

Form1.StringGrid1.Cells[2,cs]:=Form1.Edit1.Text;

Form1.StringGrid1.Cells[3,cs]:='0';

Form1.StringGrid1.Cells[4,cs]:='0';

Form1.StringGrid1.Cells[5,cs]:='0';

Form1.StringGrid1.Cells[6,cs]:='0';

Form1.StringGrid1.Cells[7,cs]:=Form1.Edit2.Text;

ynom:=StrToFloat (Form1.Edit2.Text);

Form1.Series1.AddXY (0,ynom);

cs:=cs+1;

while (cs<=steps) do

begin

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

Form1.StringGrid1.RowCount:=cs+1;

//Увеличение l

Form1.StringGrid1.Cells[0,cs]: =IntToStr (StrToInt (Form1.StringGrid1.Cells[0,cs-1])+1);

//Расчет t

Form1.StringGrid1.Cells[1,cs]: =FloatToStr (dT*StrToFloat (Form1.StringGrid1.Cells[0,cs]));

//Расчет -0,99y (l-2)

if cs<3 then

Form1.StringGrid1.Cells[3,cs]: =FloatToStr (-0.99*StrToFloat (Form1.StringGrid1.Cells[2,1]))

else

Form1.StringGrid1.Cells[3,cs]:=FloatToStr (-0.99*StrToFloat (Form1.StringGrid1.Cells[2,cs-2]));

//Расчет 1,98y (l-1)

Form1.StringGrid1.Cells[4,cs]: =FloatToStr (1.98*StrToFloat (Form1.StringGrid1.Cells[2,cs-1]));

//Расчет 0,00086V (l-5)

if cs<6 then

Form1.StringGrid1.Cells[5,cs]: =FloatToStr (0.86*StrToFloat (Form1.StringGrid1.Cells[6,1]))

else

Form1.StringGrid1.Cells[5,cs]:=FloatToStr (0.86*StrToFloat (Form1.StringGrid1.Cells[6,cs-5]));

//Расчет V (l)

if StrToFloat (Form1.StringGrid1.Cells[1,cs])<=time then

Form1.StringGrid1.Cells[6,cs]: =Form1.Edit3.Text else

Form1.StringGrid1.Cells[6,cs]: =FloatToStr (0);

//Расчет y (l)

Form1.StringGrid1.Cells[2,cs]: =FloatToStr (StrToFloat (Form1.StringGrid1.Cells[3,cs])+StrToFloat (Form1.StringGrid1.Cells[4,cs])+

StrToFloat (Form1.StringGrid1.Cells[5,cs]));

//Расчет y (l)результирующее

Form1.StringGrid1.Cells[7,cs]: =FloatToStr (StrToFloat (Form1.StringGrid1.Cells[2,cs])+ynom);

//Внесение данных на график

Form1.Series1.AddXY (StrToFloat (Form1.StringGrid1.Cells[1,cs]), StrToFloat (Form1.StringGrid1.Cells[7,cs]));

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

cs:=cs+1;

end;

end.

Показать весь текст

Заполнить форму текущей работой