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

СЛАУ методом Гаусса с поиском главного элемента - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 12, средняя оценка - 4.75
sa80956847450d
0 / 0 / 0
Регистрация: 27.05.2010
Сообщений: 9
28.05.2011, 18:11     СЛАУ методом Гаусса с поиском главного элемента #1
Вот наконец написал курсовую. Каряво, грубо, но собственными силами, чем собственно и горжусь)
Программа считает СЛАУ методом Гаусса (с поиском главного элемента по строке). Считает все отлично, находит корни и тд.

Проблема: препад захотел, что бы в програмулине был ввод с файла. Дело в том, что массивы я использую статические, а везде где не искал в основном описывают ввод/вывод динамических массивов. Дык вот я вроде смазговал что бы вводился массив статически, и программа файл вроде читает и на экран выводит, а вот дальше мой массив a[i][j] не перехватывается((

Сам листинг:

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
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
#include <iostream>
#include <conio.h>
#include <math.h>
#include <fstream>
using namespace std;
 
void main()
{
    setlocale(0,"Russian");
    double a[50][51], x[50], t;
    int k, m, n;
 
    
    cout<<"\n\tЗдравствуйте, Вас приветствует программа решения СЛАУ!!!\n\tПожалуйста, следуейте вспомогательными подсказками =) \n\n\tДля продолжения нажмите ENTER...";
    getch();
    system("cls");
    int key;
    cout<<"Для загрузки данных из файла нажмите 1\nДля ввода матрицы нажмите 2"<<endl;
    cin>>key;
 
    if (key==1)
    {
        int n=3, m=4;
        ifstream Proga("file.txt");
    //  double **a=new double *[n];
            for (int i = 0; i < n; ++i)
                for (int j = 0; j < m; ++j)
                    Proga>>a[i][j];
            for (int i = 0; i < n; ++i)
            {
                for (int j = 0; j < m; ++j)
                    //cout<<a[i][j];
                    {printf("%6.2f\t", a[i][j]);}
        printf("\n");
            }
 
 
        
 
    }
    else 
    {
        cout<<"\n\tВведите размерность Вашей матрицы СЛАУ: ";
    cin>>n;
    m=n+1;
    //a[0][0]= 1; a[0][1]= 1; a[0][2]= 1; a[0][3]= 2;
    //a[1][0]= 2; a[1][1]=2; a[1][2]=2; a[1][3]= 3;
    //a[2][0]= 1; a[2][1]= -1; a[2][2]= 2; a[2][3]= 5;
    system("cls");
    cout<<"Введите значения матрицы размерностью: "<<n<<endl;
    cout<<"Коэффициенты следует вводить построчно, как они стоят в заданной СЛАУ!"<<endl; 
        for (int i = 0; i < n; ++i)
            for (int j = 0; j < m; ++j)
               cin>> a[i][j] ;
    }
// проверка на бесконечное множество решений
    for (int j = 0; j < n; j++)
    if (a[j][j] == 0)
    {
      k = j;
      while ((a[k+1][j] == 0) && (k < n))
    k++;
      if (a[k+1][j] != 0)
    for (int i = 0; i < m; i++)
    {
      t = a[j][i];
      a[j][i] = a[k+1][i];
      a[k+1][i] = t;
    }
      else
      {
    printf("СЛАУ имеет множество решений");
    _getch();
    exit(1);
      }
    }
//проверка на отсутсвие решений СЛАУ
  for (k = 0; k < n - 1; k++)
    for (int j = k + 1; j < n; j++)
    {
      if (a[k][k] !=0)
      {
    t = a[j][k] / a[k][k];
    for (int i = 0; i < n + 1; i++)
      a[j][i] = a[k][i] * t - a[j][i];
      }
      else 
      {
        printf("СЛАУ не имеет решений\n");
        _getch();
      exit(1);
      }
    }
 
    system("pause");
    system("cls");
    //основной цикл построение диоганальной матрицы - прямой проход обнуление(исключение) элементов ниже главной диагонали
    int q;
    cout<<"Вы хотите проследить ход решения данного СЛАУ?\n\t1 - да\n\t2 - нет"<<endl;
    cin>>q;
    for(k = 0; k<n; k++)
    {
        int ch_col = k;
        //search
        for(int j = k; j<m-1; j++) // ищем главный элемент в строке k
            if(abs(a[k][j])>abs(a[k][ch_col])) ch_col = j;
 
        double bufff;
        for(int i = 0; i<n; i++)// меняем местами столбцы где главный найденный элемент и столбец которому принадлежит диагональный элемент
        {
            bufff = a[i][k];
            a[i][k] = a[i][ch_col];
            a[i][ch_col] = bufff;
        }
 
        for(int e = k+1; e<n; e++)//вычитаем из строк ниже катой строки коэф.. 
        {
            double kof = a[e][k]/a[k][k];
            for(int l = 0; l<m; l++)
                a[e][l] = a[e][l] - a[k][l]*kof; 
        }
        
        if (q==1)///смотрим или не смотрим ход решения СЛАУ
        {
        for (int j = 0; j < n; j++)
        {
            for (int i = 0; i < n + 1; i++)
            printf("%6.2f\t", a[j][i]);
            printf("\n");
        }}
        else 
        cout << endl;
        cout << endl;
    }
    cout<<"Ваша матрица приведена к диагональному виду:\n"<<endl;
    for (int j = 0; j < n; j++)// вывод диагональной матрицы
    {
        for (int i = 0; i < n + 1; i++)
            {printf("%6.2f\t", a[j][i]);}
        printf("\n");
    }
 
    for (int j = n - 1 ; j >= 0; j--) // вычисление корней
    {
        t = a[j][n];
        for (int i = j + 1; i < n; i++)
        t= t - a[j][i] * x[i];
        x[j] = t / a[j][j];
    }
 
    printf("\n Корни СЛАУ:\n\n");
    for (int i = 0; i < n; i++)
    printf("x%d = %6.2f\n", i, x[i]);
 
    _getch();
 
}
Да! В файлике у меня вот такая ересь:
1 2 4 5
7 8 9 8
6 5 3 2
от балды... Программа их все равно решает если вводить с клавиатуры

Зарание спасибо!
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
28.05.2011, 18:11     СЛАУ методом Гаусса с поиском главного элемента
Посмотрите здесь:

Метод Гаусса с выбором главного элемента C++
Решение СЛАУ методом Гаусса C++
Найти решение системы линейных уравнений методом Гаусса с выбором главного элемента в столбце C++
Метод Гаусса с выбором главного элемента по строке C++
Решение СЛОУ методом Гаусса с выбором главного элемента по строке C++
C++ Решение СЛАУ методом Гаусса
СЛАУ методом Гаусса C++
СЛАУ методом Гаусса с поиском ведущего элемента по всей матрице C++

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

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