Решение дифференциальных уравнений методом эйлера реферат

Решение дифференциальных уравнений 1 порядка методом Эйлера — реферат

1. Постановка задачи 5

2. Обзор существующих методов решения задачи 6 2.1.Метод Рунге-Кутта четвертого порядка для решения

уравнения первого порядка 6

2.2.Задача Коши 6

2.3.Метод Булирша- Штера с использованием

рациональной экстраполяции для системы уравнений 7

2.4 Метод Адамса 8

2.5. Метод Эйлера 9

3. Описание алгоритмов решения задания 13

3.1. Описание переменных 13

3.2. Блок- схема главного модуля 14

3.3. Описание алгоритма главной программы 14

3.4. Блок-схема функции “func” 15

3.5. Описание блок- схемы функции “func” 15

4. Описание программного обеспечения 16

4.1. Описание операционной системы 16

4.2. Описание языка программирования 18

4.3. Описание программы 19

5. Контрольный пример 21

6.Анализ полученных результатов 22

Список литературы 24

Уравнение называется обыкновенным дифференциальным n-го порядка, если F определена и непрерывна в некоторой области и, во всяком случае, зависит от . Его решением является любая функция u( x), которая этому уравнению удовлетворяет при всех x в определённом конечном или бесконечном интервале. Дифференциальное уравнение, разрешенное относительно старшей производной имеет вид

Решением этого уравнения на интервале I=[a,b] называется функция u( x)

Решить дифференциальное уравнение у / =f(x,y) численным методом — это значит для заданной последовательности аргументов х₀ , х₁ …, х_n и числа у₀ , не определяя функцию у=F(x), найти такие значения у₁ , у₂ ,…, у_n , что у_i =F(x_i )(i=1,2,…, n) и F(x₀ )=y₀ .

Таким образом, численные методы позволяют вместо нахождения функции y=F(x) (3) получить таблицу значений этой функции для заданной последовательности аргументов. Величина h=x_k -x_k-1 называется шагом интегрирования.

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

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

1.1. Решить приближенно дифференциальное уравнение вида методом Эйлера

1.2. Составить блок-схему алгоритма для решения данного задания.

1.3. Разработать программу на языке MicrosoftVisualC++

1.4. Протестировать программу на примере y’=2x+y (n=5, [0,1], y0=1)

1.5. Выполнить анализ результатов.

1.6. Оформить пояснительную записку с приложением.

2.Обзор методов решения задачи.

2.1. Метод Рунге-Кутта четвертого порядка для решения уравнения первого порядка.

Идея Рунге-Кута состоит в том, чтобы использовать метод неопределённых коэффициентов. Наиболее употребительным методом Рунге-Кутта решения уравнения первого порядка y’ = F(x,y) (2.1.1) является метод четвертого порядка, в котором вычисления производятся по формуле:

yk+1 = yk +(k1 +2k2 +2k3 +k4 )/6, (2.1.2)

k1 = Fk h = F(xk , yk )h

k2 = F(xk +h/2, yk +k1 /2)h

k3 = F(xk +h/2, yk +k2 /2)h

k4 = F(xk +h, yk +k3 )h,

Рассмотрим задачу Коши для уравнений первого порядка на отрезке [a,b]:

, (2.1.4)

Разобьём промежуток [a,b] на Nчастей . Обозначим , где u(x) –точное решение задачи Коши, и через значения приближенного решения в точках . Существует 2 типа численных схем :

1. явные: ) (2.2.1)

2. неявные: (2.2.2)

Здесь F некоторая функция, связывающая приближения. В явных схемах приближенное значение в точке определяется через некоторое число kуже определённых приближенных значений. В неявных схемах определяется не рекурентным способом, как в явных схемах, а для его определения возникает уравнение, поскольку равенство (2.2.2) представляет из себя именно уравнение на . Явные схемы проще, однако зачастую неявные схемы предпочтительнее.

2.3. Метод Булирша-Штера с использованием рациональной экстраполяции для системы уравнений

Метод Булирша-Штера (Bulirsch-Stoer Method) — это метод решения системы обыкновенных дифференциальных уравнений первого порядка с гладкими правыми частями. Гладкость правых частей является необходимой для работы метода. Если правые части вашей системы не являются гладкими или содержат разрывы, то лучше использовать метод Рунге-Кутта. В случае же гладкой системы метод Булирша-Штера позволяет добиться существенно большей точности, чем метод Рунге-Кутта.

Принцип работы метода

Основной идеей метода является вычисление состояния системы в точке x+h, как результата двух шагов длины h/2, четырех шагов длины h/4, восьми шагов длины h/8 и так далее с последующей экстраполяцией результатов. Метод строит рациональную интерполирующую функцию, которая в точке h/2 проходит через состояние системы после двух таких шагов, в точке h/4 проходит через состояние системы после четырех таких шагов, и т.д., а затем вычисляет значение этой функции в точке h = 0, проводя экстраполяцию.

Гладкость правых частей приводит к тому, что вычисленное при помощи экстраполяции состояние системы оказывается очень близко к действительному, а использование рациональной экстраполяции вместо полиномиальной позволяет ещё больше повысить точность.

Таким образом проводится один шаг метода, после чего принимается решение — следует ли изменять шаг, а если да — то в какую сторону. При этом используется оценка погрешности, которую мы получаем в качестве дополнительного результата при рациональной экстраполяции. Следует отметить, что алгоритм решает автономную систему, т.е. если уравнения системы содержат время, то необходимо ввести время в качестве переменной, производная от которой тождественно равна единице.

Явная схема Адамса.

Рассмотренные выше методы являются явными одношаговыми (для нахождения последующего приближения используется лишь одно предыдущее). Приведённый ниже метод является многошаговым.

Пусть задана задача Коши:

(2.4.1)

Для точного решения (которое нам не известно) выполнено:

(2.4.2)

Предположим, нам известны приближенные значения функции u ( x ) в k точках (стартовые k точек, в частности, можно найти методом Эйлера или методом Рунге-Кутта того или иного порядка), тогда функцию f ( x , u ( x )) в (2.4.2) для приближенного вычисления интеграла можно заменить на интерполяционный полином порядка k -1, построенный поk точкам , интеграл от которого считается явно и представляет собой линейную комбинацию значений cнекоторыми множителями . Таким образом, мы получаем следующую рекуррентную процедуру вычисления приближенных значений функции u ( x ) (являющимся точным решением задачи Коши) в точках :

(2.4.3)

Описанная схема является k — шаговой явной формулой Адамса.

Неявная схема Адамса.

Пусть — интерполяционный полином порядка k , построенный по k +1 значению б одно из которых, именно , мы будем считать неизвестным. Модифицируем (2.4.3), заменив в нём на полином более высокой степени , интеграл от которого выражается в виде линейной комбинации значений с некоторыми новыми коэффициентами :

(2.4.4)

Формула (2.4.4) представляет собой неявную схему Адамса и является уравнением на , которое можно решать методом последовательных приближений. Естественно, что начальное приближение , должно быть разумно выбрано. Для этого удобно объединить явную и неявную схемы Адамса в одну, называемую «методом коррекции». Именно с помощью явной схемы определяется начальное приближение (прогноз), а затем по неявной схеме оно необходимое число раз (обычно один или два) корректируется методом последовательных приближений до достижения заданной точности (коррекция).

Решить дифференциальное уравнение у / =f(x,y) численным методом — это значит для заданной последовательности аргументов х₀ , х₁ …, х_n и числа у₀ , не определяя функцию у=F(x), найти такие значения у₁ , у₂ ,…, у_n , что у_i =F(x_i )(i=1,2,…, n) и F(x₀ )=y₀ . (2.5.1)

Таким образом, численные методы позволяют вместо нахождения функции

У=F(x) получить таблицу значений этой функции для заданной последовательности аргументов. Величина h=x_k -x_k-1 называется шагом интегрирования.

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

Рассмотрим дифференциальное уравнение первого порядка (2.5.1)

с начальным условием

Требуется найти решение уравнения (2.5.1) на отрезке [а,b].

Разобьем отрезок [a, b] на n равных частей и получим последовательность х₀ , х₁ , х₂ ,…, х_n , где x_i =x₀ +ih (i=0,1,…, n), а h=(b-a)/n-шаг интегрирования.

В методе Эйлера приближенные значения у(х_i )»y_i вычисляются последовательно по формулам у_i +hf(x_i , y_i ) (i=0,1,2…).

При этом искомая интегральная кривая у=у(х), проходящая через точку М₀ (х₀ , у₀ ), заменяется ломаной М₀ М₁ М₂ … с вершинами М_i (x_i , y_i ) (i=0,1,2,…); каждое звено М_i M_{i +1} этой ломаной, называемой ломаной Эйлера, имеет направление, совпадающее с направлением той интегральной кривой уравнения (2.5.1), которая проходит через точку М_i . Если правая часть уравнения (2.5.1) в некотором прямоугольнике R<|x-x₀ |£a, |y-y₀ |£b>удовлетворяет условиям:
|f(x, y₁ )- f(x, y₂ )| £N|y₁ -y₂ | (N=const), (2.5.3)

|df/dx|=|df/dx+f(df/dy)| £ M (M=const),

то имеет место следующая оценка погрешности:

где у(х_n )-значение точного решения уравнения (2.5.1) при х=х_n , а у_n — приближенное значение, полученное на n-ом шаге.

Формула (13) имеет в основном теоретическое применение. На практике иногда оказывается более удобным двойной просчет : сначала расчет ведется с шагом h, затем шаг дробят и повторный расчет ведется с шагомh/2. Погрешность более точного значения у_n * оценивается формулой

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

Модифицированный метод Эйлера

Рассмотрим дифференциальное уравнение (2.5.1) y / =f(x,y) с начальным условием y(x₀ )=y₀ . Разобьем наш участок интегрирования на n равных частей. На малом участинтегральную кривую заменим прямой линией.

Рис.1 Метод Эйлера в графическом видa

Получаем точку N_k / . В этой точке строим следующую касательную:

Из точки М_к проводим прямую с угловым коэффициентом α_к и определяем точку пересечения этой прямой с прямой Х_к1 . Получаем точку М_к / . В качестве у_к+1 принимаем ординату точки М_к / . Тогда:

(2.5.8)-рекурентные формулы метода Эйлера.

Сначала вычисляют вспомогательные значения искомой функции у_к+1/2 в точках х_к+1/2 , затем находят значение правой части уравнения (11) в средней точке y / _k+1/2 =f(x_{k +1/2} , y_{k +1/2} ) и определяют у_к+1 .

Для оценки погрешности в точке х_к проводят вычисления у_к с шагом h, затем с шагом 2h и берут 1/3 разницы этих значений:

где у(х)-точное решение дифференциального уравнения.

Таким образом, методом Эйлера можно решать уравнения любых порядков. Например, чтобы решить уравнение второго порядка y // =f(y / ,y,x) c начальными условиями y / (x₀ )=y / ₀ , y(x₀ )=y₀ , выполняется замена:

3.Описание алгоритмов решения задачи

3.2. Блок- схема главного модуля

3.3 Описание алгоритма главной программы.

Вывод результатов на интервале X

3.4 Блок-схема функции “func”.

3.5 Описание блок- схемы функции “func”.

*Реализация алгоритма на языке программирования C++ представлена в приложении .

4.Описание программного обеспечения.

4.1 Описание операционной системы

Основное требование к операционной системе (ОС), предъявляемое поставленной задачей, это наличие ANSI или POSIX совместимого компилятора языка C++.

Для реализации задачи была выбрана последняя клиентская версия операционной системы Microsoft, основанная на ядре NT – MicrosoftWindowsXPProfessional.

Указанная операционная система обладает рядом преимуществ:

· наличие достаточного количество ANSI или POSIX совместимых компиляторов языка C++, разработанных для данной ОС, а именно –

o Microsoft C++ (version 2-6)

· достаточная управляемость, надежность и безопасность;

· широкое распространение основанных на ядре NT операционных систем Microsoft, совместимых по программному обеспечению с WindowsXPProfessional (NT/2000/XP/2003 – client & server);

· высокая скорость работы приложений, разработанных для данной ОС с использованием компиляторов C++.

Исходный код программы может быть откомпилирован и под другой операционной системой, если для таковой имеется ANSI или POSIX совместимый компилятор языка C++.

Программа была протестирована на операционной системе MicrosoftWindowsXPProfessionalSP1.

· Процессор x86 Family 15 Model 2 Stepping 7 GenuineIntel

· Версия BIOS Award Software International, Inc. F4, 06.03.2003

· Аппаратно-зависимый уровень (HAL) Версия = «5.1.2600.1106 (xpsp1.020828-1920)»

· Полный объем физической памяти 256,00 МБ

· Доступно физической памяти 29,97 МБ

· Всего виртуальной памяти 873,69 МБ

· Доступно виртуальной памяти 350,04 МБ

· Файл подкачки 618,21 МБ

4.2 Описание языка программирования

Язык программирования С++

С++ — это универсальный язык программирования, задуманный так, чтобы сделать программирование более приятным для серьезного программиста. За исключением второстепенных деталей С++ является надмножеством языка программирования C. Помимо возможностей, которые дает C, С++ предоставляет гибкие и эффективные средства определения новых типов. Используя определения новых типов, точно отвечающих концепциям приложения, программист может разделять разрабатываемую программу на лег ко поддающиеся контролю части. Такой метод построения программ часто называют абстракцией данных. Информация о типах содержится в некоторых объектах типов, определенных пользователем. Такие объекты просты и надежны в использовании в тех ситуациях, когда их тип нельзя установить на стадии компиляции. Программирование с применением таких объектов часто называют объектно-ориентированным. При правильном использовании этот метод дает более короткие, проще понимаемые и легче контролируемые программы.

В С++ нет типов данных высокого уровня и нет первичных операций высокого уровня. В нем нет, например, матричного типа с операцией обращения или типа строка с операцией конкатенации. Если пользователю понадобятся подобные типы, их можно определить в самом языке. По сути дела, основное, чем занимается программирование на С++ — это определение универсальных и специально-прикладных типов. Хорошо разработанный тип, определяемый пользователем, отличается от встроенного типа только способом определения, но не способом использования.

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

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

Компилятор MicrosoftC++ и среда разработки MicrosoftVisualStudio

В качестве компилятора для разработки приложения был выбран MicrosoftC++ по следующим причинам:

· практически полная совместимость со стандартом ANSIC++;

· наличие удобной среды разработки MicrosoftVisualStudio;

· наличие отличной документации;

· высокая скорость работы результирующих приложений;

· совместимость разработанных приложений с большим количеством широко распространенных операционных систем;

· достаточная скорость компиляции.

4.3 Описание программы

Разработанное приложение поставляется в виде 2-ух файлов:

1. methodEulera.cpp – исходный код программы на языке C++;

2. methodEulera.exe – исполняемый файл.

Для выполнения исполняемого файла необходима одна из ниже перечисленных операционных систем:

· Microsoft Windows 3.11+Win32s;

· Microsoft Windows 95/98/Me;

· MicrosoftWindowsNT/2000/XP/2003 – клиентская или серверная версия.

Программа не требует предварительной установки и может быть сразу же запущена на выполнение.

Исходный код приложения может быть откомпилирован в любом ANSI или POSIX совместимом компиляторе С++ для получения выполнимой программы. Для успешной компиляции требуется наличие стандартной библиотеки «iostream».

5. Контрольный пример

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

В результате вычислений контрольного примера вида y’=2 x+ y с интервалом [0,1],

количеством шагов равному 5 и начальным условием у равным 1 , с помощью программы, получились следующие результаты:

Рис. 2. Экран с результатами выполнения программы.

Как видно, при вычислении программа на первом шаге берёт начальные значения для вычисления, а на последующих берёт значения полученные с предыдущих шагов. Можно сделать вывод, что точность вычисления данного метода зависит от количества выбранных шагов: чем больше шагов, тем меньше фиксированное приращение , а следовательно она более точно вычисляет значение всего интервала.

По работе программы стало видно, что с её использованием намного упростилась работа пользователя. Пользователь просто вводит интервал на котором должен вычисляться пример, количество шагов и начальное значения и программа выдаёт уже готовое решение данного примера.

6.Анализ полученных результатов.

По результатам программы можно составить таблицу сравнения результатов полученных при использовании программы и результатов, полученных ручным способом:

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

Решение дифференциальных уравнений методом Эйлера можно также отобразить в графическом виде:

Рис.3.Графическое изображение решения примера y’=2x+y

Как видно из рис.3 графиком решения уравнения является кривая , форма которой зависит от количества разбиений интервала.

По результатам выполненной работы можно сделать вывод, что решение дифференциальных уравнений методом Эйлера является методом вычисления со средней точностью и точность вычисления данного метода зависит от количества разбиений интервала интегрирования. При сравнении результатов решенными разными способами можно сказать, что данный метод был верно реализован на языке программирования MicrosoftVisualC++. Полученные результаты сходятся с небольшой погрешностью.

1. Численные методы (анализ, алгебра, обыкновенные дифференциальные уравнения), Н.С. Бахвалов. Главная редакция физико- математической литературы изд-ва «Наука», М., 1975г.

2. Методы, теории обыкновенных дифференциальных уравнений. Н.И. Гаврилов . Государственное издательство «Высшая школа» Москва-1962г.

3. В.В.Пак., Ю.Л. Носенко. Высшая математика: Учебник.- Д.: Сталкер, 1997г.

4. Б. П. Демидович, И. А. Марон Основы вычислительной математике. – М., 1966

5. Загускин В. Л. – Справочник по численным методам решения уравнений. – М.: ФИЗМАТГИЗ, 1960. – 216 с.

6. Либерти, Джесс.

Освой самостоятельно С++ за 21 день, 4-е издание.:Пер с англ.-М.: Издательский дом «Вильямс», 2003.-832с.

7. П.Нортон, П.Иао «Программирование на С++ в среде Windows» («Диалектика» Киев 2003г.)

8. Янг М. Microsoft Visual C++ — М.:ЭНТРОП, 2000.

9. Марченко А.И., Марченко Л.А. – Программирование в среде

Turbo Pascal 7.0 – К.: ВЕК+, М.: Бином Универсал, 1998. – 496 с.

10. Высшая математика: Справ. материалы: Книга для учащихся .- М.:

Реферат: Решение дифференциальных уравнений 1 порядка методом Эйлера

1. Постановка задачи 5

2. Обзор существующих методов решения задачи 6 2.1.Метод Рунге-Кутта четвертого порядка для решения

уравнения первого порядка 6

2.2.Задача Коши 6

2.3.Метод Булирша- Штера с использованием

рациональной экстраполяции для системы уравнений 7

2.4 Метод Адамса 8

2.5. Метод Эйлера 9

