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

В левом маркере этого оператора должно быть определено преобразуемое выражение без знака равенства, в правом – переменная, по которой осуществляется преобразование:

Далее необходимо упростить это громоздкое выражение. Для этого переменную laplace i t ,t, s заменим на Z . Начальное условие для первой производной

После проведенных замен уравнение следует разрешить относительно переменной Z, т.е. использовать специальный оператор символьного решения уравнений solve панели Symbolic

При таких начальных условиях после применения оператора solve получаем, что

61

Теперь, чтобы получить решение дифференциального уравнения, к полученному решению-изображению нужно применить операцию обратного преобразования Лапласа. Сделать это можно при помощи оператора invlaplace панели Symbolic (Символьные), т.е.

При расчете обратного преобразования Лапласа желательно одновременно провести и упрощение получаемого громоздкого выражения, что можно сделать, задействовав одновременно с оператором invlaplace и специальный оператор Simplify (Упростить).

Найденная функция i t после применения обратного преобразования Лапласа будет иметь следующий вид

Такое решение дифференциального уравнения является пусть

изаписанным не в самой удобной форме, но полностью соответствует поставленной задаче и дает такие же результаты как

ипри применении процедуры Given odesolve, как это было сделано в пп.1.5.1 для решения дифференциального уравнения второго порядка.

После соответствующего ввода численных значений будем иметь:

E : 100 L : 0.05 C : 5 10 9 R :LC 0,5

62

t : 0, 0.000001.. 0.0002

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

Этот результат (рис.2.5) соответствует условию R 2 , где

CL .

63

2.3. Применение процедур численного дифференцирования и интегрирования к вычислениям методом интеграла Дюамеля

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

Приступая к вычислению производной какой-либо функции, следует найти соответствующий оператор. Оператор производной расположен на панели Calculus (Вычисления) и, поми-

Если в результате дифференцирования должна быть получена функция производной, следует обратиться к возможностям символьного вычисления. В качестве оператора вывода следует использовать оператор символьного вывода («→» – Evaluate Symbolically).

При символьном дифференцировании можно оперировать функциями нескольких переменных и функциями с параметрами. Оператор дифференцирования может сочетаться с любым вычислительным оператором.

Пример 2.7. Символьное дифференцирование

64

Пример 2.8. Подсчет численного значения при символьном дифференцировании

x :

Аналогично дифференцированию, интегрирование может быть как численным, так и символьным. Роль численного алгоритма выше в случае интегрирования, чем при вычислении производной. Это связано с тем, что многие функции не имеют первообразной (производную можно вычислить практически для любой функции).

Для ввода определенного интеграла следует использовать оператор панели Calculus (Вычисления) или через клавиши

[Shift] + [7]:

d

Оператор содержит четыре маркера, которые заполняются в соответствии с принятой в математике формой. Предварительно необходимо присвоить всем параметрам и переменным конкретные значения. Пределы интегрирования должны быть числами. В качестве оператора вывода следует использовать оператор численного вывода (=).

В систему Mathcad встроено несколько численных методов интегрирования. Необходимость этого связана с тем, что алгоритмы численного интегрирования не столь универсальны, как алгоритм дифференцирования. Каждый из них подходит для определенной группы функций. Так, например, такие классические методы, как метод трапеций, средних прямоугольников или Симпсона подходят лишь для гладких непрерывных функций и для небольшого интервала интегрирования. Чтобы подсчитать интеграл с бесконечным пределом или интеграл от функции, имеющей точку разрыва второго рода на границе промежутка, нужно использовать другие методы.

65

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

Лучше всего, если по умолчанию отмечен первый вариант Autoselect (Автоматический выбор). Метод интегрирования выбирается при этом автоматически. Однако, в некоторых ситуациях (прежде всего при наличии точек разрыва) система может не справиться с расчетами. В этом случае следует попробовать сменить алгоритм подсчета, сделав выбор из следующих четырех методов:

1.Romberg. Метод, применяемый для вычисления интегралов от функций, не имеющих особенностей.

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

3.Infinite Limit (Бесконечный предел). Предназначен для расчета интегралов с бесконечными пределами интегрирования. Самостоятельно выбирать данный метод не придется, поскольку система автоматически переключается на него при введении в оператор интегрирования символа бесконечности.

4.Singular Endpoint (Неопределенность в точке предела интегрирования). Данный метод предназначен для вычисления интегралов от функций, не существующих в одной или в обеих точках пределов интегрирования (так называемые несобственные интегралы второго рода).

Точность результатов интегрирования можно регулировать заданием оператора TOL перед интегралом.

По умолчанию TOL:= 10-3. Можно задавать любую величину

до 10-13.

Так, например, точное аналитическое решение интеграла

0

а при обычной точности (TOL = 10-3)

66

e x sin x dx 0,499999999986754,

0

но при TOL:=10-13 имеем, что

e x sin x dx 0.5.

0

Приведем теперь пример использования численных методов при определении отклика цепи на входное воздействие с помощью интеграла Дюамеля.

Воспользуемся условиями задачи из примера 1.4, когда на

Пример 2.9. Зададимся конкретными значениями величин,

а для интервала t t1 будем иметь отклик

Вычисляем результирующий отклик с учетом двух интервалов и строим график (рис.2.6), т.е.

u3 t : if t t1, u1 t , u2 t

67

t : 0, 0.0001..0.008

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

Вычисление отклика по формулам, полученным в примере 1.4, осуществляется в виде

(рис.2.6), рассчитанным численным методом.

68

u3 t: if t t1, u1 t , u2 t

Рис.2.9. График результата расчета с h t

В случае применения для расчетов интеграла Дюамеля на основе импульсной характеристики (1.56), переходной процесс по примеру 1.4. также вычисляется по частям на интервалах 0 t t1 и t t1 . Расчеты с применением системы Mathcad и график (рис.2.10) для этого случая показаны в листинге 2.2.

70