Форум программистов, компьютерный форум, киберфорум
Visual C++
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.73/11: Рейтинг темы: голосов - 11, средняя оценка - 4.73
0 / 0 / 0
Регистрация: 13.11.2011
Сообщений: 10

Нахождение обратной матрицы методом Крамера

13.11.2011, 20:20. Показов 2148. Ответов 3
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Нужна помощь. Есть программа,которая реализует нахождение обратной матрицы методом Крамера..Программа работает..Ищет все верно. НО одна проблема.Когда ввожу число например 0.9 программа не работает так,как надо. Когда ввожу вместо 0.9 например 1 и выше,все работает верно..Помогите исправить..Так как мне нужно ввести именно матрицу типа:
1 1 1... ...1
0.9 1 1 ... ...1
0.9 0.99 1... ...1
... ... ... ...
0.9 0.99 0.999 ...1

Код программы:
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
#include "stdafx.h"
 
#include <iostream>
#include <ctime>
#include <cmath>
 
using namespace std;
 
template <typename T> void FreeMem(T **matr, int n);
template <typename T> void PrintMtx(T **matr, int n);
template <typename T> void SetMtx(T **matr, int n);
template <typename T> void TransponMtx(T **matr, T **tMatr, int n);
void Get_matr(int **matr, int n, int **temp_matr, int indRow, int indCol);
int Det(int **matr, int n);
 
void main()
{
        srand((unsigned)time(NULL));
        setlocale(0, "");
        int n, det;
        cout << "Введите размер матрицы: ";
        cin >> n;
             cout << "Введите матрицу коэфициеэнтов" << endl;
     cout << "Вводите в строчку коэффициенты, разделяя пробелами." << endl;
        int **matr = new int * [n];
        double **obr_matr = new double * [n];
        double **tobr_matr = new double * [n];
        for(int i = 0; i < n; i++){
                matr[i] = new int[n];
                obr_matr[i] = new double[n];
                tobr_matr[i] = new double[n];
        }
        SetMtx(matr, n);
        PrintMtx(matr, n);
        det = Det(matr, n);
        cout << "Определитель матрицы = " << det << endl;
        if(det){
                for(int i = 0; i < n; i++){
                        for(int j = 0; j < n; j++){
                                int m = n - 1;
                                int **temp_matr = new int * [m];
                                for(int k = 0; k < m; k++)
                                        temp_matr[k] = new int[m];
                                Get_matr(matr, n, temp_matr, i, j);
                                obr_matr[i][j] = pow(-1.0, i + j + 2) * Det(temp_matr, m) / det;
                                FreeMem(temp_matr, m);
                        }
                }       
        }
        else
                cout << "Т.к. определитель матрицы = 0,\nто матрица вырожденная и обратной не имеет!!!" << endl;
        //Транспонирование матрицы
        TransponMtx(obr_matr, tobr_matr, n);
        //Печать обратной матрицы после транспонирования
        PrintMtx(tobr_matr, n);
        FreeMem(tobr_matr, n);
        FreeMem(matr, n);
        FreeMem(obr_matr, n);
}
//Функция транспонирования матрицы
template <typename T> void TransponMtx(T **matr, T **tMatr, int n){//
    for (int i = 0; i < n; i++) 
        for (int j = 0; j < n; j++)
            tMatr[j][i] = matr[i][j];
}
//Функция освобождения памяти
template <typename T> void FreeMem(T **matr, int n)
{
        for(int i = 0; i < n; i++)
                delete [] matr[i];
        delete [] matr;
}
 
//функция заполнения матрицы
template <typename T> void SetMtx(T **matr, int n)
{
        for (int i = 0; i < n; i++)
 
                for (int j = 0; j < n; j++)
                        cin>>matr[i][j];  
}
 
//функция печати матрицы
template <typename T> void PrintMtx(T **matr, int n)
{
        for (int i = 0; i < n; i++){
                for (int j = 0; j < n; j++)
                        cout << matr[i][j] << " ";
                cout << endl; 
        }
}
//функция вычеркивания строки и столбца
void Get_matr(int **matr, int n, int **temp_matr, int indRow, int indCol)       
{
        int ki = 0; 
        for (int i = 0; i < n; i++){
                if(i != indRow){
                        for (int j = 0, kj = 0; j < n; j++){
                                if (j != indCol){
                                        temp_matr[ki][kj] = matr[i][j];
                                        kj++;
                                }
                        }
                        ki++;           
                }
        }
}
 
//==============================================================================================================
//                                                              вычисление определителя
//==============================================================================================================
//функция вычисления определителя матрицы
int Det(int **matr, int n)      
{
        int temp = 0;   //временная переменная для хранения определителя
        int k = 1;              //степень
        if(n < 1){
                cout<<"Не верный размер матрицы!!!" << endl;
        return 0;
    }
        else if (n == 1)
                temp = matr[0][0];
        else if (n == 2)
                temp = matr[0][0] * matr[1][1] - matr[1][0] * matr[0][1];
        else{
                for(int i = 0; i < n; i++){
                        int m = n - 1; 
                        int **temp_matr = new int * [m];
                        for(int j = 0; j < m; j++)
                                temp_matr[j] = new int [m];
                        Get_matr(matr, n, temp_matr, 0, i);
                        temp = temp + k * matr[0][i] * Det(temp_matr, m);
                        k = -k;
                        FreeMem(temp_matr, m);
                }
        }
        return temp;
}
0
IT_Exp
Эксперт
34794 / 4073 / 2104
Регистрация: 17.06.2006
Сообщений: 32,602
Блог
13.11.2011, 20:20
Ответы с готовыми решениями:

Решить СЛАУ методом Крамера, методом Гаусса, с помощью обратной матрицы
Решите систему линейных уравнений а) методом Крамера; б) методом Гаусса; в) с помощью обратной матрицы:

Нахождение обратной матрицы методом LU-разложения
привет всем вот не могу найти теорию по этому вопросу может у кого-то найдется пример или алгоритм нахождения был бы очень благодарен

Нахождение обратной матрицы методом Гаусса
Помогите плиз! ни как не могу понять где ошибка: при вводе матрицы 6 16 12 0 1 3 11 2 3 6 3 2 -20 -10 5 4 получается ...

3
27 / 27 / 13
Регистрация: 07.11.2010
Сообщений: 154
13.11.2011, 21:47
У тебя в коде задана матрица целочисленного типа, исправь int на double и будет тебе счастье!
0
0 / 0 / 0
Регистрация: 13.11.2011
Сообщений: 10
13.11.2011, 23:46  [ТС]
Не получается

Добавлено через 50 минут
Запутался с типами уже...Везде конфликты вида C2664
0
0 / 0 / 0
Регистрация: 13.11.2011
Сообщений: 10
17.11.2011, 16:00  [ТС]
Так поможет кто-то ???
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
BasicMan
Эксперт
29316 / 5623 / 2384
Регистрация: 17.02.2009
Сообщений: 30,364
Блог
17.11.2011, 16:00
Помогаю со студенческими работами здесь

Нахождение обратной матрицы методом Гаусса
здравствуйте!! я чуть чуть плохо понимаю метод Гауса, хочу у вас спросить правильно ли работает программа вот сама прога, обратная...

Нахождение обратной матрицы методом Гаусса-Жордано
В данный момент проходим матрицы и у меня возникли сложности с нахождением обратной матрицы методом Гаусса-Жордано. Несколько раз уже...

Нахождение обратной матрицы методом алгебраических дополнений
#include &lt;stdio.h&gt; #include &lt;conio.h&gt; #include &lt;math.h&gt; #include &lt;stdlib.h&gt; #include &lt;time.h&gt; //int menu(); void...

Вычисление определителя и нахождение обратной матрицы методом исключения
Нужна только блок-схема, а то в схемах я пень пнем. Вот код. Function input_error() As Integer ' err – идентификатор...

Найти решение СЛАУ ( методом Крамера и обратной матрицей)
систему линейных уравнений формул Крамера и обратной матрицы: Правила форума, пункт 4.3. Создавайте темы с осмысленными и понятными...


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

Или воспользуйтесь поиском по форуму:
4
Ответ Создать тему
Новые блоги и статьи
модель ЗдравоСохранения 8. Подготовка к разному выполнению заданий
anaschu 08.04.2026
https:/ / github. com/ shumilovas/ med2. git main ветка * содержимое блока дэлэй из старой модели теперь внутри зайца новой модели 8ATzM_2aurI
Блокировка документа от изменений, если он открыт у другого пользователя
Maks 08.04.2026
Алгоритм из решения ниже реализован на примере нетипового документа, разработанного в конфигурации КА2. Задача: запретить редактирование документа, если он открыт у другого пользователя. / / . . .
Система безопасности+живучести для сервера-слоя интернета (сети). Двойная привязка.
Hrethgir 08.04.2026
Далее были размышления о системе безопасности. Сообщения с наклонным текстом - мои. А как нам будет можно проверить, что ссылка наша, а не подделана хулиганами, которая выбросит на другую ветку и. . .
Модель ЗдрввоСохранения 7: больше работников, больше ресурсов.
anaschu 08.04.2026
работников и заданий может быть сколько угодно, но настроено всё так, что используется пока что только 20% kYBz3eJf3jQ
Дальние перспективы сервера - слоя сети с космологическим дизайном интефейса карты и логики.
Hrethgir 07.04.2026
Дальнейшее ближайшее планирование вывело к размышлениям над дальними перспективами. И вот тут может быть даже будут нужны оценки специалистов, так как в дальних перспективах всё может очень сильно. . .
Горе от ума
kumehtar 07.04.2026
Эта мне ментальная установка, что вот прямо сейчас, мол, мне для полного счастья не хватает (нужное вписать), и когда я этого достигну - тогда и полный кайф. Одна из самых сильных ловушек на пути. . . .
Использование значений реквизитов справочника в документе, с определенными условиями и правами
Maks 07.04.2026
1. Контроль срока действия договора Алгоритм из решения ниже реализован на примере нетипового документа "ЗаявкаНаРаботу", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если. . .
Доступность команды формы по условию
Maks 07.04.2026
Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: сделать доступной кнопку (команда формы "ЗавершитьСписание") при. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru