Реферат: Решение системы нелинейных уравнений - текст реферата. Скачать бесплатно.
Банк рефератов, курсовых и дипломных работ. Много и бесплатно. # | Правила оформления работ | Добавить в избранное
 
 
   
Меню Меню Меню Меню Меню
   
Napishem.com Napishem.com Napishem.com

Реферат

Решение системы нелинейных уравнений

Банк рефератов / Математика

Рубрики  Рубрики реферат банка

закрыть
Категория: Реферат
Язык реферата: Русский
Дата добавления:   
 
Скачать
Архив Zip, 14 kb, скачать бесплатно
Обойти Антиплагиат
Повысьте уникальность файла до 80-100% здесь.
Промокод referatbank - cкидка 20%!

Узнайте стоимость написания уникальной работы





Теоретическая часть.


В данной расчетно-графической работе (далее РГР) требуется составить программу для решения системы нелинейных уравнений методом последовательной итерации обратной матрицы Якоби.

Суть метода в следующем:

Пусть требуется решить систему нелинейных алгебраических или трансцендентных уравнений:


F1(X1,X2,...,Xn)=0; i=1,2,...,n,


с начальным приближением к решению:


X0=(x10,x20,...xn0).


Вычислительная схема реализованного метода состоит в следующем:

В начале итерационного процесса матрица H полагается равной единичной:

H0=E.

Затем для k=0,1,...


1. Вычисляется

Pk = - Hk * F(Xk);

2. Находятся

Xk+1 = Xk + tk*Pk.


Первоначально tk=1. Затем путем последовательного деления tk на 2 находим такое tk, чтобы выполнялось неравенство:


¦ F(Xk+1) ¦ < ¦ F(Xk) ¦


Итерационный процесс заканчивается при выполнении условия:

¦ F(Xk+1) ¦ < E,

где E - заданная точность.


3. Определяется

Yk= F(Xk+1) - F(Xk)


4. Находится новое приближение матрицы:


Hk+1 = Hk - (Hk*Yk - Pk*tk) * (Pk)T * (Hk)T / ((Pk)T * Hk*Yk)


и снова повторяется вычислительный процесс с пункта 1.


Порядок работы с программой


Данная РГР представлена в виде 3 исполняемых модулей:

OBRJ.M, OBRF.M и FUN1.M. Решением поставленной задачи занимается модуль OBRF.M, а два остальных являются вспомогательными:

OBRJ.M - головной модуль, в котором вводятся входные данные и выводятся результаты вычислений, а FUN1.M - модуль, который пишет сам пользователь и который возвращает вычисленные левые части для требуемого уравнения.

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


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


В связи с тем, что данная РГР состоит из 3 частей, то опишем их по одиночке (распечатки данных модулей приведены в приложении):


1. OBRJ.M

Головной модуль

Входные данные: отсутствуют.

Выходные данные: отсутствуют.

Язык реализации: PC MathLab.

Операционная система: MS-DOS 3.30 or Higher.

Пояснения к тексту модуля:

"Стандартный" головной модуль. В данном модуле задаются начальные значения в виде вектора, например:

X0=[0.4 0.9]

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

Затем происходит вызов модуля OBRF.M с полученными входными данными. Формат вызова данного модуля описан далее (в описании самого модуля).

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


2. OBRF.M

Вычислительный модуль

Входные данные:


FunFcn - имя функции, написанной пользователем, которая вычисляет левые части для требуемой системы в определенной точке.

X0 - вектор-строка, определяющий начальные значения (начальное приближение).

E - допустимая ошибка.


Выходные данные:


Tout - Столбец итераций ("Время")

Xout - Столбцы значений вычисленных на каждом этапе для каждой итерации

DXout - Столбцы погрешностей по каждой компоненте, вычисленные на определенном этапе


Язык реализации: PC MathLab

Операционная система: MS-DOS 3.30 or Higher

Пояснения к тексту модуля:

Данный "вычислительный" модуль реализует метод последовательной итерации обратной матрицы Якоби. Общая структура вызова данного модуля:


[Tout,Xout,DXout]=OBRF(FunFcn,X0,E);


Значения каждого из параметров были описаны выше.

На начальном этапе в данном модуле инициализируются внутренние переменные (например, задается единичная матрица H, в соответствии с размерностью X0), формируются (на основе начальных значений) первичные элементы матриц Tout,Xout,DXout.

Затем данная функция, как и многие другие в численных методах,имеет вид:


While ОШИБКА > ДОПУСТИМОЙ ОШИБКИ

Оператор1

Оператор2

.........

.........

ОператорN

End


Внутри данного цикла происходят вычисления внутренней переменной Pk на каждом шаге K и, вычисляется начальное приближение Xk+1. Первоначально t=1 (Не номер итерации, а внутренний параметр!). Затем, в очередном цикле While...End в случае, если ¦F(Xk+1)¦ < ¦F(Xk)¦ t=t/2 и снова вычисляется Xk+1. Когда очередное Xk+1 найдено, вычисляется Yk, а затем и новое приближение матрицы H. Итерационный процесс заканчивается, если ¦F(Xk+1)¦ < E. Если данное условие не выполняется - итерационный процесс продолжается заново.

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


