Форум программистов, компьютерный форум, киберфорум
Наши страницы
Delphi для начинающих
Войти
Регистрация
Восстановить пароль
 
Рейтинг 4.89/27: Рейтинг темы: голосов - 27, средняя оценка - 4.89
Aleksus
0 / 0 / 0
Регистрация: 06.04.2010
Сообщений: 3
1

Метод квадратичной интерполяции-экстраполяции

06.04.2010, 21:01. Просмотров 4824. Ответов 3
Метки нет (Все метки)

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

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
program kvadr;
 
{$APPTYPE CONSOLE}
 
uses
  SysUtils;
 
var
   xm, a, b, e, h, x:real;
 
function fun(x: real): real;
{ âû÷èñëåГ*ГЁГҐ Г§Г*Г*Г·ГҐГ*ГЁГї ГґГіГ*êöèè Гў òî÷êå}
begin
   fun:=0.1*x*x*x-2*x*x+10*sqrt(x)
  end;
 
function mkv(h, x1, e: real; var xm: real): real;
var f0, f1, f2, x0, x2, c, b: real;
Begin
     xm:=x1;
     repeat
      x1:=xm;  x0:=x1-h; x2:=x1+h;
      f0:=fun(x0); f1:=fun(x1); f2:=fun(x2); c:=1/2/h/h*(f0-2*f1+f2);
      b:=((-1)*f0*(2*x1+h)+4*f1*x1-f2*(2*x1-h))/2/h/h;
      xm:=(-1)*b/2/c
     until ((xm-x1)<e);
     mkv:=fun(xm);
end;
begin
  write ('ГЌГ*Г·Г*ëüГ*ûé ГёГ*ГЈ ïîèñêГ* : '); read (h);
  write ('Òî÷Г*îñòü             : '); read (e);
  write ('ГЌГ*Г·Г*ëüГ*îå Г§Г*Г*Г·ГҐГ*ГЁГҐ X : '); read (x);
 
  h:=mkv(h,x,e,xm);
  write ('F(x)= ', h:12:5, ' x=',xm:12:5);
  readln;
end.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
06.04.2010, 21:01
Ответы с готовыми решениями:

Метод Кубической интерполяции зацикливается
unit Unit1; interface uses Windows, Messages, SysUtils, Variants,...

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

График квадратичной функции
unit Unit1; interface uses Windows, Messages, SysUtils, Variants,...

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

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

3
Inadequate
Retired
7710 / 2543 / 670
Регистрация: 17.10.2009
Сообщений: 5,100
07.04.2010, 03:25 2
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
program kvadr;
 
{$APPTYPE CONSOLE}
 
uses
  SysUtils;
 
var
   xm, a, b, e, h, x:real;
 
function fun(x: real): real;
{ вычисление значения функции в точке}
begin
   fun:=0.1*x*x*x-2*x*x+10*sqrt(x)
  end;
 
function mkv(h, x1, e: real; var xm: real): real;
var f0, f1, f2, x0, x2, c, b: real;
Begin
     xm:=x1;
     repeat
      x1:=xm;  x0:=x1-h; x2:=x1+h;
      f0:=fun(x0); f1:=fun(x1); f2:=fun(x2); c:=1/2/h/h*(f0-2*f1+f2);
      b:=((-1)*f0*(2*x1+h)+4*f1*x1-f2*(2*x1-h))/2/h/h;
      xm:=(-1)*b/2/c
     until ((xm-x1)<e);
     mkv:=fun(xm);
end;
begin
  write ('Начальный шаг поиска : '); ReadLn (h);
  write ('Точность             : '); ReadLn (e);
  write ('Начальное значение X : '); ReadLn (x);
 
  h:=mkv(h,x,e,xm);
  WriteLn ('F(x)= ', h:12:5, ' x=',xm:12:5);
  ReadLn;
end.
Не проверял результаты на правильность, но она их выводит
2
Aleksus
0 / 0 / 0
Регистрация: 06.04.2010
Сообщений: 3
11.04.2010, 20:56  [ТС] 3
Спасибо! Все понятно
0
!Master P!
0 / 0 / 0
Регистрация: 07.05.2011
Сообщений: 7
15.06.2011, 22:02 4
скиньте плз ссылку или название книги в которой этот метод описан очень прошу
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
15.06.2011, 22:02

Поиск минимума у функций методом квадратичной параболы
Реализовать поиск минимума у функций y=sin(x), y=cos(x), y=a*(x^2)+b*x+c...

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

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


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

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

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