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

Нахождение обратной матрицы (есть код! ) - C++

Восстановить пароль Регистрация
Другие темы раздела
C++ Для данной целочисленной матрицы A (aij) размером m*m найти матрицу B http://www.cyberforum.ru/cpp-beginners/thread893245.html
Будем называть соседями элемента с индексами i,j некоторой матрицы такие элементы, соответствующие индексы которых отличаются от i,j не более чем на единицу. Для данной целочисленной матрицы A (aij) размером m*m найти матрицу B, состоящую из нулей и единиц, элемент которой bij равен единице, когда: среди соседей aij есть не менее двух элементов, совпадающих с aij ...
C++ Симметричные матрицы Даны две разреженные структурно симметричные матрицы. Сложить их и результат занести в разреженную матрицу CSR. Есть код программы, считает всё нормально, но только в самом конце в массиве цифра обозначающая количество ненулевых элементов "+1" не сходиться с правильным ответом. Вот сам код: #include <cstdlib> #include <iostream> using namespace std; http://www.cyberforum.ru/cpp-beginners/thread893231.html
C++ Инвентировать биты в байте
Народ помогите пожалуйста как инвентировать биты в байте
C++ Найти положительные элементы
Нужно написать задачу в С++ с формами, надо очень срочно! В заданной целочисленной квадратной матрице Н найти в каждой строке число положительных элементов, и наибольший из положительных элементов строки. Исходную матрицу вывести в общепринятом виде. Результаты вывести в произвольной форме. Помогите пожалуйста!
C++ как найти максимальный поток в графе? http://www.cyberforum.ru/cpp-beginners/thread893213.html
как найти максимальный поток в графе?
C++ Создатть алгоритм под листинг #include <StdAfx.h> #include <iostream> #include <iomanip> using namespace std; void main() { setlocale(LC_CTYPE,"Russian"); const int N = 8; int matrix; подробнее

Показать сообщение отдельно
FireProoF
17 / 17 / 1
Регистрация: 28.10.2012
Сообщений: 190
06.06.2013, 14:11     Нахождение обратной матрицы (есть код! )
Добрый день! Имеется алгоритм нахождения обратной матрицы. Необходимо добавить в него регуляризацию, т.е. этот алгоритм может не найти обратную матрицу у несингулярной. Сразу скажу, я сам мало что понимаю)
Подскажите или объясните пожалуйста
Кликните здесь для просмотра всего текста
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
#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <math.h>
#include <conio.h>
#include <fstream>
#include <stdio.h>
#include "time.h"
#include <stdlib.h>
#include <ctime>
#include "time.h"
 
using namespace std;
 
 
 
void inversion(double **A, int N)
{
    double temp;
 
    double **E = new double *[N];
 
    for (int i = 0; i < N; i++)
        E[i] = new double [N];
 
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
        {
            E[i][j] = 0.0;
 
            if (i == j)
                E[i][j] = 1.0;
        }
 
    for (int k = 0; k < N; k++)
    {
        temp = A[k][k];
 
        for (int j = 0; j < N; j++)
        {
            A[k][j] /= temp;
            E[k][j] /= temp;
        }
 
        for (int i = k + 1; i < N; i++)
        {
            temp = A[i][k];
 
            for (int j = 0; j < N; j++)
            {
                A[i][j] -= A[k][j] * temp;
                E[i][j] -= E[k][j] * temp;
            }
        }
    }
 
    for (int k = N - 1; k > 0; k--)
    {
        for (int i = k - 1; i >= 0; i--)
        {
            temp = A[i][k];
 
            for (int j = 0; j < N; j++)
            {
                A[i][j] -= A[k][j] * temp;
                E[i][j] -= E[k][j] * temp;
            }
        }
    }
 
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
            A[i][j] = E[i][j];
 
    for (int i = 0; i < N; i++)
        delete [] E[i];
 
    delete [] E;
}
 
int main()
{
    int N;
 
    std::cout << "Enter N: ";
    std::cin >> N;
 
    double **matrix = new double *[N];
 
    for (int i = 0; i < N; i++)
        matrix[i] = new double [N];
 
    for (int i = 0; i < N; i++)
        for (int j = 0; j < N; j++)
        {
            std::cout << "Enter matrix[" << i << "][" << j << "] = ";
            std::cin >> matrix[i][j];
        }
 
    inversion(matrix, N);
 
    for (int i = 0; i < N; i++)
    {
        for (int j = 0; j < N; j++)
            std::cout << matrix[i][j] << "  ";
 
        std::cout << std::endl;
    }
 
    for (int i = 0; i < N; i++)
        delete [] matrix[i];
 
    delete [] matrix;
 
    std::cin.get();
    
    _getch();
}
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
 
Текущее время: 15:24. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru