Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

учебное_пособие_часть_2_МСПД_new

.pdf
Скачиваний:
28
Добавлен:
29.05.2015
Размер:
7.72 Mб
Скачать

Рис 2.2

Пусть задана принципиальная схема фильтра низких частот, изображенная на рис.2.3.

L1 R1

C1

R2 Uвых

Uвх

L2

Рис.2.3

103

Фильтр имеет следующие параметры:

L1=50 мГн, R1=1 Ом, С1=10 мкФ, L2=5 мГн, R2=100 Ом.

Для создания подобной схемы в среде Matlab Simulink предусмотрены модели последовательной и параллельной R-L-C цепочки, кото-

рые находятся в библиотеке SimPowerSystems / Elements / Series RLC Branch и SimPowerSystems / Elements / Раrallel RLC Branch.

В данном случае целесообразно использовать последовательную R-L-C цепочку. Если необходимо исключить какой-либо элемент из цепочки, то нужно в окне задания параметров установить «inf» (бесконечность) для исключения емкости и «0» для исключения сопротивления и индуктивности. При параллельной R-L-C цепочке для исключения емкости нужно установить «0» и «inf» для исключения сопротивления и индуктивности (рис.2.4).

Рис.2.4

Подключаем к схеме источник переменного напряжения, который находится в библиотеке SimPowerSystems / Electrical Sources / AC Voltage Sources (рис. 2.5).

104

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

лиотеке SimPowerSystems / Measurements / Voltage Measurement. Для визуализации переходных процессов нужно к вольтметру подсоединить осциллограф, который находится в библиотеке Simulink / Sinks / Scope

(рис. 2.6).

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

Рис. 2.5

105

Рис. 2.6

Рис. 2.7

106

§2. Вывод передаточных функций для заданной электрической схемы

Для заданной принципиальной схемы (рис. 2.3) по второму закону Кирхгофа запишем дифференциальные уравнения:

U

 

 

di(t )

 

 

(2.1)

ВХ (t )= L1

+i(t ) R1 +UВЫХ (t );

 

 

 

 

 

dt

 

 

 

 

 

 

 

1

 

 

 

 

 

di(t )

 

U

ВЫХ (t )=

 

i(t ) dt +i(t ) R2

+ L2

 

.

C1

 

 

 

 

 

 

 

dt

Дифференцируя второе уравнение системы дифференциальных уравнений (2.1), получим:

U

ВХ (t )= L1

di(t )

+i(t ) R1 +U

ВЫХ (t );

 

 

(2.2)

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

di(t )

 

 

d

2

i(t )

 

dUВЫХ (t )

=

i(t )+

R2

+ L2

 

.

 

 

dt

C

dt

 

dt2

 

 

 

1

 

 

 

 

 

 

 

 

 

Передаточная функция – это отношение изображения по Лапласу выходной координаты изображению к входной. Чтобы получить передаточную функцию, необходимо перейти от временной формы записи дифференциальных уравнений к операторной. Данное преобразование называется преобразованием Лапласа, при котором в искомом диффе-

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

После преобразования Лапласа система уравнений будет выглядеть следующим образом:

U

 

(p)

= L

p i(p)+i(p) R

+U

 

(p);

(2.3)

 

ВХ

 

 

1

1

 

ВЫХ

 

 

p

U

 

(p)

=

1

i(p)+i(p) R p +

L p2

i(p).

ВЫХ

 

 

 

 

 

 

C1

2

 

2

 

 

 

 

 

 

 

 

 

 

 

Выразим из второго уравнения системы (2.3) значение i(p) и подставим в первое уравнение. После несложных преобразований получим передаточную функцию заданной схемы:

107

W (p)=

 

 

 

 

 

1+ p R2 C1 + L2 C1 p2

 

 

 

 

 

(2.4)

 

L C p2

+ p R C +1

+ p R C + L C p2

 

 

 

 

1

1

 

 

1

 

1

 

2

1

2

1

 

 

 

 

следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W (p)

=

 

 

 

1+ p R2

C1 + L2 C1 p2

 

 

 

.

 

(2.5)

p2 (L C + L C )+ p (R C + R C )+1

 

 

 

 

 

 

 

1

 

1

2

1

 

1

1

2

 

1

 

 

 

 

Сделаем обозначения T 2

= L C ,

T = R C , T 2

= L C , T = R C

 

 

 

 

 

 

1

 

2

1

2

2

1

3

 

1

1

4

1

1

– постоянные времени электрической цепи.

 

 

 

 

 

 

 

 

 

 

 

W (p)=

 

 

T12 p2 +T2 p +1

 

.

 

 

(2.6)

 

 

(T32 +T12 ) p2 +(T4 +T2 ) p +1

 

 

 

 

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

Т-образную электрическую схему в общем виде можно представить, как показано на рис. 2.8.

Z1(p)

Z3(p)

 

 

 

 

Uвх(p)

Z2(p)

Uвых (p)

 

Рис.2.8

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

108

W (p)=

 

Z2 ( p)

(2.7)

Z

( p) + Z

2

( p)

 

1

 

 

 

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

W (p)=

 

 

Z21 (p) R

 

 

 

 

,

(2.8)

Z

(p) R Z

(p) Z

22

(p)+ Z

12

(p) Z

21

(p)

 

11

11

 

 

 

 

 

 

где Z21 (p)= Z2 (p), Z11 (p)= Z1 (p)+ Z2 (p), Z22 (p)= −(Z2 (p)+ Z3 (p)), Z12 (p)= −Z21 (p).

Если электрическая схема построена по П-образной схеме (рис.2.9), то передаточная функция на холостом ходу вычисляется по формуле [5]:

W (p)=

Z

(p)

(2.9)

41

 

.

 

Z31

(p)

 

Z5(p)

Uвх(p)

Z4(p)

Z6(p)

Uвых (p)

 

Рис. 2.9

Числитель передаточной функции (2.9) определяется:

Z41

(p)=

 

Z4 (p) Z6 (p)

 

.

(2.10)

Z4

(p)+ Z5 (p)+ Z6

(p)

 

 

 

 

 

109

Знаменатель передаточной функции (2.9) определяется:

Z31

(p)=

 

Z4 (p) Z6 (p)

 

.

(2.11)

Z4

(p)+ Z5 (p)+ Z6

(p)

 

 

 

 

 

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

W (p)=

 

 

 

 

Z41 (p) R

 

 

 

 

,

(2.12)

Z

31

(p) R Z

31

(p) Z

44

(p)+ Z

14

(p) Z

41

(p)

 

 

 

 

 

 

 

 

 

 

где

Z44 (p)= (Z4((p))++ Z5((p))+) Z6((p)); Z14 (p)= Z41 (p).

Z4 p Z5 p Z6 p

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

§3. Построение графика переходного процесса и частотных характеристик по передаточной функции звена

Переходная характеристика – это реакция системы на единичное входное воздействие. Объектом исследования в данном случае является передаточная функция, которая задается с помощью звена Transfer Fcn из библиотеки Simulink / Continuous. Передаточная функция исходной системы в общем виде представляется в виде:

 

ams

m

+ am1s

m1

 

(2.13)

W (s) =

 

 

+... + a1s + a0

.

 

 

 

 

 

 

b sn +b

sn1 +... +b s +b

 

n

 

n1

 

 

1

0

 

Коэффициенты аi следует вводить в поле Numerator, начиная с коэффициента аm при старшей производной. Аналогично заполняется поле знаменателя передаточной функции Denominator, начиная с коэффициента bn.

Входное воздействие на передаточную функцию подается с по-

мощью элемента Constant из библиотеки Simulink / Commonly Used

110

Blocks, а выходной сигнал регистрируется осциллографом Simulink / Sinks / Scope (рис. 2.10).

Рис.2.10

При построении переходного процесса существенно влияет на результат выбор численного метода расчета (Solver), шага расчета (Relative tolerance), начального и конечного значения времени переходного процесса (Start time, Stop time). Рекомендуется устанавливать параметры как указано на рис. 2.11.

111

Рис. 2.11

Получение графика переходного процесса в виде, удобном для обработки в графических редакторах, а также получение графиков частотных характеристик возможно при подключении к выходу пе-

редаточной функции компонента «Simulink / Sinks / To Workspace»

(рис. 2.12).

В настройках нового компонента следует обязательно указать «Array» (Массив) в поле «Save Format» (Формат сохранения) (рис. 2.13).

Повторяем расчёт переходных процессов, после чего переходим в окно диспетчера MATLAB и далее с помощью команды «Window\Workspace» переходим в окно «Workspace», где отображаются переменные текущего проекта (рис. 2.14).

112