Форум программистов, компьютерный форум CyberForum.ru

Решение СЛАУ методом простых итераций - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 39, средняя оценка - 4.87
paraphin
1 / 1 / 0
Регистрация: 09.12.2009
Сообщений: 50
26.12.2011, 23:36     Решение СЛАУ методом простых итераций #1
Нужно решить слау (4 на 4), значения матрица и вектор вводится с клавиатуры, также нужно ввести епсилон (точность, до которой необходимо делать итерации). Помогите плиз составить код, а то я скоро наверное уже моск сломаю... Только без функци и т.п., в инт мейне сразу.. В визуал 2008 делаю. Заранее спасибо.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
26.12.2011, 23:36     Решение СЛАУ методом простых итераций
Посмотрите здесь:

C++ Решение системы линейных уравнений методом простых итераций
Решение нелинейного уравнения методом простых итераций C++
C++ Решение СЛАУ методом простых итераций.
C++ Решение СЛАУ методом простых итераций, задача решена, только написать код осталось
C++ Метод простых итераций для СЛАУ
C++ Решение СЛАУ методом простых интераций
C++ Алгоритм написания для решения СЛАУ методом простых итераций
C++ Решение нелинейного уравнения методом простых итераций

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
edward_jonson
 Аватар для edward_jonson
157 / 157 / 25
Регистрация: 23.02.2011
Сообщений: 388
27.12.2011, 01:26     Решение СЛАУ методом простых итераций #2
http://ru.wikipedia.org/wiki/Метод_Я...B8.D1.82.D0.BC
case
7 / 7 / 2
Регистрация: 29.09.2010
Сообщений: 32
27.12.2011, 01:38     Решение СЛАУ методом простых итераций #3
C++
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
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
#include "stdafx.h"
#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include <math.h>
#define n 4
 
void del(double **A, double *X, double *Xk,   double *b);
 
int main()
 
{
 
    int i,j;
    double **A;     //matrix equations
    double *X;      //communications solution and iteration (k+1)
    double *Xk;     //iteration (k)
    double *b;      //free to members of
    double sumd = 0;//additional variables for computation
    double eps;
 
//--------------dynamic memory allocation------------------------------//
    A  = new double* [n];
    X  =  new double [n];
    Xk =  new double [n];
    b  =  new double [n];
    for(i = 0; i < n; i++)
        A[i] = new double [n];
//-------------input system of equations-------------------------------//
    printf("Type elements in the variables in the system: \n");
    for(i = 0; i < n; i++)
       for(j = 0; j < n; j++)
        scanf("%lf",&A[i][j]);
 
 
    printf("Members enter free system: \n");
    for(i = 0; i < n; i++)
        scanf("%lf",&b[i]);
    printf("enter the calculation accuracy: \n");
        scanf(" %lf",&eps);
//-------------test system on the diagonal elements--------------------//
    for(i = 0; i < n; i++)
       if (A[i][i] == 0) {del( A, X, Xk, b); printf("!On the main diagonal of the system is zero"); exit(0);}
//-------------checking for convergence of-----------------------------//
    for(i = 0; i < n; i++)
      {
       for(j = 0; j < n; j++)
          if (i != j)       sumd += fabs(A[i][j]);
       if (fabs(A[i][i]) < sumd)  {del( A, X, Xk, b);  printf("!Diverging methods for this system"); exit(0);}
       sumd = 0;
      }
//----- approximate solutions of the system-//
    for(i = 0; i < n; i++)
         X[i]=0.0;
    
//-------------calculation system of linear algebraic equations---------//
    int  count = 0, flag = 1;
    double x, zh = 0;
    for(i = 0; i < n; i++)
    Xk[i] = X[i];
    do
      {
       count++;
       for(i = 0; i < n; i++)
     {
       x=0;
       for(j = 0; j < n; j++)
          {
         if (i != j)  x += Xk[j] * A[i][j];
         if (i == j) zh = A[i][j];
          }
       x = (b[i] - x) / zh;
       X[i] = x;
       if ((fabs(X[i]-Xk[i]))<=eps) flag=0;
     }
     for(i = 0; i < n; i++)
        Xk[i] = X[i];
      }while(flag);
    printf("solution call method of iterations: \n");
    for(i = 0; i < n; i++)
       printf("x[%d]=%0.5lf\n",i+1,X[i]);
       printf("count=%d\n",count);
    
       del( A, X, Xk, b);
       return 0;
}
 
 
void del(double **A, double *X, double *Xk,   double *b)
  {
    for(int i = 0; i < n; i++)
      delete[]A[i];
     delete[]A;
     delete[]X;
     delete[]Xk;
     delete[]b;
  }
Yandex
Объявления
27.12.2011, 01:38     Решение СЛАУ методом простых итераций
Ответ Создать тему
Опции темы

Текущее время: 12:31. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru