Методическое пособие 370

а) в первом интервале времени от 0 до t1 (не включая,

скачка F1 )

t

б) во втором интервале времени от t1 до t2 (не вклю-

чая, скачка F2 )

t1

здесь слагаемое F1h t t1 обусловлено положительным скач-

ком воздействия в момент t1,

где слагаемое F2h t t2 обусловлено отрицательным скач-

ком воздействия в момент t2 .

Входящее в формулы (1.53)-(1.55) выражение h t

есть h t , в котором t заменено на t .

Пример 1.4. Найти напряжение на выходе цепи (рис.1.27а), если напряжение на ее входных зажимах изменяется во времени по закону

50

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

51

Возможно также определение реакции цепи на произвольное внешнее воздействие по ее импульсной характеристике.

Выражение для этого случая можно получить из формулы (1.52), используя правило интегрирования по частям и учитывая соотношение между переходной и импульсной характеристиками цепи (1.43), т.е.

52

t

Выражение (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

uвых t f1 g t d f2 t h 0 f2 g t d ;

в) в третьем интервале времени t2 t

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

2.1.Решение дифференциальных уравнений

спомощью вычислительного блока Given-Odesolve

Описание функции Odesolve можно найти в [5]. Чтобы реализовать блок Given-Odesolve, необходимо

иметь составленное заранее по электрической цепи в соответствии с правилами (пп.1.2-1.3) дифференциальное уравнение (1.8) или (1.11) и далее выполнить следующую последовательность действий:

1.Ввести вводное слово Given (дано).

2.Ниже вводного слова ввести дифференциальное уравнение типа (1.8) или (1.11). Это дифференциальное уравнение должно быть строго линейным (то есть высшая производная в нем не должна иметь никаких сомножителей или степенных показателей). Подобное ограничение связано с особенностями используемого (по умолчанию) системой при решении дифференциального уравнения в форме вычислительного блока Given-Odesolve метода Рунге-Кутта 4-го порядка точности.

3.Производные в выражении дифференциального уравнения необходимо задать при помощи специальных операторов панели Calculus (вычисления) отдельно для производной первого порядка и производных произвольного порядка.

4.Искомые функции могут быть определены произвольным образом с соблюдением общих правил синтаксиса Mathcad.

5.В качестве знака равенства в дифференциальном уравнении следует использовать логическое равенство (=) из иконки «Boolean».

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

ренциального уравнения. Так для уравнения первого порядка необходимо задать одно начальное условие при t 0, а для уравнения второго порядка еще и значение первой производной при t 0. Оно вводится через соответствующий штрих, который ставится сверху переменной

55

при помощи сочетания [Ctrl] + [F7] и при этом не требуется устанавливать значок ввода степени.

7.Затем следует задействовать специальную функцию odesolove (t, в, [step]), где t – переменная, от которой зависит искомая функция; в–правая единица интервала поиска решения, левая граница которого определяется при задании начальных условий. Обязательное условие – значение

«в» должно лежать правее.

Step – определяет число шагов, используемых численным методом Рунге-Кутта.

Задание step не является обязательным, и по умолчанию этот параметр определяется таким образом, чтобы длина шага равнялась 0,1. Нельзя задавать step равным 1, так как при этом возникает неопределенность. При уменьшении значения step увеличивается время расчетов.

Затем, перед построением графика переходного процесса, необходимо ввести изменение величины t для более полного отражения сущности возникающих отклонений тока или напряжения. Изменения вводятся по образцу:

t : 0, доля шага в .. в.

Приведем примеры решения дифференциальных уравнений первого и второго порядков с использованием блока

Given-Odesolve.

Пример 2.1. Найти решение дифференциального уравнения, составленного в примере 1.1.

Выражение для дифференциального уравнения берем из (1.8), т.е.

RC duc uc E2, dt

а начальное условие uc 0 uc 0 E1.

Теперь зададимся для определенности значениями эле-

56

Последовательность ввода и результат в виде графика (рис. 2.1) при определении переходного процесса будет иметь следующий вид.

t

Рис. 2.1. График переходного процесса

Пример 2.2. Найти решение дифференциального урав-

нения, составленного в пп.1.3. [(1.10), (1.11)].

Выражение для дифференциального уравнения берем из (1.11), т.е.

при t 0.

57