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:
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.
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
0