VPF::[Delphi] Метод Зейделя решения СЛАУ

StS

Дата 13.12.2004, 08:36 (ссылка)

(голосов: 0)

Загрузка ...

Unregistered

Помогите с методом Зейделя, нужна реализация на Delphi
Плиз!

Fedor

Дата 13.12.2004, 17:31 (ссылка)

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

Загрузка ...

Днепрянин

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

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

Ты бы хоть объяснил, что это такое...

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

Мы - Днепряне. Мы всех сильней.

ElectricalStorm

Дата 14.12.2004, 14:20 (ссылка)

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

Загрузка ...

Опытный

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

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

Вкючаем гугл и сразу попадаем на
http://www.kursovik.net/programming/180505.html

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

Нужно знать инструмент, которым пользуешься

podval

Дата 15.12.2004, 10:52 (ссылка)

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

Загрузка ...

Где я? Кто я?

Профиль
Группа: Экс. модератор
Сообщений: 3094
Регистрация: 25.3.2002
Где: СПб

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

Можешь перевести на любой язык.

Код


//////////////////////////////////////////////////////////////////////////////
//
//  Solving system of linear equations (Zaidel method)
//  (c) Johna Smith, 1996
//
//  Method description:
//    Given: system of linear equations
//       a11 x1 + a12 x2 + ... + a1n xn = b1
//       a21 x1 + a22 x2 + ... + a2n xn = b2
//       ...
//       an1 x1 + an2 x2 + ... + ann xn = bn
//
//    We use following iterative formula to solve this system:
//
//                     1    i-1                N
//     xi(j+1)=xi(j)- --- ( SUM aik xk(j+1) + SUM aik xk(j) -bi)
//                    aii   k=1               k=i
//
//     where xi(j+1) is j+1-th iteration, xi(j) is j-th iteration
//
//////////////////////////////////////////////////////////////////////////////

#include <stdio.h>
#include <math.h>

#define N               3     // size of matrix
#define N1              N+1

float matrix[N][N1]=
      {{10,1,1,12},
       {2,10,1,13},
       {2,2,10,14}
      };
float z[N]={0,0,0};  // first iteration

float epsilon=0.0001;  // required accuracy

void ShowMatrix(void)
{
  for (int i=0;i<N;i++)
  {
       for (int j=0;j<N;j++)
            printf("%+f*x%d",matrix[i][j],i+1);
       printf("=%f\n",matrix[i][N]);
  }
}

void main(void)
{
  // Variables declaration
  float x[N];
  short int i,j;
  int iterations=0,finish=0;

  // Printing given matrix
  ShowMatrix();

  while (!finish)
  {
     finish=1;
     for (i=0;i<N;i++)
    {
        x[i]=-matrix[i][N];
        for (j=0;j<N;j++) x[i]+=matrix[i][j]*z[j];
       // don't stop iterations until required accuracy is reached
       if (fabs(x[i]/matrix[i][i])>=epsilon) finish=0;
       x[i]=z[i]-x[i]/matrix[i][i];
      // next iteration
      z[i]=x[i];
    }
    iterations++;
  }

  // Printing solution
  printf("\nSolution:\n");
  for (i=0;i<N;i++)
      printf("x%d=%f\n",i+1,x[i]);
  printf("%d iterations were made",iterations);
}

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