Розробка лінійних систем автоматичного управління (САУ) з цифровою корегуючою ланкою

Вступ

Мета курсової роботи - набуття студентами навичок у розрахунку лінійних систем автоматичного управління (САУ) з цифровою корегуючою ланкою, роль якої може виконувати мікропроцесор, що управляє обчислювальною машиною, чи будь-який спеціалізований цифровий управляючий пристрій.

Відповідно до завдання необхідно розробити слідкуючу систему, що задовольняє визначеним технічним умовам. Система повинна забезпечувати синхронне і синфазне обертання двох вісей, механічно не зв'язаних між собою. Входом системи є кут повороту сельсина-датчика, а виходом - кут повороту вихідного вала редуктора, механічно зв'язаного з робочим механізмом і з ротором сельсина-приймача.

Слідкуючі системи, розглянутого типу, широко застосовуються для дистанційного управління різними механізмами, а також при побудові САУ в різних галузях промисловості.

1. Завдання на виконання курсової роботи

1.1. Мета курсової роботи - проектування слідкуючої системи, яка задовольняє заданим технічним умовам.

1.2. Вимірювальний пристрій - сельсинні пари.

1.3. Виконавчий двигун - двигун постійного струму серії МІ.

1.4. Підсилювач потужності - електромашинний підсилювач з поперечним полем.

1.5. Вихідні дані для проектування системи.

1.5.1. Статичний момент навантаження об'єкта управління - Мос, Нм

1.5.2. Момент інерції об'єкта управління - Jо, кгм2

1.5.3. Максимальна кутова швидкість об'єкта управління - о max.

1.5.4. Максимальне кутове прискорення об'єкта управління - о max

1.5.5. Вимоги, які запропоновані до якості процесу управління: максимальне перерегулювання -  мах,; час регулювання - tр ; максимальна кінетична похибка - xmax, рад.

Варіанти вихідних даних приведені в додатоку 1.

1.6. Курсова робота повинна містити такі розділи.

1.6.1. Розробка функціональної схеми.

1.6.2. Вибір елементів системи: виконавчого двигуна (серії МІ) і електромашинного підсилювача потужності (ЕМП), розрахунок пере1.6.3. Складання передаточних функцій елементів незкорегованої слідкуючої системи.

1.6.4. Побудова логарифмічних частотних характеристик (ЛАЧХ ) незкорегованої системи, бажаної системи і послідовної корегуючої ланки.

1.6.5. Побудова на ЕОМ перехідної функції H(t) і визначення по ній показників якості перехідного процесу для системи з безперервною послідовною корегуючою ланкою.

1.7. Розрахунково-пояснювальна записка повинна містити такі розділи.

1.7.1. Вступ (мета виконання роботи, опис слідкуючої системи, принцип її роботи, опис функціональної схеми системи).

1.7.2. Вихідні дані для проектування системи.

1.7.3. Розрахункова частина.

1.7.4. Висновок (основні характеристики зпроектованої системи).

1.7.5. Список літератури.

1.8. Додатком до розрахунково-пояснювальної записки повинні бути функціональна і структурна електричні схеми слідкуючої системи; ЛАЧХ і ЛФЧХ зкорегованої і незкорегованої системи, графік перехідній функції системи , таблиці і графіки, необхідні для виконання даної роботи.

1.9. Записка повинна бути оформлена відповідно до вимог ЄСКД (написана чорним чорнилом чи пастою, або надрукована на принтері на аркушах формату 11). Орієнтовний обсяг записки - 20 - 25 сторінок.

2. Порядок розрахунку слідкуючої системи

2.1 Розробка функціональної схеми

В слідкуючій системі, яка проектується, в якості виконавчого двигуна (Д) використовується двигун постійного струму серії МІ, в якості підсилювача потужності - електромашинний підсилювач з поперечним полем (ЕМП). Для вимірювального пристрою (ВП) рекомендується використовувати сельсинну пару: сельсин-датчик і сельсин-трансформатор (приймач). Оскільки вимірювальний пристрій працює на змінному струмі, а підсилювач потужності і виконавчий двигун - на постійному струмі, то після вимірювального пристрою повинний бути застосований фазовий детектор (ФД). Крім зазначених елементів у функціональну схему входять: корегувальний пристрій (КП), підсилювач напруги (П), редуктор (Р), за допомогою якого виконавчий двигун з'єднується з об'єктом управління і ротором сельсина-трансформатора, та об'єкт управління (ОУ).

