Dev - C++ - Matriz inversa de una matriz nxm

 
Vista:
Imágen de perfil de Francisco
Val: 49
Ha disminuido 1 puesto en Dev - C++ (en relación al último mes)
Gráfica de Dev - C++

Matriz inversa de una matriz nxm

Publicado por Francisco (36 intervenciones) el 12/02/2018 06:05:30
Hola, les mando código de matriz inversa y fuente de como calcularla.
La fuente es de la matriz inversa, de la matriz cofactor y la adjunta que la adjunta es la transpuesta del cofactor.
El codigo tiene errores en tiempo de pantalla, ya que muestra todos los elementos de la matriz inversa como ceros.
Y en el codigo no puedo hacer (1/det3)* ... ya que det3 es de tipo int, salvo que lo pase como float.
Les pido que me digan donde esta la falla, y si se puede mandar que sea sin hacer void, sin usar librerias de arreglos ni usar punteros, o sea todo dentro del main.


Aquí va el codigo:

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
#include<iostream>
using namespace std;
int main(){
    int n=0,m=0;
    cout<<"ingrese el tamanio de la matriz "<<endl;
    cin >>n>> m;
    int mat[n][m],cof[n][m],adj[n][m],cofaux[n][m],inv[n][m],aux_inv[n][m];
    if(n==m && n>0 && m>0){
        cout<<"es una matriz cuadrada "<<endl;
    }
    else{
        cout<<"no se puede calcular la matriz. tiene que ser cuadrada "<<endl;
    }
    int i=0,j=0,det=0,det2=0,det3=0,det4=0,aux_cof=1,d=0,k=0,e=0;
    float det5;
    cout<<"ingrese valores de la matriz "<<endl;
    for(i=0;i<n;i++){
        for(j=0;j<m;j++){
            cin >> mat[i][j];
        }
    }
    //calculamos la matriz cofactor
    for(i=0;i<n;i++){
        for(j=0;j<m;j++){
            for(d=0;d<n;d++){
            for(e=0;e<m;e++){
            for(k=0;k<m;k++){
                    if((i+k)%2==0){
                        aux_cof = 1;
                       }
                       else if((i+k)%2!=0){
                        aux_cof= -1;
                       }
                    if((d==k || d!=k) && (e==k || e==d)){
                        cofaux[i][j]= aux_cof*((mat[i][k]*mat[d][j])-(mat[i][k]*mat[e][j]));
                    }
            }
        }
    }
    }
    }
    for(i=0;i<n;i++){
        for(j=0;j<m;j++){
            cof[i][j]= cofaux[i][j];
        }
    }
    for(i=0;i<n;i++){
        for(j=0;j<m;j++){
                adj[i][j] = cof[j][i];
 
        }
    }
    //buscamos determinante
    for(i=0;i<n;i++){
        for(j=0;j<m;j++){
            for(k=0;k<n;k++){
                for(e=0;e<n;e++){
                    det=0;
                    det2=0;
                    for(d=0;d<m;d++){
                        if((d!=k || d==k) && (e==k || e==d)&& i!=j){
                            if(i%2==0 & j%2==0 ){
                                det= det + ((mat[i][d]*mat[k][j])-(mat[i][e]*mat[d][j]));
                            }
                            else if(i%2!=0 && j%2!=0 ){
                                det2= det-((mat[i][d]*mat[k][j])-(mat[i][d]*mat[e][j]));
                            }
                        }
                        if(i==j){
                                det3= det3 +(mat[i][j]*det);
                                det4= det4 +(mat[i][j]*det2);
                        }
                    }
                }
            }
        }
    }
    det5 = det3 + det4;
    if(det5==0){
        cout<<"matriz no tiene inversa"<<endl;
    }
    else{
        for(i=0;i<n;i++){
            for(j=0;j<m;j++){
                aux_inv[i][j]= (1/det5)*adj[i][j];
                inv[i][j] = aux_inv[i][j];
            }
        }
    }
    cout<<"elementos de la matriz inversa "<<endl;
    for(i=0;i<n;i++){
        for(j=0;j<m;j++){
            cout<< inv[i][j]<<" ";
        }
        cout << endl;
    }
    return 0;
}

Aqui van las fuentes:
http://mathworld.wolfram.com/images/equations/MatrixInverse/NumberedEquation4.gif
https://es.wikihow.com/invertir-una-matriz-de-3X3
http://www.academia.edu/19601954/Matriz_de_cofactores
https://es.wikipedia.org/wiki/Matriz_de_adjuntos

¿Cuáles son los errores?, sobre todo cuando hago el desarrollo por cofactores. Les agradezco.
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder

Matriz inversa de una matriz nxm

Publicado por xavi (1 intervención) el 21/09/2019 10:03:47
que pasa en la linea 34?
if((d==k || d!=k) && (e==k || e==d))
especialmente con (d==k || d!=k)
eestas dando todos los casos, solo hay dos, o es igual o no es y tu dices que si cumple uno de los dos, pues para que lo pones? o en que te equivocaste?
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar