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

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

Восстановить пароль Регистрация
Другие темы раздела
C++ Метод хорд http://www.cyberforum.ru/cpp-beginners/thread1044130.html
Есть задание:"Методом хорд найти корень уравнения x-3{x}^{3}+12=0 на отрезке . Написал,но значение выводит неправильное,т.к. отрицательным не может быть по определению,где ошибка?Помогите, пожалуйста. #include "stdafx.h" #include <stdio.h> #include <math.h> int main() { float a,b,e; printf("\nVvedite predel a:\n");
C++ задача на тему "массивы и строки " сколько раз в заданном предложении встретится одинаковая первая буква в словах. Напечатать исходный текст и число сочетаний. Помогите пожалуйста http://www.cyberforum.ru/cpp-beginners/thread1044126.html
C++ Двумерные массивы
Ребят помогите разработать программу, нужно найти самый успешный квартал для каждого года. Ре-зультат вывести в виде таблицы из двух столбцов: год, номер квартала.Если можно напишите как можно проще на уровне начинающего.Заранее спасибо!
C++ что исправить? нужно вывести значения y
#include "stdafx.h" #include "iostream" using namespace std; void main () { float x, y, r; r=1; for (x=0; x<=20; x=x+0.5) { if (x=0&&x<=2)
C++ Добавить к строке число http://www.cyberforum.ru/cpp-beginners/thread1044061.html
Как добавить к строке число? Пользователь вводит число и потом оно добавляется к строке height wchar_t height = L"height = "; чтобы получилось примерно так "height = 100"
C++ Массив как поле класса Например, имеется класс с полем массивом размерностью 10 на 10. class massiv { private: int mas; public: ... }; Конструктор класса должен позволять программисту определять реальный размер массива(например, сделать его меньше чем 10 на 10). Как это реализовать? подробнее

Показать сообщение отдельно
igorrr37
 Аватар для igorrr37
1593 / 1221 / 118
Регистрация: 21.12.2010
Сообщений: 1,868
Записей в блоге: 7
17.12.2013, 14:54     Решение СЛОУ методом Гаусса с выбором главного элемента по строке
делал давно, но, кажется, это и есть с выбором макс. элемента. Матрицу читает из файла (число строк, число столбцов и сама матрица), пример файла
Bash
1
2
3
4
5
2
3
 
1     2.2   3
3.1   5     6
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
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
#include <iostream>
#include <clocale>
#include <stdlib.h>
#include <iomanip>
#include <string>
#include <cstring>
#include <fstream>
#include <windows.h>
 
// Освобождает память из-под матрицы
template<typename T>
void FreeMemory(T** pmatrix, size_t const rows)
{
    for(size_t i = 0; i < rows; ++i)
    {
        delete[] pmatrix[i];
        pmatrix[i] = 0;
    }
    delete[] pmatrix;
    pmatrix = 0;
}
 
// Выделение памяти под матрицу
template<typename T>
T**& GetMemory(T**& pmatrix, size_t const rows, size_t const columns)
{
    size_t i = 0;
    try
    {
        pmatrix = new T*[rows];
        for(i = 0; i < rows; ++i)
            pmatrix[i] = 0;
        for(i = 0; i < rows; ++i)
            pmatrix[i] = new T[columns];
    }
    catch(std::bad_alloc const& exc)
    {
        std::cerr << "GetMemory: Exception: Неуспешное выделение памяти\n\n";
        FreeMemory(pmatrix, rows);
        throw;
    }
    return pmatrix;
}
 
// Заполняет матрицу из файла
template<typename T>
T**& FillFromFile(std::string const& filepath, T**& pmatrix, size_t& rows, size_t& columns)
{
    try
    {
        std::ifstream ifs(filepath.c_str());
        if(ifs.is_open())
        {
            ifs >> rows >> columns;
            if(!rows || !columns)
                throw "Ошибка чтения файла: !rows || !columns\n";
            GetMemory(pmatrix, rows, columns);
            for(size_t i = 0; i < rows; ++i)
            {
                for(size_t j = 0; j < columns; ++j)
                {
                    ifs >> pmatrix[i][j];
                }
            }
            if(!ifs)
                    throw "Ошибка чтения файла: !ifs\n";
            ifs.close();
        }
        else
            std::cerr << "FillFromFile: неуспешное открытие файла\n";
    }
    catch(char const* e)
    {
        std::cerr << "FillFromFile: Exception: " << e << "\n\n";
        throw;
    }
    return pmatrix;
}
 
// Вывод матрицы в поток
template<typename T>
std::ostream& MatrixToStream(T* const* const pmatrix, size_t const rows, size_t const columns, std::ostream& ost = std::cout)
{
    ost << "\n-------------------------MatrixToStream-------------------------------\n";
    ost << "Размер: " << rows << " x " << columns << "\n";
    for(size_t i = 0; i < rows; ++i)
    {
        for(size_t j = 0; j < columns; ++j)
        {
            ost << std::setw(13) << std::left << pmatrix[i][j];
        }
        ost << "\n\n";
    }
    ost << "\n----------------------------------------------------------------------\n";
    return ost;
}
 