3. Описание алгоритмов решения задания 13

3.1. Описание переменных 13

3.2. Блок- схема главного модуля 14

3.3. Описание алгоритма главной программы 14

3.4. Блок-схема функции “func” 15

3.5. Описание блок- схемы функции “func” 15

4. Описание программного обеспечения 16

4.1. Описание операционной системы 16

4.2. Описание языка программирования 18

4.3. Описание программы 19

5. Контрольный пример 21

6.Анализ полученных результатов 22

Список литературы 24

Уравнение называется обыкновенным дифференциальным n-го порядка, если F определена и непрерывна в некоторой области и, во всяком случае, зависит от . Его решением является любая функция u( x), которая этому уравнению удовлетворяет при всех x в определённом конечном или бесконечном интервале. Дифференциальное уравнение, разрешенное относительно старшей производной имеет вид

Решением этого уравнения на интервале I=[a,b] называется функция u( x)

Решить дифференциальное уравнение у / =f(x,y) численным методом — это значит для заданной последовательности аргументов х₀ , х₁ …, х_n и числа у₀ , не определяя функцию у=F(x), найти такие значения у₁ , у₂ ,…, у_n , что у_i =F(x_i )(i=1,2,…, n) и F(x₀ )=y₀ .

Таким образом, численные методы позволяют вместо нахождения функции y=F(x) (3) получить таблицу значений этой функции для заданной последовательности аргументов. Величина h=x_k -x_k-1 называется шагом интегрирования.

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

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

1.1. Решить приближенно дифференциальное уравнение вида методом Эйлера

1.2. Составить блок-схему алгоритма для решения данного задания.

1.3. Разработать программу на языке MicrosoftVisualC++

1.4. Протестировать программу на примере y’=2x+y (n=5, [0,1], y0=1)

1.5. Выполнить анализ результатов.

1.6. Оформить пояснительную записку с приложением.

2.Обзор методов решения задачи.

2.1. Метод Рунге-Кутта четвертого порядка для решения уравнения первого порядка.

Идея Рунге-Кута состоит в том, чтобы использовать метод неопределённых коэффициентов. Наиболее употребительным методом Рунге-Кутта решения уравнения первого порядка y’ = F(x,y) (2.1.1) является метод четвертого порядка, в котором вычисления производятся по формуле:

yk+1 = yk +(k1 +2k2 +2k3 +k4 )/6, (2.1.2)

k1 = Fk h = F(xk , yk )h

k2 = F(xk +h/2, yk +k1 /2)h

k3 = F(xk +h/2, yk +k2 /2)h

k4 = F(xk +h, yk +k3 )h,

Рассмотрим задачу Коши для уравнений первого порядка на отрезке [a,b]:

, (2.1.4)

Разобьём промежуток [a,b] на Nчастей . Обозначим , где u(x) –точное решение задачи Коши, и через значения приближенного решения в точках . Существует 2 типа численных схем :

1. явные: ) (2.2.1)

2. неявные: (2.2.2)

Здесь F некоторая функция, связывающая приближения. В явных схемах приближенное значение в точке определяется через некоторое число kуже определённых приближенных значений. В неявных схемах определяется не рекурентным способом, как в явных схемах, а для его определения возникает уравнение, поскольку равенство (2.2.2) представляет из себя именно уравнение на . Явные схемы проще, однако зачастую неявные схемы предпочтительнее.

2.3. Метод Булирша-Штера с использованием рациональной экстраполяции для системы уравнений

Метод Булирша-Штера (Bulirsch-Stoer Method) — это метод решения системы обыкновенных дифференциальных уравнений первого порядка с гладкими правыми частями. Гладкость правых частей является необходимой для работы метода. Если правые части вашей системы не являются гладкими или содержат разрывы, то лучше использовать метод Рунге-Кутта. В случае же гладкой системы метод Булирша-Штера позволяет добиться существенно большей точности, чем метод Рунге-Кутта.

Принцип работы метода

Основной идеей метода является вычисление состояния системы в точке x+h, как результата двух шагов длины h/2, четырех шагов длины h/4, восьми шагов длины h/8 и так далее с последующей экстраполяцией результатов. Метод строит рациональную интерполирующую функцию, которая в точке h/2 проходит через состояние системы после двух таких шагов, в точке h/4 проходит через состояние системы после четырех таких шагов, и т.д., а затем вычисляет значение этой функции в точке h = 0, проводя экстраполяцию.

Гладкость правых частей приводит к тому, что вычисленное при помощи экстраполяции состояние системы оказывается очень близко к действительному, а использование рациональной экстраполяции вместо полиномиальной позволяет ещё больше повысить точность.

Таким образом проводится один шаг метода, после чего принимается решение — следует ли изменять шаг, а если да — то в какую сторону. При этом используется оценка погрешности, которую мы получаем в качестве дополнительного результата при рациональной экстраполяции. Следует отметить, что алгоритм решает автономную систему, т.е. если уравнения системы содержат время, то необходимо ввести время в качестве переменной, производная от которой тождественно равна единице.

Явная схема Адамса.

Рассмотренные выше методы являются явными одношаговыми (для нахождения последующего приближения используется лишь одно предыдущее). Приведённый ниже метод является многошаговым.

Пусть задана задача Коши:

(2.4.1)

Для точного решения (которое нам не известно) выполнено:

(2.4.2)

Предположим, нам известны приближенные значения функции u ( x ) в k точках (стартовые k точек, в частности, можно найти методом Эйлера или методом Рунге-Кутта того или иного порядка), тогда функцию f ( x , u ( x )) в (2.4.2) для приближенного вычисления интеграла можно заменить на интерполяционный полином порядка k -1, построенный поk точкам , интеграл от которого считается явно и представляет собой линейную комбинацию значений cнекоторыми множителями . Таким образом, мы получаем следующую рекуррентную процедуру вычисления приближенных значений функции u ( x ) (являющимся точным решением задачи Коши) в точках :

(2.4.3)

Описанная схема является k — шаговой явной формулой Адамса.

Неявная схема Адамса.

Пусть — интерполяционный полином порядка k , построенный по k +1 значению б одно из которых, именно , мы будем считать неизвестным. Модифицируем (2.4.3), заменив в нём на полином более высокой степени , интеграл от которого выражается в виде линейной комбинации значений с некоторыми новыми коэффициентами :

(2.4.4)

Формула (2.4.4) представляет собой неявную схему Адамса и является уравнением на , которое можно решать методом последовательных приближений. Естественно, что начальное приближение , должно быть разумно выбрано. Для этого удобно объединить явную и неявную схемы Адамса в одну, называемую «методом коррекции». Именно с помощью явной схемы определяется начальное приближение (прогноз), а затем по неявной схеме оно необходимое число раз (обычно один или два) корректируется методом последовательных приближений до достижения заданной точности (коррекция).

Решить дифференциальное уравнение у / =f(x,y) численным методом — это значит для заданной последовательности аргументов х₀ , х₁ …, х_n и числа у₀ , не определяя функцию у=F(x), найти такие значения у₁ , у₂ ,…, у_n , что у_i =F(x_i )(i=1,2,…, n) и F(x₀ )=y₀ . (2.5.1)

Таким образом, численные методы позволяют вместо нахождения функции

У=F(x) получить таблицу значений этой функции для заданной последовательности аргументов. Величина h=x_k -x_k-1 называется шагом интегрирования.

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

Рассмотрим дифференциальное уравнение первого порядка (2.5.1)

с начальным условием

Требуется найти решение уравнения (2.5.1) на отрезке [а,b].

Разобьем отрезок [a, b] на n равных частей и получим последовательность х₀ , х₁ , х₂ ,…, х_n , где x_i =x₀ +ih (i=0,1,…, n), а h=(b-a)/n-шаг интегрирования.

В методе Эйлера приближенные значения у(х_i )»y_i вычисляются последовательно по формулам у_i +hf(x_i , y_i ) (i=0,1,2…).

При этом искомая интегральная кривая у=у(х), проходящая через точку М₀ (х₀ , у₀ ), заменяется ломаной М₀ М₁ М₂ … с вершинами М_i (x_i , y_i ) (i=0,1,2,…); каждое звено М_i M_{i +1} этой ломаной, называемой ломаной Эйлера, имеет направление, совпадающее с направлением той интегральной кривой уравнения (2.5.1), которая проходит через точку М_i . Если правая часть уравнения (2.5.1) в некотором прямоугольнике R<|x-x₀ |£a, |y-y₀ |£b>удовлетворяет условиям:
|f(x, y₁ )- f(x, y₂ )| £N|y₁ -y₂ | (N=const), (2.5.3)

|df/dx|=|df/dx+f(df/dy)| £ M (M=const),

то имеет место следующая оценка погрешности:

где у(х_n )-значение точного решения уравнения (2.5.1) при х=х_n , а у_n — приближенное значение, полученное на n-ом шаге.

Формула (13) имеет в основном теоретическое применение. На практике иногда оказывается более удобным двойной просчет : сначала расчет ведется с шагом h, затем шаг дробят и повторный расчет ведется с шагомh/2. Погрешность более точного значения у_n * оценивается формулой

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

Модифицированный метод Эйлера

Рассмотрим дифференциальное уравнение (2.5.1) y / =f(x,y) с начальным условием y(x₀ )=y₀ . Разобьем наш участок интегрирования на n равных частей. На малом участинтегральную кривую заменим прямой линией.

Рис.1 Метод Эйлера в графическом видa

Получаем точку N_k / . В этой точке строим следующую касательную:

Из точки М_к проводим прямую с угловым коэффициентом α_к и определяем точку пересечения этой прямой с прямой Х_к1 . Получаем точку М_к / . В качестве у_к+1 принимаем ординату точки М_к / . Тогда:

(2.5.8)-рекурентные формулы метода Эйлера.

Сначала вычисляют вспомогательные значения искомой функции у_к+1/2 в точках х_к+1/2 , затем находят значение правой части уравнения (11) в средней точке y / _k+1/2 =f(x_{k +1/2} , y_{k +1/2} ) и определяют у_к+1 .

Для оценки погрешности в точке х_к проводят вычисления у_к с шагом h, затем с шагом 2h и берут 1/3 разницы этих значений:

где у(х)-точное решение дифференциального уравнения.

Таким образом, методом Эйлера можно решать уравнения любых порядков. Например, чтобы решить уравнение второго порядка y // =f(y / ,y,x) c начальными условиями y / (x₀ )=y / ₀ , y(x₀ )=y₀ , выполняется замена:

3.Описание алгоритмов решения задачи

3.2. Блок- схема главного модуля

3.3 Описание алгоритма главной программы.

Вывод результатов на интервале X

3.4 Блок-схема функции “func”.

3.5 Описание блок- схемы функции “func”.

*Реализация алгоритма на языке программирования C++ представлена в приложении .

4.Описание программного обеспечения.

4.1 Описание операционной системы

Основное требование к операционной системе (ОС), предъявляемое поставленной задачей, это наличие ANSI или POSIX совместимого компилятора языка C++.

Для реализации задачи была выбрана последняя клиентская версия операционной системы Microsoft, основанная на ядре NT – MicrosoftWindowsXPProfessional.

Указанная операционная система обладает рядом преимуществ:

· наличие достаточного количество ANSI или POSIX совместимых компиляторов языка C++, разработанных для данной ОС, а именно –

o Microsoft C++ (version 2-6)

· достаточная управляемость, надежность и безопасность;

· широкое распространение основанных на ядре NT операционных систем Microsoft, совместимых по программному обеспечению с WindowsXPProfessional (NT/2000/XP/2003 – client & server);

· высокая скорость работы приложений, разработанных для данной ОС с использованием компиляторов C++.

Исходный код программы может быть откомпилирован и под другой операционной системой, если для таковой имеется ANSI или POSIX совместимый компилятор языка C++.

Программа была протестирована на операционной системе MicrosoftWindowsXPProfessionalSP1.

· Процессор x86 Family 15 Model 2 Stepping 7 GenuineIntel

· Версия BIOS Award Software International, Inc. F4, 06.03.2003

· Аппаратно-зависимый уровень (HAL) Версия = «5.1.2600.1106 (xpsp1.020828-1920)»

· Полный объем физической памяти 256,00 МБ

· Доступно физической памяти 29,97 МБ

· Всего виртуальной памяти 873,69 МБ

· Доступно виртуальной памяти 350,04 МБ

· Файл подкачки 618,21 МБ

4.2 Описание языка программирования

Язык программирования С++

С++ — это универсальный язык программирования, задуманный так, чтобы сделать программирование более приятным для серьезного программиста. За исключением второстепенных деталей С++ является надмножеством языка программирования C. Помимо возможностей, которые дает C, С++ предоставляет гибкие и эффективные средства определения новых типов. Используя определения новых типов, точно отвечающих концепциям приложения, программист может разделять разрабатываемую программу на лег ко поддающиеся контролю части. Такой метод построения программ часто называют абстракцией данных. Информация о типах содержится в некоторых объектах типов, определенных пользователем. Такие объекты просты и надежны в использовании в тех ситуациях, когда их тип нельзя установить на стадии компиляции. Программирование с применением таких объектов часто называют объектно-ориентированным. При правильном использовании этот метод дает более короткие, проще понимаемые и легче контролируемые программы.

В С++ нет типов данных высокого уровня и нет первичных операций высокого уровня. В нем нет, например, матричного типа с операцией обращения или типа строка с операцией конкатенации. Если пользователю понадобятся подобные типы, их можно определить в самом языке. По сути дела, основное, чем занимается программирование на С++ — это определение универсальных и специально-прикладных типов. Хорошо разработанный тип, определяемый пользователем, отличается от встроенного типа только способом определения, но не способом использования.

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

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

Компилятор MicrosoftC++ и среда разработки MicrosoftVisualStudio

В качестве компилятора для разработки приложения был выбран MicrosoftC++ по следующим причинам:

· практически полная совместимость со стандартом ANSIC++;

· наличие удобной среды разработки MicrosoftVisualStudio;

· наличие отличной документации;

· высокая скорость работы результирующих приложений;

· совместимость разработанных приложений с большим количеством широко распространенных операционных систем;

· достаточная скорость компиляции.

4.3 Описание программы

Разработанное приложение поставляется в виде 2-ух файлов:

1. methodEulera.cpp – исходный код программы на языке C++;

2. methodEulera.exe – исполняемый файл.

Для выполнения исполняемого файла необходима одна из ниже перечисленных операционных систем:

· Microsoft Windows 3.11+Win32s;

· Microsoft Windows 95/98/Me;

· MicrosoftWindowsNT/2000/XP/2003 – клиентская или серверная версия.

Программа не требует предварительной установки и может быть сразу же запущена на выполнение.

Исходный код приложения может быть откомпилирован в любом ANSI или POSIX совместимом компиляторе С++ для получения выполнимой программы. Для успешной компиляции требуется наличие стандартной библиотеки «iostream».

5. Контрольный пример

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

В результате вычислений контрольного примера вида y’=2 x+ y с интервалом [0,1],

количеством шагов равному 5 и начальным условием у равным 1 , с помощью программы, получились следующие результаты:

Рис. 2. Экран с результатами выполнения программы.

Как видно, при вычислении программа на первом шаге берёт начальные значения для вычисления, а на последующих берёт значения полученные с предыдущих шагов. Можно сделать вывод, что точность вычисления данного метода зависит от количества выбранных шагов: чем больше шагов, тем меньше фиксированное приращение , а следовательно она более точно вычисляет значение всего интервала.

По работе программы стало видно, что с её использованием намного упростилась работа пользователя. Пользователь просто вводит интервал на котором должен вычисляться пример, количество шагов и начальное значения и программа выдаёт уже готовое решение данного примера.

6.Анализ полученных результатов.

По результатам программы можно составить таблицу сравнения результатов полученных при использовании программы и результатов, полученных ручным способом:

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

Решение дифференциальных уравнений методом Эйлера можно также отобразить в графическом виде:

Рис.3.Графическое изображение решения примера y’=2x+y

Как видно из рис.3 графиком решения уравнения является кривая , форма которой зависит от количества разбиений интервала.

По результатам выполненной работы можно сделать вывод, что решение дифференциальных уравнений методом Эйлера является методом вычисления со средней точностью и точность вычисления данного метода зависит от количества разбиений интервала интегрирования. При сравнении результатов решенными разными способами можно сказать, что данный метод был верно реализован на языке программирования MicrosoftVisualC++. Полученные результаты сходятся с небольшой погрешностью.

1. Численные методы (анализ, алгебра, обыкновенные дифференциальные уравнения), Н.С. Бахвалов. Главная редакция физико- математической литературы изд-ва «Наука», М., 1975г.

2. Методы, теории обыкновенных дифференциальных уравнений. Н.И. Гаврилов . Государственное издательство «Высшая школа» Москва-1962г.

3. В.В.Пак., Ю.Л. Носенко. Высшая математика: Учебник.- Д.: Сталкер, 1997г.

4. Б. П. Демидович, И. А. Марон Основы вычислительной математике. – М., 1966

5. Загускин В. Л. – Справочник по численным методам решения уравнений. – М.: ФИЗМАТГИЗ, 1960. – 216 с.

6. Либерти, Джесс.

Освой самостоятельно С++ за 21 день, 4-е издание.:Пер с англ.-М.: Издательский дом «Вильямс», 2003.-832с.

7. П.Нортон, П.Иао «Программирование на С++ в среде Windows» («Диалектика» Киев 2003г.)

8. Янг М. Microsoft Visual C++ — М.:ЭНТРОП, 2000.

9. Марченко А.И., Марченко Л.А. – Программирование в среде

Turbo Pascal 7.0 – К.: ВЕК+, М.: Бином Универсал, 1998. – 496 с.

10. Высшая математика: Справ. материалы: Книга для учащихся .- М.:

Решение дифференциальных уравнений по методу Эйлера

Все приложения, графические материалы, формулы, таблицы и рисунки работы на тему: Решение дифференциальных уравнений по методу Эйлера (предмет: Математика) находятся в архиве, который можно скачать с нашего сайта. Приступая к прочтению данного произведения (перемещая полосу прокрутки браузера вниз), Вы соглашаетесь с условиями открытой лицензии Creative Commons «Attribution» («Атрибуция») 4.0 Всемирная (CC BY 4.0) .

Министерство образования Республики Беларусь

БЕЛОРУССКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

ИНФОРМАТИКИ И РАДИОЭЛЕКТРОНИКИ

Кафедра вычислительной математики и программирования

к курсовому проекту

«Решение дифференциальных уравнений по методу Эйлера»

2. Математическое объяснение метода

2.1 Метод Эйлера

2.2 Исправленный метод Эйлера

2.3 Модифицированный метод Эйлера

