Dada una matriz de nxm hallar su inversa
Publicado por Francisco (36 intervenciones) el 19/02/2018 17:26:58
Hola, ¿Cómo están?, ¿Este código está bien?, ¿O hay que corregirlo?, ¿Tiene alguna falla?, por favor, díganme.
Lo hice mediante desarrollo por cofactores.
Aquí va el código.
Lo hice mediante desarrollo por cofactores.
Aquí va el código.
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;
}Valora esta pregunta
0