por favor me pueden ayudar con este ejercicio
Publicado por isak (10 intervenciones) el 18/07/2020 02:54:25
realizar un código de programación lineal que resuelva ejercicios con el método simplex dual
aqui ya tenia algo pero tiene algunos errores
aqui ya tenia algo pero tiene algunos errores
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
#include <iostream>
#include <stdlib.h>
#include <stdio.h>
using namespace std;
double funcion[100], restricciones[100][100];
double restricc[100][100];
double independiente[100];
int main(){
int rest, var, ni ;
int i,j, k, columna, renglon=0, count, cuenta;
double maximo = 0, menor, aux, temp;
int c;
cout<<" Funcion a Optimizar \nNumero de Variables: ";
cin>>var;
cout<<"\nFUNCION OBJETIVO"<<endl;
for(c=0;c<var;c++){
cin>>funcion[c];
}
cout<<"\nNumero de Restricciones: ";
cin>>rest;
for(i=0;i<rest;i++){
cout<<"Restriccion "<<i+1<<endl;
for(c=0;c<var;c++){
cin>>restricc[i][c];
}
cout<<"Termino independiente ";
cin>>independiente[i];
}
for( i=0;i<var;i++)
for( j=0;j<rest;j++)
if(restricc[j][i]!=0)
restricciones[j][i]=restricc[j][i]*-1;
else
restricciones[j][i]=restricc[j][i];
//Crear tabla inicial simplex
for(i=0;i<rest;i++){
for(j=0; j<rest; j++)
if (i==j)
restricciones[i][var+j]=1;
else
restricciones[i][var+j]=0;
}
for(i=0; i<rest;i++)
restricciones[i][var+rest] = independiente[i]*-1;
for(j=0; j<var; j++)
restricciones[rest][j] = funcion[j]*-1;
for(j=var; j<rest+var; j++)
restricciones[rest][j] = 0;
ni=1;
cout<<"\n-----------------------------------------------";
cout<<"\n ITERACION "<<ni;
cout<<"\n-----------------------------------------------\n";
ni++;
//rest cantidad de restricciones para la tabla
for(i=0; i<=rest; i++) {
for(j=0; j<=rest+var; j++)
cout<<" \t"<<restricciones[i][j];
cout<<endl;
}
//Encontrar la variable de decision que entra a la base
//y la de holgura que saldra de la base
do{
///Encontrar la variable de decision que entrara a la base
maximo = abs(restricciones[rest][0]);
columna = 0;
count = 0;
aux = 120 ;
menor =120;
for(i=0; i<rest; i++) {
if(restricciones[i][columna] < 0)
aux = restricciones[i][rest+var] / restricciones[i][columna];
else
count++;
if( menor < aux ) {
menor = aux;
renglon = i;
}
else
continue;
}
for(j=0; j<=rest+var; j++)
if( restricciones[rest][j] < 0 ){
temp = -1*restricciones[rest][j];
if ( maximo < temp){
maximo = temp;
columna = j;
}
else
continue;
}
for(i=0; i<rest; i++) {
if(restricciones[i][columna] < 0)
aux = restricciones[i][rest+var] / restricciones[i][columna];
else
count++;
if( menor > aux ) {
menor = aux;
renglon = i;
}
}
if(count==var) {
cout<<"Solucion no acotada";
exit(1);
}
cout<<"\nPivote:"<<restricciones[renglon][columna];
//Encontrar los coeficientes de la nueva tabla
aux = restricciones[renglon][columna];
for(j=0; j<=(rest+var); j++)
restricciones[renglon][j] = restricciones[renglon][j] / aux;
for(i=0; i<=rest; i++) {
if (i == renglon ){
continue;}
else {
aux = restricciones[i][columna]*-1;
for(k=0; k <= rest+var; k++) {
restricciones[i][k] = (aux *restricciones[renglon][k]) + restricciones[i][k];
}
}
}
cout<<"\n-----------------------------------------------";
cout<<"\n ITERACION "<<ni;
cout<<"\n-----------------------------------------------\n";
ni++;
///la nueva tabla
for(i=0; i<=rest; i++) {
for(j=0; j<=rest+var; j++)
cout<<" \t"<<restricciones[i][j];
cout<<"\n";
}
cuenta=0;
for(j=0;j<rest+var;j++)
if(restricciones[var][j] >= 0)
cuenta++;
if (cuenta == rest+var)
break;
}while(cuenta);
return 0;
}
Valora esta pregunta
0