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

СЛАУ на C++ с использованием ООП - C++

Восстановить пароль Регистрация
 
Mirazhana
1 / 1 / 0
Регистрация: 04.09.2012
Сообщений: 94
21.06.2013, 20:58     СЛАУ на C++ с использованием ООП #1
Помогите пожалуйста, у самой ничего не получается. Нужно написать программу по подсчету системы линейных уравнений третьего порядка с использованием ООП (классы и пр.), метод решения любой (Гаусс/Крамер).
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
21.06.2013, 20:58     СЛАУ на C++ с использованием ООП
Посмотрите здесь:

Метод Гаусса для решения СЛАУ с использованием одномерного массива C++
C++ Решение слау.
C++ игра BlackJAck с использованием ООП
Переписать реализацию стека с использованием ООП C++
C++ Выполнить задания с использованием инструкции выбора if и с использованием условного оператора
C++ СЛАУ
Написать программу для вычисления СЛАУ методом простой итерации с использованием динамических массивов C++
Написать программу с использованием функции и с использованием макроса C++

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
igorrr37
 Аватар для igorrr37
1593 / 1221 / 118
Регистрация: 21.12.2010
Сообщений: 1,868
Записей в блоге: 7
21.06.2013, 21:17     СЛАУ на C++ с использованием ООП #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
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
#include <iostream>
#include <clocale>
#include <stdlib.h>
#include <iomanip>
#include <string>
#include <cstring>
#include <fstream>
 
// Освобождает память из-под матрицы
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("Hilbert_matrix.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;
}
Mirazhana
1 / 1 / 0
Регистрация: 04.09.2012
Сообщений: 94
21.06.2013, 23:10  [ТС]     СЛАУ на C++ с использованием ООП #3
У меня очень большие проблемы с классами и пр. сколько бы я ни читала, никак не могу усвоить и понять как их использовать
Yandex
Объявления
21.06.2013, 23:10     СЛАУ на C++ с использованием ООП
Ответ Создать тему
Опции темы

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