Dev - C++ - AYUDA MÉTODO DE GAUSS SIMPLE,

 
Vista:
sin imagen de perfil
Val: 22
Ha disminuido su posición en 2 puestos en Dev - C++ (en relación al último mes)
Gráfica de Dev - C++

AYUDA MÉTODO DE GAUSS SIMPLE,

Publicado por Bryan Jeronimo (13 intervenciones) el 02/03/2020 06:34:43
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
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
////LO QUE PASA ES QUE SACA BIEN LAS SOLUCIONES EL PROGRAMA Y LOS CEROS DEBAJO DE LA DIAGONAL, PERO QUISIERA PEDIR AYUDA PARA QUE EN LA MATRIZ FINAL IMPRIMA LOS 1 EN DIAGONAL OBVIAMENTE CON LAS SOLUCIONES CORRECTAS CON EL METODO DE GAUSS SIMPLE, OJO EN GAUSS-JORDAN ES DISTINTO EL MÉTODO, EN GAUSS SIMPLE TIENE QUE SALIR EL ESCALONADO DE CEROS SOLO ABAJO DE LA DIAGONAL DE LOS UNOS (EN GAUSS JORDAN ES TANTO ARRIBA COMO ABAJO) , EL CÓDIGO:
 
 
#include <iostream>
#include <cmath>
#include <iomanip>
#include <graphics.h>
#include <conio.h>
#include <windows.h>
#include <math.h>
#include <stdio.h>
#include <fstream>
#include <string>
#include <ctime>
#include <time.h>
#include <vector>
 
 
using namespace std;
 
void leerDatos( double **M, int fil, int col, double B[] );
void gaussEliminacion( double **M, int fil, int col, double v[] );
void intercambiarFilas( double **M, int col, int k, int max );
void imprimirMatriz( double **M, int fil, int col, double v[] );
void imprimirSolucion( double v[], int fil );
 
int main()
 
{
	char repetir= 'S';
    int opcion;
    while (repetir != 'N')
	{
	  setlocale(LC_CTYPE, "Spanish");
	  int fil, col;
    //gotoxy (30,1);
	cout<<"METODO DE GAUSS SIMPLE"<<endl;
	//gotoxy (1,2);
	cout<<"Este programa muestra las soluciones de sistemas de ecuaciones"<<endl;
	//gotoxy (1,3);
	cout<<"por el método de Gauss Simple"<<endl;
	//gotoxy (1,4);
    cout << "\Digite el número de variables a encontrar del sistema de ecuaciones: ";
    cin >> fil;
    col = fil;
 
    double** A = new double*[fil];
 
    for( int i = 0; i < fil; i++ )
        A[i] = new double[col];
        double B[fil];
 
    leerDatos( A, fil, col, B );
    imprimirMatriz( A, fil, col, B );
    gaussEliminacion( A, fil, col, B );
 
	cout<<endl;
    cout<<"1) Otro calculo  2) Salir" <<endl;
    cin>>opcion;
 
    switch (opcion)
	{
    	case 1:
    		cout<<"..."<<endl;
    		break;
    	case 2:
    		repetir= 'N';
    		break;
    	default:
    		cout<<"Opcion invalida"<<endl;
    		break;
	}
	if (repetir != 'N')
	{
		repetir = ' ';
	}
	while (repetir != 'S' && repetir != 'N')
	{
			cout<<" Ejecutar de nuevo el programa: S/N: " ;
			cin>>repetir;
			repetir=toupper(repetir);
			system ("cls");
	}
    }
 
 
}
 
 
 
void leerDatos( double **A, int fil, int col, double B[] )
 
{
 
    cout << "\nDigite coeficientes o elementos que conformarán a la matriz:\n";
 
    for( int i = 0; i < fil; i++ ){
 
        cout << "\nFILA " << i + 1 << "\n";
 
        for( int j = 0; j < col; j++ ) {
 
            cout << "Columna " << j + 1 << ": ";
 
            cin >> A[i][j];
 
        }
 
    }
 
    cout << "\nIntroduzca cada valor de la igualdad de la ecuación:\n";
 
    for(int i = 0; i < fil; i++ ){
 
        cout << "\nFila " << i + 1 << ": ";
 
        cin >> B[i];
 
    }
 
 
 
}
 
 
 
void gaussEliminacion( double **A, int fil, int col, double B[] )
 
{
    int max;
    int N = fil;
    double factor, temp;
 
    for( int k = 0; k < N; k++ ){
        max = k;
 
        for( int i = k + 1; i < N; i++ )
            if( abs( A[i][k] ) > abs( A[max][k] ) )
                max = i;
 
        intercambiarFilas( A, col, k, max );
        temp = B[k];
        B[k] = B[max];
        B[max] = temp;
 
        for( int i = k + 1; i < N; i++ ) {
            factor = A[i][k] / A[k][k];
            B[i] -= factor * B[k];
 
            for( int j = k; j < N; j++ )
                A[i][j] -= factor * A[k][j];
        }
 
    }
 
    imprimirMatriz( A, fil, col, B );
 
    double solucion[fil] = { 0.0 };
    double suma;
 
    for (int i = N - 1; i >= 0; i-- ){
 
        suma = 0.0;
 
        for( int j = i + 1; j < N; j++ )
            suma += A[i][j] * solucion[j];
 
        solucion[i] = (B[i] - suma) / A[i][i];
    }
 
    imprimirSolucion( solucion, fil );
 
}
 
 
 
void intercambiarFilas( double **M, int col, int k, int max )
 
{
    double temp[col];
 
    for( int i = k; i < col; i++ ) {
            temp[i] = M[k][i];
            M[k][i] = M[max][i];
            M[max][i] = temp[i];
    }
 
}
 
 
 
void imprimirMatriz( double **M, int fil, int col, double v[] )
 
{
 
    cout << setprecision(3) << fixed;
 
    cout << "\nMatriz digitada/final:\n";
 
    for( int i = 0; i < fil; i++){
 
        for( int j = 0; j < col; j++){
 
            cout << M[i][j] << "\t";
 
        }
 
        cout << "| " << v[i] << "\n";
 
    }
 
    cout << "\n";
 
}
 
 
 
void imprimirSolucion( double v[], int fil )
 
{
 
    cout << "\nSolucion: ";
 
    for( int i = 0; i < fil; i++ )
 
        cout << "\nX" << i + 1 << ": " << v[i];
 
 
 
    cout << endl;
}



AYUDA
ayuda2
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
Imágen de perfil de Alfil
Val: 4.344
Oro
Ha mantenido su posición en Dev - C++ (en relación al último mes)
Gráfica de Dev - C++

AYUDA MÉTODO DE GAUSS SIMPLE,

Publicado por Alfil (1444 intervenciones) el 02/03/2020 18:00:15
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
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
#include <iostream>
#include <iomanip>
 
using namespace std;
 
void LeerMatriz( float **M, int fil, int col );
void ImprimirFase( float **M, int fil, int col );
void Gauss( float **M, int fil, int col );
void ImprimirResultado( float **M, int fil, int col );
 
int main()
{
    int col, fil;
 
    cout << "# filas de la matriz: ";
    cin >> fil;
    cout << "# columnas de la matriz: ";
    cin >> col;
 
    float** M = new float*[fil];
    for( int i = 0; i < fil; i++ )
        M[i] = new float[col];
 
    LeerMatriz( M, fil, col );
 
    cout << "\nMatriz inicial:";
    ImprimirFase( M, fil, col );
 
    Gauss( M, fil, col );
    ImprimirResultado( M, fil, col );
 
    return 0;
}
 
void LeerMatriz( float **M, int fil, int col )
{
    cout << "\nDatos de la matriz:\n";
    for( int i = 0; i < fil; i++ ){
        cout << "Fila: " << i + 1 << "\n";
        for( int j = 0; j < col; j++ ){
            cout << "Columna " << j + 1 << ": ";
            cin >> M[i][j];
        }
    }
}
 
void ImprimirFase( float **M, int fil, int col )
{
    cout << "\n\n";
    for( int i = 0; i < fil; i++ ){
        for( int j = 0; j < col; j++ ){
            cout << setprecision(2) << fixed << setw(6)
                 << M[i][j] << " ";
        }
        cout << "\n";
    }
}
 
void Gauss( float **M, int fil, int col )
{
    int k = 0;
    float temp;
 
    while( k < col - 1 ){
        for( int i = k; i < fil; i++ ){
            for( int j = col - 1; j >= k; j-- ){
                if( i == k ){
                    temp = 1 / M[k][k];
                    M[i][j] = temp * M[i][j];
                }
                if( i > k ){
                    temp = (-M[i][k]);
                    M[i][j] = ( temp * (M[k][j])) + M[i][j];
                }
            }
        }
        k++;
    }
    cout << "\nFase 1:";
    ImprimirFase( M, fil, col);
 
    for( int i = 0; i < fil - 1; i++ ){
        for( int j = col - 1; j >= col - 2; j-- ){
            temp = (-M[i][col - 2]);
            M[i][j] = (temp * M[fil - 1][j])+ M[i][j];
        }
    }
    cout << "\nFase 2:";
    ImprimirFase( M, fil, col);
 
    k = 0;
    int c = 1;
    while( k < col - 1 ){
        for( int i = k; i < fil - 2; i++ ){
            for( int j = col - 1; j >= k + c; j-- ){
                temp = (-M[i][k + c]);
                M[i][j] = temp * M[i + 1][j] + M[i][j];
            }
        }
        k++;
        c++;
    }
    cout << "\nFase 3:";
    ImprimirFase( M, fil, col );
 
}
 
void ImprimirResultado( float **M, int fil, int col )
{
    cout << "\nImprimir resultado:\n\n";
    for( int i = 0; i < fil; i++ ){
        cout << setprecision(0) << fixed << setw(3)
             << char(120 + i) << " = "
             << setw(3) << setiosflags( ios::right )
             << M[i][col - 1] << endl;
    }
}
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
2
Comentar