Вот наконец написал курсовую. Каряво, грубо, но собственными силами, чем собственно и горжусь)
Программа считает СЛАУ методом Гаусса (с поиском главного элемента по строке). Считает все отлично, находит корни и тд.
Проблема: препад захотел, что бы в програмулине был ввод с файла. Дело в том, что массивы я использую статические, а везде где не искал в основном описывают ввод/вывод динамических массивов. Дык вот я вроде смазговал что бы вводился массив статически, и программа файл вроде читает и на экран выводит, а вот дальше мой массив 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
от балды... Программа их все равно решает если вводить с клавиатуры
Зарание спасибо!