Программные математические комплексы. Практикум. Учебное пособие

8.

9.

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

41

Решить систему линейных уравнений численно, методом Гаусса, матричным методом, по формулам Крамера и с помощью символьных преобразований:

42

7. Решить однородную систему линейных уравнений:

ì2x1 + x2 + 4x3 = 0, 2), 8), 14) ïí3x1 + 5x2 + 7x3 = 0,

ïî4x1 + 5x2 + 6x3 = 0;

ì2x1 + x2 + 5x3 + 7x4 = 0, 3), 9), 15) ïí4x1 + 2x2 + 7x3 + 5x4 = 0,

ïî2x1 + x2 + x3 + 5x4 = 0;

43

ïî2x1 + 5x3 - 3x4 + 4x5 = 0; ì2x1 + 3x2 - x3 + 5x4 = 0,

ï3x - 2x + x = 0,

ï 1 2 4

6), 12), 18) í

ï5x1 + x2 - x3 + 6x4 = 0, ïîx1 - 3x3 - x4 = 0.

Практическая работа № 5. Решение дифференциальных уравнений

Цель работы: изучить способы решения дифференциальных уравнений и систем.

Теоретические сведения

Дифференциальные уравнения - это уравнения, в которых неизвестными являются не переменные(т. е. числа), а функции одной или нескольких переменных. Эти уравнения (или системы) включают соотношения между искомыми функциями и их производными. Если в уравнения входят производные только по одной переменной, то они называются обыкновенными дифференциальными уравнениями (далее чаще используется сокращение ОДУ). В противном случае говорят об уравнениях в частных производных. Таким образом, решить (иногда говорят проинтегрировать) дифференциальное уравнение - значит определить неизвестную функцию на определенном интервале изменения ее переменных. Как известно, одно обыкновенное дифференциальное уравнение или система ОДУ имеет единственное решение, если помимо уравнения определенным образом заданы начальные или граничные условия. В соответствующих курсах высшей

44

математики доказываются теоремы о существовании и единственности решения в зависимости от тех или иных условий.

1. Решение ОДУ с помощью группы решения

Дифференциальное уравнение первого порядка может по определению содержать помимо самой искомой функцииy(t) только ее первую производнуюy'(t). В подавляющем большинстве случаев дифференциальное уравнение можно записать в стандартной форме (форме Коши): у' (t)=f (y(t),t).

Вычислительный блок для решения одного ОДУ, реализующий численный метод Рунге - Кутты, состоит из трех частей:

Given - ключевое слово;

ОДУ и начальное условие, записанное с помощью логических операторов, причем начальное условие должно быть в фор-

ме y(t0)=b;

odesolve(t,t1) - встроенная функция для решения ОДУ относительно переменной t на интервале (t0,t1).

Например:

y := odesolve (t ,10)

Чтобы визуализировать полученную функцию можно -по строить ее график:

45

Обыкновенное дифференциальное уравнение с неизвестной функцией y(t), в которое входят производные этой функции вплоть до y(N)(t), называется ОДУ N-го порядка. Если имеется та-

кое уравнение, то для корректной постановки задачи Коши требуется задать N начальных условий на саму функциюy(t) и ее производные от первого до(N - l)-ro порядка включительно. В MathCAD можно решать ОДУ высших порядков как с помощью вычислительного блока Given/odesolve, так и путем сведения их к системам уравнений первого порядка. Внутри вычислительного блока:

-ОДУ должно быть линейно относительно старшей производной, т. е. фактически должно быть поставлено в стандартной форме;

-начальные условия должны иметь формуy(t)=b или

у(t)=b.

46

Например:

2. Решение систем дифференциальных уравнений

Систему дифференциальных уравнений можно решить, используя тот же блок Given/odesolve. При этом количество начальных условий равно сумме порядков всех уравнений, входящих в систему.

3. Решение дифференциальных уравнений без помощи блока решений

MathCAD содержит следующие встроенные функции для решения дифференциальных уравнений:

- rkfixed – метод Рунге - Кутта с фиксированным шагом интегрирования;

47

-Rkadapt – метод Рунге - Кутта с переменным шагом интегрирования (величина шага адаптируется к скорости изменения функции решения);

-Bulstoer – метод Булирша - Штера (более эффективен в случае если решение является плавной функцией).

Эти функции имеют сходные параметры (y, a, b, n, D):

-y – вектор начальных значений неизвестных функций, входящих в систему (в случае одного уравнений и одной неизвестной функции – это просто число);

-a – начало отрезка, на котором ищется решение (именно в

этой точке находятся значения вектора );

-b – конец отрезка интегрирования;

-n – количество частей, на которые разбивается отрезок интегрирования;

-D(x, y) – векторная функция, элементы которой содержат правые части уравнений системы в нормальной записи(когда левые части – производные от соответствующих функций, а в правых частях производные отсутствуют).

Например:

æ 1 ö y := ç ÷

è 0 ø

48

Рассмотрим решение системы дифференциальных уравнений на примере системы ВольтерраЛотки. Эта система описывает динамику численности хищников и жертв на замкнутой территории и задается следующими уравнениями:

где N1 – количество хищников;

N2 – количество жертв;

e1 , e 2 , g1 , g 2 - параметры системы.

Запишем данные уравнения, параметры системы и начальные значения количества хищников и жертв в группу решения:

49

ORIGIN:= 1

Z := Rkadapt( N, 0, t , n , D)

k := 1.. n

Контрольные задания

1. Решить дифференциальное уравнение с помощью функций odesolve, Rkadapt, Bulstoer. Для сравнения решений построить все 3 графика полученных функций на одной координатной плоскости.

1.

3.

4.

5.

6.

7.

8.

9.

50