3. Блок-схема алгоритма программы

4. Описание программы

Список использованной литературы

Приложение 1 (Текст программы)

Приложение 2 (Результаты работы программы)

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

Решение обыкновенных дифференциальных уравнений (нелинейных) первого порядка с начальными данными (задача Коши) — классическая область применения численных методов. Имеется много разностных методов, часть из которых возникла в домашинную эпоху и оказалось пригодным для современных ЭВМ.

В этой программе использовался метод Эйлера, один из самых старых и широко известных методов численного интегрирования дифференциальных уравнений. Этот метод имеет довольно большую ошибку; кроме того, он очень часто оказывается неустойчивым — малая начальная ошибка быстро увеличивается с ростом Х. Поэтому чаще используют более точные методы, такие как: исправленный метод Эйлера и модифицированный метод Эйлера. Нужно, однако, заметить, что метод Эйлера является методом Рунге — Кутта первого порядка.

2. Математическое объяснение метода

2.1 Метод Эйлера

Решить дифференциальное уравнение у / =f(x,y) численным методом — это значит для заданной последовательности аргументов х₀, х₁…, х_n и числа у₀, не определяя функцию у=F(x), найти такие значения у₁, у₂,…, у_n, что у_i=F(x_i)(i=1,2,…, n) и F(x₀)=y₀.

Таким образом, численные методы позволяют вместо нахождения функции

У=F(x) получить таблицу значений этой функции для заданной последовательности аргументов. Величина h=x_k-x_k-1 называется шагом интегрирования.

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

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

с начальным условием

Требуется найти решение уравнения (1) на отрезке [а,b].

Разобьем отрезок [a, b] на n равных частей и получим последовательность х₀, х₁, х₂,…, х_n, где x_i=x₀+ih (i=0,1,…, n), а h=(b-a)/n-шаг интегрирования.

В методе Эйлера приближенные значения у(х_i)y_i вычисляются последовательно по формулам у_i+hf(x_i, y_i) (i=0,1,2…).

При этом искомая интегральная кривая у=у(х), проходящая через точку М₀(х₀, у₀), заменяется ломаной М₀М₁М₂… с вершинами М_i(x_i, y_i) (i=0,1,2,…); каждое звено М_iM_i+1 этой ломаной, называемой ломаной Эйлера, имеет направление, совпадающее с направлением той интегральной кривой уравнения (1), которая проходит через точку М_i, смотри рисунок 1.

Если правая часть уравнения (1) в некотором прямоугольнике R<|x-x₀|a, |y-y₀|b>удовлетворяет условиям:

|df/dx|=|df/dx+f(df/dy)| M (M=const),

то имеет место следующая оценка погрешности:

где у(х_n)-значение точного решения уравнения(1) при х=х_n, а у_n— приближенное значение, полученное на n-ом шаге.

Формула (3) имеет в основном теоретическое применение. На практике иногда оказывается более удобным двойной просчет: сначала расчет ведется с шагом h, затем шаг дробят и повторный расчет ведется с шагом h/2. Погрешность более точного значения у_n * оценивается формулой

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

2.2 Исправленный метод Эйлера

В исправленном методе Эйлера мы находим средний тангенс наклона касательнй для двух точек: x_m, y_m и x_m+h, y_m+hy’_m. Последняя точка есть та самая, которая в простом методе обозначалась x_m+1, y_m+1. Геометрический процесс нахождения точки x_m+1, y_m+1 можно проследить по рисунку 2. С помощью метода Эйлера находится точка x_m+h, y_m+hy’_m, лежащая на прямой L_1. В этой точке снова вычисляется тангенс угла наклона касательной, на рисунке этому значению соответствует прямая L₂. Усреднение двух тангенсов дает прямую L’₃. Наконец, через точку x_m, y_m мы проводим прямую L₃ параллельную L’₃. Точка, в которой прямая L₃ пересечется с ординатой, восстановленной из x= x_m+1=x_m+h, и будет искомой точкой y= y_m+1= y_m+hy’_m. Тангенс угла наклона L₃ равен:

Уравнение линии L₃ при всём этом записывается в виде:

Соотношения 5, 6, 7 и 8 описывают исправленный метод Эйлера. (рис. 2)

2.3 Модифицированный метод Эйлера

Этот метод более точен. Рассмотрим дифференциальное уравнение (1) с начальным условием y(x₀)=y₀. Разобьем наш участок интегрирования на n равных частей. На малом участке [x₀,x₀+h] интегральную кривую заменим прямой линией. Получаем точку М_к(х_к,у_к). (рис. 3)

Получаем точку N_k / . В этой точке строим следующую касательную:

Из точки М_к проводим прямую с угловым коэффициентом б_к и определяем точку пересечения этой прямой с прямой x_к1. Получаем точку М_к / . В качестве у_к+1 принимаем ординату точки М_к / . Тогда:

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

Сначала вычисляют вспомогательные значения искомой функции y_к+1/2 в точках x_к+1/2, затем находят значение правой части уравнения (1) в средней точке

3. Блок-схема алгоритма

где A — начальное значение x, B — конечное значение x, F(x) — значение функции в точке x_n, N — количество промежутков, st — выбор операции, C1,C2,C3 — константы для формул, nom — сохраняет номер используемой функции.

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

Подсчитывая каждый раз новое значение уравнения F(x), получаем последовательность значений x_n y_n, n=1,2,…

По этим значениям строим график.

4. Описание программы

Программа весьма проста. В ней много предусмотрено моментов неправильного ввода данных, о которых программа предупреждает пользователя и сразу же просит повторно ввести данные.

С самого начала программа предоставляет пользователю меню выполняемых функций, которые выделяются при помощи стрелок ^ и v выбор клавишей Enter:

После запуска программы нужно выбрать Formyla -> Enter, эта опция позволит из предложенного списка формул выбрать одну, по которой компьютер будет производить расчет и строить график. Все предложенные формулы имеют номерацию; чтобы выбрать интересующий вас пример нажмите на цифру равную номеру примера, и сразу же появится новое окно, в котором сверху будет записан ваш пример. Также в окне будет этот же пример но с нулями на месте констант. Под примером будет высвечена большая буква С, это используется для ввода констант. Для этого вам нужно нажать номер константы, он появится, и после знака равно запишите чему она равна (вводятся целые и вещественные значения). По окончании набора нажать Enter. Операцию повторять пока не будут введены все числа. По окончании нажать Esc. После появится строчка «уточните границы изменения Х, от A= до B= » здесь нужно занести данные на каком промежутке абсциссы будет рассматриваться функция. Следующая строчка попросит ввести начальные данные y(A)=. Следующей строчкой будет вопрос: «сохранить данные в файле? Да/Нет» ответить на этот вопрос с помощью клавиш Д и Н (рус), после чего программа вернется в первоначальное меню. Если данные были сохранены (в папке с программой появляется файл form.txt), то в следующий раз чтобы не набирать снова выберите в меню опцию Formyla -> Open in fails и на экране появятся введенные данные с пометкой снизу, сообщая что данные были прочитаны из файла.

Следующая опция Reshenie. После нажатия в окне просят ввести N(целое число) — число промежутков, на которые разделится рассматриваемый участок (ось ОХ). После появится таблица рассчитанных данных (номер точки, значение абсциссы, значение ординаты). При нажатии любой клавиши произойдет переход в меню.

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

Опция Exit применяется для выхода из программы.

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

Данная программа решает заданную пользователем дифференциальное уравнение за минимальный промежуток времени. При этом пользователю предоставляется возможность визуально оценить решение, рассматривая график полученного решения.

К достоинствам программы можно отнести также удобный пользовательский интерфейс, возможность ввода пользовательских дифференциальных уравнений, а также давольно высокая стабильность работы. При этом имеются и некоторые недостатки. К недостаткам программы можно отнести: критичность к вводимым пользователем урававней, отсутствие обработки исключительных событий. Это, естественно, ограничивает возможности программы.

Список использованной литературы

1. Д. Мак-Кракен, У. Дорн. Численные методы и программирование на фортране. -М.: Мир,1977.-389,396-408 с.

2. А.А. Самарский. ВВЕДЕНИЕ в численные методы. — М.:Наука,1987.-176 с.

3. Алгоритмы вычислительной математики: Лабораторный практикум по курсу «Программирование» для студентов 1 — 2-го курсов всех специальностей БГУИР/А.К. Синицын, А.А. Навроцкий.- Мн.: БГУИР, 2002.- 65-69 с.

4. ГОСТ 2.105-95. Общие требования к текстовым документам.

5. ГОСТ 7.32-91. Система стандартов по информации, библиотечному и издательскому делу. Отчет о НИР. Структура и правила оформления.

Приложение 1. Текст программы.

void formyl(int p)

if(p==1) printf(«\n 1. C1*y’ = C2*y + C3*x + C4*x*y»);

else if(p==2) printf(«\n 2. y’/(C1-100) = C2*y + C3*x + (C4+x)*y»);

else if(p==3) printf(«\n 3. pow(e,C1)*y’ = C2*y + C3*cos(x) + (C4+x+y)»);

else if(p==4) printf(«\n 4. C1*sin(x)*y’ = e*C2*y + C3*arcsin(x) + C4*y/x»);

else if(p==5) printf(«\n 5. C1*y’ = sin(C2)*y + tg(C3*x) + C4*ln(x)*y»);

else if(p==6) printf(«\n 6. C1*y’ = y*C2 + C3*sin(x) + C4*cos(x)*y»);

else if(p==7) printf(«\n 7. (C1+C2+C3+C4)*y’ = C2*y + (C3-x) + lg(C4*x)*y»);

else if(p==8) printf(«\n 8. y’/C1 = y/C2 + C3*sin(x) + C4*x*y»);

else if(p==9) printf(«\n 9. sin(C1)*y’ = C2*y + |C3|*x + x*y/C4»);

void formyl2(int p,double C1,double C2,double C3,double C4)

else if(p==2) printf(«y’/(%.2f-100)=%.2f*y+%.2f*x+(%.2f+x)*y»,C1,C2,C3,C4);else if(p==3) printf(«pow(e,%.2f)*y’=%.2f*y+%.2f*cos(x)+(%.2f+x+y)»,C1,C2,C3,C4);

else if(p==4) printf(«%.2f*sin(x)*y’=e*%.2f*y+%.2f*arcsin(x)+%.2f*y/x»,C1,C2,C3,C4);

else if(p==5) printf(«%.2f*y’=sin(%.2f)*y+tg(%.2f*x)+%.2f*ln(x)*y»,C1,C2,C3,C4);

else if(p==6) printf(«%.2f*y’=y*%.2f+%.2f*sin(x)+%.2f*cos(x)*y»,C1,C2,C3,C4);

else if(p==8) printf(«y’/%.2f=y/%.2f+%.2f*sin(x)+%.2f*x*y»,C1,C2,C3,C4);

else if(p==9) printf(«sin(%.2f)*y’=%.2f*y+|%.2f|*x+x*y/%.2f»,C1,C2,C3,C4);

double formyl3(int p,double h,double x,double y,double C1,double C2,double C3,double C4)

else if(p==2) Y=h*(C1-100)*(y*C2+C3*x+(C4+x)*y)+y;

else if(p==3) Y=h*(C2*y+C3*cos(x)+C4+x+y)/exp(C1)+y;

else if(p==4) Y=h*(exp(1)*C2*y+C3*asin(x)+C4*y/x)/(C1*sin(x))+y;

else if(p==5) Y=h*(sin(C2)*y+tan(C3*x)+C4*log10(x)*y)/C1+y;

else if(p==6) Y=h*(y*C2+C3*sin(x)+C4*cos(x)*y)/C1+y;

else if(p==7) Y=h*(C2*y+(C3-x)+log10(C4*x)*y)/(C1+C2+C3+C4)+y;

else if(p==8) Y=h*(y/C2+C3*sin(x)+C4*x*y)*C1+y;

else if(p==9) Y=h*(C2*y+abs(C3)*x+x*y/C4)/sin(C1)+y;

int vv=0,vv1=0; // руководит операциями

int N=0,W; // кол промежутков

int i,j,k; // используются во всех «for»

int nom; // номер примера

int st=4,vst=0; // строчка в меню

double C1,C2,C3,C4; // константы

double M; // масштаб

double xtoch,ytoch; // считает y(x) по графику

double A=0,B=0,ii,jj,kk; // пределы интегрирования

double x[102],y[102]; // главные переменные x,y

int g_driver=9,g_mode=2, g_error;

printf(«\n error=%d, reason=%s\n», g_error, grapherrormsg(g_error));