РефератыИнформатика, программированиеМеМетод наискорейшего спуска

Метод наискорейшего спуска

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ


КУРСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
Кафедра ПО ВТ

Курсовая работа


по дисциплине «Методы и алгоритмы обработки изображений»


тема: Метод наискорейшего спуска


2004


1. Техническое задание




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


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


Общую постановку задачи оптимизации можно описать следующим образом:


Задачей многомерной оптимизации является минимизация функции U=f(x1,x2,..xm) от m переменных x1,x2,…xm. Если нет ограничений на параметры x1,x2,…xm, то говорят о глобальной минимизации, если ограничения есть, то говорят об условной минимизации. Вектор g(x) называется градиентом функции f(x) и обозначается g(x)=gradf(x)=(df/dx1)*e1+…+(df/dxm)*em.


Вектор g(x) указывает направление наискорейшего возрастания функции, а вектор –g(x) называется антиградиентом и указывает направление наискорейшего убывания функции.


Метод наискорейшего спуска.


Из начальной точки x0 перпендикулярно линии уровня то есть в направлении градиента (антиградиента), двигаются до тех пор, пока функция возрастает (убывает), то есть решают одномерную задачу минимизации для функции F(a)=f(x0-a*gradf(x0)), где a выступает в качестве параметра.


В результате находится значение x1=x0-a* gradf(x0) соответствующее максимальному (минимальному) значению функции на выбранной прямой. Затем вычислительный процесс повторяется для точки x1 и так далее. Критерием окончания является условие f’
(xm)=|df/dxi|<=ε, i=1,…,m.


Данное программное изделие (далее по тексту именуемое «программа») носит название «Метод наискорейшего спуска» – это одна из задач оптимизации.



Основания для разработки


Разработка ведется на основании задания по курсовой работе.



Назначение разработки


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



Функциональное и эксплуатационное назначение изделия



Перечень требований пользователя к программному изделию


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



Рассмотренные альтернативы


При постановке задачи на разработку данного программного изделия рассматривались следующие альтернативы:


1. Программная платформа.
Язык Delphi доступен в операционных системах Windows 95/98/Me, WindowsNT/2000/XP. Выбрана была система Windows 2000/NT в связи с высоким качеством программного кода, надёжностью и наличием расширенных возможностей разработки программного обеспечения под эту платформу.


2. Средства разработки.
Любой компилятор для системы Windows 2000 мог быть использован для написания данного программного продукта. Из наиболее распространенных, таких как MicrosoftVisual С++, WebGainVisualCafe, Borland C++, Borland C++ Builder, BorlandDelphi, PowerBuilder, BorlandJBuilder был выбран компилятор BorlandDelphi. Выбор обоснован широким распространением, удобством использования, высоким качеством генерируемого кода данной среды.



Требования к программе или программному изделию



Стандарты


Разработка программной документации и программного продукта должна производится согласно ГОСТ 19.701-90, ГОСТ 2.304-88.



Информационная совместимость


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



Программные ограничения, совместимость


Программное изделие должно быть написано на языке
Delphi 5.0, и работать под управлением операционных систем MicrosoftWindows 98, Windows 2000.



Требования к составу и параметрам технических средств (аппаратные ограничения)


Для работы с программным изделием должно требоваться:


· процессор, аналогичный IntelCeleron 900;


· оперативная память объемом 64 Мб;


· накопитель на жестком диске с объемом 1 Гб;


· видеоадаптер SVGA с цветным дисплеем;


· клавиатура и манипулятор типа “мышь”.



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


Входные данные должны вводится вручную пользователем с помощью клавиатуры.



Процессы обработки


Процессы обработки должны разделяться на следующие этапы:


· решение задачи минимизации методом наискорейшего спуска;



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


Выходными данными программы должно являться:


· вывод на дисплей, полученного решения задачи минимизации методом наискорейшего спуска;



Результирующие компоненты изделия


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


· исходный текст программы на языке
Borland Delphi 5.0;


· исполняемый файл;


· программная документация на изделие
.



Носители информации


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


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


· используемый в процессе работы программы объем оперативного запоминающего устройства 64 Мб,


· используемый в процессе работы программы объем накопителя на жестком диске
1 Гб.