3. FUN1.M

Модуль, вычисляющий левые части

Входные данные:

X - вектор-строка, задающий точки для вычислений по каждой компоненте.

Выходные данные:

FF - вектор-строка, возвращающий значения каждой компоненты в определенной точке

Язык реализации: PC MathLab

Операционная система: MS-DOS 3.30 or Higher

Пояснения к тексту модуля:

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


Сравнительный анализ и

оценка быстродействия.



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


1. Метод последовательной итерации обратной матрицы Якоби Число операций: порядка 682

Быстродействие: порядка 0.11 секунды


2. Метод Ньютона дискретный

Число операций: порядка 990

Быстродействие: порядка 0.22 секунды


Как видно из вышеприведенных данных, эти два метода очень близки между собой, но метод Ньютона дискретный более сложен в реализации, однако обладает лучшей сходимостью, например при начальных значениях X0=[2.0 2.0]; метод последовательной итерации обратной матрицы Якоби уже не справляется, в то время как дискретный метод Ньютона продолжает неплохо работать. Однако метод Ньютона требует больших затрат машинного времени и поэтому при выборе метода необходимо исходить их конкретных условий задачи и если известно довольно точное приближение и требуется быстрота вычислений, то к таким условиям отлично подходит разработанный метод последовательной итерации обратной матрицы Якоби.


Выводы


В данной РГР был разработан и реализован метод последовательной итерации обратной матрицы Якоби, предназначенный для решения системы нелинейных уравнений. Программа, реализованная на языке PC MathLab хотя и не является оптимальной, однако выполняет поставленную задачу и решает системы уравнений. Реализованный метод не отличается повышенной сходимостью и требует довольно точного начального приближения, однако довольно быстро сходится к точному решению, то есть его можно порекомендовать для вычисления непростых систем нелинейных уравнений при наличии довольно точного начального приближения и наличия временных ограничений.


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


1. О.М.Сарычева. "Численные методы в экономике. Конспект лекций", Новосибирский государственный технический университет, Новосибирск 1995г.

2. Д.Мак-Кракен, У.Дорн. "Численные методы и программирование на Фортране", Издательство "Мир", М. 1977г.

3. Н.С.Бахвалов. "Численные методы", Издательство "Наука", М. 1975г.

1Авиация и космонавтика
2Архитектура и строительство
3Астрономия
 
4Безопасность жизнедеятельности
5Биология
 
6Военная кафедра, гражданская оборона
 
7География, экономическая география
8Геология и геодезия
9Государственное регулирование и налоги
 
10Естествознание
 
11Журналистика
 
12Законодательство и право
13Адвокатура
14Административное право
15Арбитражное процессуальное право
16Банковское право
17Государство и право
18Гражданское право и процесс
19Жилищное право
20Законодательство зарубежных стран
21Земельное право
22Конституционное право
23Конституционное право зарубежных стран
24Международное право
25Муниципальное право
26Налоговое право
27Римское право
28Семейное право
29Таможенное право
30Трудовое право
31Уголовное право и процесс
32Финансовое право
33Хозяйственное право
34Экологическое право
35Юриспруденция
36Иностранные языки
37Информатика, информационные технологии
38Базы данных
39Компьютерные сети
40Программирование
41Искусство и культура
42Краеведение
43Культурология
44Музыка
45История
46Биографии
47Историческая личность
 
48Литература
 
49Маркетинг и реклама
50Математика
51Медицина и здоровье
52Менеджмент
53Антикризисное управление
54Делопроизводство и документооборот
55Логистика
 
56Педагогика
57Политология
58Правоохранительные органы
59Криминалистика и криминология
60Прочее
61Психология
62Юридическая психология
 
63Радиоэлектроника
64Религия
 
65Сельское хозяйство и землепользование
66Социология
67Страхование
 
68Технологии
69Материаловедение
70Машиностроение
71Металлургия
72Транспорт
73Туризм
 
74Физика
75Физкультура и спорт
76Философия
 
77Химия
 
78Экология, охрана природы
79Экономика и финансы
80Анализ хозяйственной деятельности
81Банковское дело и кредитование
82Биржевое дело
83Бухгалтерский учет и аудит
84История экономических учений
85Международные отношения
86Предпринимательство, бизнес, микроэкономика
87Финансы
88Ценные бумаги и фондовый рынок
89Экономика предприятия
90Экономико-математическое моделирование
91Экономическая теория

 Анекдоты - это почти как рефераты, только короткие и смешные Следующий
А что если снег все-таки выпадает, но его по ночам собирают и отвозят в Сочи?..
Anekdot.ru

Узнайте стоимость курсовой, диплома, реферата на заказ.


Банк рефератов - РефератБанк.ру
© РефератБанк, 2002 - 2017
Рейтинг@Mail.ru