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

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

По таблице 3.1 строим годограф Михайлова

Рис. 3.1 Годограф Михайлова

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

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

Критерий Г. Найквиста позволяет по амплитудно-фазовой частотной характеристике (АФЧХ) разомкнутой системы оценить устойчивость замкнутой системы с отрицательной обратной связью.

АФЧХ разомкнутой системы можно построить на комплексной плоскости [+1; j], откладывая вещественную UРАЗ(щ) и мнимую VРАЗ(щ) составляющую частотной передаточной функции разомкнутой системы WРАЗ(jщ), или в полярной системе координат, откладывая угол фазы ц() разомкнутой системы относительно оси +1 и в этом направлении откладывать вектор амплитуды разомкнутой системы А().

Вещественную UРАЗ(щ) и мнимую VРАЗ(щ) функции WРАЗ(jщ) можно определить по АЧХ АРАЗ() и ФЧХ цРАЗ() разомкнутой системы:

;

.

АЧХ передаточной функции разомкнутой системы АРАЗ() равна произведению АЧХ отдельных звеньев, а ФЧХ цРАЗ() - сумме ФЧХ звеньев:

;

.

Найти АЧХ А() и ФЧХ ц() можно по вещественной U(щ) и мнимой V(щ) составляющим частотной передаточной функции W(jщ) звена.

Амплитуда А() и фаза ц() частотной передаточной функции W(jщ):

;

.

Формулировка критерия Найквиста

Если разомкнутая система устойчива, то для устойчивости системы в замкнутом состоянии необходимо и достаточно, чтобы АФЧХ разомкнутой системы на комплексной плоскости [+1; j] при изменении частоты от 0 до не охватывала точку с координатами (-1; j0). Если АФЧХ разомкнутой системы проходит через точку с координатами (-1; j0), то система находится на границе устойчивости.

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

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

;

.

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

;

.

;

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

.

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

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

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

АЧХ А()

ФЧХ ц()

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

Частота

щ

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

Фазы звеньев

Амплитуда

WРАЗ(jщ)

Фаза

WРАЗ(jщ)

Действит.

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

Мнимая

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

A1(щ)

A2(щ)

ц1(щ)

ц2(щ)

AРАЗ(щ)

цРАЗ(щ)

UРАЗ(щ)

VРАЗ(щ)

0,0

8,190

1,000

0,000

0,000

8,190

0,000

8,19000

0,00000

1,0

7,496

0,893

-0,415

-0,655

6,692

-1,069

3,21843

-5,86721

2,0

6,148

0,679

-0,722

-1,185

4,175

-1,906

-1,37434

-3,94283

4,0

4,046

0,352

-1,054

1,292

1,426

0,238

1,38568

0,33590

6,0

2,901

0,197

-1,209

0,937

0,571

-0,272

0,55031

-0,15333

8,0

2,238

0,122

-1,294

0,728

0,273

-0,566

0,23028

-0,14639

10,0

1,815

0,082

-1,347

0,592

0,149

-0,755

0,10824

-0,10183

12,0

1,524

0,058

-1,384

0,499

0,089

-0,885

0,05640

-0,06892

14,0

1,312

0,044

-1,410

-2,712

0,057

-4,121

-0,03192

0,04758

16,0

1,152

0,034

-1,430

-2,764

0,039

-4,194

-0,01931

0,03381

18,0

1,026

0,027

-1,445

-2,805

0,028

-4,250

-0,01232

0,02471

20,0

0,925

0,022

-1,458

-2,838

0,020

-4,296

-0,00821

0,01854

22,0

0,842

0,018

-1,468

-2,865

0,015

-4,333

-0,00567

0,01422

24,0

0,772

0,015

-1,476

-2,888

0,012

-4,364

-0,00404

0,01113

26,0

0,713

0,013

-1,484

-2,907

0,009

-4,391

-0,00295

0,00886

28,0

0,663

0,011

-1,490

-2,924

0,007

-4,414

-0,00221

0,00717

30,0

0,619

0,010

-1,495

-2,938

0,006

-4,433

-0,00168

0,00587

32,0

0,580

0,009

-1,500

-2,951

0,005

-4,451

-0,00130

0,00487

34,0

0,546

0,008

-1,504

-2,962

0,004

-4,466

-0,00103

0,00408

36,0

0,516

0,007

-1,508

-2,972

0,004

-4,480

-0,00082

0,00345

38,0

0,489

0,006

-1,511

-2,981

0,003

-4,492

-0,00066

0,00295

40,0

0,465

0,006

-1,514

-2,989

0,003

-4,503

-0,00054

0,00254

42,0

0,443

0,005

-1,517

-2,996

0,002

-4,513

-0,00044

0,00220

44,0

0,422

0,005

-1,519

-3,003

0,002

-4,522

-0,00037

0,00192

46,0

0,404

0,004

-1,521

-3,009

0,002

-4,530

-0,00031

0,00168

48,0

0,387

0,004

-1,523

-3,014

0,002

-4,538

-0,00026

0,00148

50,0

0,372

0,004

-1,525

-3,019

0,001

-4,545

-0,00022

0,00131

52,0

0,358

0,003

-1,527

-3,024

0,001

-4,551

-0,00019

0,00117

54,0

0,344

0,003

-1,529

-3,028

0,001

-4,557

-0,00016

0,00104

56,0

0,332

0,003

-1,530

-3,032

0,001

-4,563

-0,00014

0,00094

58,0

0,321

0,003

-1,532

-3,036

0,001

-4,568

-0,00012

0,00084

60,0

0,310

0,002

-1,533

-3,040

0,001

-4,573

-0,00011

0,00076

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

Рис. 3.2.Годограф Найквиста (АФЧХ)

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

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

Амплитудная частотная характеристика (АЧХ) - это зависимость отношения амплитуды сигнала на выходе звена к амплитуде на входе A = ymax/xmax от частоты входного сигнала .

Фазовая частотная характеристика (ФЧХ) - зависимость угла сдвига по фазе сигнала на выходе звена от частоты входного сигнала .

Логарифмические частотные характеристики - АЧХ и ФЧХ, представленные в логарифмическом масштабе.

Параметры L и lgщ определяются следующим образом:

где ;

.

Ордината ЛАХ L измеряется в децибелах [дБ].

Децибел - логарифмическая единица уровней, затуханий и усилений.

1 Б = 10 дБ - это увеличение мощности сигнала в 10 раз.

2 Б = 20 дБ - это увеличение мощности сигнала в 100 раз;

3 Б = 30 дБ - это увеличение мощности сигнала в 1000 раз и т.д.

Абсцисса ЛАХ и ЛФХ lgщ измеряется в декадах [дек].

Декада - логарифмическая единица частот, соответствующая изменению частоты щ в 10 раз.

Логарифмические критерии, так же как и критерий Найквиста, позволяют судить об устойчивости замкнутой системы по виду ЛАХ и ЛФХ разомкнутой системы.

Формулировка логарифмического критерия

Если разомкнутая система устойчива, то для устойчивости системы в замкнутом состоянии, ЛАХ разомкнутой системы должна пересечь ось абсцисс раньше, чем ЛФХ, спадая окончательно, перейдёт через значение - р. То есть, на частоте среза щср величина фазы ц должна быть меньше значения | - р |.

Запас устойчивости по амплитуде ДL - это величина допустимого увеличения общего коэффициента усиления разомкнутой системы, при котором замкнутая система окажется на границе устойчивости .

Запас устойчивости по фазе Дц - это величина допустимого увеличения запаздывания по фазе разомкнутой системы на частоте среза щср, при котором замкнутая система окажется на границе устойчивости .

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

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

;

.

Таблица 3.4 Асимптотические логарифмические частотные характеристики типовых динамических звеньев

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

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

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

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

Фаза ц()

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

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

Частота

щ

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

Общая ЛАХ

LРАЗ(щ)

Фазы звеньев

Общая ЛФХ

цРАЗ(щ)

L1(щ)

L2(щ)

ц1(щ)

ц2(щ)

0,100

18,257

-5,194

-0,001

13,062

-0,215

-0,020

0,500

18,060

-5,236

-0,013

12,811

-0,747

-0,100

0,900

17,633

-5,331

-0,042

12,260

-0,889

-0,178

1,300

17,036

-5,477

-0,088

11,472

