Банк рефератов содержит более 364 тысяч рефератов, курсовых и дипломных работ, шпаргалок и докладов по различным дисциплинам: истории, психологии, экономике, менеджменту, философии, праву, экологии. А также изложения, сочинения по литературе, отчеты по практике, топики по английскому.
Полнотекстовый поиск
Всего работ:
364150
Теги названий
Разделы
Авиация и космонавтика (304)
Административное право (123)
Арбитражный процесс (23)
Архитектура (113)
Астрология (4)
Астрономия (4814)
Банковское дело (5227)
Безопасность жизнедеятельности (2616)
Биографии (3423)
Биология (4214)
Биология и химия (1518)
Биржевое дело (68)
Ботаника и сельское хоз-во (2836)
Бухгалтерский учет и аудит (8269)
Валютные отношения (50)
Ветеринария (50)
Военная кафедра (762)
ГДЗ (2)
География (5275)
Геодезия (30)
Геология (1222)
Геополитика (43)
Государство и право (20403)
Гражданское право и процесс (465)
Делопроизводство (19)
Деньги и кредит (108)
ЕГЭ (173)
Естествознание (96)
Журналистика (899)
ЗНО (54)
Зоология (34)
Издательское дело и полиграфия (476)
Инвестиции (106)
Иностранный язык (62792)
Информатика (3562)
Информатика, программирование (6444)
Исторические личности (2165)
История (21320)
История техники (766)
Кибернетика (64)
Коммуникации и связь (3145)
Компьютерные науки (60)
Косметология (17)
Краеведение и этнография (588)
Краткое содержание произведений (1000)
Криминалистика (106)
Криминология (48)
Криптология (3)
Кулинария (1167)
Культура и искусство (8485)
Культурология (537)
Литература : зарубежная (2044)
Литература и русский язык (11657)
Логика (532)
Логистика (21)
Маркетинг (7985)
Математика (3721)
Медицина, здоровье (10549)
Медицинские науки (88)
Международное публичное право (58)
Международное частное право (36)
Международные отношения (2257)
Менеджмент (12491)
Металлургия (91)
Москвоведение (797)
Музыка (1338)
Муниципальное право (24)
Налоги, налогообложение (214)
Наука и техника (1141)
Начертательная геометрия (3)
Оккультизм и уфология (8)
Остальные рефераты (21697)
Педагогика (7850)
Политология (3801)
Право (682)
Право, юриспруденция (2881)
Предпринимательство (475)
Прикладные науки (1)
Промышленность, производство (7100)
Психология (8694)
психология, педагогика (4121)
Радиоэлектроника (443)
Реклама (952)
Религия и мифология (2967)
Риторика (23)
Сексология (748)
Социология (4876)
Статистика (95)
Страхование (107)
Строительные науки (7)
Строительство (2004)
Схемотехника (15)
Таможенная система (663)
Теория государства и права (240)
Теория организации (39)
Теплотехника (25)
Технология (624)
Товароведение (16)
Транспорт (2652)
Трудовое право (136)
Туризм (90)
Уголовное право и процесс (406)
Управление (95)
Управленческие науки (24)
Физика (3463)
Физкультура и спорт (4482)
Философия (7216)
Финансовые науки (4592)
Финансы (5386)
Фотография (3)
Химия (2244)
Хозяйственное право (23)
Цифровые устройства (29)
Экологическое право (35)
Экология (4517)
Экономика (20645)
Экономико-математическое моделирование (666)
Экономическая география (119)
Экономическая теория (2573)
Этика (889)
Юриспруденция (288)
Языковедение (148)
Языкознание, филология (1140)

Курсовая работа: Визуализация численных методов

Название: Визуализация численных методов
Раздел: Рефераты по математике
Тип: курсовая работа Добавлен 09:18:05 30 мая 2010 Похожие работы
Просмотров: 221 Комментариев: 2 Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать

Сибирский государственный университет телекоммуникации и информатики

Уральский технический институт связи и информатики

Факультет телекоммуникации, информатики и управления

Кафедра организации управления связи

По курсу: “Информатика”

По теме: “Визуализация численных методов”

Написал:

Плишкин М. Ю

группа МЕ-72

Преподаватель:

Кандидат технических наук , доцент

Е.Е.Минина

г. Екатеринбург. 2010 г.


Содержание

Введение

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

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

1.2 Метод Рунге – Кутта

2. Блок-схемы

3. Виды, формы

3.1 Начальная форма

3.2 Конечная форма

4. Программа для решения дифференциального уравнения в Visual Basic

Заключение


Введение

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

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

Чтобы решить ОДУ, необходимо знать значение зависимой переменной и (или) её производные при некоторых значениях независимой переменной. Если эти дополнительные условия задаются при одном значении независимой переменной, то такая задача называется задачей с начальными условиями, или задачей Коши.

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

Одношаговые методы, в которых для нахождения следующей точки на кривой y=f(x) требуется информация лишь об одном предыдущем шаге. Одношаговыми являются метод Эйлера и методы Рунге – Кутта.

