Курсовая работа: Анализ и синтез системы автоматического регулирования

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Коэффициенты характеристического полинома замкнутой скорректированной системы имеют следующие значения: d0 = 0,002625, d1 = 0,089175, d2 = 0,91025, d3 = 3,01497, d4 = 1,5499.

Матрица Гурвица имеет вид:

Диагональные миноры матрицы Гурвица:

Д1 = 0,089175 > 0;

;

.

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

4.3 Исследование скорректированной АСР по критерию Михайлова

В пункте 4.2 был получен характеристический полином такой замкнутой скорректированной системы DЗАМ(p):

Тогда получим частотный годограф Михайлова путем перевода характеристического полинома (4.15) в частотную область:

Получим действительную и мнимую часть частотного годографа Михайлова (4.16), возведя частотный оператор j в соответствующую степень:

Вычисляем значения действительной и мнимой части частотного годографа Михайлова при изменении часты от 0 до требуемого значения, при котором можно сделать вывод об устойчивости системы. Шаг изменения частоты принимаем произвольный с учетом удобства восприятия графика. Вычисления проводим в MS Excel.

Таблица 4.1 Значения действительной и мнимой части частотного годографа Михайлова

Частота сигнала

Действительная часть UD()

Мнимая часть VD()

0

1,55

0

2

-2,049

5,31654

4

-12,342

6,35268

6

-27,817

-1,17198

8

-45,954

-21,5378

10

-63,225

-59,0253

12

-75,094

-117,915

14

-76,017

-202,487

16

-59,442

-317,021

18

-17,809

-465,799

20

57,45

-653,101

По таблице 4.1 строим годограф Михайлова (рис. 4.1).

Рис. 4.1 Годограф Михайлова скорректированной системы

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

4.4 Исследование скорректированной АСР по критерию Найквиста

Для построения АФЧХ разомкнутой скорректированной системы необходимо перейти к частотной форме записи передаточной функции.

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

Проведем преобразования, выделив общий коэффициент усиления разомкнутой системы kРАЗ:

Передаточная функция разомкнутой системы в частотной форме (p = jщ):

Запишем выражения амплитуд А() и фаз ц() частотных передаточных функции для каждого звена (табл. 4.2). Вычисляем значения амплитуд А() и фаз ц() при изменении часты от 0 до значения, при котором UРАЗ(щ) и VРАЗ(щ) станут равны нулю. Также вычислим амплитуду АРАЗ() и фазу цРАЗ() передаточной функции разомкнутой системы, а также вещественную UРАЗ(щ) и мнимую VРАЗ(щ) составляющую частотной передаточной функции разомкнутой системы WРАЗ(jщ) (3.18 - 3.23). Шаг изменения частоты принимаем произвольный с учетом удобства восприятия графика. Вычисления проводим в MS Excel (табл. 4.3).

Таблица 4.2 Выражения амплитуд А() и фаз ц() передаточных функции звеньев системы

Частотная передаточная функция

Амплитуда А()

Фаза ц()

Таблица 4.3 Значения амплитуд А() и фаз ц() передаточных функции звеньев системы

щ

Амплитуды звеньев

Фазы звеньев

Амплитуда

WРАЗ(jщ)

Фаза

WРАЗ(jщ)

Действит.

часть WРАЗ(jщ)

Мнимая

часть WРАЗ(jщ)

A1(щ)

A2(щ)

A3(щ)

ц1(щ)

ц2(щ)

ц3(щ)

AРАЗ(щ)

цРАЗ(щ)

UРАЗ(щ)

VРАЗ(щ)

0,0

1,000

0,550

1,000

0,000

0,000

0,000

0,550

0,000

0,55000

0,00000

0,1

0,971

0,550

1,000

-0,215

-0,020

-0,015

0,534

-0,250

0,51709

-0,13203

0,2

0,896

0,550

1,000

-0,404

-0,040

-0,030

0,492

-0,474

0,43810

-0,22458

0,3

0,803

0,549

0,999

-0,554

-0,060

-0,045

0,441

-0,659

0,34855

-0,26976

0,4

0,712

0,548

0,999

-0,666

-0,080

-0,060

0,390

-0,806

0,27015

-0,28139

0,5

0,632

0,547

0,999

-0,747

-0,100

-0,075

0,345

-0,922

0,20863

-0,27500

0,6

0,564

0,546

0,998

-0,805

-0,119

-0,090

0,307

-1,014

0,16229

