а) в первом интервале времени от 0 до t1 (не включая,
скачка F1 )
t
u t |
|
(1.53) |
|
f1 0 h t f1 |
h t d , |
||
|
0 |
|
|
б) во втором интервале времени от t1 до t2 (не вклю-
чая, скачка F2 )
t1
|
|
|
u t f1 0 h t f1 |
h t dt F1h t t1 |
|
0 |
|
|
t |
|
|
f2 h t d , |
(1.54) |
|
t1 |
|
|
здесь слагаемое F1h t t1 обусловлено положительным скач-
ком воздействия в момент t1,
в) в третьем интервале времени от |
t2 до |
|
с учетом |
|
(1.53) и (1.54) |
|
|
|
|
t1 |
|
t2 |
|
|
|
f |
|
h t d |
|
u t f1 0 h t f1 |
h t d F1h t t1 |
2 |
||
0 |
|
t1 |
|
|
t |
|
|
|
|
F2h t t2 f3 h t d , |
|
|
(1.55) |
|
t2 |
|
|
|
|
где слагаемое F2h t t2 обусловлено отрицательным скач-
ком воздействия в момент t2 .
Входящее в формулы (1.53)-(1.55) выражение h t
есть h t , в котором t заменено на t .
Пример 1.4. Найти напряжение на выходе цепи (рис.1.27а), если напряжение на ее входных зажимах изменяется во времени по закону
50
|
|
0 |
приt 0 |
u t |
U e t |
при0 t t |
|
1 |
|
|
1 |
|
|
0 |
приt t , |
|
|
|
1 |
Rt
а переходная характеристика цепи h t e L .
График входного воздействия u1 t на электрическую цепь показан на рис. 1.27б.
Рис. 1.27. RL-цепь (а) и входное воздействие (б)
Решение. В интервале времени 0 t t1 воспользуемся формулой интеграла Дюамеля (1.53)
t
f t f1 0 h t f1 h t d .
0
Подставляя соответствующие значения, получаем, что
t
u2 t u1 0 h t du1 h t d d
|
|
|
|
0 |
|
|
|
|
|
R |
t |
t |
d Uea |
|
R t /L |
|
|
|
|
|
|
|||||
U e L |
|
e |
|
d |
||||
d |
|
|
||||||
|
|
|
0 |
|
|
|
|
|
51
|
|
|
|
|
|
|
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
Rt |
|
|
t |
R |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
U e L U a e L e |
|
L d |
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Rt |
|
|
|
R |
|
t |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
t |
|
U a e L |
|
|
a |
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
U e L |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
e |
L |
|
|
|
|
||||||||||||||||||||
|
|
a R L |
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Rt |
|
|
R |
|
|
|
|
|
|
|
|
|
|
Rt |
|
|
|
|||||||||||
|
R |
t |
|
U a e L |
a |
|
|
t |
|
|
|
|
|
aeat Re L /L |
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
Ue L |
|
|
|
|
|
|
|
|
|
|
e |
L |
1 |
|
|
|
|
|
|
|
U. |
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
a R/L |
|
||||||||||||||||||||||||||||
|
|
|
|
|
a |
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
В интервале времени t t1 |
|
|
получаем с учетом скачка в |
||||||||||||||||||||||||||||||||||||||
момент t1, что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
u t u 0 h t t1 |
du1 |
h t d u t h t t |
|||||||||||||||||||||||||||||||||||||||
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
2 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
|
|
|
1 |
|
1 |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
Rt |
|
|
|
|
|
|
|
|
|
|
|
|
|
Rt |
|
|
e a R/L t1 |
1 Ueat1 e R t t1 L |
||||||||||||||||||
|
Ue |
|
U a e L |
|
||||||||||||||||||||||||||||||||||||||
|
L |
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a R/ L |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
U |
|
|
R |
e Rt /L 1 e a R/L t1 . |
|
|
|||||||||||||||||||||||||||||||||
|
a R/L |
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
Результирующий отклик u2 t определяется последова- |
|||||||||||||||||||||||||||||||||||||||||
тельной |
|
стыковкой |
выходного |
|
|
воздействия на |
интервале |
|||||||||||||||||||||||||||||||||||
0 t t1 |
и интервале t t1. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
Возможно также определение реакции цепи на произвольное внешнее воздействие по ее импульсной характеристике.
Выражение для этого случая можно получить из формулы (1.52), используя правило интегрирования по частям и учитывая соотношение между переходной и импульсной характеристиками цепи (1.43), т.е.
52
t
uвых t u t h 0 u g t d , |
(1.56) |
|
где h 0 |
0 |
|
- значение h t цепи при t 0. |
|
|
Выражение (1.56) можно использовать для определения реакции цепи и в том случае, когда внешнее воздействие на цепь описывается кусочно-непрерывной функцией. При этом интервал интегрирования разбивается на несколько промежутков в соответствии с интервалами непрерывности входного сигнала u t .
Так, если на цепь воздействует колебание сложной формы f t (рис. 1.26), то в случае применения импульсной характеристики, отклик цепи uвых t необходимо определять с учетом различных интервалов времени:
а) в первом интервале времени 0 t t1
t
uвых t f1 t h 0 f1 g t d ;
0
б) во втором интервале времени t1 t t2
t1 |
t |
uвых t f1 g t d f2 t h 0 f2 g t d ;
0 |
t1 |
в) в третьем интервале времени t2 t
t1 |
t2 |
uвых t f1 g t d f2 g t d f3 t h 0
0 t1
t
f3 g t d .
t2
53
2. ПРИМЕНЕНИЕ ЭВМ ДЛЯ РАСЧЕТА ПЕРЕХОДНЫХ ПРОЦЕССОВ В ЛИНЕЙНЫХ ЦЕПЯХ
Для расчета переходных процессов в линейных электрических цепях достаточно успешно может быть использован программный продукт компьютерной математики Mathcad
последних версий: Mathcad 8/2000, Mathcad 2001 [3, 4].
Применение пакета Mathcad позволяет, используя стандартные программы, получить численное решение дифференциального уравнения цепи с определением переходного процесса и построения его графика.
Здесь имеются в виду так называемые обыкновенные дифференциальные уравнения (ОДУ). Это уравнения, в которые входят производные только по одной переменной.
Mathcad, к сожалению, не позволяет непосредственно находить аналитические решения таких дифференциальных уравнений, однако реально упрощает выполнение различных вычислительных этапов при поиске указанных решений традиционными методами .
Переходные процессы, отражаемые ОДУ, обычно описывают изменение анализируемой характеристики с течением времени, а начальные условия должны определяться в момент времени t 0.
Заданные таким образом дифференциальные уравнения называются задачами Коши [2].
Методы численного решения ОДУ в форме задачи Коши разработаны весьма досконально. Самыми популярными из них заслуженно являются алгоритмы Рунге-Кутта [3], успешно используемые для решения подавляющего большинства дифференциальных уравнений.
Решение одного линейного дифференциального уравнения в Mathcad может быть реализовано в двух формах: в виде вычислительного блока или в виде встроенных функций. Первая форма предпочтительнее с точки зрения наглядности представления решения и технической простоты.
54