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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 12, средняя оценка - 4.75
sa80956847450d
0 / 0 / 0
Регистрация: 27.05.2010
Сообщений: 9
#1

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

28.05.2011, 18:11. Просмотров 1477. Ответов 0
Метки нет (Все метки)

Вот наконец написал курсовую. Каряво, грубо, но собственными силами, чем собственно и горжусь)
Программа считает СЛАУ методом Гаусса (с поиском главного элемента по строке). Считает все отлично, находит корни и тд.

Проблема: препад захотел, что бы в програмулине был ввод с файла. Дело в том, что массивы я использую статические, а везде где не искал в основном описывают ввод/вывод динамических массивов. Дык вот я вроде смазговал что бы вводился массив статически, и программа файл вроде читает и на экран выводит, а вот дальше мой массив 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++
СЛАУ методом Гаусса с выбором главного элемента по столбцу.Помогите пожалуйста с кодом программы,ошибки выдает #include &lt;locale&gt; ...

СЛАУ методом Гаусса с поиском ведущего элемента по всей матрице - C++
Методом Гаусса все легко получилось, но вот только задание у меня включает поиск ведущего элемента (самого большого в матрице). ...

Решение СЛОУ методом Гаусса с выбором главного элемента по строке - C++
Задание:Решение СЛОУ методом Гаусса с выбором главного элемента по строке.Для тестирования решить систему. Запускается,но вывод...

Решение систем линейных алгебраических уравнений методом Гаусса с выбором главного элемента - C++
Помогите пожалуйста найти решение системы линейных уравнений с матричными элементами. Условия задания матрицы в файле.

Метод Гаусса с выбором главного элемента - C++
Нужно написать программу, решающую СЛАУ методом Гаусса, так же надо учесть, что матрица может быть вырожденной, первый элемент может быть...

Метод Гаусса с выбором главного элемента по строке - C++
Метод Гаусса с выбором главного элемента по строке. Помогите разобраться в чем косяки программы. #include &lt;iostream&gt; #include...

СЛАУ методом Жордана-Гаусса - C++
Нужна программа для решения СЛАУ методом Жордана-Гаусса с выбором главного элемента.

Реализовать программу, решающую СЛАУ методом Гаусса - C++
Реализовать программу, решающую СЛАУ методом Гаусса. Требования: 1. Программа должна выдавать решение для СЛАУ, имеющих единственное...

Нужен рабочий пример решения СЛАУ методом Гаусса - C++
Нужен рабочий (ну или более менее рабочий) код для решения СЛАУ методом Гауса на С++. Уже замучился искать

Квадратичная аппроксимация, ошибка при решении СЛАУ методом Гаусса - C++
всем привет, возникла проблема при квадратичной аппроксимации функции, в результате левые данные... сама функция имеет вид у = х^0.25, а...


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

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

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