Форум по Delphi программированию

Delphi Sources



Вернуться   Форум по Delphi программированию > Все о Delphi > Синтаксис
Ник
Пароль
Регистрация <<         Правила форума         >> FAQ Пользователи Календарь Поиск Сообщения за сегодня Все разделы прочитаны

Ответ
 
Опции темы Поиск в этой теме Опции просмотра
  #1  
Старый 11.12.2009, 16:49
vbs vbs вне форума
Прохожий
 
Регистрация: 11.12.2009
Сообщения: 1
Репутация: 10
По умолчанию Ошибки вычитания

Здравствуйте, уважаемые форумчане!
В общем, название темы само за себя говорит, а теперь более подробно.
Есть программа, решающая задачу Коши методом Эйлера и методом двойного пересчета. Так как метод двойного пересчета, по-сути, усовершенствованный метод Эйлера, то они для удобства объединены в одну процедуру, которая за каждый проход цикла считает и выводит в таблицу: Х, Y, Y2, Е. Е - это погрешность, т.е. разность между Y2 (значение y методом двойного пересчета) и Y (значение у методом Эйлера). Так вот, проблема в том и заключается, что по какой-то неведомой мне причине вычитаение вычисляется неправильно. Надеюсь на то, что вы мне поможете разобраться с этим.

Итак, код классов:
Код:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
type
  ADouble=array[1..2] of double;
 
  //Класс - родитель, для вывода результатов и хранения заданной функции
  EchoMethods=class(TObject)
    //Функция
    Function MyFunction(x,dy:double):ADouble;
    //Процедура вывода в таблицу
    Procedure EchoTab(x,y,y2:double);
  end;
 
  //Класс - наследник, для хранения методов счета
  DifferentialMethod2=class(EchoMethods)
    //Метод Эйлера/Двойного пересчета
    Procedure Euler(x0,x1,y0,dy0,h:double);
  end;
 
  ...

Коды Метода Эйлера, заданной функции и вывода:
Код:
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
//Заданная функция
Function EchoMethods.MyFunction(x,dy:double):ADouble;
begin
  MyFunction[1]:=dy;
  MyFunction[2]:=-0.01*exp(0.4*x);
end;
 
//Метод Эйлера
Procedure DifferentialMethod2.Euler(x0,x1,y0,dy0,h:double);
var
  F:ADouble;
  x,x2:double;
  y,y2:double;
  dy,dy2:double;
 
begin
  x:=x0;
  x2:=x0;
  y:=y0;
  y2:=y0;
  dy:=dy0;
  dy2:=dy0;
  Form1.Series1.AddXY(x0,y0);
  Form1.Series2.AddXY(x0,y0);
  repeat
    EchoTab(x,y,y2);
    //Стандартный метод Эйлера
    F:=MyFunction(x,dy);
    x:=x+h;
    y:=y+F[1]*h;
    dy:=F[1]+F[2]*h;
    Form1.Series1.AddXY(x, y);
    //Метод двойного пересчета
    repeat
      F:=MyFunction(x2,dy2);
      x2:=x2+h/2;
      y2:=y2+F[1]*(h/2);
      dy2:=dy2+F[2]*(h/2);
      if (Round(x2*100) mod 2 = 0) then Form1.Series2.AddXY(x2,y2);
    until (Round(x2*100) mod 2 = 0);
  until (x>=x1);
end;
 
//Процедура вывода в таблицу
Procedure EchoMethods.EchoTab(x,y,y2:double);
var
  i:integer;
 
begin
  i:=Round(x*10);
  Form1.StringGrid1.Cells[0,i+1]:=FloatToStr(x);
  Form1.StringGrid1.Cells[1,i+1]:=FloatToStr(y);
  Form1.StringGrid1.Cells[2,i+1]:=FloatToStr(y2);
  <font color="Red">Form1.StringGrid1.Cells[3,i+1]:=FloatToStr(y2-y);</font>
end;

Собственно, в выделенной строке вся засада и кроется. Вычитание производится совершенно неверно. Хотя, конечно, может и не в выделенной...
В общем, такие пироги. Хотелось бы с вашей помощью понять, в чем причина неправильного счета.
Ответить с цитированием
Ответ


Delphi Sources

Опции темы Поиск в этой теме
Поиск в этой теме:

Расширенный поиск
Опции просмотра

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

BB-коды Вкл.
Смайлы Вкл.
[IMG] код Вкл.
HTML код Выкл.
Быстрый переход


Часовой пояс GMT +3, время: 04:59.


 

Сайт

Форум

FAQ

Соглашения

Прочее

 

Copyright © Форум "Delphi Sources" by BrokenByte Software, 2004-2025