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

Обратная мартрица - C++

Восстановить пароль Регистрация
 
lenalena
39 / 2 / 1
Регистрация: 04.06.2012
Сообщений: 19
08.11.2012, 16:08     Обратная мартрица #1
Ребят, пожалуйста помогите исправить код. Препод сказал что в конце неправильно!!!

Для заданной матрицы A(m,n) найти обратную http://www.cyberforum.ru/cgi-bin/latex.cgi?{A}^{-1} , используя интерационную формулу: http://www.cyberforum.ru/cgi-bin/latex.cgi?{A}^{-1}_{k}={A}^{-1}_{k-1}*(2E-A*{A}^{-1}_{k-1}) , где E -единичная матрица; http://www.cyberforum.ru/cgi-bin/latex.cgi?{A}^{-1}_{0}=E. Интерационный процесс заканчивается, если для заданной погрешности http://www.cyberforum.ru/cgi-bin/latex.cgi?\varepsilon справедливо |det(A*http://www.cyberforum.ru/cgi-bin/latex.cgi?{A}^{-1}_{k})-1<=E[/QUOTE]

#include <iostream.h>
#include <conio.h>
#include <math.h>
int main()
{
const int a=3, b=3;
int arr[a][b];
float arr2[a][b];
int x=0, y=0;
float det, tmp;
cout<<"\nMatrix:\n";
for (x=0; a>x; x++)
{
for (y=0; b>y; y++)
{
cout<<"\n"<<x+1<<":"<<y+1<<" ";
cin>>arr[x][y];
}
}
cout<<"\n";
for (x=0; a>x; x++)
{
for (y=0; b>y; y++)
{
cout<<"\t"<<arr[x][y];
}
cout<<"\n";
}
cout<<"\n";
det=(arr[0][0]*arr[1][1]*arr[2][2])+(arr[0][1]*arr[1][2]*arr[2][0])+(arr[0][2]*arr[1][0]*arr[2][1])-(arr[0][2]*arr[1][1]*arr[2][0])-(arr[0][1]*arr[1][0]*arr[2][2])-(arr[0][0]*arr[1][2]*arr[2][1]);
arr2[0][0]=((arr[1][1]*arr[2][2])-(arr[1][2]*arr[2][1]))/det;
arr2[0][1]=((arr[1][0]*arr[2][2])-(arr[1][2]*arr[2][0]))/det;
arr2[0][2]=((arr[1][0]*arr[2][1])-(arr[1][1]*arr[2][0]))/det;
arr2[1][0]=((arr[0][1]*arr[2][2])-(arr[0][2]*arr[2][1]))/det;
arr2[1][1]=((arr[0][0]*arr[2][2])-(arr[0][2]*arr[2][0]))/det;
arr2[1][2]=((arr[0][0]*arr[2][1])-(arr[0][1]*arr[2][0]))/det;
arr2[2][0]=((arr[0][1]*arr[1][2])-(arr[0][2]*arr[1][1]))/det;
arr2[2][1]=((arr[0][0]*arr[1][2])-(arr[0][2]*arr[1][0]))/det;
arr2[2][2]=((arr[0][0]*arr[1][1])-(arr[0][1]*arr[1][0]))/det;
от сюда уже что то неправильно((
for (x = 1; a > x; x++)
for ( y = 0; y < x; y++)
{
tmp = arr2[x][y];
arr2[x][y] = arr2[y][x];
arr2[y][x] = tmp;
}
for (y=0; b>y; y++)
{
for (x=0; a>x; x++)
{
cout<<"\t"<<arr2[y][x];
}
cout<<"\n";
}
return 0;
}
Помогите пожалуйста, буду очень признательна))
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
08.11.2012, 16:08     Обратная мартрица
Посмотрите здесь:

C++ Обратная матрица на С
C++ Обратная матрица
C++ Обратная матрица
C++ Обратная формула
C++ Обратная матрица
Обратная матрица C++
C++ Обобщенная обратная матрица
C++ Обратная матрица

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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Simply me
 Аватар для Simply me
235 / 31 / 2
Регистрация: 05.05.2012
Сообщений: 292
08.11.2012, 16:14     Обратная мартрица #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
void inversion(float A[4][4], int N)         
{                                                
    float temp;
    float B[4][4];
    float res[4][4];
    int i,j,k;
 
    float **E = new float *[N];
 
    for (int i = 0; i < N; i++)
        E[i] = new float [N];
        for (i=0; i<N; i++)
        for (j=0; j<N; j++)
        B[i][j]=A[i][j];
 
    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;
}
В качестве параметров передаешь матрицу и ее размерность.
lenalena
39 / 2 / 1
Регистрация: 04.06.2012
Сообщений: 19
08.11.2012, 16:54  [ТС]     Обратная мартрица #3
Спасибо, но только мне нужно чтобы был такой же код как у меня, но с исправленной концовкой((

Добавлено через 37 минут
Ребят, ну помогите мне пожалуйста(( завтра сдать нужно
Yandex
Объявления
08.11.2012, 16:54     Обратная мартрица
Ответ Создать тему
Опции темы

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