Методы прогноза и коррекции (многошаговые), в которых для отыскивания следующей точки кривой y=f(x) требуется информация более чем об одной из предыдущих точек. Чтобы получить достаточно точное численное значение, часто прибегают к итерации. К числу таких методов относятся методы Милны, Адамса – Башфорта и Хемминга.

Явные методы, в которых функция Ф в выражении (1) не зависит от yn+1.

Неявные методы, в которых функция Ф зависит от yn+1.

В данной курсовой работе будут рассматриваться два одношаговых метода: метод Эйлера первого порядка точности и Рунге – Кутта четвёртого порядка точности.


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

В данной курсовой работе необходимо решить ОДУ вида y` = 4y/x с заданными начальными значениями x0=1, xk=1.4, y0=2, h=0.05. Для проверки точности результатов дано общее решение данного уравнения y=x^4с. Требуется решить уравнение двумя методами: Эйлера модифицированного и Рунге – Кутта четвёртого порядка, сравнить результаты и сделать вывод какой метод эффективнее использовать, построить графики.

Численное решение задачи Коши сводится к табулированию искомой функции.

График решения дифференциального уравнения называется интегральной кривой.

Геометрический смысл задачи:

y`=f(x,y) – тангенс угла наклона касательной к графику решения в точке (x,y) к оси OX (угловой коэффициент (в общей формуле прямой,

y=k*x+b,

обозначается как “k”)(рис 1).

Рисунок 1. Геометрический смысл задачи Коши

Существующие решения:

Если правая часть f(x,y) непрерывная в некоторой области R, определяемой неравенствами |x – x0| < a; |y – y0| > b, то существует, по меньшей мере, одно решение y=y(x), определённое в окрестности |x – x0| < h, где h > 0.

При использовании численных методов выполняется замена отрезка [x0,X] – области непрерывного изменения аргумента x множеством wh, состоящего из конечного числа точек x0<x1<...<xn=X – сеткой.

При этом xi называют узлами решётки.

Задача Коши, определённая ранее на непрерывном отрезке [x0,X], заменяется её дискретным аналогом – системой уравнений, решая которую можно последовательно найти значения y1,y2,...,yn – приближённые значения функции в узлах сетки.

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

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

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

y`=f(x,y)

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

y(x0)=y0.

Выберем шаг h и введём обозначения:

xi=x0+i*h и yi=y(xi), где i=0,1,2,...,


xi- узлы сетки,

yi- значение интегральной функции в узлах.

Иллюстрации к решению приведены на рисунке 2.

Проведём прямую АВ через точку (xi,yi) под углом α. При этом

tgα=f(xi,yi) (1)

В соответствии с геометрическим смыслом задачи, прямая АВ является касательной к интегральной функции. Произведём замену точки интегральной функции точкой, лежащей на касательной АВ.

Тогда

y i+1=yi+Δy (2).

Из прямоугольного треугольника АВС

tgα= Δy/h (3).

Приравниваем правые части (1) и (3). Получим

Δy/h= f(xi,yi).

Отсюда

Δy= f(xi,yi)*h.

Подставим в это выражение формулу (2), а затем преобразуем его. В результате получаем формулу расчёта очередной точки интегральной функции:


y i+1=yi+ h*f (xi,yi) (4).

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

Рисунок 2. Метод Эйлера

Метод Эйлера – один из простейших методов численного решения ОДУ. Но существенным его недостатком является большая погрешность вычислений. На рисунке 2 погрешность вычислений дляi-го шага обозначена ε. С каждым шагом погрешность вычислений увеличивается.

1.2 Метод Рунге – Кутта

Пусть дано дифференциальное уравнение первого порядка y`=f(x,y) с начальным условием y(x0)=y0.

Выберем шаг h и введём обозначения:

xi=x0+i*h и yi=y(xi), где i=0,1,2,...,

xi- узлы сетки,

yi- значение интегральной функции в узлах.

Проведём решение в несколько этапов.

Обозначим точки: A(xi,yi), B(xi+1,yi+1), C, D, E.

Через точку А проведём прямую под углом α, где tg α = f(xi,yi).

На прямой (1) найдём точку С. Через точку С проведём прямую под углом α1, где

tg α1 = f(xi+h/4, yi+h/4*f(xi,yi).

Через точку А проведём прямую параллельную последней прямой.

Найдём точку D на прямой (2) и через неё проведём прямую под углом α2, где

tg α2 = f(xi+h/2, yi+h/2*f(xi,yi)).

Через точку А проведём прямую параллельную последней прямой.

По примеру, описанному выше, построим прямую, которая пересечётся с прямой x = xi+1. Эта точка и будет решением дифференциального уравнения при x = xi+1.

Согласно методу Рунге – Кутта четвёртого порядка, последовательные значения yi искомой функции y определяется по формуле:

y i+1=yi+Δy,

где

Δy=(k1+2*k2+2*k3+k4)/6, i=0,1,2,...


а числа k1(i),k2(i),k3(i),k4(i) на каждом шаге вычисляются по формулам:

k1=h*f(xi,yi)

k2 =h*f(xi+h/2,yi+k1/2)

k3=h*f(xi+h/2,yi+k2/2)

k4 =h*f(xi+h,yi+k3)

Это явный четырёхэтапный метод четвёртого порядка точности.


2. Блок-схемы





3. Виды, формы

3.1 Начальная форма

3.2 Конечная форма


4. Программа для решения дифференциального уравнения в Visual Basic

Dim y(9) As Single

Dim YE(9) As Single

Dim YR(9) As Single

Dim YT(9) As Single

Dim l(9) As Single

Private x0 As Single

Private Function fun(a As Single, b As Single) As Single

f = (b + 2) / (a + 1)

fun = f

End Function

Private Sub Command1_Click()

x0 = Val(Text1.Text)

xk = Val(Text2.Text)

y0 = Val(Text3.Text)

h = Val(Text4.Text)

N = (xk - x0) / h

MSFlexGrid1.Rows = N + 2

MSFlexGrid1.TextMatrix(0, 0) = "x"

MSFlexGrid1.TextMatrix(0, 1) = "YE"

MSFlexGrid1.TextMatrix(0, 2) = "YR"

MSFlexGrid1.TextMatrix(0, 3) = "YT"

Min = y0

Max = y0

l(0) = x0

y(0) = y0

YE(0) = y0

YR(0) = y0

YT(0) = y0

For i = 0 To N

l(i + 1) = x0 + i * h

k1 = h * fun(l(i), YR(i))

k2 = h * fun(l(i) + h / 2, YR(i) + k1 / 2)

k3 = h * fun(l(i) + h / 2, YR(i) + k2 / 2)

k4 = h * fun(l(i) + h, YR(i) + k3)

k = (k1 + 2 * k2 + 2 * k3 + k4) / 6

YR(i + 1) = YR(i) + k

YE(i + 1) = YE(i) + h * fun(l(i), YE(i))

YT(i) = (l(i + 1) + 1) * 2 - 2

MSFlexGrid1.TextMatrix(i + 1, 0) = l(i + 1)

MSFlexGrid1.TextMatrix(i + 1, 1) = YE(i)

MSFlexGrid1.TextMatrix(i + 1, 2) = YR(i)

MSFlexGrid1.TextMatrix(i + 1, 3) = YT(i)

If YE(i) > Max Then max1 = YE(i)

If YE(i) < Min Then min1 = YE(i)

If YR(i) > Max Then Max2 = YR(i)

If YR(i) < Min Then Min2 = YR(i)

If YT(i) > Max Then Max3 = YT(i)

If YT(i) < Min Then Min3 = YT(i)

Next i

For i = 0 To N - 1

px = (5415 / (xk - x0))

py = (6705 / (max1 - min1))

u1 = (l(i) - x0) * px + 600

u2 = 7440 - (YE(i) - min1) * py

u3 = (l(i + 1) - x0) * px + 600

u4 = 7440 - (YE(i + 1) - min1) * py

Picture1.Line (u1, u2)-(u3, u4)

u5 = (l(i) - x0) * px + 600

u6 = 7440 - (YR(i) - min1) * py

u7 = (l(i + 1) - x0) * px + 600

u8 = 7440 - (YR(i + 1) - min1) * py

Picture1.Line (u5, u6)-(u7, u8)

u9 = (l(i) - x0) * px + 600

u10 = 7440 - (YT(i) - min1) * py

u11 = (l(i + 1) - x0) * px + 600

u12 = 7440 - (YT(i + 1) - min1) * py

Picture1.Line (u9, u10)-(u11, u12)

Next i

End Sub


Заключение

В данной курсовой рассматривались два метода решения ОДУ с начальными условиями, то есть задачи Коши: метод Эйлера и метод Рунге – Кутта четвёртого порядка.

Данные полученные этими методами идентичны друг другу, но с точки зрения простоты использования метод Эйлера гораздо проще в описании, чем метод Рунге – Кутта четвертого порядка.

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

Оценить/Добавить комментарий
Имя
Оценка
Комментарии:
Где скачать еще рефератов? Здесь: letsdoit777.blogspot.com
Евгений07:17:33 19 марта 2016
Кто еще хочет зарабатывать от 9000 рублей в день "Чистых Денег"? Узнайте как: business1777.blogspot.com ! Cпециально для студентов!
14:22:30 25 ноября 2015

Работы, похожие на Курсовая работа: Визуализация численных методов

Назад
Меню
Главная
Рефераты
Благодарности
Опрос
Станете ли вы заказывать работу за деньги, если не найдете ее в Интернете?

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



Результаты(150738)
Комментарии (1839)
Copyright © 2005-2016 BestReferat.ru bestreferat@mail.ru       реклама на сайте

Рейтинг@Mail.ru