VPF::Метод наименьших квадратов - Форум программистов

Axenia

Дата 22.5.2013, 22:27 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 5
Регистрация: 18.5.2013

Репутация: нет
Всего: нет

Помогите, пожалуйста. Почему дают Error код Delphi для Метод наименьших квадратов

Код


Procedure RegresieLiniara(x, y: array of extended;var a1, b: extended);
var n, i: integer;
    Sx, Sxy, Sy, Sxx: extended;

begin
  n := length(x);
  Sx := 0; Sxy := 0; Sy := 0; Sxx := 0;

  for i:=0 to n-1 do
  begin
   Sx := Sx + x[i];
   Sxy := Sxy + (x[i] * y[i]);
  Sy := Sy + y[i];
  Sxx := Sxx + (x[i] * x[i]);
  end;

  a1:=((n * Sxy) - (Sx * Sy))/((n * Sxx) - (Sx * Sx));
  b :=((Sxx * Sy) - (Sx * Sxy))/((n * Sxx) - (Sx * Sx));

 end;

M
Poseidon

При вставке фрагмента кода пользуйтесь кнопкой "Код"

Это сообщение отредактировал(а) Poseidon - 23.5.2013, 10:39

northener

Дата 23.5.2013, 01:26 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Завсегдатай
Сообщений: 1361
Регистрация: 2.9.2010

Репутация: 12
Всего: 20

Цитата(Axenia @ 22.5.2013, 22:27

)

Почему дают Error код Delphi для Метод наименьших квадратов

А где вы взяли сей код? И правильно ли вы его скопировали?

--------------------

Но только лошади летают вдохновенно.
Иначе лошади разбились бы мгновенно!

Axenia

Дата 23.5.2013, 06:28 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 5
Регистрация: 18.5.2013

Репутация: нет
Всего: нет

из теория http://ru.scribd.com/doc/96403336/metode

northener

Дата 24.5.2013, 01:39 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Завсегдатай
Сообщений: 1361
Регистрация: 2.9.2010

Репутация: 12
Всего: 20

Цитата(Axenia @ 23.5.2013, 06:28

)

из теория http://ru.scribd.com/doc/96403336/metode

Лень было искать в "той теория" решение МНК для линейной функции.
Так что либо в "той теория" ошибки, либо вы неправильно что-то поняли.

Вот система двух уравнений для линейной функции (Даю уравнения с использованием ваших имен переменных):

Код


  Sxy = a1 * Sxx + b * Sx
  Sy = a1 * Sx + n * b

Надеюсь решать систему двух уравнений вы умеете?

Это сообщение отредактировал(а) northener - 24.5.2013, 02:53

--------------------

Но только лошади летают вдохновенно.
Иначе лошади разбились бы мгновенно!

david_yu

Дата 24.5.2013, 12:12 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 3
Регистрация: 13.10.2006

Репутация: нет
Всего: нет

Код


Procedure RegresieLiniara(x, y: array of extended;var a1, b: extended);
var 
  iLow, iHigh, i, nx, ny: integer;
  Sx, Sxy, Sy, Sxx: extended;
begin
  nx := Length(x);
  ny := Length(y);
  
  if nx <> ny then begin
    a1 := 0;
    b := 0;
    Exit;
  end;
  
  iLow := Low(x);
  iHigh := High(x);
  
  Sx := 0; 
  Sxy := 0; 
  Sy := 0; 
  Sxx := 0;
     
  for i:= iLow to iHigh do begin
    Sx := Sx + x[i];
    Sxy := Sxy + (x[i] * y[i]);
    Sy := Sy + y[i];
    Sxx := Sxx + (x[i] * x[i]);
  end;
  
  if (Sxx = 0) xor (Sx = 0) then begin
   a1 := ((n * Sxy) - (Sx * Sy))/((n * Sxx) - (Sx * Sx));
   b  := ((Sxx * Sy) - (Sx * Sxy))/((n * Sxx) - (Sx * Sx));
  end;
end;

Этот ответ добавлен с нового Винграда - http://vingrad.com

!
Poseidon

Используйте кнопку "Код"! На новом Винграде она тоже есть.

Это сообщение отредактировал(а) Poseidon - 24.5.2013, 13:50

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Delphi: Для новичков \| Следующая тема »