Ф ункціональна схема цифрової слідкуючої системи, наведена на рис.1. Рис.1. Функціональна схема цифрової слідкуючої системи

Рис.1. Функціональна схема слідкуючої системи

слідкуюча цифрова корегуюча ланка

2.2 Вибір виконавчого двигуна

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

Необхідна потужність, Вт:

,

де >р>ККД редуктора, >р = >0,72 .

З каталогу (додаток 2) вибираємо найближчий двигун більшої потужності Р>н > >Р>тр >і виписуємо його паспортні дані:

Р>н> - номінальна потужність (Вт);

n>н >- номінальна швидкість обертання (об/хв);

U>н - >номінальна напруга (В);

I>н> - номінальний струм якоря (А);

R>д>- опір ланцюга обмотки якоря (Ом);

J>д> - момент інерції якоря (кг.м2);

>д>- ККД двигуна.

Потім послідовно визначаємо такі величини:

номінальна кутова швидкість двигуна >н >(с-1) -

>н> = n>н>/30 ;

номінальний момент двигуна М>н> (Н.м) -

М>н> = 9,55Р>н>/n>н ;>

оптимальне передаточне число редуктора i>р >-

,

де J>р> = 1.10-4 кг.м2 - момент інерції редуктора.

Визначаємо необхідний момент на валу двигуна:

.

Для обраного двигуна необхідно перевірити вимогу чи задовольняє він моменту і швидкості у відповідності до таких умов:

М>тр>/М>н>  >o max >i>р / >>н > a ;

