Размещено на http://www.allbest.ru/

Министерство образования и науки РФ

Федеральное Государственное Бюджетное Образовательное Учреждение Высшего Профессионального Образования

Тульский Государственный Университет

Кафедра ССМиК

Отчет по курсовой работе

на тему «Численные методы решения нелинейных уравнений»

по дисциплине «Информатика»

Выполнил ст. гр. 321731

Храмеев Максим Александрович

Проверил к.т.н., доц. Теличко В.Г.

Тула 2013

• Содержание
• Цель работы
• Решение нелинейного уравнения
• Использование программы Excel для решения нелинейных уравнений. Отделение корней
• Уточнение корней методом хорд
• Уточнение корней методом касательных (Ньютона)
• Уточнение корней методом простой итерации
• Циклические ссылки
• Подбор параметра
• Поиск решения
• Решение нелинейного уравнения методом «Подбор параметра»
• Список использованных источников информации
• Цель работы
• нелинейный уравнение корень подбор
• Целью заключается освоение методики решения нелинейных уравнений на компьютере в программе Excel.

Решение нелинейного уравнения

Найти корень нелинейного уравнения f1(x)=f2(x) на заданном отрезке [a,b] средствами Excel тремя возможными способами:

1. методом касательных, либо методом простой итерации с использованием циклических ссылок;

2. с помощью средства Подбор параметра

3. используя возможности Поиска решения при ограничениях корень>a и корень <b.

2x+x5-1=0, x(0;1]

Использование программы Excel для решения нелинейных уравнений. Отделение корней

В общем случае отделение корней уравнения f(x)=0 базируется на известной теореме, утверждающей, что если непрерывная функция f(x) на концах отрезка [a,b] имеет значения разных знаков, в общем случае выбирают некоторый диапазон, где могут обнаружиться корни с выбранным шагом h для обнаружения перемены знаков f(x).

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

Уточнение корней методом хорд

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

Рис.2. Метод хорд

Здесь вычисляются значения функции на концах отрезка, и строится «хорда», соединяющая точки (a,f(a)) и (b,f(b)). Точка пересечения ее с осью абсцисс

принимается за очередное приближение к корню. Анализируя знак f(z) в сопоставлении со знаком f(x) на концах отрезка, сужаем интервал до [a,z] или [z,b] и продолжаем процесс построения хорд до тех пор, пока разница между очередными приближениями не окажется достаточно малой (в пределах допустимой погрешности) |Zn-Zn-1|<e--.

Можно доказать, что истинная погрешность найденного приближения:

,

где X* - корень уравнения, Zn и Zn+1 - очередные приближения, m и M - наименьшее и наибольшее значения f(x) на интервале [a,b].

Алгоритм метода хорд

Уточнение корней методом касательных (Ньютона)

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

Наиболее популярным из итерационных методов является метод Ньютона (метод касательных).

Рис. 3. Метод касательных

Пусть известно некоторое приближенное значение Zn корня X*. Применяя формулу Тейлора и ограничиваясь в ней двумя членами, имеем

,

откуда

.

Геометрически этот метод предлагает построить касательную к кривой y=f(x) в выбранной точке x=Zn, найти точку пересечения её с осью абсцисс и принять эту точку за очередное приближение к корню (3).

Размещено на http://www.allbest.ru/

Рис. 4. Расходящийся процесс

Размещено на http://www.allbest.ru/

Рис.5. Приближение к другому корню

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

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

Если производная функции мало изменяется в окрестности корня, то можно использовать видоизменение метода

.

Существуют и другие модификации метода Ньютона.

Уточнение корней методом простой итерации

Другим представителем итерационных методов является метод простой итерации.

Здесь уравнение f(x)=0 заменяется равносильным уравнением и строится последовательность значений

.

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Рис. 6 Геометрическая интерпретация метода простой итерации

Отсюда может быть построен итерационный процесс

.

Возьмем для примера уравнение x3 + x -1000 = 0. Очевидно, что корень данного уравнения несколько меньше 10. Если переписать это уравнение в виде x =1000 - x3 и начать итерационный процесс при x0=10, то из первых же приближений очевидна его расходимость. Если же учесть f '(x)=3x2+1>0 и принять за приближенное значение максимума f '(x) M=300, то можно построить сходящийся итерационный процесс на основе представления

.

Можно и искусственно подобрать подходящую форму уравнения, например:

или .

Циклические ссылки

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

Рассмотрим задачу нахождения корня уравнения методом Ньютона с использованием циклических ссылок. Возьмем для примера квадратное уравнение: х2 - 5х + 6=0, графическое представление которого приведено на 7 Найти корень этого (и любого другого) уравнения можно, используя всего одну ячейку Excel.

