обратная матрица

[Vitek.i]

Да, это басик ))). Просто сам допендрил и решил сдесь написать. Вдруг когда пригодится кому нить))

[vitaska]

Может можно и проще, но я решил таким образом, так что не обессудьте:

#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;
int **matr = new int * [n];
for(int i = 0; i < n; i++)
matr[i] = new int[n];
SetMtx(matr, n);
PrintMtx(matr, n);
det = Det(matr, n);
cout << "Определитель матрицы = " << det << endl;

if(det){
double **obr_matr = new double * [n];
double **tobr_matr = new double * [n];
for(int i = 0; i < n; i++){
obr_matr[i] = new double[n];
tobr_matr[i] = new double[n];
}
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);
}
}
TransponMtx(obr_matr, tobr_matr, n);
PrintMtx(tobr_matr, n);
FreeMem(tobr_matr, n);
FreeMem(obr_matr, n);	
}
else cout << "Т.к. определитель матрицы = 0,\nто матрица вырожденная и обратной не имеет!!!" << endl;
FreeMem(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)
{
if(&matr != NULL){
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++)
matr[i][j] = rand()%9 + 1; 
}
Функция печати матрицы
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)
         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;
}