Как привести матрицу к унитреульной матрице или верхней треугольной матрице?

@Naatikin · Регистрация: 01.11.2010

Как привести матрицу к унитреульной матрице или верхней треугольной матрице?

Добавлено через 14 часов 4 минуты

C

#include <stdio.h>
#include <conio.h>
#include <math.h>
#include <iostream>
#define N 5 
void main (){
int mt, nt, i,j,k,c;
int A[N+3][N]={
1,1,0,0,0,
1,1,0,1,0,
1,0,1,0,0,
1,0,1,1,1,
1,0,0,0,0,
0,0,1,0,0,
0,1,0,0,1,
0,1,0,1,0}, B[N+3]={0,1,2,3,4,5,6,7};
mt=8; nt=5;
for(int u=0; u<mt; u++){
    for(int y=0; y<nt; y++){
        printf("%d\t",A[u][y]);
    if(y==(nt-1))printf("\n");}
}
printf("\n");
 
 
k=0;
while(k<mt){
j=k;
if(A[B[j]][k]==0 && j<mt){j++;}
if(j>(mt-1)){k++;}
if(j>(k-1)){c=B[j];B[j]=B[k];B[k]=c;}
i=k;
while(i<(mt-1)){
    i++;
if(A[B[i]][k]==1){
    for(j=0;j<nt;j++){
        A[B[i]][j]=(A[B[i]][j]+A[B[k]][j])%2;}
}
}
for(int u=0; u<N+3; u++){
    for(int y=0; y<N; y++){
        printf("%d\t",A[u][y]);
    if(y==(N-1))printf("\n");}
}
printf("\n");
k++;
}
}

написал код, но работает не так. может кто посмотрит

@A555 · 09.12.2011, 13:12

а вот лови прямо здесь на форуме лежит я правдо немного отредактировал что бы сдать преподу мог вот функция

C++

#define FOR_k  for (int k = 0; k < dim; ++k)
#define FOR_j  for (int j = 0; j < dim; ++j)
#define FOR_i  for (int i = k + 1; i < dim; ++i)
void swap(double &a, double &b)
{
    double temp = a;
    a = b;
    b = temp;                         
 
}
 
void Gauss(double **A, double *B, int dim)
{
    FOR_k
    {
        if (A[k][k] == 0.0)
        {
           FOR_i
            {
                if (A[i][k] != 0.0)
                {
                    FOR_j
                        swap(A[i][j], A[k][j]);
 
                    swap(B[i], B[k]);
 
                    break;
                }
            }
        }
 
        double div = A[k][k];
 
        FOR_j
            A[k][j] /= div;
 
        B[k] /= div;
 
       FOR_i
        {
            double multi = A[i][k];
 
            for (int j = 0; j < dim; ++j)
                A[i][j] -= multi * A[k][j];
 
            B[i] -= multi * B[k];
        }
    }
 
    for (int k = dim - 1; k > 0; --k)
    {
        for (int i = k - 1; i + 1 > 0; --i)
        {
            double multi = A[i][k];
 
            for (int j = 0; j < dim; ++j)
                A[i][j] -= multi * A[k][j];
 
            B[i] -= multi * B[k];
        }
    }                                                    
}

решение матриц методом Гауса

@Naatikin · 09.12.2011, 13:29 **[ТС]**

dim что за переменная? размерность?

@A555 · 09.12.2011, 13:37

Сообщение от Naatikin

dim что за переменная? размерность?

да размерность

@Naatikin · 09.12.2011, 13:51 **[ТС]**

у меня есть матрица из нулей и единиц, нужно получить унитреугольную матрицу. что тогда я должен в качестве матрицы B[] подавать?

@A555 · 09.12.2011, 13:53

C++

int A[N+3][N]={
1,1,0,0,0,
1,1,0,1,0,
1,0,1,0,0,
1,0,1,1,1,
1,0,0,0,0,
0,0,1,0,0,
0,1,0,0,1,
0,1,0,1,0}, B[N+3]={0,1,2,3,4,5,6,7};

я так понимаю это матрица и столбец свободных членов ?? так если так то программа сделает еденичную матрицу треугольную

@Naatikin · 09.12.2011, 13:55 **[ТС]**

нет мне нужно только совершить эквивалентные преобразования над матрицей А, матрица В тут не при чём

@A555 · 09.12.2011, 14:21

Сообщение от Naatikin

нет мне нужно только совершить эквивалентные преобразования над матрицей А, матрица В тут не при чём

матрица) Б у тебя там массив когда вижу такое сразу понимаю что это матрица со столбоцом свободных членов

@Naatikin 4 / 4 / 0 Регистрация: 01.11.2010 Сообщений: 97
	09.12.2011, 13:29 [ТС]	3
	dim что за переменная? размерность? 0

@A555 51 / 51 / 4 Регистрация: 04.04.2011 Сообщений: 209
	09.12.2011, 13:37	4
	Сообщение от Naatikin dim что за переменная? размерность? да размерность 0

@Naatikin 4 / 4 / 0 Регистрация: 01.11.2010 Сообщений: 97
	09.12.2011, 13:51 [ТС]	5
	у меня есть матрица из нулей и единиц, нужно получить унитреугольную матрицу. что тогда я должен в качестве матрицы B[] подавать? 0

@Naatikin 4 / 4 / 0 Регистрация: 01.11.2010 Сообщений: 97
	09.12.2011, 13:55 [ТС]	7
	нет мне нужно только совершить эквивалентные преобразования над матрицей А, матрица В тут не при чём 0

@A555 51 / 51 / 4 Регистрация: 04.04.2011 Сообщений: 209
	09.12.2011, 14:21	8
	Сообщение от Naatikin нет мне нужно только совершить эквивалентные преобразования над матрицей А, матрица В тут не при чём матрица) Б у тебя там массив когда вижу такое сразу понимаю что это матрица со столбоцом свободных членов 0

Опции темы

Как привести матрицу к унитреульной матрице или верхней треугольной матрице?

Решение