Форум программистов, компьютерный форум, киберфорум
Наши страницы
С++ для начинающих
Войти
Регистрация
Восстановить пароль
 
bobiko
0 / 0 / 0
Регистрация: 16.11.2015
Сообщений: 59
1

Метод Гаусса посредством использования valarray

14.12.2017, 00:47. Просмотров 127. Ответов 1

Здравствуйте, дорогие форумчане!
Задался идеей реализовать метод гаусса посредством использования valarray, ну уже пожалел об этом. Идея в чем: хотел использовать 1 vallarray для всех вычислений, но не вижу это возможным, поэтому использую 1 valarray для хранения всех матричных значений + ответы и 2 для хранения временных значений, чтобы при преобразовании их не потерять. Для простоты решил задать все значения сразу, прямой ход делал без хранения ответов в виде нулей и все получалось, а как только добавил их, сразу несуразные неопределенности в выводе. Как я понимаю, где-то оно дальше считает и не видит выхода, но хоть убей, не могу найти ошибку. Потом еще нужен обратный ход же сделать. Вопрос: я правильно начал делать его, или нет?
Исходник
Кликните здесь для просмотра всего текста
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
#include "stdafx.h"
#include <iostream>
#include <stdlib.h>
#include <valarray>
 
using namespace std;
 
int main()
{
    int i, j, n, m, a, b, c, d, k, l;
    cout << "Number of unknown quantities: " << endl;
    cin >> m;
    n = (m + 1) *m;
    valarray<float> gauss = { 2, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 0, 0, 0 }; //там, где нули будут сохранятся ответы
    valarray<float> tmp1(1);
    valarray<float> tmp2(1);
    l = gauss.size();
 
    for (i = 0; i < m; i++) {
        for (j = 0; j < m - i; j++) {
                a = (m + 1) * i;            //значимый диагональный элемент строки
                b = a + ((m + 1)*(j + 1));  // индекс соответствующей ему неизвестной в след.строке
                d = a + ((m + 1)*(j + 1));  //индексы строки
                    tmp1 = gauss[slice(a + i, 1, 1)]; //сохраняем а, чтобы не потерять при изменениях
                    tmp2 = gauss[slice(b + i, 1, 1)]; //также с b
                    for (k = 0; k < (m + 1) - i; k++) { gauss[slice(a + i + k, 1, 1)] /= tmp1; }  //делим первую строку на а
 
                    for (k = 0; k < (m + 1) - i; k++) { gauss[slice(a + i + k, 1, 1)] *= tmp2; } //умножаем первую строку на b
 
                    for (k = 0; k < (m + 1) - i; k++) { gauss[slice(d + i + k, 1, 1)] -= gauss[slice(a + i + k, 1, 1)]; } // отнимает от второй строки преобразованную первую
 
                    for (k = 0; k < (m + 1) - i; k++) { gauss[slice(a + i + k, 1, 1)] /= tmp2; } // возвращаем вид первой в исходное
 
 
                    for (size_t i = 0; i < 4; i++) { cout << gauss[i] << ' '; }cout << endl;
                    for (size_t i = 4; i < 8; i++) { cout << gauss[i] << ' '; }cout << endl;
                    for (size_t i = 8; i < 12; i++) { cout << gauss[i] << ' '; } cout << endl << endl;
                }
 
    }
    //якобы обратный ход
    for (i = 0; i <= m; i++) {
        tmp1 = gauss[slice(n - 1, 1, 1)];
        for (j = 1; j <=i; j++) {
            a = (n - 1) - ((m + 1) * i);
            tmp1 = gauss[slice(a, 1, 1)];
            tmp2 = gauss[slice(a - j, 1, 1)];
            tmp2 *= gauss[slice(l - (j - 1), 1, 1)];
            tmp1 -= tmp2;
        }
        gauss[slice(l - i, 1, 1)] = tmp1;
    }
 
    for (size_t i = 0; i < 4; i++) { cout << gauss[i] << ' '; }cout << endl;
    for (size_t i = 4; i < 8; i++) { cout << gauss[i] << ' '; }cout << endl;
    for (size_t i = 8; i < 12; i++) { cout << gauss[i] << ' '; } cout << endl << endl;
    for (size_t i = l - m; i < l; i++) {
        cout << ' ' << gauss[i] << endl;
    }
    //check --- http://matematikam.ru/solve-equations/sistema-gaus.php
    return 0;
}