-0,900

-0,254

1,700

16,336

-5,668

-0,149

10,519

-0,868

-0,328

2,100

15,585

-5,898

-0,226

9,461

-0,820

-0,398

2,500

14,822

-6,162

-0,318

8,341

-0,769

-0,464

2,900

14,069

-6,452

-0,425

7,192

-0,718

-0,526

3,300

13,340

-6,763

-0,546

6,031

-0,670

-0,583

3,700

12,642

-7,089

-0,680

4,873

-0,626

-0,637

4,100

11,977

-7,426

-0,826

3,725

-0,585

-0,687

4,500

11,346

-7,770

-0,984

2,592

-0,549

-0,733

4,900

10,746

-8,116

-1,153

1,477

-0,516

-0,775

5,300

10,178

-8,463

-1,332

0,383

-0,486

-0,815

5,700

9,638

-8,809

-1,520

-0,691

-0,459

-0,851

6,100

9,125

-9,152

-1,716

-1,743

-0,435

-0,884

6,500

8,637

-9,490

-1,921

-2,773

-0,412

-0,915

6,900

8,172

-9,823

-2,132

-3,783

-0,392

-0,944

7,300

7,728

-10,150

-2,350

-4,772

-0,374

-0,970

7,700

7,304

-10,471

-2,573

-5,740

-0,357

-0,995

8,100

6,898

-10,785

-2,801

-6,688

-0,341

-1,018

8,500

6,508

-11,092

-3,034

-7,618

-0,327

-1,039

8,900

6,134

-11,392

-3,270

-8,528

-0,313

-1,059

9,300

5,775

-11,686

-3,510

-9,421

-0,301

-1,077

9,700

5,429

-11,972

-3,753

-10,296

-0,290

-1,095

10,100

5,096

-12,252

-3,998

-11,154

-0,279

-1,111

12,100

3,590

-13,554

-5,248

-15,211

-0,236

-1,179

14,100

2,301

-14,712

-6,508

-18,919

-0,204

-1,230

16,100

1,174

-15,750

-7,754

-22,330

-0,179

-1,270

18,100

0,175

-16,686

-8,971

-25,482

-0,160

-1,301

20,100

-0,722

-17,538

-10,148

-28,409

-0,144

-1,327

22,100

-1,537

-18,318

-11,283

-31,138

-0,132

-1,348

24,100

-2,282

-19,037

-12,372

-33,691

-0,121

-1,366

26,100

-2,969

-19,703

-13,417

-36,089

-0,112

-1,382

30,000

-4,171

-20,875

-15,331

-40,376

-0,097

-1,406

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

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

3.5 Проверка статистической ошибки

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

Раздел 4. Анализ автоматической системы регулирования, c корректирующим звеном.

4.1 Расчет передаточной функции корректирующего звена

Коррекция динамических свойств линейных систем осуществляется для выполнения рассмотренных ранее требований по точности, устойчивости и качеству переходных процессов. Осуществляется коррекция путем введения в систему специальных корректирующих звеньев с особо подобранной передаточной функцией WКЗ(p).

В данном случае корректирующее звено включается в систему относительно прямой ветви последовательно, образуя новую структуру системы (рис. 4.1).

Рис. 4.1 Блок-схема стабилизации угловой скорости ДПТ с корректирующим звеном

Корректирующие звенья аналоговых линейных электрических систем обычно выполняются в виде пассивных или активных четырехполюсников (АЧ и ПЧ) постоянного тока.

Рассмотрим принципиальную электрическую схему заданного корректирующего звена в виде пассивного четырехполюсника (рис. 4.2).

Рис.4.2.Пассивное корректирующее звено.

Параметры корректирующей цепи:

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

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

Где:

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

4.2 Исследование скорректированной АСР по критерию Гурвица

Для исследования устойчивости замкнутой АСР по критерию устойчивости Гурвица понадобится характеристический полином замкнутой скорректированной системы.

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

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

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

; (4.9)

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

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

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

%Вводим в командную строку Matlab:%

>> syms p;

expand((2.5*p+1)*(0.2*p+1)*(0.035*0.15*p^2+0.15*p+1)+0.235*(0.3*p+1)*10*0.5*1.3*0.36)

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