де коефіцієнт припустимого перевантаження двигуна по моменту (для двигуна постійного струму a - коефіцієнт припустимого короткочасного збільшення швидкості двигуна понад номінальний, за звичай a = 1,2 - 1,5.

Приклад. Вибрати виконавчий двигун слідкуючої системи, якщо відповідно до технічного завдання: Jо=100кг.м2; М>ос>= 120 Нм; >омах>= 0,7 з-1; >omax>= 0,44 з-2 ; >р> = 0,72.

Необхідна потужність двигуна

Р>тр> = 2(120 + 100. 0,44). 0,7/0,72 = 319 Вт .

Вибираємо двигун типу МІ-22 (додаток2) з такими параметрами: Р>н> = 370Вт; U>н> = 110 В; n>н> = 3000 об/хв; J>д> = 0,004 кг.м2.

Задамо момент інерції редуктора, який підведений до вала двигуна J>р>=1.10-4 кг.м2, і знайдемо передаточне число редуктора

.

Необхідний момент

.

Номінальний момент обраного двигуна дорівнює:

.

Проведемо перевірку двигуна: по моменту - 1,28/1,18< 10 ;

по швидкості - 0,7.355/314 = 0,79 < 1,5 .

У результаті перевірок двигуна по моменту і швидкості видно, що він не перевантажений. Отже, двигун МІ-22 обраний правильно.

2.3 Вибір підсилювача потужності

В якості підсилювача потужності використовуємо ЕМП з поперечним полем. При виборі підсилювача необхідно дотримуватися таких умов.

1. Номінальна потужність підсилювача повинна задовольняти нерівності

Рпн  Р>н> / >д> ,

де >д> - ККД двигуна.

2. Номінальна напруга підсилювача повинна бути не менше номінальної напруги виконавчого двигуна.

3. Номінальний струм підсилювача повинний бути не менше, ніж номінальний струм двигуна.

З огляду на зазначені умови, вибираємо тип ЕМП (див. додаток 3).

2.4 Складання передаточних функцій елементів слідкуючої системи

2.4.1 Виконавчий двигун

Передаточна функція виконавчого двигуна по куту повороту має вид (якщо знехтувати індуктивністю ланцюга якоря)

де К>д> - коефіцієнт підсилення двигуна, рад/В.с : К>д> = >н / >U>н ;>

Т>д - >електромеханічна постійна часу:

В останній формулі =1,2 - постійний коефіцієнт;

J>c> - сумарний момент інерції, приведений до вала двигуна

;

2.4.2 Електромашинний підсилювач (ЕМП)

Передаточна функція ЕМП

К>ЕМП> - коефіцієнт підсилення ЕМП по напрузі

де U>ЕМП> - напруга на виході ЕМП ;

- напруга обмотки управління ЕМП ;

Р>у>, R>у> - відповідно потужність і опір обмотки управління ЕМП,

Ту, Т>кз> - постійні часу обмотки управління і короткозамкненої обмотки

якоря ЕМП.

2.4.3 Підсилювач

Підсилювач (на функціональній схемі - П) слугує для узгодження вихідного сигналу ЦАП із вхідним опором обмотки управління ЕМП. Його можна вважати безінерційною ланкою з передаточною функцією W>п>(Р) = К>п>.

У розрахунках прийняти К>п>=> >1> .>

2.4.4 Фазовий детектор

Передаточна функція фазового детектора W>фд >(P) = К>фд>,

де К>фд>=1 - коефіцієнт підсилення фазового детектора.

2.4.5 Вимірювальний пристрій

Передаточна функція вимірювального пристрою W>вп>(Р)=К>вп>> >,

де К>вп>=1 - коефіцієнт підсилення вимірювального пристрою.

2.4.6 Редуктор

Передаточна функція редуктора W>ред>(Р)=К>ред>=1/i>р> .

С
труктурна схема незкорегованої слідкуючої системи, наведена на рис.2.

Рис.2. Структурна схема незкорегованої слідкуючої системи

2.5 Розрахунок послідовної безперервної корегуючої ланки методом ЛАЧХ

Побудова логарифмічної амплітудно-частотної характеристики (ЛАЧХ) послідовної корегуючої ланки проводиться в такій послідовності.

1. Будується ЛАЧХ заданої (незкорегованої) системи.

2. Будується бажана ЛАЧХ за заданими показниками якості перехідного процесу.

3. Будується ЛАЧХ послідовної корегуючої ланки шляхом графічного вирахування ЛАЧХ заданої системи з бажаної ЛАЧХ.

4. По виду ЛАЧХ корегуючої ланки визначається його передаточна функція (безперервна).

2.5.1 Побудова ЛАЧХ заданої системи по виду передаточної функції

Передаточну функцію розімкненої системи потрібно представити у виді добутку передаточних функцій типових динамічних ланок (обмежимося випадком, коли в системі відсутні коливальні ланки і ланки із запізненням). Наприклад, нехай передаточна функція розімкненої системи має вид:

.

Побудову зручно проводити в такій послідовності.

1. Визначити частоти, при яких відбувається зміна нахилу ЛАЧХ і які розраховуються як >i>=1/Т>i>>,> та відкласти їх по вісі абсцис у логарифмічному масштабі (Т>i> - постійні часу передаточної функції К>з>(Р) ).

2. Відкласти точку A1 з координатами >A1>=1c-1 і L(>A1>)=20lgК>з> ( див. рис. 3). Через точку A1 провести пряму з нахилом -20 дБ/дек. Побудована в такий спосіб пряма лінія збігається з ЛАЧХ при частотах, менших першої зупрягаючої частоти (один по одному їхнього розташування на вісі частот ліворуч-праворуч).

3. На частоті сполучення >i >характеристика змінює свій нахил: або на +20 дБ/дек, якщо постійна часу Т>i=>1/>i >знаходиться в чисельнику вихідної передаточної функції, або на -20 дБ/дек, якщо постійна часу Т>i >знаходиться в знаменнику передаточної функції.

Для дослідження системи на стійкість за амплітудно-фазовим критерієм стійкості за допомогою логарифмічних частотних характеристик, необхідно крім ЛАЧХ побудувати ще логарифмічну фазочастотну характеристику (ЛФЧХ). По вісі абсцис відкладається частота в логарифмічному масштабі (використовують ту ж вісь частот, що і для побудови ЛАЧХ), а по вісі ординат відкладають аргумент амплітудно-фазової характеристики чи у градусах чи у радіанах в лінійному масштабі. Для розглянутого приклада  розраховується по формулі

arctgT>у>arctgT>кз>arctgT>д>[град.].

2.5.2 Побудова бажаної ЛАЧХ

При побудові бажаної ЛАЧХ виділяють три області: область низьких частот, область середніх частот і область високих частот. Вид ЛАЧХ у кожній з областей по різному впливає на якість системи. В області низьких частот вид ЛАЧХ визначає точність роботи системи в сталих режимах. Область середніх частот визначає динамічні властивості системи (швидкодія, коливальність). Вид ЛАЧХ в області високих частот практично мало впливає на якість системи.

Побудову бажаної ЛАЧХ зручно починати з області середніх частот у такій послідовності.

По заданих величинах >max > і t>р> визначаємо за допомогою табл. 1 частоту зрізу >с.>

Наприклад, нехай задано, що >max>=30  t>р>=1,5 с. З табл.1 для >max >=30% визначаємо t>р>>с>/=3,2, звідки випливає:

>с>>=> 3,26c-1 .

Таблиця 1

>max >, %

10

15

20

25

30

35

40

t>р>>с>/

5

4,4

4

3,6

3,2

3

2,8

L>1>, дБ

18

15

13,5

12

11

10,5

10

, гр

85

80

65

55

45

40

35

Р
ис.3. а-а-а - ЛАЧХ заданої (незмінної частини) системи;

б-б-б - ЛАЧХ бажана;

c-c-c - ЛАЧХ послідовної корегуючої ланки.

Наносимо на вісь абсцис частоту зрізу >с >і проводимо через неї пряму лінію з нахилом -20 дБ/дек (рис.3).

Частота >2>, яка обмежує область середніх частот бажаної ЛАЧХ ліворуч, визначається величиною відрізка L>1>, що може бути знайдена в залежності від заданої величини >max> з табл.1. Частота >3>, яка обмежує область середніх частот праворуч, визначається величиною відрізка L>2>, при цьому .

В області високих частот бажану ЛАЧХ потрібно будувати у виді прямолінійних відрізків з нахилом, кратним -20 дБ/дек. (тобто -40, - 60, -80 і т.д.), таким чином, щоб різниця характеристик бажаної і заданої (у межі при складала пряму лінію, рівнобіжну вісі частот.

В області низьких частот бажана ЛАЧХ будується в такий спосіб. По заданій величині коефіцієнта підсилення системи К>с>=>о>>m>>ах>/x>max >визначаємо величину L>A2>=20LgK>c> і відзначаємо на кресленні точку A>2> c координатами >A2>=1 c-1 і L>A2> (див. рис.3). Через точку A>2> проводимо пряму лінію з нахилом -20 дБ/дек.

Від точки М, що обмежує область середніх частот ліворуч, проводимо пряму лінію з нахилом -40дБ/дек до перетинання з низькочастотною частиною бажаної ЛАЧХ.

Якщо в завданні на розробку слідкуючої системи, зазначена максимально припустима похибка спостереження Хмах за умови, що вхідний сигнал може змінюватися з максимальною кутовою і з максимальним кутовим прискоренням max , то для виконання цих вимог необхідно, щоб бажана ЛАЧХ не попадала у заборонену область.

Заборонена область будується в такий спосіб. Визначаємо на кресленні точку В з координатами:

.

Від точки B вправо проводимо пряму лінію з нахилом -40 дБ/дек, а вліво - пряму лінію з нахилом -20 дБ/дек.

Якщо ЛАЧХ, побудована по заданому коефіцієнту Кс, попадає в заборонену область, то це означає, що при даному коефіцієнті Кс задана точність спостереження не може бути забезпечена і коефіцієнт Кс потрібно збільшити, (бажану ЛАЧХ зрушувати вправо), щоб бажана ЛАЧХ не попадала в заборонену область.

По виду ЛАЧХ бажаної можна записати передаточну функцію безперервної зкорегованої (бажаної) системи. Для розглянутого приклада (крива б-б-б... рис.3) передаточна функція має вид:

. (1)

Для визначення передаточної функції бажаної системи можна скористатися програмою, що наведена в додатку 4. Програма написана мовою BASIC і дозволяє знайти постійні часу T1, Т2, Т3 бажаної передаточної функції по показниках якості mах і tp. При цьому передаточна функція записується у виді

. (2)

Типова бажана ЛАЧХ, по якій записана передаточна функція (2), показана на рис. 4.

Показник ступеня k визначається нахилом заданої ЛАЧХ в області високих частот.

Рис 4. Типова ЛАЧХ бажаної слідкуючої системи

Для забезпечення заданих показників якості перехідного процесу скорегована система повинна мати визначений запас стійкості по фазі. Необхідна величина запасу сталості по фазі (для заданої величини максимального перерегулювання зазначена в табл.1. Після побудови бажаної ЛАЧХ потрібно розрахувати і побудувати ЛФЧХ зкорегованої системи і визначити 

2.5.3 Розрахунок послідовної корегуючої ланки

ЛАЧХ послідовного безперервної корегуючої ланки будується шляхом графічного вирахування з ЛАЧХ бажаної ЛАЧХ заданої частини системи ( на рис.3 - лінія с-с-с).

По виду ЛАЧХ потрібно записати передаточну функцію безперервної послідовної корегуючої ланки. Для ЛАЧХ с-с-с на рис.3 передаточна функція має вид:-

. (2)

Коефіцієнт К>кор> визначається зі співвідношення: 20lgК>кор>=L>3> .

Цю ж передаточну функцію можна одержати, якщо бажану передаточну функцію К>Б>(Р) поділити на передаточну функцію заданої частини системи К>з>(Р).

2.6 Моделювання слідкуючої системи з безперервною послідовною корегуючою ланкою

Для того щоб переконатися, що корекція системи проведена правильно і скорегована система має показники якості перехідного процесу не гірше заданих, потрібно провести моделювання. Рекомендується моделювати скореговану систему на ПЕОМ, використовуючи пакет Matlab simulink.

Якщо виявиться, що показники якості скорегованої системи гірше заданих, то потрібно вносити в корекцію відповідні зміни. Так, якщо перерегулювання виявилося більше заданого, то потрібно збільшувати відрізки L1 і L2. Якщо час регулювання виявилося більше заданого, потрібно збільшувати частоту зрізу бажаної ЛАЧХ.

Література

1. Бесекерский В.А., Попов Е.П. Теория систем автоматического регулирования. М.: Наука, 2001. 766 с.

2. Батоврин А.А., Дашевский П.Г. и др. Цифровые следящие системы судовой автоматики. Л.: Судостроение, 1999. 445 с.

Додаток 1

Варіанти вихідних даних для проектування слідкуючої системи

Вариант

Статичний момент навантаження

Мос, Нм

Момент інерції об'єкта управління

Jо, кгм2

Максим. кутова швидкість

о max, с-1

Максим. кутове прискорення

о max, с-2

Максим. кинетич. похибка

Xmax, рад

Максим. пере -регулювання

max, %

1

2

3

4

5

6

7

1

50

20

0,5

0,02

0,01

20

2

100

70

1,0

0,06

0,02

25

3

200

100

1,5

0,06

0,01

25

4

300

90

2,0

0,08

0,02

30

5

350

120

2,5

0,08

0,02

30

6

50

30

3,0

0,05

0,03

35

7

100

60

0,5

0,02

0,01

20

8

150

90

1,0

0,06

0,02

20

9

200

80

1,5

0,05

0,02

25

10

250

120

2,0

0,08

0,03

25

11

300

100

0,5

0,02

0,01

30

12

350

100

1,0

0,04

0,02

30

13

50

50

0,5

0,05

0,01

35

14

100

100

0,2

0,08

0,02

35

15

150

120

0,5

0,02

0,05

35

16

200

150

1,0

0,03

0,02

20

17

250

180

1,0

0,04

0,01

25

18

300

150

1,5

0,02

0,02

25

19

350

150

0,5

0,02

0,05

30

20

50

70

2,5

0,05

0,01

30

21

100

50

3,0

0,05

0,01

20

1

2

3

4

5

6

7

22

150

100

2,5

0,08

0,05

25

23

200

120

1,5

0,05

0,02

30

24

70

25

0,75

0,025

0,01

20

25

120

50

0,8

0,01

0,02

20

26

110

60

0,5

0,04

0,01

25

27

300

120

0,5

0,01

0,02

25

28

150

100

1,0

0,02

0,05

30

29

100

70

1.2

0,05

0,02

20

30

75

50

1,0

0,01

0,05

25

Учбові групи

Час регулювання, tp, с

АТ - хх1

2

АТ - хх2

2,5

АТ - хх3

3

Додаток 2

Технічні дані двигунів серії МІ

Тип

двигуна

Потужність

Швидкість обертання

Напруга живлення

Струм

якоря

Опір

ланцюга

якоря

ККД

Момент

інерц.

Рн,

кВт

nн,

об/хв

Uн,

В

Iн,

А

Rд,

Ом

д,

%

Jд,

кгм2

1

2

3

4

5

6

7

8

МІ-11

0,12

3000

60

2,87

0,46

62

0,0015

0,1

2000

60

2,27

0,94

63

0,0015

0,12

3000

100

1,53

1,48

62

0,0015

0,1

2000

110

1,22

3,60

63

0,0015

МІ-12

0,2

3000

60

4,57

0,23

66

0,002

0,12

2000

60

2,72

0,52

64

0,002

0,2

3000

110

2,46

0,765

66

0,002

0,12

2000

110

1,46

1,74

64

0,002

МІ-21

0,25

3000

60

5,6

0,284

67

0,0035

0,2

2000

60

4,3

0,645

68

0,0035

0,25

3000

110

3,05

0,945

67

0,0035

0,2

2000

110

2,33

2,20

68

0,0035

МІ-22

0,37

3000

60

8,2

0,195

71

0,04

0,25

2000

60

5,5

0,360

75

0,004

0,12

1000

60

2,6

1,44

64

0,004

0,37

3000

110

4,4

0,546

72

0,004

0,25

2000

110

2,9

1,29

70

0,004

0,12

1000

110

1,4

4,58

64

0,004

МІ-31

0,45

3000

60

10,3

0,204

68

0,009

0,37

2000

60

8,2

0,405

70

0,009

0,45

3000

110

5,6

0,585

68

0,009

0,2

1000

60

4,4

1,32

66

0,009

0,37

2000

110

4,4

1,16

70

0,009

0,2

1000

110

2,4

3,9

66

0,009

МІ-32

0,76

2500

110

8,2

0,368

80

0,0132

0,45

1500

110

5,0

0,975

75

0,0132

0,37

1000

110

4,2

2,21

73

0,0132

0,76

2500

220

4,1

1,36

80

0,0132

0,37

1000

220

2,1

8,37

73

0,0132

1,6

2500

110

19,5

0,249

73

0,035

1,1

1500

110

13,0

0,67

74

0,035

0,76

1000

110

9,0

1,3

72

0,035

1,6

2500

220

9,5

0,93

73

0,035

1,1

1500

220

6,4

2,63

75

0,035

0,76

1000

220

4,5

5,32

72

0,035

МІ-42

3,2

2500

110

36,3

0,1

78

0,065

1,6

1500

110

18,2

0,32

78

0,065

1,1

1000

110

12,6

0,75

75

0,065

3,2

2500

220

18,0

0,376

79

0,065

1,6

1500

220

9,1

1,28

78

0,065

1,1

1600

220

6,3

2,95

75

0,065

МІ-51

5,5

2500

220

27,2

0,164

82

0,125

3,2

1500

220

17,1

0,46

82

0,125

1,6

1000

220

8,7

1,1

79

0,125

МІ-52

7,0

2500

220

37,0

0,088

84

0,15

4,5

1500

220

23,3

0,26

85

0,15

2,5

1000

220

13,1

0,569

82

0,15

Додаток 3

Технічні дані ЕМП

Тип ЕМП

Потужність ЕМП

Потужність

управління

Напруга

Струм якоря

Опір

обмотки

управління

Постійні часу

Ту

Ткз

кВт

Вт

В

А

Ом

с

с

ЕМП-3А3

0,2

0,4

115

1,75

1000

0,005

0,018

ЕМП-5А3

0,5

0,4

115

4,35

1000

0,01

0,033

ЕМП-12А3

1,0

0,4

115

8,7

2200

0,015

0,06

ЕМП-25А3

2,0

0,4

230

9,1

1500

0,02

0,1

ЕМП-50А3

4,0

0,5

230

17,4

2200

0,03

0,17

ЕМП-70А3

6,0

0,5

230

26

1500

0,04

0,22

ЕМП-100А3

8,5

0,5

230

37

1000

0,06

0,28

Програма LOGHAR.

Визначення постійних часу передаточної функції бажаної системи

10 PRINT "Визначення пост. часу передат. функції бажаної системи"

20 PRINT " Перед. ф-ція, яка відшукується, має вид:"

30 PRINT " Kc(T(2)*P+1)/((T(1)*P+1)*(T(3)*P+1)^(n- m)*(T(4)P+1)*..(T(m+3)P+1)*P)"

40 PRINT "n-ступінь полінома Q(P) - знаменника перед. ф-ції заданої сист."

50 PRINT "m-кількість пост. часу Q(P), менших, за Т(3)"

60 PRINT "Q(P)=(Tз(1)*P+1)*(Tз(2)*P+1)*...*(Tз(N)*P+1)"

70 PRINT "Введіть порядок полінома знаменника Q(P) заданої перед.ф-ції"

80 PRINT "N="

90 INPUT N

100 DIM T3(5), T(8)

110 FOR I = 1 TO 8

120 T(I) = 0

130 NEXT I

140 PRINT "Введіть пост. часу знаменника заданої перед. ф-ції Q(P)"

150 FOR I = 1 TO N

160 PRINT "Tз("; I; ")="

170 INPUT T3(I)

180 NEXT I

190 PRINT "Введіть величину макс. перерегулювання"

200 PRINT "Сигма макс. ,% ="

210 INPUT SM

220 IF SM = 10 THEN C = 5: L1 = 18

230 IF SM = 15 THEN C = 4.4: L1 = 15

240 IF SM = 20 THEN C = 4: L1 = 13.5

250 IF SM = 25 THEN C = 3.6: L1 = 12

260 IF SM = 30 THEN C = 3.2: L1 = 11

270 IF SM = 35 THEN C = 3: L1 = 10.5

280 IF SM = 40 THEN C = 2.8: L1 = 10

290 PRINT "Порядок астатизму NU="

300 INPUT NU

310 PRINT "Коефіцієнт підсилення бажаної системи Кс="

320 INPUT KC

330 PRINT "Час регулювання TR="

340 INPUT TR

350 M1 = .434

360 OC = C * 3.14 / TR

370 XC = LOG(OC) * M1

380 B = 20 * XC

390 X2 = (B - L1) / 20

400 T(2) = 1 / (10 ^ X2)

410 A = 20 * LOG(KC) * M1

420 X1 = (L1 + 40 * X2 - A) / 20

430 T(1) = 1 / (10 ^ X1)

440 X3 = (L1 + 20 * XC) / 20

450 T(3) = 1 / 10 ^ X3

460 IF NU = 1 THEN GOTO 490

470 X1 = (40 * X2 + L1 - A) / 40

480 T(1) = 1 / (10 ^ X1)

490 M = 0

510 FOR I = 1 TO N

520 IF T3(I) <= T(3) * .75 THEN M = M + 1: T(3 + M) = T3(I)

540 NEXT I

550 S = N - M

560 IF S = 2 THEN T(3) = .8 * T(3)

570 IF S > 2 THEN T(3) = .7 * T(3)

580 X3 = LOG(1 / T(3))

590 PRINT "T(1)="; T(1); "T(2)="; T(2); "T(3)="; T(3)

600 PRINT "Вираз у дужках знаменника (Т(3)*Р+1) повинен зводитися в ступінь (n-m)"

610 PRINT "n-m="; S

620 FOR I = 1 TO M

630 K = 3 + M

640 PRINT "T("; K; ")="; T(K)

650 NEXT I

660 PRINT "Wc="; OC; "LgWc="; XC

670 PRINT "Lg(1/T1)="; X1; "Lg(1/T2)="; X2; "Lg(1/T3)="; X3

680 END