Безопасность и секретность


Информация, содержащаяся в программном изделии, не является секретной, поэтому ограничение доступа к ней не предусматривается.



Требования к надежности



Требований к надежности предъявлено не было.


Рестарт


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



Удобство эксплуатации


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



Мобильность


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



Требования к программной документации


Программная документация должна соответствовать ГОСТ 19.701-90, ГОСТ 2.304-88 и включать следующие документы:


· техническое задание
;


· рабочий проект.


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



Технико-экономические показатели


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



Стадии и этапы разработки


Выполнение разработки должно включать две стадии:


· техническое задание
;


· рабочий проект.


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


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



Порядок контроля и приемки


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



2. Рабочий проект


Примером задачи оптимизации является минимизация функции U=f(x1,x2) от переменных x1,x2.


Метод наискорейшего спуска.


Из начальной точки x1 перпендикулярно линии уровня то есть в направлении градиента (антиградиента), двигаются до тех пор, пока функция возрастает (убывает), то есть решают одномерную задачу минимизации для функции F(a)=f(x0-a*gradf(x0)),


где a выступает в качестве параметра.


В результате находится значение x1=x0-a* gradf(x0) соответствующее максимальному (минимальному) значению функции на выбранной прямой. Затем вычислительный процесс повторяется для точки x1 и так далее. Критерием окончания является условие f’
(xm)=|df/dxi|<=ε, i=1,…,m.


Общие сведения


Данное программное изделие носит название «Метод наискорейшего спуска» и предназначено для решения одного из видов задач оптимизации.


Назначение разработки


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


Функциональные и эксплуатационные характеристики изделия


Функциональные возможности программного изделия


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


Рассмотренные альтернативы


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


Требования к программе или программному изделию


Стандарты


Разработка программной документации и программного продукта производилась согласно ГОСТ 19.701-90, ГОСТ 2.304-88.


Информационная совместимость


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


Программные ограничения, совместимость


Программное изделие написано на языке Delphi 5.0, и работает под управлением операционных систем MicrosoftWindows 95/98/Me, WindowsNT/2000/XP. Для использования исходных текстов, поставляемых с программным изделием, требуется система разработки программного обеспечения BorlandDelphi 5.0 или выше и минимальные аппаратные ограничения, накладываемые этой системой.


Требования к составу и параметрам технических средств (аппаратные ограничения)


Для работы с программным изделием требуется:


· процессор, аналогичный Intel 486, или выше;


· оперативная память объемом не менее 64 Мб;


· накопитель на жестком диске со свободным местом не менее 1 Гб;


· видеоадаптер SVGA с цветным дисплеем;


· клавиатура или манипулятор типа “мышь”.


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


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


Процессы обработки


Процессы обработки разделяются на 2 группы:


· решение задачи минимизации методом наискорейшего спуска;


· отображение полученного результата на экране;


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


Выходными данными программы являются:


· отображение на экране результата решения задачи минимизации методом наискорейшего спуска.


Результирующие компоненты изделия


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


a) исходный текст программы на языке Delphi 5.0, включающий файлы:


· файл реализации модуля Unit1.– текстовый файл с кодом алгоритма решения задачи динамического программирования,


· файл опций проекта Project1.dpr –файл c установками опций проекта,


· исполняемый файл Project1.exe – автономный исполняемый файл (законченное приложение Windows).


b) программная документация на изделие (в том числе в цифровом виде в файле Курсовик.doc).


Носители информации


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


Программное изделие использует следующие носители информации:


• используемый в процессе работы программы объем оперативного запоминающего устройства не превышает 64 Мб,


• используемый в процессе работы программы объем накопителя на жестком диске не превышает 1 Гб.


Безопасность и секретность


Информация, содержащаяся в программном изделии, не является секретной, поэтому ограничение доступа к ней не предусмотрено.


Требования к надежности


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


Рестарт


В случаях возникновения аппаратного или программного сбоя работа программы возобновляется повторным запуском.


Удобство эксплуатации


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


Мобильность


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



Состав проекта


Рабочий проект состоит
из исполняемого файла Project1.exe.



Описание модул
я


Модуль
Unit1


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


Литература


1. А.Я. Архангельский Программирование в Delphi 5.0. – М.: ЗАО “Издательство БИНОМ”, 2000 г.


2. Боресков А.В., Шикин Е.В., Шикина Г.Е. Компьютерная графика: первое знакомство; Под ред. Е.В. Шикина. - М.: Финансы и статистика, 1996. - 176 с.: ил. - (Диалог с компьютером).


3. Краснов М.В. OpenGL графика в проектах Delphi, БХВ-Санкт-Петербург, 2000. – 352 с.


Приложение А




Листинг программы




Файл Unit1.pas


unit Unit1;


interface


uses


Windows, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs,


StdCtrls, ExtCtrls,OpenGL, Grids, ComCtrls;


type


TForm1 = class(TForm)


Button1: TButton;


Label1: TLabel;


Panel1: TPanel;


Label2: TLabel;


StringGrid1: TStringGrid;


RadioButton1: TRadioButton;


Button2: TButton;


Label3: TLabel;


Label4: TLabel;


Edit1: TEdit;


Edit2: TEdit;


Label5: TLabel;


Label6: TLabel;


Edit3: TEdit;


Label7: TLabel;


UpDown1: TUpDown;


Edit4: TEdit;


ComboBox1: TComboBox;


Edit5: TEdit;


Label8: TLabel;


Edit6: TEdit;


Label9: TLabel;


procedure Button1Click(Sender: TObject);


procedure FormCreate(Sender: TObject);


procedure FormDestroy(Sender: TObject);


procedure FormResize(Sender: TObject);


procedure FormPaint(Sender: TObject);


procedure Panel1MouseDown(Sender: TObject; Button: TMouseButton;


Shift: TShiftState; X, Y: Integer);


procedure Panel1MouseUp(Sender: TObject; Button: TMouseButton;


Shift: TShiftState; X, Y: Integer);


procedure Panel1MouseMove(Sender: TObject; Shift: TShiftState; X,


Y: Integer);


procedure Button2Click(Sender: TObject);


procedure Edit3Change(Sender: TObject);


procedure Edit1Change(Sender: TObject);


procedure Edit2Change(Sender: TObject);


procedure Edit4Change(Sender: TObject);


procedure Edit5Change(Sender: TObject);


procedure Edit6Change(Sender: TObject);


private


DC : HDC;


hrc: HGLRC;


{ Private declarations }


end;


type TE = array [0..2,0..2] of glfloat;


type TM = array [0..2] of glfloat;


var


Form1: TForm1;


x,y,z,m,a,b,c,d,masht : glfloat;


mx,my,mz: TE;


mm : TM;


xmaus,ymaus,rt,verx: integer;


implementation


{$R *.DFM}


{Формат пикселя}


procedure SetDCPixelFormat (hdc : HDC);


var


pfd : TPixelFormatDescriptor;


nPixelFormat : Integer;


begin


FillChar (pfd, SizeOf (pfd), 0);


pfd.dwFlags := PFD_DRAW_TO_WINDOW or PFD_SUPPORT_OPENGL or PFD_DOUBLEBUFFER;


nPixelFormat := ChoosePixelFormat (hdc, @pfd);


SetPixelFormat (hdc, nPixelFormat, @pfd);


end;


{====================================================


Создание формы}


procedure TForm1.FormCreate(S

ender: TObject);


begin


DC := GetDC (Form1.Panel1.Handle);


SetDCPixelFormat(DC);


hrc := wglCreateContext(DC);


wglMakeCurrent(DC, hrc);


glClearColor (0.5, 0.5, 0.75, 1.0); // цвет фона


glColor3f (1.0, 0.0, 0.5); // текущий цвет примитивов


glEnable(GL_DEPTH_TEST);


x:=0.0;


y:=0.0;


z:=0.0;


m:=1.0;


a:=0.0;


b:=0.0;


c:=0.0;


xmaus:= 0;


ymaus:=0 ;


masht:=0;


rt:= 0 ;


verx:= 0;


StringGrid1.Cells[0,0]:='1';


tringGrid1.Cells[0,1]:='0';


tringGrid1.Cells[0,2]:='0';


tringGrid1.Cells[0,3]:='0';