Для включения режима циклических вычислений в меню Сервис/Параметры/вкладка Вычисления включаем флажок Итерации, при необходимости изменяем число повторений цикла в поле Предельное число итераций и точность вычислений в поле Относительная погрешность (по умолчанию их значения равны 100 и 0,0001 соответственно). Кроме этих установок выбираем вариант ведения вычислений: автоматически или вручную. При автоматическом вычислении Excel выдает сразу конечный результат, при вычислениях, производимых вручную, можно наблюдать результат каждой итерации.

Рис. 7. График функции

Выберем произвольную ячейку, присвоим ей новое имя, скажем - Х, и введем в нее рекуррентную формулу, задающую вычисления по методу Ньютона:

,

где F и F1 задают соответственно выражения для вычисления значений функции и ее производной. Для нашего квадратного уравнения после ввода формулы в ячейке появится значение 2, соответствующее одному из корней уравнения (рис. 8). В нашем случае начальное приближение не задавалось, итерационный вычислительный процесс начинался со значения, по умолчанию хранимого в ячейке Х и равного нулю. А как получить второй корень? Обычно это можно сделать изменением начального приближения. Решать проблему задания начальных установок в каждом случае можно по-разному. Мы продемонстрируем один прием, основанный на использовании функции ЕСЛИ. С целью повышения наглядности вычислений ячейкам были присвоены содержательные имена (рис 8).

· В ячейку Хнач (В4) заносим начальное приближение - 5.

· В ячейку Хтекущ (С4) записываем формулу:

· =ЕСЛИ(Хтекущ=0;Хнач; Хтекущ-(Хтекущ^2-5*Хтекущ+6)/(2*Хтекущ-5)).

· В ячейку D4 помещаем формулу, задающую вычисление значения функции в точке Хтекущ, что позволит следить за процессом решения.

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

· Чтобы сменить начальное приближение, недостаточно изменить содержимое ячейки Хнач и запустить процесс вычислений. В этом случае вычисления будут продолжены, начиная с последнего вычисленного значения.

Размещено на http://www.allbest.ru/

Рис. 8. Определение начальных установок

Чтобы обнулить значение, хранящееся в ячейке Хтекущ, нужно заново записать туда формулу. Для этого достаточно для редактирования выбрать ячейку, содержащую формулу, дважды щелкнув мышью на ней (при этом содержимое ячейки отобразится в строке формул). Щелчок по кнопке (нажатие клавиши) Enter запустит вычисления с новым начальным приближением.

Подбор параметра

Когда желаемый результат вычислений по формуле известен, но неизвестны значения, необходимые для получения этого результата, можно воспользоваться средством Подбор параметра, выбрав команду Подбор параметра в меню Сервис. При подборе параметра Excel изменяет значение в одной конкретной ячейке до тех пор, пока вычисления по формуле, ссылающейся на эту ячейку, не дадут нужного результата. Возьмем в качестве примера все то же квадратное уравнение х2-5х+6=0. Для нахождения корней уравнения выполним следующие действия:

* В ячейку С3 (рис. 10) введем формулу для вычисления значения функции, стоящей в уравнении слева от знака равенства. В качестве аргумента используем ссылку на ячейку С2, т.е. =С2^2-5*C2+6.

Размещено на http://www.allbest.ru/

* В окне диалога Подбор параметра (рис. 10) в поле Установить в ячейке введем ссылку на ячейку с формулой, в поле Значение - ожидаемый результат, в поле Изменяя значения ячейки - ссылку на ячейку, в которой будет храниться значение подбираемого параметра (содержимое этой ячейки не может быть формулой).

* После нажатия на кнопку Ok Excel выведет окно диалога Результат подбора параметра. Если подобранное значение необходимо сохранить, то нажмите на Оk, и результат будет сохранен в ячейке, заданной ранее в поле Изменяя значения ячейки. Для восстановления значения, которое было в ячейке С2 до использования команды Подбор параметра, нажмите кнопку Отмена. При подборе параметра Excel использует итерационный (циклический) процесс. Количество итераций и точность устанавливаются в меню Сервис/Параметры/вкладка Вычисления. Если Excel выполняет сложную задачу подбора параметра, можно нажать кнопку Пауза в окне диалога Результат подбора параметра и прервать вычисление, а затем нажать кнопку Шаг, чтобы выполнить очередную итерацию и просмотреть результат. При решении задачи в пошаговом режиме появляется кнопка Продолжить - для возврата в обычный режим подбора параметра.

Вернемся к примеру. Опять возникает вопрос: как получить второй корень? Как и в предыдущем случае необходимо задать начальное приближение. Это можно сделать следующим образом (рис. 11,а):

А)

Б)

Рис. 11. Поиск второго корня

* В ячейку Х (С2) вводим начальное приближение.

* В ячейку Хi (С3) вводим формулу для вычисления очередного приближения к корню, т.е. =X-(X^2-5*X+6)/(2*X-5).

* В ячейку С4 поместим формулу, задающую вычисление значения функции, стоящей в левой части исходного уравнения, в точке Хi.

* После этого выбираем команду Подбор параметра, где в качестве изменяемой ячейки принимаем ячейку С2Результат вычислений изображен на рис. 11,б (в ячейке С2 - конечное значение, а в ячейке С3 предыдущее). Однако все это можно сделать и несколько проще. Для того чтобы найти второй корень, достаточно в качестве начального приближения (рис. 10) в ячейку C2 поместить константу 5 и после этого запустить процесс Подбор параметра.

Поиск решения

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

Искомые переменные - ячейки рабочего листа Excel - называются регулируемыми ячейками. Целевая функция F(х1, х2, … , хn), называемая иногда просто целью, должна задаваться в виде формулы в ячейке рабочего листа. Эта формула может содержать функции, определенные пользователем, и должна зависеть (ссылаться) от регулируемых ячеек. В момент постановки задачи определяется, что делать с целевой функцией. Возможен выбор одного из вариантов:

· найти максимум целевой функции F(х1, х2, … , хn);

· найти минимум целевой функции F(х1, х2, … , хn);

· добиться того, чтобы целевая функция F(х1, х2, … , хn) имела фиксированное значение: F(х1, х2, … , хn) = a.

Функции G(х1, х2, … , хn) называются ограничениями. Их можно задать как в виде равенств, так и неравенств. На регулируемые ячейки можно наложить дополнительные ограничения: неотрицательности и/или целочисленности, тогда искомое решение ищется в области положительных и/или целых чисел.

Под эту постановку попадает самый широкий круг задач оптимизации, в том числе решение различных уравнений и систем уравнений, задачи линейного и нелинейного программирования. Такие задачи обычно проще сформулировать, чем решать. И тогда для решения конкретной оптимизационной задачи требуется специально для нее сконструированный метод. Решатель имеет в своем арсенале мощные средства решения подобных задач: метод обобщенного градиента, симплекс-метод, метод ветвей и границ.

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

Размещено на http://www.allbest.ru/

Рис. 11. Окно диалога Поиск решения

После открытия диалога Поиск решения (11) необходимо выполнить следующие действия:

1. в поле Установить целевую ячейку ввести адрес ячейки, содержащей формулу для вычисления значений оптимизируемой функции, в нашем примере целевая ячейка - это С4, а формула в ней имеет вид: = C3^2 - 5*C3 + 6;

2. для максимизации значения целевой ячейки, установить переключатель максимальному значению в положение ?, для минимизации используется переключатель минимальному значению, в нашем случае устанавливаем переключатель в положение значению и вводим значение 0;

3. в поле Изменяя ячейки ввести адреса изменяемых ячеек, т.е. аргументов целевой функции (С3), разделяя их знаком «;» (или щелкая мышью при нажатой клавише Сtrl на соответствующих ячейках), для автоматического поиска всех влияющих на решение ячеек используется кнопка Предположить;

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

5. для запуска процесса поиска решения нажать кнопку Выполнить.

Для сохранения полученного решения необходимо использовать переключатель Сохранить найденное решение в открывшемся окне диалога Результаты поиска решения. После чего рабочий лист примет вид, представленный на рис.12. Полученное решение зависит от выбора начального приближения, которое задается в ячейке С4 (аргумент функции). Если в качестве начального приближения в ячейку С4 ввести значение, равное 1,0, то с помощью Поиска решения найдем второй корень, равный 2,0.

Размещено на http://www.allbest.ru/

Рис. 12. Результаты поиска

Опции, управляющие работой Поиска решения, задаваемые в окне Параметры (окно появляется, если нажать на кнопку Параметры окна Поиск решения), следующие (рис. 13):



Размещено на http://www.allbest.ru/

Рис. 13. Настройка параметров Решателя

· Максимальное время - ограничивает время, отведенное на процесс поиска решения (по умолчанию задано 100 секунд, что достаточно для задач, имеющих около 10 ограничений, если задача большой размерности, то время необходимо увеличить).

· Предельное число итераций - еще один способ ограничения времени поиска путем задания максимального числа итераций. По умолчанию задано 100, и, чаще всего, если решение не получено за 100 итераций, то при увеличении их количества (в поле можно ввести время, не превышающее 32767 секунд) вероятность получить результат мала. Лучше попытаться изменить начальное приближение и запустить процесс поиска заново.

· Относительная погрешность - задает точность, с которой определяется соответствие ячейки целевому значению или приближение к указанным ограничениям (десятичная дробь от 0 до 1).

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

· Сходимость - когда относительное изменение значения в целевой ячейке за последние пять итераций становится меньше числа (дробь из интервала от 0 до 1), указанного в данном параметре, поиск прекращается.

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

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

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

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

· Оценки - эта группа служит для указания метода экстраполяции - линейная или квадратичная, - используемого для получения исходных оценок значений переменных в каждом одномерном поиске. Линейная служит для использования линейной экстраполяции вдоль касательного вектора. Квадратичная служит для использования квадратичной экстраполяции, которая дает лучшие результаты при решении нелинейных задач.

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

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

Сохранить модель поиска решения можно следующими способами

1. при сохранении книги Excel после поиска решения все значения, введенные в окнах диалога Поиск решения, сохраняются вместе с данными рабочего листа. С каждым рабочим листом в рабочей книге можно сохранить один набор значений параметров Поиска решения;

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

3. еще один способ сохранения параметров поиска - сохранение их в виде именованных сценариев. Для этого необходимо нажать на кнопку Сохранить сценарий диалогового окна Результаты поиска решений.

Кроме вставки оптимальных значений в изменяемые ячейки Поиск решения позволяет представлять результаты в виде трех отчетов: Результаты, Устойчивость и Пределы. Для генерации одного или нескольких отчетов необходимо выделить их названия в окне диалога Результаты поиска решения. Рассмотрим более подробно каждый из них.

Рис. 14. Отчет по устойчивости

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

Рис. 15. Отчет по результатам

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

Решение нелинейного уравнения методом «Подбор параметра»

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

Подбор параметра - это итерационный процесс. Для ячейки, которая задана в поле Изменяя значение ячейки, Excel поочередно проверяет различные значения до тех пор, пока не отыщет наиболее подходящее. По умолчанию программа производит от 100 итераций или продолжает вычислять значения до достижения относительной погрешности 0,001.

Дано: 2x+x5-1=0, x(0;1]

Найти корни уравнения с помощью метода Excel «Подбор параметра».

1) Построим график данной функции

2) В ячейку G2 введем формулу для вычисления значения функции. В качестве аргумента используем ссылку на ячейку F5, т.е. = 2*F5+F5^5-1

3) В окне диалога Подбор параметра в поле Установить в ячейке ввел ссылку на ячейку с формулой (ячейку L6), в поле Значение - ожидаемый результат (в данном примере ожидаемый результат = 0), в поле Изменяя значение ячейки - ссылку на ячейку $L$7, в которой будет храниться значение подбираемого параметра.

4) После нажатия на кнопку Ok, получаю результат подбора параметра.

Ответ:

Вывод

Я ознакомился со способами решения нелинейных уравнений методом Excel и «Подбор параметра». Это очень удобный метод нахождения неизвестных параметров.

Список использованных источников информации

1. Додж М. Эффективная работа в Microsoft Excel 2000/ М. Додж, К.Стинсон.- СПб.: Питер, 2002. - 1056с.

2. Биллиг В.А. VBA и Office 97. Офисное программирование/ В.А.Биллиг, М.И. Дехтярь.- М.: Издательский отдел «Русская Редакция» ТОО «Channel Trading Ltd.», 1998.- 720 с.

3. Курицкий Б. Поиск оптимальных решений средствами Excel 7.0.- СПб.: BHV - Санкт-Петербург, 1997.- 384с.

4. Гетц К. Программирование в Microsoft Office. Полное руководство по VBA/ К. Гетц, М. Джилберт.- Киев: Издательская группа BHV, 1999.-768 с.

5. Гарнаев А.Ю. Excel, VBA, Internet в экономике и финансах.- СПб.: BHV - Санкт-Петербург, 2001.- 816 с.

6. Петруцос Э. Visual Basic 6 и VBA для профессионалов/ Э. Петруцос, К. Хау.- СПб.: Питер, 2000.- 432 с.

7. Гусева О.Л. Практикум по Excel/ О.Л. Гусева, Н.Н. Миронова.- М.: Финансы и статистика, 1997.- 160 с.

8. Турчак Л.И. Основы численных методов: Учеб. пособие.- М.: Наука; Гл. ред. физ.-мат. лит., 1987.- 320 с.

9. Мудров А.Е. Численные методы для ПЭВМ на языках Бейсик, Фортран и Паскаль.- Томск: МП «Раско», 1991.- 272 с.

10. http://www.exponenta.ru

11. http://necomod.narod.ru

Размещено на Allbest.ru