Министерство  образования Ставропольского края Государственно-образовательного учреждения                             Среднего профессионального образования  ГРК «Интеграл» 
 
 

              КУРСОВОЙ ПРОЕКТ

по дисциплине: «Математические методы»

тема: «Классические способы определения экстремумов, функций нескольких переменных» 
 

     Разработал                                                                           Руководитель

                                                                                                Капаильская С.А.

Группа                                                                                     Дата___________                                                              Дата____________                                                             Оценка_________

Подпись_________                                                             Подпись________ 
 
 

                                               Георгиевск

                                                     2011г.

      Содержание  

Содержание 2

Введение 3

Классические методы поиска экстремума  функции одной    переменной 3

Определение глобального максимума или минимума функции     одной переменной 6

Выпуклые и вогнутые функции 6

Методы исключения интервалов 10

Правило исключения интервалов 11

Поиск экстремумов функции нескольких переменных 15

Заключение 18

Литература 19

 

      Введение

 

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

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

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

 

      

      1.Пояснительная записка

      Определение глобального максимума  или минимума функции  одной переменной

      Пусть требуется максимизировать f(x) при ограничениях a<=x<=b, где a и b – установленные границы измерений переменной x. (Очевидно в этом случае проверку наличия локального оптимума необходимо проводить не только в стационарных точках, но и в граничных точках интервала). Алгоритм следущий:

      Шаг 1: приравнять df/dx=0 и найти все стационарные точки.

      Шаг 2: выбрать все стационарные точки, которые расположены а интервале [a,b]. Обозначим эти точки через x₁,x_2,…,x_n. Проверку наличия локального оптимума следует проводить только на множестве указанных точек, дополненном точками a и b.

      Шаг 3:найти наибольшее значение f(x) из множества f(a),f(b),f(x₁),…,f(x_n). Это значение соответствует глобальному максимуму.

      Выпуклые и вогнутые функции

      Это важный класс унимодальных функций. Введем обозначение: x=(x₁,x₂,…,x_n)-n-мерный вектор.

      Определение: Функция n мерных f(x), определенная на выпуклом множестве D, называется выпуклой функцией тогда и только тогда, когда для любых двух точек x⁽¹⁾ и x⁽²⁾ принадлежащих D, и любого числа L (0<=L<=1) выполняется неравенство: f(Lx⁽¹⁾ +(1-L)x⁽²⁾)<=Lf(x⁽¹⁾)+(1-L)f(x⁽²⁾).

 

      Свойства  выпуклых функций:

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

      2.Выпуклая  функция лежит над своими касательными

      3.Тангенс  угла наклона касательной, или  первая производная f(x), возрастает или по крайней мере не убывает при увеличении x.

      4.Вторая  производная f(x) всегда не отрицательна на рассматриваемом интервале.

      5.Для  выпуклой функции локальный минимум  всегда является глобальным минимумом.

      Определение: Градиент функции f(x₁,x₂,…,x_n) определяется как вектор

      grad  f(x₁,…,x_n)=(df/dx₁,df/dx₂,…,df/dx_n)^T.

      Определение:  Матрица Гессе (гессиан) для функции  f(x₁,…,x_n) есть симметрическая матрица порядка n*n: H_f(x₁,…,x_n)=[d²f/dx_idx_j]=    H(f).

      Проверка  функции на выпуклость: Функция f(x₁,…,x_n) выпуклая, если ее матрица Гессе положительно определена или положительно полуопределена для всех значений x₁,x₂,…,x_n .

      Для функции одной переменной: функция  f(x) выпуклая, если ее вторая производная неотрицательна для всех значений x: d²f/dx²=>0, для всех x.

      Если  матрица Гессе H_f – положительно определенная матрица, то f называется строго выпуклой функцией и обладает единственной точкой минимума.

      Проверка  матриц на положительную определенность:  

      1) Все диагональные элементы должны  быть положительными.

      2) Все ведущие главные определители  должны быть положительными.

      Проверка  матриц на положительную полуопределенность:

      1) Все диагональные элементы неотрицательны.

      2) Все главные определители неотрицательны.

      Замечание: Чтобы установить, что данная матрица  является отрицательно определенной (полуотрицательно определенной), следует умножить ее на -1 и проверить полученную матрицу  на положительную определенность (полуположительную  определенность).

      Вогнутая  функция. Функция f(x₁,…,x_n) является вогнутой функцией на множестве D тогда и только тогда, когда –f(x) есть выпуклая функция на D.

      Проверка  функции на вогнутость. Функция f(x₁,…,x_n) вогнутая, если ее матрица Гессе отрицательно определена, или отрицательно полуопределена для всех значений x₁,…,x_n.

      Пример: Исследовать функцию на выпуклость.

      f(x₁,x₂,x₃)=3x₁² +2x₂² +x₃² –2x₁x₂ –2x₁x₃ +2x₂x₃ –6x₁ –4x₂ –2x₃

                                  6x₁ –2x₂ –2x₃ –6

      grad(x₁,x₂,x₃)= 4x₂ –2x₁ +2x₃ –4

                                2x₃ –2x₁ =2x₂ –2

      H_f (x₁,x2,x3)=

      Исследуем H_f.

1. H_f –симметрическая матрица.
2. Все диагональные элементы H_f положительны.
3. Ведущие главные определители Н равны:

      |6|>0   

      Следовательно, H_f – положительно определенная матрица. Отсюда следует, что f-выпуклая функция. Более того, f строго выпуклая функция и обладает единственной точкой минимума.

 

      Методы  исключения интервалов

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

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

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

      Для унимодальной функции f(x) сравнение значений f(x) в двух различных точках интервала поиска позволяет определить, в каком из заданных этими двумя точками подинтервалов точка оптимума отсутствует.

 

      Правило исключения интервалов

      Пусть функция f унимодальна на интервале a£x£b, а ее минимум достигается в точке x*.

      Рассмотрим  точки x₁ и x₂, расположенные в интервале таким образом, что a<x₁<x₂<b. Сравнивая значения функции в точках x₁ и x₂, можно сделать следующие выводы:

1. Если f(x₁)>f(x₂), то точка минимума f(x) не лежит в интервале (a,x₁), т.е. x*Î(x₁,b)

      2. Если f(x₁)<f(x₂), то точка минимума не лежит в интервале (x₂,b), т.е. x*Î(a,x₂)

       3. Если f(x₁)=f(x₂), то можно исключить оба крайних интервала (a,x₁) и (x₂,b), при этом x*Î(x₁,x₂).

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

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

      Достоинства этих методов:

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

      (при  этом не требуется, чтобы исследуемые  функции были дифференцируемы).

 

      Метод золотого сечения 

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