Метод итераций и метод Зейделя

[Sazary]

омг...
Внимательно смотрим на первую строку:

int a1,b1,x1;

теперь на вторую

int a1[N][N],b1[N],x1[N];

Что вы изменили? По сути - ничего. У вас по-прежнему по два раза объявлены a,b и x.

Удалите вообще первую строку!

Я просто не понимаю, как это можно не видеть.

[prikolist]

Хорошо,спасибо, смотрите, написав вот так:

int main()
{

double a1[N][N],b1[N],x[N];

int i;

a1[0][0] = -1.01;    a1[0][1] = 0.02;  a1[0][2] = -0.62;   a1[0][3] = 0.08;     b1[0] = 1.3;
a1[1][0] = -1.03;    a1[1][1] = 0.28;  a1[1][2] = 0.33;    a1[1][3] = -0.07;    b1[1] = -1.1;
a1[2][0] = -1.09;    a1[2][1] = 0.13;  a1[2][2] = 0.42;    a1[2][3] = 0.28;     b1[2] = 1.7;
a1[3][0] = -1.19;    a1[3][1] = -0.23; a1[3][2] = 0.08;    a1[3][3] = 0.37;     b1[3] = 1.5;


iteration(a,b,x);
for(i=0; i<N; i++)
 cout<<"Rezyltat metoda iterasiy = "<<x[i]<<endl;
}

Выдаёт ошибку a1 Undeclared Identifier, я же объявил как вы сказали.

[Sazary]

А вы внимательно посмотрите на те переменные, которые у вас объявлены, и те, которые вы посылаете в функцию. Не замечаете нестыковку?

Сразу добавлю, что массив x у вас непроинициализирован. А там, вроде как, нужно задать начальное приближение.

[prikolist]

Вы про эту строку имеете ввиду:

iteration(a,b,x);

Ну, если я изменю её вот так,

iteration(a1,b1,x1);

вновь выдаст ошибку.

Тоесть если напишу вот так:

int main()
{

double a1[N][N],b1[N],x1[N];

int i;

a1[0][0] = -1.01;    a1[0][1] = 0.02;  a1[0][2] = -0.62;   a1[0][3] = 0.08;     b1[0] = 1.3;
a1[1][0] = -1.03;    a1[1][1] = 0.28;  a1[1][2] = 0.33;    a1[1][3] = -0.07;    b1[1] = -1.1;
a1[2][0] = -1.09;    a1[2][1] = 0.13;  a1[2][2] = 0.42;    a1[2][3] = 0.28;     b1[2] = 1.7;
a1[3][0] = -1.19;    a1[3][1] = -0.23; a1[3][2] = 0.08;    a1[3][3] = 0.37;     b1[3] = 1.5;


iteration(a1,b1,x);
for(i=0; i<N; i++)
 cout<<"Rezyltat metoda iterasiy = "<<x[i]<<endl;
}

[Sazary]

Сообщение от prikolist

Вы про эту строку имеете ввиду:

iteration(a,b,x);

Ну, если я изменю её вот так,

iteration(a1,b1,x1);

вновь выдаст ошибку.

А что, у вас где-то существует x1? Я вот что-то только x вижу.

Тоесть если напишу вот так:

int main()
{

double a1[N][N],b1[N],x1[N];

int i;

a1[0][0] = -1.01;    a1[0][1] = 0.02;  a1[0][2] = -0.62;   a1[0][3] = 0.08;     b1[0] = 1.3;
a1[1][0] = -1.03;    a1[1][1] = 0.28;  a1[1][2] = 0.33;    a1[1][3] = -0.07;    b1[1] = -1.1;
a1[2][0] = -1.09;    a1[2][1] = 0.13;  a1[2][2] = 0.42;    a1[2][3] = 0.28;     b1[2] = 1.7;
a1[3][0] = -1.19;    a1[3][1] = -0.23; a1[3][2] = 0.08;    a1[3][3] = 0.37;     b1[3] = 1.5;


iteration(a1,b1,x);
for(i=0; i<N; i++)
 cout<<"Rezyltat metoda iterasiy = "<<x[i]<<endl;
}

Не знаю, как вам еще объяснить, что у одной и той же переменной должно быть одно и то же имя.

[prikolist]

Спасибо,теперь скопилировалось, когда написал во таткой код:

#include "stdafx.h"
#include <iostream>
#include <cmath>
#define eps 0.0001 //требуемая точность вычислений
#define N 4
using namespace std;


//----------------------------------------------------------------------------------
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 ;
}


int main()
{

double a1[N][N],b1[N],x[N];

int i;

a1[0][0] = -1.01;    a1[0][1] = 0.02;  a1[0][2] = -0.62;   a1[0][3] = 0.08;     b1[0] = 1.3;
a1[1][0] = -1.03;    a1[1][1] = 0.28;  a1[1][2] = 0.33;    a1[1][3] = -0.07;    b1[1] = -1.1;
a1[2][0] = -1.09;    a1[2][1] = 0.13;  a1[2][2] = 0.42;    a1[2][3] = 0.28;     b1[2] = 1.7;
a1[3][0] = -1.19;    a1[3][1] = -0.23; a1[3][2] = 0.08;    a1[3][3] = 0.37;     b1[3] = 1.5;


iteration(a1,b1,x);
for(i=0; i<N; i++)
 cout<<"Rezyltat metoda iterasiy = "<<x[i]<<endl;
}

Но дело в том,что выдало всё по -1, вот фото:

[Sazary]

Вы что-то в функции изменили. Если взять функцию из первого поста, то что-то считает.

[prikolist]

Попробуйте пожалуйста,я честно говорю, выдаёт все -1, я ничего не изменял.

[Sazary]

Взял функцию из первого поста.

#include <iostream>
#include <cmath>
#define eps 0.0001 //требуемая точность вычислений
#define N 4
using namespace std;


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

 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];
  }
}

int main()
{

double a[N][N],b[N],x[N];
int i;

a[0][0] = -1.01; a[0][1] = 0.02; a[0][2] = -0.62; a[0][3] = 0.08; b[0] = 1.3;
a[1][0] = -1.03; a[1][1] = 0.28; a[1][2] = 0.33; a[1][3] = -0.07; b[1] = -1.1;
a[2][0] = -1.09; a[2][1] = 0.13; a[2][2] = 0.42; a[2][3] = 0.28; b[2] = 1.7;
a[3][0] = -1.19; a[3][1] = -0.23; a[3][2] = 0.08; a[3][3] = 0.37; b[3] = 1.5;

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

iteration(a,b,x);
for(i=0; i<N; i++)
 cout<<x[i]<<endl;
 
cin.get(); 
}

Выдает:

-1.28713
-8.66337
3.38873
-6.20366

Насчет начального приближения точно не знаю. Поэтому просто нулями заполнил.

[prikolist]

Спасибо большое, скажите, а почему при запуске, когда я вставил данный код, выдаёт такое фото ниже, но когда нажать на кнопку "пропустить ", только тогда выдаёт результат.