end;


procedure TForm1.Button1Click(Sender: TObject);


var


hx,hy,x0,y0,f0,gx,gy,e,f1: real;


i : integer;


begin


StringGrid1.ColCount:=1;


hx:=0.001;


hy:=0.002;


e:=0.002;


x0:=StrToInt(Edit1.Text);


y0:=StrToInt(Edit2.Text);


if(ComboBox1.ItemIndex=0) then begin


f0:=(exp(x0)+exp(y0))/(exp(x0*x0)+exp(y0*y0));


gx:=exp(x0)/(exp(x0*x0)+exp(y0*y0))-exp(x0*x0)*x0*(2*(exp(x0)+exp(y0))/((exp(x0*x0)+exp(y0*y0))*(exp(x0*x0)+exp(y0*y0))));


gy:=exp(y0)/(exp(y0*y0)+exp(x0*x0))-exp(y0*y0)*y0*(2*(exp(y0)+exp(x0))/((exp(y0*y0)+exp(x0*x0))*(exp(y0*y0)+exp(x0*x0))));


end;


if(ComboBox1.ItemIndex=1) then begin


f0:=sin(x0)*sin(x0)/(sin(y0)*sin(y0)+1);


gx:=2*sin(x0)*cos(x0)/(sin(y0)*sin(y0)+1);


gy:=-2*sin(x0)*sin(x0)*cos(y0)*sin(y0)/((sin(y0)*sin(y0)+1)*(sin(y0)*sin(y0)+1));


end;


if(ComboBox1.ItemIndex=2) then begin


f0:=sin(x0)*cos(y0*y0);


gx:= cos(x0)*cos(y0*y0);


gy:=-2*sin(x0)*sin(y0*y0)*y0;


end;


x0:=x0+hx*gx;


y0:=y0+hy*gy;


if(ComboBox1.ItemIndex=0) then


f1:=(exp(x0)+exp(y0))/(exp(x0*x0)+exp(y0*y0));


if(ComboBox1.ItemIndex=1) then


f1:=sin(x0)*sin(x0)/(sin(y0)*sin(y0)+1);


if(ComboBox1.ItemIndex=2) then


f1:=sin(x0)*cos(y0*y0);


StringGrid1.Cells[0,0]:='1';


StringGrid1.Cells[0,1]:=FloatToStr(x0);


StringGrid1.Cells[0,2]:=FloatToStr(y0);


StringGrid1.Cells[0,3]:=FloatToStr(f1);


if f1<f0 then begin


hx:=hx/2;


hy:=hy/2;


end;


i:=1;


while not((abs(gx)<e) and (abs(gy)<e))do


begin


f0:=f1;


if(ComboBox1.ItemIndex=0) then begin


gx:=exp(x0)/(exp(x0*x0)+exp(y0*y0))-exp(x0*x0)*x0*(2*(exp(x0)+exp(y0))/((exp(x0*x0)+exp(y0*y0))*(exp(x0*x0)+exp(y0*y0))));


gy:=exp(y0)/(exp(y0*y0)+exp(x0*x0))-exp(y0*y0)*y0*(2*(exp(y0)+exp(x0))/((exp(y0*y0)+exp(x0*x0))*(exp(y0*y0)+exp(x0*x0))));


end;


if(ComboBox1.ItemIndex=1) then begin


gx:=2*sin(x0)*cos(x0)/(sin(y0)*sin(y0)+1);


gy:=-2*sin(x0)*sin(x0)*cos(y0)*sin(y0)/((sin(y0)*sin(y0)+1)*(sin(y0)*sin(y0)+1));


end;


if(ComboBox1.ItemIndex=2) then begin


gx:= cos(x0)*cos(y0*y0);


gy:=-2*sin(x0)*sin(y0*y0)*y0;


end;


x0:=x0+hx*gx;


y0:=y0+hy*gy;


if(ComboBox1.ItemIndex=0) then


f1:=(exp(x0)+exp(y0))/(exp(x0*x0)+exp(y0*y0));


if(ComboBox1.ItemIndex=1) then


f1:=sin(x0)*sin(x0)/(sin(y0)*sin(y0)+1);


