VPF::[С] Алгебраические дополнения - Форум программистов

IronM

Дата 14.2.2010, 16:15 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 1
Регистрация: 14.2.2010

Репутация: нет
Всего: нет

Здравствуйте, помогите пожалуйста с алгоритмом нахождения алгебраических дополнений матрицы. Матрица максимальным размером 5x5.
Пример:
дана матрица A:
1 2 3
2 1 3
3 2 1
Необходимо получить матрицу B в которой будут алгебраические дополнения:
-5 7 1
4 -8 4
3 3 -3

И если не трудно, помогите еще с поиском определителя матрицы, но только чтоб без указателей (их нам еще нельзя использовать), а во всех примерах работа идет через указатели и переделать у меня чего-то не получается.

Luyan

Дата 14.2.2010, 20:03 (ссылка)

(нет голосов)

Загрузка ...

Бывалый

Профиль
Группа: Участник
Сообщений: 180
Регистрация: 3.12.2008

Репутация: 5
Всего: 7

Сначала нахожу детерминант матрицы A, потом заполняю матрицу В алгебраическими дополнениями. Сделал как у вас в примере:

Код


#include <stdio.h>
#include <stdlib.h>

#define N 3

void FindAlgDop(int A[N][N], int size, int B[N][N])
{
    int i, j;
    
    // находим определитель матрицы A
    int det = FindDeterminant(A, size);
    
    if(det > 0) // это для знака алгебраического дополнения
        det = -1;
    else
        det = 1;

    int minor[size-1][size-1];
    
    for(j = 0; j < size; j++)
    {
        for(i = 0;i < size; i++)
        {
            // получаем алгебраическое дополнение
            GetMinor(A, minor, j, i, size);
            if( (i+j)%2 == 0)
                B[j][i] = -det*FindDeterminant(minor, size-1);
            else
                B[j][i] = det*FindDeterminant(minor, size-1);
        }
    }
}

int GetMinor(int A[N][N], int B[N][N], int x, int y, int size)
{
    int xCount = 0, yCount = 0;
    int i, j;
    for(i = 0; i < size; i++ )
    {
        if( i != x )
        {
            yCount = 0;
            for(j = 0; j < size; j++)
            {
                if( j != y )
                {
                    B[xCount][yCount] = A[i][j];
                    yCount++;
                }
            }
            xCount++;
        }
    }
    return 0;
}

int FindDeterminant(int A[N][N], int size)
{
    // останавливаем рекурсию, если матрица
    // состоит из одного элемента
    if(size == 1)
    {
        return A[0][0];
    }
    else
    {
        int det = 0;
        int i;
        int Minor[size-1][size-1];
        for(i = 0; i < size; i++ )
        {
            GetMinor( A, Minor, 0, i , size);
            // Рекурсия
            det += pow( -1, i ) * A[0][i] * FindDeterminant( Minor, size-1 );
        }
        return det;
    }
}

int main()
{
    int i, j;

    int A[N][N] = { {1, 2, 3},
                    {2, 1, 3},
                    {3, 2, 1} };
                        
    int B[N][N];
    
    printf("matr A:\n");
    for (i = 0; i < N; i++)
    {
        for (j = 0; j < N; j++)
            printf("%i ", A[i][j]);
        printf("\n");
    }
    
    // находим определитель матрицы matr1
    int determinant = FindDeterminant(A, N);
    printf("Determinant matr A: %i\n", determinant);
    
    
    // находим алгебраические дополнения
    // матрицы A и записываем их в B
    FindAlgDop(A, N, B);
    
    printf("matr B:\n");
    for (i = 0; i < N; i++)
    {
        for (j = 0; j < N; j++)
            printf("%i ", B[i][j]);
        printf("\n");
    }
    
    _getch();
    return 0;
}

1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »