Форум программистов, компьютерный форум, киберфорум
Наши страницы
Delphi для начинающих
Войти
Регистрация
Восстановить пароль
 
SkN[i]lpl[e][r]
3 / 3 / 2
Регистрация: 12.05.2011
Сообщений: 248
#1

нужно немного подправить программу. метод наискорейшего спуска - Delphi

18.05.2013, 09:35. Просмотров 793. Ответов 2
Метки нет (Все метки)

Нужно немного подкорректировать программу метод наискорейшего спуска. Работает для нахождения минимума функции из 3-х переменных. Хотел переделать для 2-х переменных, но что то не получилось, программа зависает... может кто поможет переделать чтоб адекватно работала...
Вот пример функции: Z=9*x2+16*y2-90*x-128*y.
Ответ: zmin=481 в точке (5;4).

Листинг программы и архив с файлами программы ниже.
листинг:
Delphi
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
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
unit Unit1;
 
interface
 
uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, Grids;
 
type
  TForm1 = class(TForm)
    Edit1: TEdit;
    Label1: TLabel;
    Label2: TLabel;
    Edit3: TEdit;
    Label3: TLabel;
    Label4: TLabel;
    Edit4: TEdit;
    Edit5: TEdit;
    Label5: TLabel;
    Button1: TButton;
    StringGrid1: TStringGrid;
    StringGrid2: TStringGrid;
    Label6: TLabel;
   function f(a: array of double ): double;
    function sgn(a, b: double): integer;
    function pr(x: array of double): double;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;
 
var
  Form1: TForm1;
    k:integer;
  type
  float = double;
 
 
implementation
 
{$R *.dfm}
 
 
function TForm1.f(a: array of float): float;
begin
  inc(k);
  f := sqr(a[0]-1)+sqr(a[1]-3)+4*sqr(a[2]+5);
end;
 
function TForm1.sgn(a, b: float): integer;
begin
  if (a - b) >= 0 then sgn := 1
  else sgn := -1;
end;
 
 
const
  max_n = 75;
var
  x, y, g, d, l, ff: array[0 .. pred(max_n)] of float;
  n: integer;
 
 
function TForm1.pr(x: array of float): float;
var
  p: float;
  i: integer;
begin
  p := 0;
  g[0]:=2*(x[0]);
  g[1]:=2*(x[1]);
  g[2]:=8*(x[2]+5);
  for i := 0 to pred(n) do p := p + sqr(g[i]);
  pr := sqrt(p);
end;
 
procedure TForm1.Button1Click(Sender: TObject);
var
 
  h, eps, z, g0, g2, dn: float;
  vrem, v: float;
  i, j, m, s1, s2, s3: integer;
begin
k:=0;
n:=strtoint(Edit1.text);
for i := 0 to pred(n) do begin
    x[i]:=strtofloat(stringgrid1.Cells[i,0]);
    end;
 h:=strtofloat(Edit3.Text);
 eps:=strtofloat(edit4.Text);
 
 for i := 0 to pred(n) do begin
    y[i] := x[i];
  end;
  z := f(x);
  g0 := pr(x);
 
  while g0 > eps do begin
    for i := 0 to pred(n) do d[i] := -g[i] / g0;
    l[0] := 0;
    ff[0] := z;
    repeat
      m := 0; l[2] := h;
      for i := 0 to pred(n) do x[i] := y[i]+l[2]*d[i];
      z := f(x);
      ff[2] := z;
      g0 := pr(x);
      g2 := 0;
      for i := 0 to pred(n) do g2 := g2 + g[i]*d[i];
      if (ff[2]>=ff[0]) or (g2>=0) then m := 0
      else begin
        { Удвоить длину шага, чтобы накрыть min }
        h := h * 2;
        m := 1;
      end;
    until m = 0;
 
    l[1] := h/2;
    for i := 0 to pred(n) do x[i] := y[i]+l[1]*d[i];
    z := f(x);
    ff[1] := z;
 
    { Выполнить первую квадратичную интерполяцию }
    l[3]:=h*(ff[1]-0.75*ff[0]-0.25*ff[2]);
    l[3]:=h*(ff[1]-0.75*ff[0]-0.25*ff[2]);
    l[3]:= l[3] / (2*ff[1]-ff[0]-ff[2]);
 
    for i := 0 to pred(n) do x[i] := y[i]+l[3]*d[i];
    z := f(x);
    ff[3] := z;
    { Имеем 4 значения L и 4 значения функции, упорядочим их в порядке убывания}
    repeat
      m := 0;
      for i := 0 to pred(n) do
        for j := i + 1 to n do
          if ff[i] > ff[j] then begin
            vrem:=l[i]; l[i]:=l[j]; l[j]:=vrem;
            vrem:=ff[i]; ff[i]:=ff[j]; ff[j]:=vrem;
          end;
 
      { Закончить поиск в данном направлении если точность достигнута }
      if abs(l[0]-l[1])<eps*50 then m := 0
      else begin
        s1:=sgn(l[1],l[0]);
        s2:=sgn(l[2],l[0]);
        s3:=sgn(l[3],l[0]);
        if (s1 = s2) and (s1=-s3) then begin
          l[2]:=l[3]; ff[2]:=ff[3];
        end;
        dn:=(l[1]-l[2])*ff[0]+(l[2]-l[0])*ff[2]+(l[0]-l[1])*ff[1];
        v:=(ff[0]-ff[1])/(2*dn);
        v := v * (l[1]-l[2])*(l[2]-l[0]);
        l[3]:=(l[0]+l[1])/2+v;
        for i := 0 to pred(n) do x[i]:=y[i]+l[3]*d[i];
        z:=f(x);
        ff[3]:=z;
        m:=1;
      end
    until m = 0;
 
    for i := 0 to pred(n) do begin
      x[i]:=y[i]+l[0]*d[i];
      y[i]:=x[i];
 
    end;
    z:=f(x);
    g0:=pr(x);
 
    h := h / 2;
  end;
 
  for i := 0 to pred(n) do begin
  stringgrid2.Cells[i,0]:=floattostr(x[i]);
 
  end;
 edit5.Text:=floattostr(z);
 
 
end;
 
end.

http://www.cyberforum.ru/delphi-beginners/thread760822.html
0
Вложения
Тип файла: rar Naisk_spusk.rar (194.8 Кб, 28 просмотров)
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
18.05.2013, 09:35
Я подобрал для вас темы с готовыми решениями и ответами на вопрос нужно немного подправить программу. метод наискорейшего спуска (Delphi):

Из метода градиента сделать метод наискорейшего спуска
Всем добра) Сижу вот с рабочим кодом метода градиента, требуется переделать...

Зависает программа расчета методом наискорейшего спуска
делаю программу для расчета методом наискорейшего спуска функцию...

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

нужно немного дополнить рабочую программу
unit Unit1; interface uses Windows, Messages, SysUtils, Variants,...

Метод градиентного спуска
unit Unit1; interface uses Windows, Messages, SysUtils, Variants,...

2
SkN[i]lpl[e][r]
3 / 3 / 2
Регистрация: 12.05.2011
Сообщений: 248
24.05.2013, 13:42  [ТС] #2
Никто не шарит чтоли?)
0
mss
2631 / 2256 / 275
Регистрация: 24.12.2010
Сообщений: 13,725
24.05.2013, 13:45 #3
Отладчик "шарит", у него и спроси.
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
24.05.2013, 13:45
Привет! Вот еще темы с решениями:

Помогите сделать метод скорейшего градиентного спуска
Имеется код НАискорейшего градиентного спуска,но мне нужен СКорейший...

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

Файлы. Подправить программу
program katalog5; {$APPTYPE CONSOLE} uses Windows,system; const...

Немного изменить программу
как сделать что бы программа СЧИТАЛА мои баллы??


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
3
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Рейтинг@Mail.ru