if(ComboBox1.ItemIndex=2) then


f1:=sin(x0)*cos(y0*y0);


StringGrid1.Cells[i,0]:=IntToStr(i+1);


StringGrid1.Cells[i,1]:=FloatToStr(x0);


StringGrid1.Cells[i,2]:=FloatToStr(y0);


StringGrid1.Cells[i,3]:=FloatToStr(f1);


i:=i+1;


StringGrid1.ColCount:=i;


if f1<f0 then begin


hx:=hx/2;


hy:=hy/2;


end;


end;


Label1.Caption:=FloatToStr(f1);


refresh;


end;


procedure TForm1.FormDestroy(Sender: TObject);


begin


{====================================================


Конец работы приложения}


wglMakeCurrent(0, 0);


wglDeleteContext(hrc);


ReleaseDC (Panel1.Handle, DC);


DeleteDC (DC);


end;


procedure TForm1.FormResize(Sender: TObject);


begin


glViewport(0, 0, Panel1.ClientWidth, Panel1.ClientHeight);


glLoadIdentity;


glFrustum (-1, 1, -1, 1, 3, 10); // задаем перспективу


// этот фрагмент нужен для придания трёхмерности


glTranslatef(0.0, 0.0, -5); // перенос объекта - ось Z


InvalidateRect(Panel1.Handle, nil, False);


end;


procedure TForm1.FormPaint(Sender: TObject);


var


i,j: real;


zz,xl,yl,zl,zlm,stg,h: glfloat;


fy,fx,clr,clb,obl : real;


it,wn : integer;


mm1,mm2,mm3,mm4 : array [0..5] of TM;


begin


glClear (GL_COLOR_BUFFER_BIT or GL_DEPTH_BUFFER_BIT); // очистка буфера цвета


xl:=0;


yl:=0;


zl:=0;


zlm:=StrToFloat(StringGrid1.cells[0,3]);


stg:=StringGrid1.ColCount;


clr:=0;


clb:=1.0;


h:=StrToFloat(Edit3.text);


wn:=StrToInt(Edit5.text);


obl:= StrToFloat(Edit6.Text);


//ox oy oz


glscalef(masht,masht,masht);


glLineWidth(3);


glcolor3f(1.0,1.0,1.0);


glBegin(GL_Lines);


glVertex3f(-1,-1,1);


glVertex3f(1,-1,1);


glVertex3f(-1,-1,1);


glVertex3f(-1,1,1);


glVertex3f(-1,-1,1);


glVertex3f(-1,-1,-1);


glend;


glLineWidth(1);


glcolor3f(1.0,1.0,0.0);


glBegin(GL_Lines);


//x


glVertex3f(1.1-0.03,-1+0.1,1);


glVertex3f(1.1+0.03,-1-0.1,1);


glVertex3f(1.1+0.03,-1+0.1,1);


glVertex3f(1.1-0.03,-1-0.1,1);


//z


glVertex3f(-1-0.03,1.1,1);


glVertex3f(-1+0.03,1.1,1);


glVertex3f(-1-0.03,1.1,1);


glVertex3f(-1+0.03,1.2,1);


glVertex3f(-1-0.03,1.2,1);


glVertex3f(-1+0.03,1.2,1);


//y


glVertex3f(-1,-1,-1-0.1);


glVertex3f(-1,-1+0.1,-1-0.1);


glVertex3f(-1,-1+0.1,-1-0.1);


glVertex3f(-1-0.03,-1+0.2,-1-0.1);


glVertex3f(-1,-1+0.1,-1-0.1);


glVertex3f(-1+0.03,-1+0.2,-1-0.1);


glend;


// if(ComboBox1.ItemIndex=1) then begin


/f1:=sin(x0)*sin(x0)/(sin(y0)*sin(y0)+1);


/f1:=sin(x0)/(x0+y0*y0+1);


i:=-wn;


j:=-wn;


while j<=wn do begin


i:=-wn;


yl:=j;


//glBegin(GL_points);


glBegin(GL_Lines);


glcolor3f(clr,0.0,clb);


while i<=wn do begin


xl:=i; //f(x) ==


if(ComboBox1.ItemIndex=0) then


zl:=(exp(xl)+exp(yl))/(exp(xl*xl)+exp(yl*yl));


if(ComboBox1.ItemIndex=1) then


zl:=sin(xl)*sin(xl)/(sin(yl)*sin(yl)+1);


if(ComboBox1.ItemIndex=2) then


zl:=sin(xl)*cos(yl*yl);


if zl>=zlm then begin


zlm:=zl;


clr:=clr+0.01;


clb:=clb-0.01;


glcolor3f(clr,0.0,clb);


end;


glVertex3f(xl*obl,zl*obl,yl*obl); //+++


i:=i+0.1;


xl:=i;


if(ComboBox1.ItemIndex=0) then


zl:=(exp(xl)+exp(yl))/(exp(xl*xl)+exp(yl*yl));


if(ComboBox1.ItemIndex=1) then


zl:=sin(xl)*sin(xl)/(sin(yl)*sin(yl)+1);


if(ComboBox1.ItemIndex=2) then


zl:=sin(xl)*cos(yl*yl);


glVertex3f(xl*obl,zl*obl,yl*obl); //++++


end;


j:=j+h;


end;


glend;


//----------------------


i:=-wn;


j:=-wn;


clr:=0;


clb:=1.0;


zlm:=StrToFloat(StringGrid1.cells[0,3]);


while j<=wn do begin


i:=-wn;


xl:=j;


glBegin(GL_Lines);


glcolor3f(clr,0.0,clb);


while i<=wn do begin


yl:=i; //f(x) ||


if(ComboBox1.ItemIndex=0) then


zl:=(exp(xl)+exp(yl))/(exp(xl*xl)+exp(yl*yl));


if(ComboBox1.ItemIndex=1) then


zl:=sin(xl)*sin(xl)/(sin(yl)*sin(yl)+1);


if(ComboBox1.ItemIndex=2) then


zl:=sin(xl)*cos(yl*yl);


if zl>=zlm then begin


zlm:=zl;


clr:=clr+0.01;


clb:=clb-0.01;


glcolor3f(clr,0.0,clb);


end;


glVertex3f(xl*obl,zl*obl,yl*obl); //+++


i:=i+0.1;


yl:=i;


if(ComboBox1.ItemIndex=0) then


zl:=(exp(xl)+exp(yl))/(exp(xl*xl)+exp(yl*yl));


if(ComboBox1.ItemIndex=1) then


zl:=sin(xl)*sin(xl)/(sin(yl)*sin(yl)+1);


if(ComboBox1.ItemIndex=2) then


zl:=sin(xl)*cos(yl*yl);


glVertex3f(xl*obl,zl*obl,yl*obl); //+++


end;


j:=j+h;


end;


glend;


//------


it:=0;


glLineWidth(2);


glBegin(GL_Line_strip);


glcolor3f(0.0,1.0,0.0);


while it<StringGrid1.ColCount do begin


xl:=StrToFloat(StringGrid1.cells[it,1]); //+


yl:=StrToFloat(StringGrid1.cells[it,2]); //+


zl:=StrToFloat(StringGrid1.cells[it,3]); //+


glVertex3f(xl*obl,zl*obl,yl*obl);


it:=it+1;


end;


glend;


glLineWidth(1);


//tmax


glPointSize(6);


glBegin(GL_points);


glcolor3f(1.0,0.0,0.0);


glVertex3f(xl*obl,zl*obl,yl*obl);


glcolor3f(0.0,0.0,1.0);


glVertex3f(StrToFloat(StringGrid1.cells[0,1])*obl,StrToFloat(StringGrid1.cells[0,3])*obl,StrToFloat(StringGrid1.cells[0,2])*obl);


glend;


glscalef(1.0/masht,1.0/masht,1.0/masht);


glRotatef(verx,1,0,0);


glRotatef(rt,0,1,0);


SwapBuffers(DC);


Label2.Caption:=FloatToStr(zlm);


end;


procedure TForm1.Panel1MouseDown(Sender: TObject; Button: TMouseButton;


Shift: TShiftState; X, Y: Integer);


begin


if(Button=mbleft) then begin


radiobutton1.Checked:=true;


end;


end;


procedure TForm1.Panel1MouseUp(Sender: TObject; Button: TMouseButton;


Shift: TShiftState; X, Y: Integer);


begin


radiobutton1.Checked:=false;


end;


procedure TForm1.Panel1MouseMove(Sender: TObject; Shift: TShiftState; X,Y: Integer);


begin


rt:=0;


verx:=0;


if radiobutton1.Checked=true then begin


rt:=x- xmaus;


verx:=y-ymaus;


refresh;


end;


xmaus:=x;


ymaus:=y;


end;


procedure TForm1.Button2Click(Sender: TObject);


var


hx,hy,h,x0,y0,f0,gx,gy,g1x,g1y,e,f1: real;


i : integer;


begin


StringGrid1.ColCount:=1;


h :=0.001;


hx:=0.001;


hy:=0.001;


e:=0.002;


x0:=StrToInt(Edit1.Text);


y0:=StrToInt(Edit2.Text);


if(ComboBox1.ItemIndex=0) then begin


f0:=(exp(x0)+exp(y0))/(exp(x0*x0)+exp(y0*y0));


gx:=exp(x0)/(exp(x0*x0)+exp(y0*y0))-exp(x0*x0)*x0*(2*(exp(x0)+exp(y0))/((exp(x0*x0)+exp(y0*y0))*(exp(x0*x0)+exp(y0*y0))));


gy:=exp(y0)/(exp(y0*y0)+exp(x0*x0))-exp(y0*y0)*y0*(2*(exp(y0)+exp(x0))/((exp(y0*y0)+exp(x0*x0))*(exp(y0*y0)+exp(x0*x0))));


end;


if(ComboBox1.ItemIndex=1) then begin


f0:=sin(x0)*sin(x0)/(sin(y0)*sin(y0)+1);


gx:=2*sin(x0)*cos(x0)/(sin(y0)*sin(y0)+1);


gy:=-2*sin(x0)*sin(x0)*cos(y0)*sin(y0)/((sin(y0)*sin(y0)+1)*(sin(y0)*sin(y0)+1));


end;


if(ComboBox1.ItemIndex=2) then begin


f0:=sin(x0)*cos(y0*y0);


gx:= cos(x0)*cos(y0*y0);


gy:=-2*sin(x0)*sin(y0*y0)*y0;


end;


g1x:=gx;


g1y:=gy;


x0:=x0+hx*gx;


y0:=y0+hy*gy;


if(ComboBox1.ItemIndex=0) then


f1:=(exp(x0)+exp(y0))/(exp(x0*x0)+exp(y0*y0));


if(ComboBox1.ItemIndex=1) then


f1:=sin(x0)*sin(x0)/(sin(y0)*sin(y0)+1);


if(ComboBox1.ItemIndex=2) then


f1:=sin(x0)*cos(y0*y0);


StringGrid1.Cells[0,0]:='1';


StringGrid1.Cells[0,1]:=FloatToStr(x0);


StringGrid1.Cells[0,2]:=FloatToStr(y0);


StringGrid1.Cells[0,3]:=FloatToStr(f1);


if f1<f0 then begin


hx:=h;


hy:=h;


end


else begin


hx:=hx*2;


hy:=hy*2;


end;


i:=1;


while not((abs(g1x)<e) and (abs(g1y)<e))do


begin


f0:=f1;


x0:=x0+hx*gx;


y0:=y0+hy*gy;


if(ComboBox1.ItemIndex=0) then begin


g1x:=exp(x0)/(exp(x0*x0)+exp(y0*y0))-exp(x0*x0)*x0*(2*(exp(x0)+exp(y0))/((exp(x0*x0)+exp(y0*y0))*(exp(x0*x0)+exp(y0*y0))));


g1y:=exp(y0)/(exp(y0*y0)+exp(x0*x0))-exp(y0*y0)*y0*(2*(exp(y0)+exp(x0))/((exp(y0*y0)+exp(x0*x0))*(exp(y0*y0)+exp(x0*x0))));


f1:=(exp(x0)+exp(y0))/(exp(x0*x0)+exp(y0*y0));


end;


if(ComboBox1.ItemIndex=1) then begin


g1x:=2*sin(x0)*cos(x0)/(sin(y0)*sin(y0)+1);


g1y:=-2*sin(x0)*sin(x0)*cos(y0)*sin(y0)/((sin(y0)*sin(y0)+1)*(sin(y0)*sin(y0)+1));


f1:=sin(x0)*sin(x0)/(sin(y0)*sin(y0)+1);


end;


if(ComboBox1.ItemIndex=2) then begin


f1:=sin(x0)*cos(y0*y0);


g1x:= cos(x0)*cos(y0*y0);


g1y:=-2*sin(x0)*sin(y0*y0)*y0;


end;


StringGrid1.Cells[i,0]:=IntToStr(i+1);


StringGrid1.Cells[i,1]:=FloatToStr(x0);


StringGrid1.Cells[i,2]:=FloatToStr(y0);


StringGrid1.Cells[i,3]:=FloatToStr(f1);


i:=i+1;


StringGrid1.ColCount:=i;


if f1<f0 then begin


if(ComboBox1.ItemIndex=0) then begin


gx:=exp(x0)/(exp(x0*x0)+exp(y0*y0))-exp(x0*x0)*x0*(2*(exp(x0)+exp(y0))/((exp(x0*x0)+exp(y0*y0))*(exp(x0*x0)+exp(y0*y0))));


gy:=exp(y0)/(exp(y0*y0)+exp(x0*x0))-exp(y0*y0)*y0*(2*(exp(y0)+exp(x0))/((exp(y0*y0)+exp(x0*x0))*(exp(y0*y0)+exp(x0*x0))));


end;


if(ComboBox1.ItemIndex=1) then begin


gx:=2*sin(x0)*cos(x0)/(sin(y0)*sin(y0)+1);


gy:=-2*sin(x0)*sin(x0)*cos(y0)*sin(y0)/((sin(y0)*sin(y0)+1)*(sin(y0)*sin(y0)+1));


end;


if(ComboBox1.ItemIndex=2) then begin


gx:= cos(x0)*cos(y0*y0);


gy:=-2*sin(x0)*sin(y0*y0)*y0;


end;


hx:=h;


hy:=h;


end


else begin


hx:=hx*2;


hy:=hy*2;


end;


end;


Label1.Caption:=FloatToStr(f1);


refresh;


end;


procedure TForm1.Edit3Change(Sender: TObject);


var


c : real;


begin


try


c:=StrToFloat(Edit3.text);


except


Edit3.text:='0,05';


end;


end;


procedure TForm1.Edit1Change(Sender: TObject);


var


c : integer;


begin


try


c:=StrToInt(Edit1.text);


except


Edit1.text:='-1';


end;


end;


procedure TForm1.Edit2Change(Sender: TObject);


var


c : integer;


begin


try


c:=StrToInt(Edit2.text);


except


Edit2.text:='-1';


end;


end;


procedure TForm1.Edit4Change(Sender: TObject);


begin


masht:=StrToFloat(Edit4.text)/10;


refresh;


end;


procedure TForm1.Edit5Change(Sender: TObject);


var


c : integer;


begin


try


c:=StrToInt(Edit5.text);


except


Edit5.text:='1';


end;


end;


procedure TForm1.Edit6Change(Sender: TObject);


var


c : real;


begin


try


c:=StrToFloat(Edit6.text);


except


Edit6.text:='1,0';


end;


end;


end.


Приложение Б


Тест 1:


Найдём максимум функции f(x,y)=(exp(x)+exp(y))/(exp(x^2)+exp(y^2))


Простой градиентный метод



программный изделие технический средство


Метод наискорейшего спуска



Тест 2:


Найдём максимум функции f(x,y)=sin(x0)*sin(x0)/(sin(y0)*sin(y0)+1)


Простой градиентный метод



Метод наискорейшего спуска



Тест 3:


Найдём максимум функции f(x,y)=sin(x0)*cos(y0*y0)


Простой градиентный метод



Метод наискорейшего спуска



При изменении начальных точек, может меняться и положение максимума



Сохранить в соц. сетях:
Обсуждение:
comments powered by Disqus

Название реферата: Метод наискорейшего спуска

Слов:2611
Символов:36707
Размер:71.69 Кб.