ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.03.2024
Просмотров: 135
Скачиваний: 0
ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
СОДЕРЖАНИЕ
4) Определяем напряжение на индуктивности U(t):
u(t)Т/—i(t) u(t) = p-A-ep 1 = -^-L— -ep 1 = -EeP 1 dt L R1
t
Puc. 4.11
Пример-2.
Ищем решения в виде:
u(t) = ucB(t) + unp = Aef 1 + unp
-
inp определяет принуждённую составляющую схеме после коммутации : цпр :=Е unp 20 unp := F -
из ИНУ определяет константу интегрирования А в схеме до коммутации : u(-0) = и(0) = 0 = А + unp Aj- -imp А - -20 -
Корень характеристического уравнения через входное сопротивление в схеме после коммутации :
Z=RI+ — = Ср:- — Р--1.667Х103
р-С ** RI-C
-
Записываем окончательное решение и строим график i(t):
uft) Л-еР 1 + unp т-Д- Т“6х 10 4 t0, т-0.5.. 4 т
|р|
up) :=-F. ep * + Е
4) Определяем напряжение на индуктивности U(t):
С = 6х1(Г5 i(t) := С —u(t) i(t) = С-р Л-еР ‘ СЕер ‘ = — -еР ’ dt CRI R1
и » “"Е p-t
t
| | | | | R2 | Ri L | |
| "Гример- Гщем рс i(t) = icB(t
inp. RI-
A e1 t:= O,tOJ = 0 0.0033 0.0067 0.01 0.0133 0.0167 0.02 0.0233 0.0267 | 1 ШС ) + pc; E 4- R У ) = b X fyr R2 ЫВ< ?t j..4 | ПИЯ в в inp = A- 1СЛЯСТ - inp - апреле; E R2+ — 2 арактс] ации: + p-L = гем око hinP £ 1(0 = 0.4 0.505 0.569 0.607 0.631 0.645 0.653 0.659 0.662 | Puc. 4.15 иде: p t e + inp 1ринуждённую составляющую схеме после коммутации: 0.667 1яет константу интегрирования А в схеме до коммутации: •— = А + А : — inp А - -0.267 2 RI + R2 RI 2 R2+ — 2 мистического уравнения через входное сопротивление схеме по- 0 RI + R2+x-Lsolve,х -+-150 -150 нчатсльное решение и строим график i(t): :=-Дг т = 6.667х ю“3 |р| |
Пример- 4.
Дано:
R1 := 10 R2 := 20 L := 0.2 Е 20 С •- 60-10" 6
ЛАДДА ААААА ААА ДАД АЛА/
Ищем решения в виде:
р-t i( t) = icb( t) + inp = Л-е + inp
-
inp определяет принуждённую составляющую схеме после коммутации:
-
определяет ННУ в схеме до коммутации: ilX-O) = ЩО) = о
-
определяет ЗНУ в схеме после коммутации (Рис. 4.19):
i( 0 + ■) = — io—-— io - 0.667 RI + R2 R1 + R2
А := -inp + io Л = -0.133 АЛАА •
-
Корень характеристического уравнения через входное сопротивление в схеме
после коммутации:
„ RIR2 , Л 1 RI R2
Z — + р-1. = 0 п
RI + R2 L RI + R2
р - -33.333
RI-R2
RI + R2
- 6.667
-
Записываем окончательное решение и строим график i(t):
t =
(t) =
ДО —-еР 1 + —
м RI R1
0
0.667
0.015
0.719
0.03
0.751
0.045
0.77
0.06
0.782
0.075
0.789
0.09
0.793
0.105
0.796
0.12
0.798
-
Unp определяет принуждённую составляющую схеме после коммутации
-
(рис. 4.22):
Unp:-Г* Unp-20
2) Определяет ИНУ в схеме до коммутации (рис. 4.23):
що>.iS.fi.„Г' js,iS.fi.R'rl
2 \ 2 ) *** 2 I 2 )
-
Определяет ЗНУ в схеме после коммутации (рис. 4.24):
U(0) = R (J - Щ0)) Uo:-R (J-io) Uo - 26.667
А -Unp + Uo А - 6.667 АЛЛА I
-
Корень характеристического уравнения через входное сопротивление в схеме после коммутации (рис. 4.22):
Z=y+y + R+p-L=0 р = -200
RIR2
RI + R2
- 6.667
5) Записываем окончательное решение и строим график i(t):
U(i) := Л-ер 1 + Unp г:-—!— т - 0.005 t0, т-0.5.. 4-т
|р|
t -
U(l) =
0
26.67
0.0025
24.04
0.005
22.45
0.0075
21.49
0.01
20.9
0.0125
20.55
0.015
20.33
0.0175
20.2
0.02
20.12
Rj- 20 £/-=40-106 ^:=2
Ищем решения в виде:
U(t) = Ucb(i) + inp = ЛеР 1 + Unp
1) Unp определяет принуждённую составляющую схеме после коммутации (рис. 4.27):
Unp J R-3 Unp - 120
АААААЛА •
-
определяет ННУ в схеме до коммутации (рис. 4. 28):
А/с. 4.2Я Рис. 4.29
-
определяет ЗНУ в схеме после коммутации (рис. 4.29 и 4.30):
*:= М- W = J (Re + R) + Ее
J (Re R) + Ее Uo - 93.333
-Unp + Uo А = -26.667
-
Корень характеристического уравнения через входное сопротивление в схеме после коммутации:
-
Записываем окончательное решение и строим график i(t):
т:=-гт т=2.4х IO-3 t 0,т-0.5.. 4-т
1р|
t- U(t)= UJt)А еР ’ + Unp
0
26.67
0.0012
25.24
0.0024
24.13
0.0036
23.25
0.0048
22.55
0.006
22.01
0.0072
21.58
0.0084
21.24
0.0096
20.98
Лекция № 10
§4.3 Метод расчета переходных процессов в цепи переменный тока
Для расчета переходный процессов в цепи переменного тока используют символический метод
Пример: Определи 1ь ток источника напряжения если:
/?1=20Ом, R2 = 10 Ом, L = 0,2 Гн, е(/) - 20 sin( со/ - 60° )В, f = 50 Гц, со-27и/’-313рад/с.
Решение: Находим индуктивное сопротивление и комплекс напряжения XL - (&L - 62,80м, £ - 20е”/6° В.
Ищем решение в виде i(/)-iCB(/) + /np(/)- Аер1 + /(/).
1. Определяем принужденную составляющую в цепи после коммутации, используя символический метод
Z(p) -Rt + R2 + pL -0 p _-R'+R1- _-150c’1.
L
-
Определяем независимые начальные условия, i£(0) используя символический метод.
£0
I ‘ 2/?|+у%£ J
D
! = -0,118 - /О,156 A, i, (0) =
2Л,+^£ J ,/Л'
-
= -0,156 А.
ie(0+) = fr(0)=-0,156A
-
Определяем зависимые начальные условия в схеме после коммутации, заменяя индуктивность/, источником тока равным <(0).Определяем константу интегрирования А
/,.(0+) - А + ги/,(0) -» А -/„(0+) -/„„(0) -0,081 А.
Записываем решение и строим график.
i(t) -гсв(/) + Гпр(О-^' + гпр(/) -0,081е’150' +0,287sin(w/-124,466°)A . 6
1 ... 10 11 12 13 14 15 16 17 ... 21
inp определяет принуждённую составляющую схеме после коммутации:
определяет ННУ в схеме до коммутации: ilX-O) = ЩО) = о
определяет ЗНУ в схеме после коммутации (Рис. 4.19):
Корень характеристического уравнения через входное сопротивление в схеме
-
Записываем окончательное решение и строим график i(t):
t = | (t) = | ДО —-еР 1 + — м RI R1 | |
0 | | 0.667 | |
0.015 | | 0.719 | |
0.03 | | 0.751 | |
0.045 | | 0.77 | |
0.06 | | 0.782 | |
0.075 | | 0.789 | |
0.09 | | 0.793 | |
0.105 | | 0.796 | |
0.12 | | 0.798 | |
Unp определяет принуждённую составляющую схеме после коммутации
(рис. 4.22):
Unp:-Г* Unp-20
2) Определяет ИНУ в схеме до коммутации (рис. 4.23):
що>.iS.fi.„Г' js,iS.fi.R'rl
2 \ 2 ) *** 2 I 2 )
-
Определяет ЗНУ в схеме после коммутации (рис. 4.24):
U(0) = R (J - Щ0)) Uo:-R (J-io) Uo - 26.667
А -Unp + Uo А - 6.667 АЛЛА I
-
Корень характеристического уравнения через входное сопротивление в схеме после коммутации (рис. 4.22):
Z=y+y + R+p-L=0 р = -200
RIR2
RI + R2
- 6.667
5) Записываем окончательное решение и строим график i(t):
U(i) := Л-ер 1 + Unp г:-—!— т - 0.005 t0, т-0.5.. 4-т
|р|
t - | | U(l) = |
0 | | 26.67 |
0.0025 | 24.04 | |
0.005 | 22.45 | |
0.0075 | 21.49 | |
0.01 | 20.9 | |
0.0125 | 20.55 | |
0.015 | 20.33 | |
0.0175 | 20.2 | |
0.02 | 20.12 |
Rj- 20 £/-=40-106 ^:=2
Ищем решения в виде:
U(t) = Ucb(i) + inp = ЛеР 1 + Unp
1) Unp определяет принуждённую составляющую схеме после коммутации (рис. 4.27):
Unp J R-3 Unp - 120
АААААЛА •
-
определяет ННУ в схеме до коммутации (рис. 4. 28):
А/с. 4.2Я Рис. 4.29
-
определяет ЗНУ в схеме после коммутации (рис. 4.29 и 4.30):
*:= М- W = J (Re + R) + Ее
J (Re R) + Ее Uo - 93.333
-Unp + Uo А = -26.667
-
Корень характеристического уравнения через входное сопротивление в схеме после коммутации: -
Записываем окончательное решение и строим график i(t):
т:=-гт т=2.4х IO-3 t 0,т-0.5.. 4-т
1р|
t- U(t)= UJt)А еР ’ + Unp
0 | | 26.67 |
0.0012 | 25.24 | |
0.0024 | 24.13 | |
0.0036 | 23.25 | |
0.0048 | 22.55 | |
0.006 | 22.01 | |
0.0072 | 21.58 | |
0.0084 | 21.24 | |
0.0096 | 20.98 |
Лекция № 10
§4.3 Метод расчета переходных процессов в цепи переменный тока
Для расчета переходный процессов в цепи переменного тока используют символический метод
Пример: Определи 1ь ток источника напряжения если:
/?1=20Ом, R2 = 10 Ом, L = 0,2 Гн, е(/) - 20 sin( со/ - 60° )В, f = 50 Гц, со-27и/’-313рад/с.
Решение: Находим индуктивное сопротивление и комплекс напряжения XL - (&L - 62,80м, £ - 20е”/6° В.
Ищем решение в виде i(/)-iCB(/) + /np(/)- Аер1 + /(/).
1. Определяем принужденную составляющую в цепи после коммутации, используя символический метод
Z(p) -Rt + R2 + pL -0 p _-R'+R1- _-150c’1.
L
-
Определяем независимые начальные условия, i£(0) используя символический метод.
£0
I ‘ 2/?|+у%£ J
D
! = -0,118 - /О,156 A, i, (0) =
2Л,+^£ J ,/Л'
-
= -0,156 А.
ie(0+) = fr(0)=-0,156A
-
Определяем зависимые начальные условия в схеме после коммутации, заменяя индуктивность/, источником тока равным <(0).Определяем константу интегрирования А
/,.(0+) - А + ги/,(0) -» А -/„(0+) -/„„(0) -0,081 А.
Записываем решение и строим график.
i(t) -гсв(/) + Гпр(О-^' + гпр(/) -0,081е’150' +0,287sin(w/-124,466°)A . 6
t =
0
0 001
0.002
0 003
0.004
0 005
0 006
0 007
0.008
0 009
0.01001
0.01101
0.01201
0.01301
001401
0.01501
П) =
-0.156
-0.206
-0.227
-0.219
-0.183
-0.124
-0 049
0.036
012
0.196
0.255
0.291
03
0.282
0.237
0.171
Рис. 4.36
§4.4 Переходные процессы в цепи второго порядка
Рассмотрим цепь второго порядка представленную на рис. 4.37 с параметрами: Е = 50В, R - 1 ООм, L = 0.1 Гн, С - 40мкФ.
Записываем уравнения по второму закону Кирхгофа, в результате получаем систему дифференциальных уравнений:
+ис+/•_/? =£* —> L — + Ri + uc=E. i = C—;
Первое слагаемое это исв = 4 -ехр(/>|/) + А2 ехр(р2Осв°бодная составляющая. Она зависит только от параметров схемы, а также от начальных и конечных запасов энергии. Эта составляющая решения не зависит от формы воздействующего напряжения.
Второе слагаемое это ww/, = uc(oo) принуждённая составляющая. Она зависит от внешнего воздействия и имеет форму этого воздействия. Очевидно, что в нашем случае она определяется как ulip = i/c(oo) - Е.
Постоянные интегрирования определяются из начальных условий, отражающих невозможность мгновенного изменения начальных запасов энергии в конденсаторе и в катушке.
Для определения констант интегрирования используем независимые начальные условия ис(0) = 0, iL(0) = 0.
м^(0) — 0 — /4| + А2 + Е\ / к zA(O)-f^(O+)-C^-O = C'(^lp1 + Л2р,).
dt
Откуда следует, что
. (4)
P2
Pl P2
Pl
Теперь можно записать окончательное решение
иг(0 = —ехр(р/)+ Р,£ exp(p2t)+E = ———(-р2ер1 + р,еА' )■
Р2-Р1 Р2-Р1 “ Р2
Р\
Определим корни характеристического уравнения входящие в решение uc\t) Р\^ Р2 через входное сопротивление схемы.
1 CL • р
+ RC• р +1 А 2 п/> 1 Л /с\
pL + н R — — 0 —> CL - р 4- RC • р 4" 1 — 0. (5)
Ср Ср
В результате решения уравнения получаются корни:
-h±4o -RC±J(RC)2-4LC R If R f Г
1,2 2а 2CL 2L LC (6)
г- о 1
Где о = — показатель затухания контура, соо - - .— угловая часто-
та незатухающих колебаний, при выполнении условия (Oq > б2 имеем
Здесь соГ6, частота свободных колебаний,
Корни уравнения определяются параметрами цепи и могут принимать следующие возможные значения (рис. 4.38).
Дискриминант равен нулю. Кони вещественные, отрицательные и тг с R
кратные. Критический режим р^ 2 =“'О = - —
uc(/) = £(l + 8/)e"5' + Е.
• Дискриминант положительный. Корни вещественные отрицательные и неравные. Апериодический режим р} 2 = -6± Jb2 -со^ ;
Uc (0 = —/ Е (Р\ еР:‘ - Р-> qP}> ) + С.
• Дискриминант отрицательный. Корни комплексно-сопряжённые, с отрицательной вещественной частью. Колебательный режим
=-d±j(0„.
uc(t)-Ee-6'
cos(w,,60 + Sin(wrsr) + E .
p1=p2=-6
+1
“ ^св
Puc. 4.38. Расположение корней на комплексной плоскости.
Примеры определения корней характеристического уравнения в
Mathcad
-456.2341361360570100)
-182.6547527528318788^
R 10 С:=60-10 6 L:=0.2
(R+Lp)-2R I , p : + + R solve ,p —►
R+Lp+2R Cp
(R+ Lp)-2 R J_
R + E-р + 2-R {"'’P
2 2
5 R С р + 3 R Ср L + 3R+ Lp
(3-R + Lp)Cp
Р“
-456.234")
-182.655)
С- 100 ю 6
L:=0.1
Л (R+LP)R I
P := +
R+LP+R CP
R solve, Р —>
(-275) - 156.12494995995995515
.(-275) + 156.1249499599599551Н
Р-
-275- 156.125
-275 + 156.125iJ
(R+LP)R I
+ + R
R+LP+R CP
3R2CP + 2R-C P2L+ 2-R + LP
(2R+ LP)CP
Примеры определения корней характеристического уравнения и
зависимых и независимых начальных условий
Пример: Определить независимые/} (0),(/с(0) и зависимые начальные
условия (/£(0+),/с(0+). Определить корень характеристического урав-
22Я+Д
2
(УС(О) = /£(О)Я = -.
1 .Определяем независимые начальные условия /}(0)5С/с(0).
2.
Определяем зависимые начальные условия (7Z (0+),/’с(0+) из схемы после коммутации (см. документ Mathcad).
3.
Определяем корень характеристического уравнения из схемы после коммутации (см. документ Mathcad).
Документ Mathcad
ORIGIN:= 1
АААЛАЛААЛЛАА
Определи гь напряжение па конденсаторе.
Е?= 100 R 10 L:- 0.1 Г;-50106
Е 2R(Lp + R) I
UnD := - Z р) := — + R +
пр 3 З-R + bp Ср
Корни характеристического уравнения
solve,р Г(—416.62)-162.451'
р •- Z(p)
float, 5
(-416.61) + 162.45i_
' -416.67- 162.45^
ч -416.67+ 162.451J
Независимые начальные условия Е п Е 'Lo 1R Ьс0‘ 2
Зависимые начальные условия
Е“ иС0+ *Lo*r ’Lo"5 UC0"50 ’R0“
iCO:",RO",Lc
Пос гоя иные интегрирования
- I ( иС0" ипр
<АП
Л)
р.ЗЗЗ-81.2211^
t 8.333+ 8l.22liJ
со - 162.45 Unp - 33.333
ч । y 1
*Lo *np
UL0 ।
< L J
цсо+ 'LoR 2R
UL0^ F“ ‘R0R
pit pj-t
i(t):=Are + Aye ‘
t = 3.O31x 10 3 t:-0,0.0l-T..T-IC
I
|Rc(P)|
-0.417 + 0.716i>
-0.417-0.7l6iJ
§4.5 Операторный метол расчёта переходных процессов
Операторный метод (преобразование Лапласа) расчета переходных процессов используется для того, чтобы обыкновенные дифференциальные уравнения с постоянными коэффициентами (в пространстве оригиналов) преобразовать в алгебраические (в пространстве изображений). Очевидно, что алгебраические уравнения решаются проще. После решения алгебраического уравнения над полученной функцией (изображением) производится обратное преобразование Лапласа, получается оригинал. Полученный оригинал это функция, которая и буде1 решением дифференциального уравнения.
Любой функции можно сопоставить её преобразование Лапласа
F{p)^]f{t)ep,dt, (7)
О
здесь F(p) изображение, /'(/) оригинал. Выражении (7) записывают ещё и в операторной формеF(p) = L[/(/)] •
Приведём изображение нескольких часто встречающихся функций. Определим изображение константы - /’(/) = A (const):
00 j
F(p) = A\e-t>'dt=-— |;=-.
о Р Р
F(p) = f e(l'ep‘dt
О
Найдем изображение экспоненциальной функции - f(t) = eal:
р - а ° р - а
Изображение экспоненциальной функции поможет нам найти изображения синусоидальной косинусной функций- sin(cof), cos(cof). Для этого запишем эти функции через формулу Эйлера. Далее осуществляем следующую цепочку преобразований:
sin( cat) =
2J
cos( cot)
if I I I _ I \ p + jeo - p + jeo co
2j\p-jco /?+ jeo] 2j\ p2 +co2 J р2+ш2'
Определим изображение производной функции /(f), имеющей dt
изображение F(p)
1 ... 11 12 13 14 15 16 17 18 ... 21
t =
0
0 001
0.002
0 003
0.004
0 005
0 006
0 007
0.008
0 009
0.01001
0.01101
0.01201
0.01301
001401
0.01501
П) =
-0.156
-0.206
-0.227
-0.219
-0.183
-0.124
-0 049
0.036
012
0.196
0.255
0.291
03
0.282
0.237
0.171
Рис. 4.36
§4.4 Переходные процессы в цепи второго порядка
Рассмотрим цепь второго порядка представленную на рис. 4.37 с параметрами: Е = 50В, R - 1 ООм, L = 0.1 Гн, С - 40мкФ.
Записываем уравнения по второму закону Кирхгофа, в результате получаем систему дифференциальных уравнений:
+ис+/•_/? =£* —> L — + Ri + uc=E. i = C—;
Первое слагаемое это исв = 4 -ехр(/>|/) + А2 ехр(р2Осв°бодная составляющая. Она зависит только от параметров схемы, а также от начальных и конечных запасов энергии. Эта составляющая решения не зависит от формы воздействующего напряжения.
Второе слагаемое это ww/, = uc(oo) принуждённая составляющая. Она зависит от внешнего воздействия и имеет форму этого воздействия. Очевидно, что в нашем случае она определяется как ulip = i/c(oo) - Е.
Постоянные интегрирования определяются из начальных условий, отражающих невозможность мгновенного изменения начальных запасов энергии в конденсаторе и в катушке.
Для определения констант интегрирования используем независимые начальные условия ис(0) = 0, iL(0) = 0.
м^(0) — 0 — /4| + А2 + Е\ / к zA(O)-f^(O+)-C^-O = C'(^lp1 + Л2р,).
dt
Откуда следует, что
. (4)
P2
Pl P2
Pl
Теперь можно записать окончательное решение
иг(0 = —ехр(р/)+ Р,£ exp(p2t)+E = ———(-р2ер1 + р,еА' )■
Р2-Р1 Р2-Р1 “ Р2
Р\
Определим корни характеристического уравнения входящие в решение uc\t) Р\^ Р2 через входное сопротивление схемы.
1 CL • р
+ RC• р +1 А 2 п/> 1 Л /с\
pL + н R — — 0 —> CL - р 4- RC • р 4" 1 — 0. (5)
Ср Ср
В результате решения уравнения получаются корни:
-h±4o -RC±J(RC)2-4LC R If R f Г
1,2 2а 2CL 2L LC (6)
г- о 1
Где о = — показатель затухания контура, соо - - .— угловая часто-
та незатухающих колебаний, при выполнении условия (Oq > б2 имеем
Здесь соГ6, частота свободных колебаний,
Корни уравнения определяются параметрами цепи и могут принимать следующие возможные значения (рис. 4.38).
Дискриминант равен нулю. Кони вещественные, отрицательные и тг с R
кратные. Критический режим р^ 2 =“'О = - —
uc(/) = £(l + 8/)e"5' + Е.
• Дискриминант положительный. Корни вещественные отрицательные и неравные. Апериодический режим р} 2 = -6± Jb2 -со^ ;
Uc (0 = —/ Е (Р\ еР:‘ - Р-> qP}> ) + С.
• Дискриминант отрицательный. Корни комплексно-сопряжённые, с отрицательной вещественной частью. Колебательный режим
=-d±j(0„.
uc(t)-Ee-6'
cos(w,,60 + Sin(wrsr) + E .
p1=p2=-6
+1
“ ^св
Puc. 4.38. Расположение корней на комплексной плоскости.
Примеры определения корней характеристического уравнения в
Mathcad
-456.2341361360570100)
-182.6547527528318788^
R 10 С:=60-10 6 L:=0.2
(R+Lp)-2R I , p : + + R solve ,p —►
R+Lp+2R Cp
(R+ Lp)-2 R J_
R + E-р + 2-R {"'’P
2 2
5 R С р + 3 R Ср L + 3R+ Lp
(3-R + Lp)Cp
Р“
-456.234")
-182.655)
С- 100 ю 6
L:=0.1
Л (R+LP)R I
P := +
R+LP+R CP
R solve, Р —>
(-275) - 156.12494995995995515
.(-275) + 156.1249499599599551Н
Р-
-275- 156.125
-275 + 156.125iJ
(R+LP)R I
+ + R
R+LP+R CP
3R2CP + 2R-C P2L+ 2-R + LP
(2R+ LP)CP
Примеры определения корней характеристического уравнения и
зависимых и независимых начальных условий
Пример: Определить независимые/} (0),(/с(0) и зависимые начальные
условия (/£(0+),/с(0+). Определить корень характеристического урав-
22Я+Д
2
(УС(О) = /£(О)Я = -.
1 .Определяем независимые начальные условия /}(0)5С/с(0).
2.
Определяем зависимые начальные условия (7Z (0+),/’с(0+) из схемы после коммутации (см. документ Mathcad).
3.
Определяем корень характеристического уравнения из схемы после коммутации (см. документ Mathcad).
Документ Mathcad
ORIGIN:= 1
АААЛАЛААЛЛАА
Определи гь напряжение па конденсаторе.
Е?= 100 R 10 L:- 0.1 Г;-50106
Е 2R(Lp + R) I
UnD := - Z р) := — + R +
пр 3 З-R + bp Ср
Корни характеристического уравнения
solve,р Г(—416.62)-162.451'
р •- Z(p)
float, 5
(-416.61) + 162.45i_
' -416.67- 162.45^
ч -416.67+ 162.451J
Независимые начальные условия Е п Е 'Lo 1R Ьс0‘ 2
Зависимые начальные условия
Е“ иС0+ *Lo*r ’Lo"5 UC0"50 ’R0“
iCO:",RO",Lc
Пос гоя иные интегрирования
- I ( иС0" ипр
<АП
Л)
р.ЗЗЗ-81.2211^
t 8.333+ 8l.22liJ
со - 162.45 Unp - 33.333
ч । y 1
*Lo *np
UL0 ।
< L J
цсо+ 'LoR 2R
UL0^ F“ ‘R0R
pit pj-t
i(t):=Are + Aye ‘
t = 3.O31x 10 3 t:-0,0.0l-T..T-IC
I
|Rc(P)|
-0.417 + 0.716i>
-0.417-0.7l6iJ
§4.5 Операторный метол расчёта переходных процессов
Операторный метод (преобразование Лапласа) расчета переходных процессов используется для того, чтобы обыкновенные дифференциальные уравнения с постоянными коэффициентами (в пространстве оригиналов) преобразовать в алгебраические (в пространстве изображений). Очевидно, что алгебраические уравнения решаются проще. После решения алгебраического уравнения над полученной функцией (изображением) производится обратное преобразование Лапласа, получается оригинал. Полученный оригинал это функция, которая и буде1 решением дифференциального уравнения.
Любой функции можно сопоставить её преобразование Лапласа
F{p)^]f{t)ep,dt, (7)
О
здесь F(p) изображение, /'(/) оригинал. Выражении (7) записывают ещё и в операторной формеF(p) = L[/(/)] •
Приведём изображение нескольких часто встречающихся функций. Определим изображение константы - /’(/) = A (const):
00 j
F(p) = A\e-t>'dt=-— |;=-.
о Р Р
F(p) = f e(l'ep‘dt
О
Найдем изображение экспоненциальной функции - f(t) = eal:
р - а ° р - а
Изображение экспоненциальной функции поможет нам найти изображения синусоидальной косинусной функций- sin(cof), cos(cof). Для этого запишем эти функции через формулу Эйлера. Далее осуществляем следующую цепочку преобразований:
sin( cat) =
2J
cos( cot)
if I I I _ I \ p + jeo - p + jeo co
2j\p-jco /?+ jeo] 2j\ p2 +co2 J р2+ш2'
Определим изображение производной функции /(f), имеющей dt
изображение F(p)
1 ... 11 12 13 14 15 16 17 18 ... 21
t =
0
0 001
0.002
0 003
0.004
0 005
0 006
0 007
0.008
0 009
0.01001
0.01101
0.01201
0.01301
001401
0.01501
П) =
-0.156
-0.206
-0.227
-0.219
-0.183
-0.124
-0 049
0.036
012
0.196
0.255
0.291
03
0.282
0.237
0.171
Рис. 4.36
§4.4 Переходные процессы в цепи второго порядка
Рассмотрим цепь второго порядка представленную на рис. 4.37 с параметрами: Е = 50В, R - 1 ООм, L = 0.1 Гн, С - 40мкФ.
Записываем уравнения по второму закону Кирхгофа, в результате получаем систему дифференциальных уравнений:
+ис+/•_/? =£* —> L — + Ri + uc=E. i = C—;
Первое слагаемое это исв = 4 -ехр(/>|/) + А2 ехр(р2Осв°бодная составляющая. Она зависит только от параметров схемы, а также от начальных и конечных запасов энергии. Эта составляющая решения не зависит от формы воздействующего напряжения.
Второе слагаемое это ww/, = uc(oo) принуждённая составляющая. Она зависит от внешнего воздействия и имеет форму этого воздействия. Очевидно, что в нашем случае она определяется как ulip = i/c(oo) - Е.
Постоянные интегрирования определяются из начальных условий, отражающих невозможность мгновенного изменения начальных запасов энергии в конденсаторе и в катушке.
Для определения констант интегрирования используем независимые начальные условия ис(0) = 0, iL(0) = 0.
м^(0) — 0 — /4| + А2 + Е\ / к zA(O)-f^(O+)-C^-O = C'(^lp1 + Л2р,).
dt
Откуда следует, что
. (4)
P2
Pl P2
Pl
Теперь можно записать окончательное решение
иг(0 = —ехр(р/)+ Р,£ exp(p2t)+E = ———(-р2ер1 + р,еА' )■
Р2-Р1 Р2-Р1 “ Р2
Р\
Определим корни характеристического уравнения входящие в решение uc\t) Р\^ Р2 через входное сопротивление схемы.
1 CL • р
+ RC• р +1 А 2 п/> 1 Л /с\
t =
0
0 001
0.002
0 003
0.004
0 005
0 006
0 007
0.008
0 009
0.01001
0.01101
0.01201
0.01301
001401
0.01501
П) =
-0.156
-0.206
-0.227
-0.219
-0.183
-0.124
-0 049
0.036
012
0.196
0.255
0.291
03
0.282
0.237
0.171
Рис. 4.36
§4.4 Переходные процессы в цепи второго порядка
Рассмотрим цепь второго порядка представленную на рис. 4.37 с параметрами: Е = 50В, R - 1 ООм, L = 0.1 Гн, С - 40мкФ.
Записываем уравнения по второму закону Кирхгофа, в результате получаем систему дифференциальных уравнений:
+ис+/•_/? =£* —> L — + Ri + uc=E. i = C—;
Первое слагаемое это исв = 4 -ехр(/>|/) + А2 ехр(р2Осв°бодная составляющая. Она зависит только от параметров схемы, а также от начальных и конечных запасов энергии. Эта составляющая решения не зависит от формы воздействующего напряжения.
Второе слагаемое это ww/, = uc(oo) принуждённая составляющая. Она зависит от внешнего воздействия и имеет форму этого воздействия. Очевидно, что в нашем случае она определяется как ulip = i/c(oo) - Е.
Постоянные интегрирования определяются из начальных условий, отражающих невозможность мгновенного изменения начальных запасов энергии в конденсаторе и в катушке.
Для определения констант интегрирования используем независимые начальные условия ис(0) = 0, iL(0) = 0.
м^(0) — 0 — /4| + А2 + Е\ / к zA(O)-f^(O+)-C^-O = C'(^lp1 + Л2р,).
dt
Откуда следует, что
. (4)
P2
Pl P2
Plt =
0
0 001
0.002
0 003
0.004
0 005
0 006
0 007
0.008
0 009
0.01001
0.01101
0.01201
0.01301
001401
0.01501
П) =
-0.156
-0.206
-0.227
-0.219
-0.183
-0.124
-0 049
0.036
012
0.196
0.255
0.291
03
0.282
0.237
0.171
Рис. 4.36
§4.4 Переходные процессы в цепи второго порядка
Рассмотрим цепь второго порядка представленную на рис. 4.37 с параметрами: Е = 50В, R - 1 ООм, L = 0.1 Гн, С - 40мкФ.
Записываем уравнения по второму закону Кирхгофа, в результате получаем систему дифференциальных уравнений:
+ис+/•_/? =£* —> L — + Ri + uc=E. i = C—;
Первое слагаемое это исв = 4 -ехр(/>|/) + А2 ехр(р2Осв°бодная составляющая. Она зависит только от параметров схемы, а также от начальных и конечных запасов энергии. Эта составляющая решения не зависит от формы воздействующего напряжения.
Второе слагаемое это ww/, = uc(oo) принуждённая составляющая. Она зависит от внешнего воздействия и имеет форму этого воздействия. Очевидно, что в нашем случае она определяется как ulip = i/c(oo) - Е.
Постоянные интегрирования определяются из начальных условий, отражающих невозможность мгновенного изменения начальных запасов энергии в конденсаторе и в катушке.
Для определения констант интегрирования используем независимые начальные условия ис(0) = 0, iL(0) = 0.
м^(0) — 0 — /4| + А2 + Е\ / к zA(O)-f^(O+)-C^-O = C'(^lp1 + Л2р,).
dt
Откуда следует, что
. (4)
P2
Теперь можно записать окончательное решение
иг(0 = —ехр(р/)+ Р,£ exp(p2t)+E = ———(-р2ер1 + р,еА' )■
Р2-Р1 Р2-Р1 “ Р2
pL + н R — — 0 —> CL - р 4- RC • р 4" 1 — 0. (5)
Ср Ср
В результате решения уравнения получаются корни:
-h±4o -RC±J(RC)2-4LC R If R f Г
1,2 2а 2CL 2L LC (6)
г- о 1
Где о = — показатель затухания контура, соо - - .— угловая часто-
та незатухающих колебаний, при выполнении условия (Oq > б2 имеем
Здесь соГ6, частота свободных колебаний,
Корни уравнения определяются параметрами цепи и могут принимать следующие возможные значения (рис. 4.38).
Дискриминант равен нулю. Кони вещественные, отрицательные и тг с R
кратные. Критический режим р^ 2 =“'О = - —
uc(/) = £(l + 8/)e"5' + Е.
• Дискриминант положительный. Корни вещественные отрицательные и неравные. Апериодический режим р} 2 = -6± Jb2 -со^ ;
Uc (0 = —/ Е (Р\ еР:‘ - Р-> qP}> ) + С.
• Дискриминант отрицательный. Корни комплексно-сопряжённые, с отрицательной вещественной частью. Колебательный режим
=-d±j(0„.
uc(t)-Ee-6'
cos(w,,60 + Sin(wrsr) + E .
p1=p2=-6
+1
“ ^св
Puc. 4.38. Расположение корней на комплексной плоскости.
Примеры определения корней характеристического уравнения в
Mathcad
-456.2341361360570100)
-182.6547527528318788^
R 10 С:=60-10 6 L:=0.2
(R+Lp)-2R I , p : + + R solve ,p —►
R+Lp+2R Cp
(R+ Lp)-2 R J_
R + E-р + 2-R {"'’P
2 2
5 R С р + 3 R Ср L + 3R+ Lp
(3-R + Lp)Cp
Р“
-456.234")
-182.655)
С- 100 ю 6
L:=0.1
Л (R+LP)R I
P := +
R+LP+R CP
R solve, Р —>
(-275) - 156.12494995995995515
.(-275) + 156.1249499599599551Н
Р-
-275- 156.125
-275 + 156.125iJ
(R+LP)R I
+ + R
R+LP+R CP
3R2CP + 2R-C P2L+ 2-R + LP
(2R+ LP)CP
Примеры определения корней характеристического уравнения и
зависимых и независимых начальных условий
Пример: Определить независимые/} (0),(/с(0) и зависимые начальные
условия (/£(0+),/с(0+). Определить корень характеристического урав-
22Я+Д
2
(УС(О) = /£(О)Я = -.
1 .Определяем независимые начальные условия /}(0)5С/с(0).
2.
Определяем зависимые начальные условия (7Z (0+),/’с(0+) из схемы после коммутации (см. документ Mathcad).
3.
Определяем корень характеристического уравнения из схемы после коммутации (см. документ Mathcad).
Документ Mathcad
ORIGIN:= 1
АААЛАЛААЛЛАА
Определи гь напряжение па конденсаторе.
Е?= 100 R 10 L:- 0.1 Г;-50106
Е 2R(Lp + R) I
UnD := - Z р) := — + R +
пр 3 З-R + bp Ср
Корни характеристического уравнения
solve,р Г(—416.62)-162.451'
р •- Z(p)
float, 5
(-416.61) + 162.45i_
' -416.67- 162.45^
ч -416.67+ 162.451J
Независимые начальные условия Е п Е 'Lo 1R Ьс0‘ 2
Зависимые начальные условия
Е“ иС0+ *Lo*r ’Lo"5 UC0"50 ’R0“
- I ( иС0" ипр
<АП
Л)
р.ЗЗЗ-81.2211^
t 8.333+ 8l.22liJ
со - 162.45 Unp - 33.333
| ч । y 1 | *Lo *np |
| | UL0 । < L J |
цсо+ 'LoR 2R
UL0^ F“ ‘R0R
pit pj-t
i(t):=Are + Aye ‘
t = 3.O31x 10 3 t:-0,0.0l-T..T-IC
I
|Rc(P)|
-0.417 + 0.716i>
-0.417-0.7l6iJ
§4.5 Операторный метол расчёта переходных процессов
Операторный метод (преобразование Лапласа) расчета переходных процессов используется для того, чтобы обыкновенные дифференциальные уравнения с постоянными коэффициентами (в пространстве оригиналов) преобразовать в алгебраические (в пространстве изображений). Очевидно, что алгебраические уравнения решаются проще. После решения алгебраического уравнения над полученной функцией (изображением) производится обратное преобразование Лапласа, получается оригинал. Полученный оригинал это функция, которая и буде1 решением дифференциального уравнения.
Любой функции можно сопоставить её преобразование Лапласа
F{p)^]f{t)ep,dt, (7)
О
здесь F(p) изображение, /'(/) оригинал. Выражении (7) записывают ещё и в операторной формеF(p) = L[/(/)] •
Приведём изображение нескольких часто встречающихся функций. Определим изображение константы - /’(/) = A (const):
00 j
F(p) = A\e-t>'dt=-— |;=-.
о Р Р
F(p) = f e(l'ep‘dt
О
Найдем изображение экспоненциальной функции - f(t) = eal:
р - а ° р - а
Изображение экспоненциальной функции поможет нам найти изображения синусоидальной косинусной функций- sin(cof), cos(cof). Для этого запишем эти функции через формулу Эйлера. Далее осуществляем следующую цепочку преобразований:
sin( cat) =
2J
cos( cot)
if I I I _ I \ p + jeo - p + jeo co
2j\p-jco /?+ jeo] 2j\ p2 +co2 J р2+ш2'
Определим изображение производной функции /(f), имеющей dt
изображение F(p)
ОС 7 / V СО 00
j ^c-^dl -J-.1\1)е
р' f +pj f{l)e
p,dt = -/(0) + pF(p).
0^0 0
И, наконец, определим изображение интегрального выражения j/'(/)J/ о
> e'-j/M*' j/Ж'Л
J //«■>*' е-А --JI ./(е-) — |>л—
око ) ” О \ О 7 ” !' ”
Таблица преобразований Лапласа
/(г)-оригинал
F(p)- изображение
1
l/P
еа>
V(p-a)
ёш
l/(p + a)
sin (со/)
оУ(/?-со2)
cos(cof)
p/(p2+(O2)
-/(0) + pF(p)
0
F(p)
P
Вернемся теперь к переходным процессам.
Итак, мы будем сопоставлять каждой функции его изображение. Например /(/)—>/(/>), u(t) —>U(p). С учётом полученной таблицы можно сопоставить каждому элементу его изображение:
u£(f) = I^y^->l(p/(p)-/£(O))^-pL I(p)- —J J1'1—> i?P“
dt F ifojL
-II >’ ню-”
с о Р PC с i/cp ис(оур
схема изображений, в пространстве изображений находятся желаемые токи и напряжения с использованием известных методов расчета (МКирхгофа, МУП, МКТ и т.д.). Для перехода от изображения к оригиналу (к временной зависимости) необходимо использовать теорему
разложения:
М(р)
W)
где /?к корни уравнения N(p) - 0.
p-N(p) МО)
где рк корни уравнения N(p) = 0.
Пример: Определить ток источника напряжения если
Е = 50В,/? = 10Ом,Л = 0,4Гн.
Рис. 4.49
-
Определим независимые начальные условия i£(0) iL(0) = E/ R = 50/\() = 5А.
-
Изображаем операторную схему замещения после коммутации и находим изображение тока
E/p+it(0)L = Е + iL(O)Lp = М(р)
Р 2R + Lp p(2R + Lp) pN(p)'
где: М(р) -Е + iL(0)Lp - 50 + 2р, N(p) -(2R + Lp} -20 +0,4р. Находим корень знаменателя и его производную
JV(p) = 20 + 0,4p=0 —>р = -2/?/£ = -20/0,4 =-50с-1,
N\p) = L = ^4.
-
Для определения оригинала /(/) используем теорему разложения
1(р} MW > ,(0-^(0), MW ср; 50 1 50-100^, p-N(p) А(0) pJV'(p) 20 -50 0,4
-2,5 + 2,5e-50/A.
-
Переходные процессы в электрических цепях при воздействии импульсного напряжения. (Метод пространства состояний)
-
Интегрирующие и дифференцирующие цепи (дифференцирование и интегрирование как операции фильтрации сигналов). Частотные характеристики
-
Интеграл Дюамеля - аналитический метод расчета переходных процессов при импульсном воздействии
L:-0.l2.‘ -Л--0.07
IpI
:=0..N Dl(t,x):-
Puc. 4.52
R:- 10( C:-700 10 6 p —— R-C
-1 F(t) 2
D(t,x):- x + — N:- 10 -4
RC RC —
Puc. 4.53
140 := 1Пt < T,f(t),E(t -1)) t := 0.0.01-T..3 I
Д>1 f(t):= 4-t-l if0
E(t)
RC
Dc(t,x) := Ac x + Bc(t)
al' T
DL(l,x) = AL x + BL(t) х^:= rkfixed(0,0,4 T,N,DiJ
AL^ BLa>:-^
DL(t,x) Al x + BL(t) x^ = rkfixed(0,0,4 T.N.DjJ
Ас:=^
Dc(t,x) Ac x + Bc(t)
xc := ricfixedi0.0,4 T.N,DCI t := xq<0>
W. (jm) =
L j
ПЛ Referenced ‘.Documents and Settings\Yusup\Pa6o4Hii стол\Аним_^Анимации\Уч_Пособ xmcd Rp IC 100- IO- 6 L^:- 0.02 RL1 Of
co.o*- н
arg(z) °’1 deg ci,o*-Rc(z> Cj ]<-lw(z)
Puc. 4.59
<0.928 -2I.8O1A mo - 400 com( a) -
V 0.862 -0.345 )
ДО2 sin|
FXD
f:-5f (o:-2-T[-f сп - 314.159 T — t := 0,0.01 Т..2 Т W (WA
(О
Al A
0.01 \ 0.02 I 0.03 \ 0.04
V/v ч
, . (0.954 -17.441) (0.847 -32.142) , . (0.623-51.48$
conl(a|) ^091 -0.286 J 10.717 -0.45 J [0.388 -0.487J
(0.414-65.54)
[o.!71-O.377J
El(t) :- 2- |a| I sin^t-cu — 9O.deg + argfa|)j + 0.2'|a->| sin(o)'2-t + 30-deg + arg(a^)j ...
+ 0.3-|ад| sin^o-4-t + 45 deg + arg^a^jj + 0.1- |a-y| -sin(nj-7-t - 60 deg + argfa-yjj
Puc. 4.62
Лекция №11
§4.6 Интеграл Дюамеля
g(0
Прежде всего, уместно ввести понятие переходная функция. Переходная функция это отклик системы на единичное воздействие. При известной переходной функции g(t) для заданной схемы можно найти ток в цепи
i(t) = g(t)UQ
Здесь (70 постоянное внешнее воздействие. Для того чтобы
Определить ток при произвольном внешнем воздействии U(t), разобьем функцию С;(г)на прямоугольники как показано на рисунке 4.64. Полный ток в момент t получаем, используя метод наложения. Просуммируем все частичные токи от отдельных скачков и прибавим их к току «(0)g(0:
КО = u(0)g(/) + £u'(T)g(/ -т)Дт
Число членов суммы равно числу ступенек напряжения. Очевидно, что ступенчатая кривая тем лучше заменяет плавную кривую, чем больше число ступенек. С этой целью заменим конечный интервал времени Ат на бесконечно малый di и перейдем от суммы к интегралу:
/(/) - i/(0)g(/) + J u'(z)g(/ - x)di,
о
или /
K0-«(0g(0) + f*K'0g'('-'0 о
Пример:
R IR°
— + RI + L-x solve, x —► -1200 RI + R2
-1200
Uo ?= 200 т :=-r—r
Ipi
t - 8.333x 10 4
0 otherwise
RI- 200 R2-50 L - 0.2
RI + R2 RI + R2 RI
R2 + — о
Л - -0.003
0.007
RI
R2+ —
2
RI -200
R2=50 L=0.2
Находим переходную проводимость i(t):
e(t) .•= 0.0Лер1 + ! т - 8.333x 10“ 4
** R|
Находим ток на первом интервале i(t) 0 < t < т :
ОС 7 / V СО 00
j ^c-^dl -J-.1\1)е
р' f +pj f{l)e
p,dt = -/(0) + pF(p).
0^0 0
И, наконец, определим изображение интегрального выражения j/'(/)J/ о
> e'-j/M*' j/Ж'Л
J //«■>*' е-А --JI ./(е-) — |>л—
око ) ” О \ О 7 ” !' ”
Таблица преобразований Лапласа
/(г)-оригинал
F(p)- изображение
1
l/P
еа>
V(p-a)
ёш
l/(p + a)
sin (со/)
оУ(/?-со2)
cos(cof)
p/(p2+(O2)
-/(0) + pF(p)
0
F(p)
P
Вернемся теперь к переходным процессам.
Итак, мы будем сопоставлять каждой функции его изображение. Например /(/)—>/(/>), u(t) —>U(p). С учётом полученной таблицы можно сопоставить каждому элементу его изображение:
u£(f) = I^y^->l(p/(p)-/£(O))^-pL I(p)- —J J1'1—> i?P“
dt F ifojL
-II >’ ню-”
с о Р PC с i/cp ис(оур
схема изображений, в пространстве изображений находятся желаемые токи и напряжения с использованием известных методов расчета (МКирхгофа, МУП, МКТ и т.д.). Для перехода от изображения к оригиналу (к временной зависимости) необходимо использовать теорему
разложения:
М(р)
W)
где /?к корни уравнения N(p) - 0.
p-N(p) МО)
где рк корни уравнения N(p) = 0.
Пример: Определить ток источника напряжения если
Е = 50В,/? = 10Ом,Л = 0,4Гн.
Рис. 4.49
-
Определим независимые начальные условия i£(0) iL(0) = E/ R = 50/\() = 5А.
-
Изображаем операторную схему замещения после коммутации и находим изображение тока
E/p+it(0)L = Е + iL(O)Lp = М(р)
Р 2R + Lp p(2R + Lp) pN(p)'
где: М(р) -Е + iL(0)Lp - 50 + 2р, N(p) -(2R + Lp} -20 +0,4р. Находим корень знаменателя и его производную
JV(p) = 20 + 0,4p=0 —>р = -2/?/£ = -20/0,4 =-50с-1,
N\p) = L = ^4.
-
Для определения оригинала /(/) используем теорему разложения
1(р} MW > ,(0-^(0), MW ср; 50 1 50-100^, p-N(p) А(0) pJV'(p) 20 -50 0,4
-2,5 + 2,5e-50/A.
-
Переходные процессы в электрических цепях при воздействии импульсного напряжения. (Метод пространства состояний)
-
Интегрирующие и дифференцирующие цепи (дифференцирование и интегрирование как операции фильтрации сигналов). Частотные характеристики
-
Интеграл Дюамеля - аналитический метод расчета переходных процессов при импульсном воздействии
L:-0.l2.‘ -Л--0.07
IpI
:=0..N Dl(t,x):-
Puc. 4.52
R:- 10( C:-700 10 6 p —— R-C
-1 F(t) 2
D(t,x):- x + — N:- 10 -4
RC RC —
Puc. 4.53
140 := 1Пt < T,f(t),E(t -1)) t := 0.0.01-T..3 I
Д>1 f(t):= 4-t-l if0
E(t)
RC
Dc(t,x) := Ac x + Bc(t)
al' T
DL(l,x) = AL x + BL(t) х^:= rkfixed(0,0,4 T,N,DiJ
AL^ BLa>:-^
DL(t,x) Al x + BL(t) x^ = rkfixed(0,0,4 T.N.DjJ
Ас:=^
Dc(t,x) Ac x + Bc(t)
xc := ricfixedi0.0,4 T.N,DCI t := xq<0>
W. (jm) =
L j
ПЛ Referenced ‘.Documents and Settings\Yusup\Pa6o4Hii стол\Аним_^Анимации\Уч_Пособ xmcd Rp IC 100- IO- 6 L^:- 0.02 RL1 Of
co.o*- н
arg(z) °’1 deg ci,o*-Rc(z> Cj ]<-lw(z)
Puc. 4.59
<0.928 -2I.8O1A mo - 400 com( a) -
V 0.862 -0.345 )
ДО2 sin|
FXD
f:-5f (o:-2-T[-f сп - 314.159 T — t := 0,0.01 Т..2 Т W (WA
(О
Al A
0.01 \ 0.02 I 0.03 \ 0.04
V/v ч
, . (0.954 -17.441) (0.847 -32.142) , . (0.623-51.48$
conl(a|) ^091 -0.286 J 10.717 -0.45 J [0.388 -0.487J
(0.414-65.54)
[o.!71-O.377J
El(t) :- 2- |a| I sin^t-cu — 9O.deg + argfa|)j + 0.2'|a->| sin(o)'2-t + 30-deg + arg(a^)j ...
+ 0.3-|ад| sin^o-4-t + 45 deg + arg^a^jj + 0.1- |a-y| -sin(nj-7-t - 60 deg + argfa-yjj
Puc. 4.62
Лекция №11
§4.6 Интеграл Дюамеля
g(0
Прежде всего, уместно ввести понятие переходная функция. Переходная функция это отклик системы на единичное воздействие. При известной переходной функции g(t) для заданной схемы можно найти ток в цепи
i(t) = g(t)UQ
Здесь (70 постоянное внешнее воздействие. Для того чтобы
Определить ток при произвольном внешнем воздействии U(t), разобьем функцию С;(г)на прямоугольники как показано на рисунке 4.64. Полный ток в момент t получаем, используя метод наложения. Просуммируем все частичные токи от отдельных скачков и прибавим их к току «(0)g(0:
КО = u(0)g(/) + £u'(T)g(/ -т)Дт
Число членов суммы равно числу ступенек напряжения. Очевидно, что ступенчатая кривая тем лучше заменяет плавную кривую, чем больше число ступенек. С этой целью заменим конечный интервал времени Ат на бесконечно малый di и перейдем от суммы к интегралу:
/(/) - i/(0)g(/) + J u'(z)g(/ - x)di,
о
или /
K0-«(0g(0) + f*K'0g'('-'0 о
Пример:
R IR°
— + RI + L-x solve, x —► -1200 RI + R2
-1200
Uo ?= 200 т :=-r—r
Ipi
t - 8.333x 10 4
0 otherwise
RI- 200 R2-50 L - 0.2
RI + R2 RI + R2 RI
R2 + — о
Л - -0.003
0.007
RI
R2+ —
2
RI -200
R2=50 L=0.2
Находим переходную проводимость i(t):
e(t) .•= 0.0Лер1 + ! т - 8.333x 10“ 4
** R|
Находим ток на первом интервале i(t) 0 < t < т :
ОС 7 / V СО 00
j ^c-^dl -J-.1\1)е
р' f +pj f{l)e
/(г)-оригинал | F(p)- изображение |
1 | l/P |
еа> | V(p-a) |
ёш | l/(p + a) |
sin (со/) | оУ(/?-со2) |
cos(cof) | p/(p2+(O2) |
| -/(0) + pF(p) |
0 | F(p) P |
где /?к корни уравнения N(p) - 0.
p-N(p) МО)
где рк корни уравнения N(p) = 0.
Пример: Определить ток источника напряжения если
Е = 50В,/? = 10Ом,Л = 0,4Гн.
Рис. 4.49
-
Определим независимые начальные условия i£(0) iL(0) = E/ R = 50/\() = 5А. -
Изображаем операторную схему замещения после коммутации и находим изображение тока
E/p+it(0)L = Е + iL(O)Lp = М(р)
Р 2R + Lp p(2R + Lp) pN(p)'
где: М(р) -Е + iL(0)Lp - 50 + 2р, N(p) -(2R + Lp} -20 +0,4р. Находим корень знаменателя и его производную
JV(p) = 20 + 0,4p=0 —>р = -2/?/£ = -20/0,4 =-50с-1,
N\p) = L = ^4.
-
Для определения оригинала /(/) используем теорему разложения
1(р} MW > ,(0-^(0), MW ср; 50 1 50-100^, p-N(p) А(0) pJV'(p) 20 -50 0,4
-2,5 + 2,5e-50/A.
-
Переходные процессы в электрических цепях при воздействии импульсного напряжения. (Метод пространства состояний) -
Интегрирующие и дифференцирующие цепи (дифференцирование и интегрирование как операции фильтрации сигналов). Частотные характеристики -
Интеграл Дюамеля - аналитический метод расчета переходных процессов при импульсном воздействии
L:-0.l2.‘ -Л--0.07
IpI
:=0..N Dl(t,x):-
Puc. 4.52
R:- 10( C:-700 10 6 p —— R-C
-1 F(t) 2
D(t,x):- x + — N:- 10 -4
RC RC —
140 := 1Пt < T,f(t),E(t -1)) t := 0.0.01-T..3 I
Д>1 f(t):= 4-t-l if0
E(t)
RC
Dc(t,x) := Ac x + Bc(t)
al' T
DL(l,x) = AL x + BL(t) х^:= rkfixed(0,0,4 T,N,DiJ
AL^ BLa>:-^
DL(t,x) Al x + BL(t) x^ = rkfixed(0,0,4 T.N.DjJ
Ас:=^
Dc(t,x) Ac x + Bc(t)
xc := ricfixedi0.0,4 T.N,DCI t := xq<0>
ДО2 sin|
FXD
f:-5f (o:-2-T[-f сп - 314.159 T — t := 0,0.01 Т..2 Т W (WA
(О
0.01 \ 0.02 I 0.03 \ 0.04
V/v ч
fx i(t) = g(t)-U(0) + g(t- | r) — U(t) dr il(t) :-g(0-IJ(0) + g(t - r)-Ud(r)dr |
| dr Jo |
il(z) float,4 -+(-6.00)e< l2°°')Z + 13.33 + 3200,-z
(- PCX) i t
il(t) --б.ООе' ' + 13.33+ 3200.t
N:-5f k:-0.. N At:-— tAtk 1, t- 8.333x10”4
**** jq к к W
§6.1.1. Переходная характеристика (иди переходная функция)
Дельта функция Дирака б(х-х0) и G(x-xQ)-ступенчатая функция
Хевисайда
Свойство дельта функции Дирака:
0>
1 ... 13 14 15 16 17 18 19 20 21