Добавлено через 3 часа 35 минут
Все, тему можно закрывать, уже все решил.) Волшебный форум какой-то, ведь как только написал, так сразу все замечать начал)
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
14.12.2017, 00:47
Ответы с готовыми решениями:

СЛАУ. Метод обратной матрицы, метод Гаусса, метод Крамера, метод Зейделя
Помогите ребят. Не могу построить алгоритмы для этих методов Язык C++

Метод Гаусса
Всем привет. Ребят, что скажите по поводу такой реализации метода Гаусса: ...

метод гаусса
найти метод гаусса n+1 значение определителя в n+1 точке, а затем решить...

Метод Гаусса в С++
Подскажите please, как запустить эту программу, через C++??? Решение системы...

Метод Гаусса
Как решить квадратную матрицу методом Гаусса, я ее привел к треугольному виду,...

1
AndrSlav
65 / 53 / 14
Регистрация: 20.12.2013
Сообщений: 459
15.12.2017, 22:22 2
Понимаю, что уже не надо, но я как раз пытаюсь изучать библиотеку, так что попробовал решить тоже СЛАУ методом Гаусса. Может кто что подскажет насчет скрытых ошибок. Так вроде решается.
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
#pragma hdrstop
#pragma argsused
 
#ifdef _WIN32
#include <tchar.h>
#else
  typedef char _TCHAR;
  #define _tmain main
#endif
 
#include <stdio.h>
#include <iostream>
#include <vector>
#include <algorithm>
#include <functional>
#include <string>
#include <memory>
#include <math.h>
#include <valarray>
#include <iterator>
   using std::vector;
   using std::valarray;
   using std::slice;
   using std::slice_array;
 
   int amax_val(const valarray<double>& v, vector<bool> trig)
   {int m= 0;
    for(int i=0; i<v.size(); ++i, fabs(v[i])>fabs(v[m]) && trig[i] ? m=i : m);
    return m;
   };
 
 int _tmain(int argc, _TCHAR* argv[])
{
  int dim= 3;
  vector<int> ind(3);
  vector<bool> trig;
  for(int i=0; i<dim; ++i, trig.push_back(true));
 
  std::cout<<"Initial matrix:";
  double mas[]= {1,5,4,2,9,15,7,2,4,3,8,2};
  valarray<double> matr(mas, 12);    //Матрица коэффициентов
  valarray<double> col(dim);         //Столбец
  valarray<double> rr(dim+1);        //Строка
  vector<double> rez;
  int m;
 
          for(int i=0; i<matr.size(); ++i)
          {
           if(!(i%(dim+1))) std::cout<<'\n';
           std::cout<<matr[i]<<"  ";
          };
 
  for(int i=0; i<dim; ++i)
  {
   col= matr[slice(i, dim, dim+1)];
   m= amax_val(col, trig);
   ind[i]= m;
   trig[m]= false;
   const slice_array<double>& row= matr[slice(m*(dim+1), dim+1, 1)];
   row/= (rr= col[m]);
 
        for(int j=0; j<dim; ++j)
        {
             if(j!=m)
             {
              const slice_array<double>& row2= matr[slice(j*(dim+1), dim+1, 1)];
              row2/= (rr= col[j]);
              row2-= matr[slice(m*(dim+1), dim+1, 1)];
             };
        };
  };
 
  for(int i=0; i<dim; ++i)
  {
   const slice_array<double>& row= matr[slice(ind[i]*(dim+1), dim+1, 1)];
   row/= (rr= matr[i+ ind[i]*(dim+1)]);
   rez.push_back(matr[(ind[i]+1)*(dim+1)- 1]);
  };
 
 
  std::cout<<'\n'<<'\n'<<"Rezult: ";
  for(int i=0; i<matr.size(); ++i)
  {
   if(!(i%(dim+1))) std::cout<<'\n';
   std::cout<<matr[i]<<"  ";
  };
 
  std::cout<<'\n'<<'\n';
  std::cout<<rez[0]<<" "<<rez[1]<<" "<<rez[2]<<'\n';
 
  char ch;
  std::cin>>ch;
 
    return 0;
}
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
15.12.2017, 22:22

Метод Гаусса
#include&lt;stdio.h&gt; #include&lt;iostream.h&gt; #include&lt;conio.h&gt; #include&lt;stdlib.h&gt;...

Метод Гаусса
Суть в том что бы СЛАУ методом Гаусс с последовательным исключениям...

Метод Гаусса
Добрый вечер, можете написать программу для решения СЛАУ методом гаусса на Си,...


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

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Рейтинг@Mail.ru