-0,26074

0,7

0,507

0,545

0,997

-0,845

-0,139

-0,105

0,275

-1,089

0,12767

-0,24394

0,8

0,460

0,543

0,996

-0,872

-0,159

-0,120

0,249

-1,150

0,10162

-0,22712

0,9

0,421

0,541

0,995

-0,889

-0,178

-0,135

0,227

-1,202

0,08175

-0,21135

1,0

0,388

0,539

0,994

-0,899

-0,197

-0,150

0,208

-1,246

0,06636

-0,19700

2,0

0,229

0,511

0,977

-0,833

-0,381

-0,297

0,114

-1,511

0,00683

-0,11386

3,0

0,178

0,472

0,949

-0,705

-0,540

-0,441

0,080

-1,687

-0,00924

-0,07905

4,0

0,155

0,429

0,913

-0,595

-0,675

-0,580

0,061

-1,850

-0,01678

-0,05861

5,0

0,144

0,389

0,871

-0,508

-0,785

-0,712

0,049

-2,006

-0,02053

-0,04418

6,0

0,137

0,352

0,825

-0,441

-0,876

-0,837

0,040

-2,154

-0,02192

-0,03323

7,0

0,133

0,320

0,778

-0,387

-0,951

-0,955

0,033

-2,293

-0,02180

-0,02475

8,0

0,130

0,291

0,729

-0,345

-1,012

-1,065

0,028

-2,422

-0,02076

-0,01818

9,0

0,128

0,267

0,682

-0,310

-1,064

-1,168

0,023

-2,542

-0,01922

-0,01313

10,0

0,126

0,246

0,636

-0,282

-1,107

-1,264

0,020

-2,653

-0,01745

-0,00927

12,0

0,124

0,212

0,551

-0,238

-1,176

-1,436

0,014

-2,850

-0,01388

-0,00417

14,0

0,123

0,185

0,476

-0,205

-1,228

-1,585

0,011

-3,018

-0,01078

-0,00134

16,0

0,123

0,164

0,412

-0,180

-1,268

-1,713

0,008

-3,161

-0,00829

0,00016

18,0

0,122

0,147

0,358

-0,161

-1,300

-1,825

0,006

-3,286

-0,00637

0,00092

20,0

0,122

0,133

0,313

-0,145

-1,326

-1,922

0,005

-3,393

-0,00492

0,00126

22,0

0,121

0,122

0,275

-0,132

-1,347

-2,008

0,004

-3,487

-0,00382

0,00138

24,0

0,121

0,112

0,242

-0,121

-1,365

-2,083

0,003

-3,570

-0,00299

0,00137

26,0

0,121

0,104

0,215

-0,112

-1,381

-2,150

0,003

-3,643

-0,00237

0,00130

28,0

0,121

0,097

0,191

-0,104

-1,394

-2,209

0,002

-3,707

-0,00189

0,00120

30,0

0,121

0,090

0,171

-0,097

-1,406

-2,262

0,002

-3,765

-0,00152

0,00109

32,0

0,121

0,085

0,154

-0,091

-1,416

-2,310

0,002

-3,817

-0,00123

0,00099

34,0

0,121

0,080

0,139

-0,086

-1,425

-2,353

0,001

-3,864

-0,00101

0,00089

36,0

0,121

0,076

0,126

-0,081

-1,433

-2,392

0,001

-3,906

-0,00083

0,00080

38,0

0,120

0,072

0,115

-0,077

-1,440

-2,428

0,001

-3,945

-0,00069

0,00071

40,0

0,120

0,068

0,105

-0,073

-1,446

-2,460

0,001

-3,980

-0,00058

0,00064

По значениям вещественной UРАЗ(щ) и мнимой VРАЗ(щ) части частотной передаточной функции разомкнутой скорректированной системы из таблицы 4.3 строим годограф Найквиста (рис. 4.2).

Рис. 4.2. Годограф Найквиста (АФЧХ) скорректированной системы

Вывод: замкнутая скорректированная система устойчива, так как АФЧХ разомкнутой системы не охватывает точку с координатами (-1; j0).

4.5 Исследование скорректированной АСР по логарифмическому критерию

Для построения ЛАХ и ЛФХ системы необходимо разложить передаточную функцию разомкнутой системы на элементарные звенья, амплитуды А() и фазы ц() которых приведены в пункте 4.3.

Затем на плоскости строятся ЛАХ и ЛФХ каждого отдельного звена и методом графического суммирования находятся результирующие ЛАХ и ЛФХ разомкнутой системы [**]:

Зная выражения амплитуд А() и фаз ц() частотных передаточных функции для каждого звена (табл. 4.2), представим их в логарифмическом масштабе (табл. 4.4).

Таблица 4.4 Выражения логарифмических амплитуд L() и фаз ц() передаточных функции звеньев

Частотная передаточная функция

Амплитуда L()

Фаза ц()

Вычисляем значения логарифмических амплитуд L() и фаз ц() для каждого звена при изменении часты от 0 до значения, при котором результирующая ЛФХ пересекает значение - р. Также вычислим амплитуду LРАЗ() и фазу цРАЗ() передаточной функции разомкнутой системы WРАЗ(jщ). Шаг изменения частоты принимаем произвольный с учетом удобства восприятия графика. Вычисления проводим в MS Excel (табл. 4.5).

Таблица 4.5 Значения логарифмических амплитуд L() и фаз ц() передаточных функции звеньев

Частота щ

Логарифмические амплитуды

Общая ЛАХ

LРАЗ(щ)

Фазы звеньев

Общая ЛФХ

цРАЗ(щ)

L1(щ)

L2(щ)

L3(щ)

ц1(щ)

ц2(щ)

ц3(щ)

0,100

-0,259

-5,194

-0,001

-5,454

-0,215

-0,020

-0,015

-0,250

0,500

-3,990

-5,236

-0,013

-9,239

-0,747

-0,100

-0,075

-0,922

0,900

-7,521

-5,331

-0,042

-12,894

-0,889

-0,178

-0,135

-1,202

1,300

-10,016

-5,477

-0,088

-15,580

-0,900

-0,254

-0,194

-1,349

1,700

-11,798

-5,668

-0,149

-17,615

-0,868

-0,328

-0,253

-1,449

2,100

-13,106

-5,898

-0,226

-19,231

-0,820

-0,398

-0,312

-1,530

2,900

-14,841

-6,452

-0,425

-21,718

-0,718

-0,526

-0,427

-1,670

3,300

-15,426

-6,763

-0,546

-22,735

-0,670

-0,583

-0,483

-1,737

4,100

-16,254

-7,426

-0,826

-24,506

-0,585

-0,687

-0,593

-1,866

4,500

-16,551

-7,770

-0,984

-25,304

-0,549

-0,733

-0,647

-1,929

5,300

-16,994

-8,463

-1,332

-26,789

-0,486

-0,815

-0,750

-2,051

5,700

-17,160

-8,809

-1,520

-27,489

-0,459

-0,851

-0,801

-2,110

6,100

-17,300

-9,152

-1,716

-28,169

-0,435

-0,884

-0,849

-2,168

6,900

-17,520

-9,823

-2,132

-29,475

-0,392

-0,944

-0,944

-2,280

7,300

-17,607

-10,150

-2,350

-30,107

-0,374

-0,970

-0,989

-2,333

8,100

-17,748

-10,785

-2,801

-31,334

-0,341

-1,018

-1,076

-2,435

8,500

-17,805

-11,092

-3,034

-31,931

-0,327

-1,039

-1,118

-2,484

9,300

-17,900

-11,686

-3,510

-33,096

-0,301

-1,077

-1,198

-2,576

9,700

-17,939

-11,972

-3,753

-33,664

-0,290

-1,095

-1,236

-2,621

10,100

-17,974

-12,252

-3,998

-34,224

-0,279

-1,111

-1,273

-2,664

12,100

-18,103

-13,554

-5,248

-36,905

-0,236

-1,179

-1,444

-2,859

14,100

-18,184

-14,712

-6,508

-39,404

-0,204

-1,230

-1,591

-3,025

16,100

-18,237

-15,750

-7,754

-41,741

-0,179

-1,270

-1,719

-3,168

18,100

-18,274

-16,686

-8,971

-43,930

-0,160

-1,301

-1,830

-3,291

По таблице 4.5 строим логарифмические частотные характеристики системы (рис. 4.3).

Рис. 4.3 Логарифмические частотные характеристики разомкнутой системы (ЛАХ и ЛФХ)

По частотным характеристикам (рис. 4.3) графически определим запасы устойчивости системы: ДL = **** дБ и Дц = **** рад.

Вывод: замкнутая скорректированная система устойчива, так как ЛАХ разомкнутой системы пересекает ось абсцисс раньше, чем ЛФХ переходит через значение - р.

Раздел 5. Исследование АСР в среде Simulink

Программный комплекс MATLAB поставляется вместе с пакетом расширения Simulink, предназначенным для моделирования динамических систем, модели которых составляются из отдельных блоков (компонентов) [**].

Для исследования скорректированной системы в среде Simulink составим структурную схему модели АСР по управляющему воздействию (рис. 5.1).

Для получения линейных характеристик системы необходимо задать точки входа и выхода для линеаризации: «Input point» и «Output point».

Для получения переходных и частотных характеристик в пакете Simulink с помощью LTI Viewer необходимо выполнить следующие действия:

1. Выполнить команду меню Tools > Control Design > Linear Analysis, после чего появится окно менеджера Control and Estimation Tools Manager, в котором показываются активные точки, кнопка запуска линеаризации «Linearize Model» и окно выбора типа характеристик.

2. Выбрать нужную характеристику и нажать кнопку «Linearize Model», после чего на экран выводится окно LTI Viewer с построенным графиком.

При получении частотных характеристик необходимо поставить галочку «Open Loop» на активной точке выхода «Output point», так как частотные характеристики строятся только для разомкнутых систем.

На рисунке 5.2 представлена переходная характеристика АСР стабилизации угловой скорости ДПТ с корректирующим звеном, полученная в Simulink.

На рисунках 5.3 и 5.4 представлены логарифмические частотные характеристика и АФЧХ АСР угловой скорости ДПТ с корректирующим звеном.

Таким образом, …………..(вывод).

Рис. 3.3 Логарифмические частотные характеристики разомкнутой системы (ЛАХ и ЛФХ)

По частотным характеристикам (рис. 3.3) графически определим запасы устойчивости системы: ДL = **** дБ и Дц = **** рад.

Вывод: замкнутая система устойчива, так как ЛАХ разомкнутой системы пересекает ось абсцисс раньше, чем ЛФХ переходит через значение - р.

3.2 Корректирующее звено:

Передаточная функция корректирующего звена:

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

Где:

Находим такое значение k0, чтобы статическая ошибка соответствовала разрешенной:

3.3 Критерий устойчивости Гурвица:

Составим характеристическое уравнение:

Определитель Гурвица:

Вывод: Система устойчива, так как не имеет отрицательных миноров.

3.4 Критерий устойчивости Михайлова:

?

Ud(?)

j*Vd(?)

?

Ud(?)

j*Vd(?)

0

1

0

11

-12,13026

-40,1687

0,5

0,973845

4,7369125

12

-14,62832

-73,8216

1

0,89334

9,3923

13

-17,34378

-115,3009

2

0,56928

18,1324

14

-20,27664

-165,2588

3

0,02782

25,5681

15

-23,4269

-224,3475

4

-0,73104

31,0472

16

-26,79456

-293,2192

5

-1,7073

33,9175

17

-30,37962

-372,5261

6

-2,90096

33,5268

18

-34,18208

-462,9204

7

-4,31202

29,2229

19

-38,20194

-565,0543

8

-5,94048

20,3536

20

-42,4392

-679,58

9

-7,78634

6,2667

21

-46,89386

-807,1497

10

-9,8496

-13,69

22

-51,56592

-948,4156

Список использованной литературы

1. Бесекерский, В. А. Теория автоматического управления [Текст]: учебное пособие / В.А. Бесекерский, Е.П. Попов. - СПб.: Изд-во «Профессия», 2003. - 751 c.

2. Мирошник, И. В. Теория автоматического управления. Линейные системы [Текст]: учебное пособие / И. В. Мирошник. - СПб.: Изд-во «Питер», 2005. - 336 с.

3. Щербаков, В.С. Теория автоматического управления. Линейные непрерывные системы [Текст]: учебное пособие / В.С. Щербаков, И.В. Лазута. - Омск: СибАДИ, 2013. - 142 с.

4. Щербаков, В.С. Основы моделирования систем автоматического регулирования и электротехнических систем в среде MatLab и Simulink [Текст]: учебное пособие / В.С. Щербаков, А.А. Руппель, В.А. Глушец. - Омск: Изд-во СибАДИ, 2003. - 160 с.

5. Щербаков В.С. Теория автоматического управления [Текст]: методические указания по выполнению курсовых работ (для студентов специальностей 220301, 140607) / В.С. Щербаков, Р.Ю. Сухарев. - Омск: СибАДИ, 2011. - 36 с.