// Умножение массива на число
double* MultiplicationArrayToNumber(double* parray, size_t const columns, double const number)
{
    for(size_t j = 0; j < columns; ++j)
    {
        parray[j] = parray[j] * number;
    }
    return parray;
}
 
// Вычитание массива 2 из массива 1 (они равной длины)
double* SubtractionOfArrays(double* parray1, double* parray2, size_t const columns)
{
    for(size_t j = 0; j < columns; ++j)
    {
        parray1[j] = parray1[j] - parray2[j];
    }
    return parray1;
}
 
// Обмен двух массивов
void SwapArrays(double* parray1, double* parray2, size_t const columns)
{
    double* ptmparray = new double[columns];
    memcpy(ptmparray, parray1, columns * sizeof(double));
    memcpy(parray1, parray2, columns * sizeof(double));
    memcpy(parray2, ptmparray, columns * sizeof(double));
    delete[] ptmparray;
    ptmparray = 0;
}
 
// Поиск строки среди текущей и нижележащих строк, содержащей в заданном столбце максимальный элемент
size_t MaxRowIndex(double** pdmatrix, size_t const rows, size_t const columnindex)
{
    size_t maxrowindex = columnindex, i, j;
    for(i = columnindex + 1; i < rows; ++i)
    {
        if(pdmatrix[i][columnindex] > pdmatrix[columnindex][columnindex])
            maxrowindex = i;
    }
    return maxrowindex;
}
 
// Приведение матрицы к верхнетреугольному виду
size_t MakeTriangularMatrix(double** pdmatrix, size_t const rows, size_t const columns)
{
    size_t rootnumber = rows;
    size_t i = 0, j = 0, k = 0, maxrowindex = 0;
    double* prowcopy = new double[columns];
    for(j = 0; j < rows - 1; ++j)
    {
        maxrowindex = MaxRowIndex(pdmatrix, rows, j);
        if(maxrowindex != j)
            SwapArrays(pdmatrix[maxrowindex], pdmatrix[j], columns);
        if(double(0) != pdmatrix[j][j])
        {
            MultiplicationArrayToNumber(pdmatrix[j], columns, 1 / pdmatrix[j][j]);
            for(i = j + 1; i < rows; ++i)
            {
                memcpy(prowcopy, pdmatrix[j], columns * sizeof(double));
                MultiplicationArrayToNumber(prowcopy, columns, pdmatrix[i][j]);
                SubtractionOfArrays(pdmatrix[i], prowcopy, columns);
            }
        }
    }
    delete[] prowcopy;
    prowcopy = 0;
    // Проверка системы на вырожденность
    for(i = 0; i < rows; ++i)
    {
        bool zerorow = true;
        for(j = 0; j < columns - 1; ++j)
        {
            if(pdmatrix[i][j] != double(0))
            {
                zerorow = false;
                break;
            }
        }
        if(zerorow && pdmatrix[i][columns - 1] != double(0))
        {
            std::cout << "\nCистема НЕСОВМЕСТНА\n\n";
            rootnumber = 0;
            break;
        }
        else if(zerorow && pdmatrix[i][columns - 1] == double(0))
        {
            std::cout << "\nCистема имеет БЕСКОНЕЧНОЕ МНОЖЕСТВО РЕШЕНИЙ\n\n";
            rootnumber = rows + 5;
            break;
        }
    }
    return rootnumber; // Возвращает число корней
}
 
//Вычисляет корни СЛАУ по верхнетреугольной матрице
double* GetRoots(double** pdmatrix, size_t const rows, double* pdroots)
{
    int i = 0, j = 0;
    for(i = rows - 1; i >= 0; --i)
    {
        double sum = 0;
        for(j = rows - 1; j > i; --j)
            sum += pdroots[j] * pdmatrix[i][j];
        pdroots[i] = (pdmatrix[i][rows] - sum) / pdmatrix[i][i];
    }
    return pdroots;
}
 
int main()
{
    setlocale(LC_ALL, "rus");
    std::string const filepath("in.txt"); // файл с матрицей
    size_t rows = 0, columns = 0, i = 0, j = 0;
    double** pdmatrix = 0;
    try
    {
        FillFromFile(filepath, pdmatrix, rows, columns);
        if(rows + 1 == columns)
        {
            std::cout << "\n\nИсходная матрица:\n";
            MatrixToStream(pdmatrix, rows, columns);
            size_t rootnumber = MakeTriangularMatrix(pdmatrix, rows, columns);
            std::cout << "\n\nТреугольная матрица:\n";
            MatrixToStream(pdmatrix, rows, columns);
            if(rows == rootnumber)
            {
                double* pdroots = new double[rows];
                GetRoots(pdmatrix, rows, pdroots);
                std::cout << "\n\nКорни системы:\n\n";
                for(i = 0; i < rows; ++i)
                    std::cout << "X" << i << ": " << pdroots[i] << '\n';
                delete[] pdroots;
                pdroots = 0;
            }
            else
                std::cout << "Матрица вырождена\n";
        }
        else
            std::cerr << "Не соблюдено: rows + 1 == columns\n";
        FreeMemory(pdmatrix, rows);
    }
    catch(...){}
    return 0;
}
 
Текущее время: 22:52. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru