Минимизация функции методом ломаных

@Rainbow In Dark · Регистрация: 13.04.2014

Студворк — интернет-сервис помощи студентам

Доброго времени суток!
Помогите, пожалуйста, с программой.
Мне нужно реализовать метод ломаных.
Я написала программу, но она вылетает. Не знаю почему. Может я не правильно в until условие прописала? Хотя что я бы там не прописывала всегда не правильно работает.
Может я в методе не правильно разобралась. Или руки слишком кудрявые.

Вот код:

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
program Project2;
 
{$APPTYPE CONSOLE}
 
uses
  SysUtils,
  Math;
 
function f(xx: real):real;
begin
  f:=xx*xx+3*xx*(ln(xx)-1);
end;
 
var
a, b,  e,n : real;
i, j, k, t, l:integer;
x:  array [0..999999] of real;
p:  array [0..999999] of real;
 
begin
a:=0.5;
b:=1;
l:=2;
i:=0;
j:=0;
t:=0;
k:=0;
 
x[k]:=(-f(a)+l*a+l*b+f(b))/2/l;
p[k]:=(f(a)+f(b)+l*(a-b))/2;
 
repeat
E:=abs(f(x[k])-p[k]);
n:=f(x[k]-p[k])/2/l;
i:=1+i;
x[i]:=x[ k]-n;
i:=i+1;
x[i]:=x[ k]+n;
j:=1+j;
p[j]:=(f(x[ k ])+p[ k ])/2;
k:=k+1;
 
until abs(p[k]-p[k-1])<0.0001;
for i:=0 to j do
writeln(p[i]);
readln;
readln;
readln;
end.

@ZX Spectrum-128 · 02.11.2014, 21:37

Pascal
1
E:=abs(f(x[k])-p[k]);

Уже на втором шаге разность x[k])-p[k] отрицательна. Что не нравится функции ln.
Надо вкуривать алгоритм.
Или поискать решение на форуме.

Добавлено через 1 час 4 минуты
То бишь не разность, а значение (f(x[k])

@Rainbow In Dark · 02.11.2014, 21:42 **[ТС]**

Спасибо, только тогда не знаю как с этим бороться, ведь в итерационном процессе получилось отрицательное значение. Что делать? Брать модуль от х?

@ZX Spectrum-128 · 02.11.2014, 21:47

Rainbow In Dark, если вам сдать и забыть, то поищите готовое решение на форуме. А если вы хотите разобраться, то нужно думать алгоритм.

@Rainbow In Dark · 02.11.2014, 22:01 **[ТС]**

Мне проще будет самой написать, чем разобраться в чужом коде. Попробую что-то придумать! =)
Спасибо за помощь!

Новые блоги и статьи Все статьи Все блоги /
Использование SDL3-callbacks вместо функции main() на Android, Desktop и WebAssembly 8Observer8 24.01.2026 Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а привычная функция main(). . .	моя боль iceja 24.01.2026 Выложила интерполяцию кубическими сплайнами www. iceja. net REST сервисы временно не работают, только через Web. Написала за 56 рабочих часов этот сайт с нуля. При помощи perplexity. ai PRO , при. . .	Модель сукцессии микоризы anaschu 24.01.2026 Решили писать научную статью с неким РОманом	http://iceja.net/ математические сервисы iceja 20.01.2026 Обновила свой сайт http:/ / iceja. net/ , приделала Fast Fourier Transform экстраполяцию сигналов. Однако предсказывает далеко не каждый сигнал (см ограничения http:/ / iceja. net/ fourier/ docs ). Также. . .
http://iceja.net/ сервер решения полиномов iceja 18.01.2026 Выкатила http:/ / iceja. net/ сервер решения полиномов (находит действительные корни полиномов методом Штурма). На сайте документация по API, но скажу прямо VPS слабенький и 200 000 полиномов. . .	Расчёт переходных процессов в цепи постоянного тока igorrr37 16.01.2026 / * Дана цепь(не выше 3-го порядка) постоянного тока с элементами R, L, C, k(ключ), U, E, J. Программа находит переходные токи и напряжения на элементах схемы классическим методом(1 и 2 з-ны. . .	Восстановить юзерскрипты Greasemonkey из бэкапа браузера damix 15.01.2026 Если восстановить из бэкапа профиль Firefox после переустановки винды, то список юзерскриптов в Greasemonkey будет пустым. Но восстановить их можно так. Для этого понадобится консольная утилита. . .	Сукцессия микоризы: основная теория в виде двух уравнений. anaschu 11.01.2026 https:/ / rutube. ru/ video/ 7a537f578d808e67a3c6fd818a44a5c4/

@Rainbow In Dark 16 / 0 / 0 Регистрация: 13.04.2014 Сообщений: 45
	02.11.2014, 21:42 [ТС]
	Спасибо, только тогда не знаю как с этим бороться, ведь в итерационном процессе получилось отрицательное значение. Что делать? Брать модуль от х? 0

@ZX Spectrum-128 6812 / 4568 / 4820 Регистрация: 05.06.2014 Сообщений: 22,434
	02.11.2014, 21:47
	Rainbow In Dark, если вам сдать и забыть, то поищите готовое решение на форуме. А если вы хотите разобраться, то нужно думать алгоритм. 0

@Rainbow In Dark 16 / 0 / 0 Регистрация: 13.04.2014 Сообщений: 45
	02.11.2014, 22:01 [ТС]
	Мне проще будет самой написать, чем разобраться в чужом коде. Попробую что-то придумать! =) Спасибо за помощь! 0