VPF::[Алгоритм] Метод Зейделя и метод итерации

prikolist

Дата 17.6.2009, 22:46 (ссылка)

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

Загрузка ...

Шустрый

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

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

Здравсвуйте программисты! Подскажите пожалуйста,как решить метод Зейделя и методом Итерации, вот эти уравнения. Вот код, как эти уравнения правильно использовать в данном коде.

Вот метод итераций, как его подстроить, чтоб решило моё.

Метод итерации

Код


#include "stdafx.h"
#define eps 0.0001 //требуемая точность вычислений
void iteration(double a[N][N],double b[N],double x[N])
{
 int i,j;
 double norma; //чебышевская норма вектора
 double xn[N]={0};//вектор для текущей итерации, начальное значение
       //должно быть равно начальному приближению
 


 do{
  norma=0.0;
  for(i=0;i < N;i++)
  {
   xn[i]=-b[i];

   for(j=0;j < N;j++)
   {
    if(i!=j)
     xn[i]+=a[i][j]*x[j];
   }

   xn[i]/=-a[i][i];
  }

  for(i=0;i < N;i++)
  {
   if(fabs(x[i]-xn[i]) > norma)
    norma=fabs(x[i]-xn[i]); //Вычисление нормы вектора
   x[i]=xn[i];
  }
  
 }
 while(norma > eps); //проверка на необходимую точность вычислений


 return;
}

Метод Зейделя

Код


 for(i=0;i < N;i++)
  {
   x[i]=-b[i];

   for(j=0;j < N;j++)
   {
    if(i!=j)
     x[i]+=a[i][j]*x[j];
   }

   x[i]/=-a[i][i];
  }

  for(i=0;i < N;i++)
  {
   if(fabs(x[i]-xn[i]) > norma)
    norma=fabs(x[i]-xn[i]);
   xn[i]=x[i];
  }

[COLOR="Blue"]Подскажите пожалуйста, завтра или в пятницу хочу сдать этот код, очень нужно. Как подстроить его под мои уравнения?[/COLOR]

Вот фото уравнений

Присоединённый файл ( Кол-во скачиваний: 12 )

111.jpg 11,27 Kb

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