martes, 20 de marzo de 2012

Python

Python
Es un lenguaje de programación de alto nivel. Python fue creado a finales de los ochenta por  Guido Van Rossum en el Centro para las Matemáticas y la Informática (CWI, Centrum Wiskunde & Informatica), en los Países Bajos, como un sucesor del lenguaje de programación ABC. 


¿Por qué usar Python?
Por que es facil de usar,  es de código abierto

Ejemplos de python
 Un ejemplo sencillo de como generar números aleatorios y guardarlos en  un archivo. También un ejemplo de como trabajar con importat gnuplot en python.





#!/usr/bin/python
import random
import Gnuplot

f = open("archivo.txt","w")
for i in range(0, 10):
  a=random.uniform(1,1000)
  f.write("%fl\n" %a)
  print(a)
f.close()

b=[[1,2],[3,4],[5,6],]
gp = Gnuplot.Gnuplot(persist=1)
gp.title('Grafica sencilla ')
gp.xlabel('primera entrada')
gp.ylabel('segunda entrada')
gp.plot(b)

 Un ejemplo de como ingresar datos desde teclado  y como poner mensajes de error







valor=input('Dame un numero: ')
try:
    temp=int(valor)
    print 'el valor es: ', valor
except ValueError as variable:
    print 'Debes ingresar un numero entero'



Python y Cplex 
Es posible trabajar con Cplex desde pyhton  en el enlace 
se puede encontrar todo la información y de como trabajar con python y Cplex.

A continuación se presenta un demo

#!/usr/bin/python
# diet.py -- A diet problem
# Descripción del problema
# Minimizar el costo de la dieta sujeto a las restricciones de los nutrientes.
# To run from within the python interpreter, use
# >>> import diet


import cplex
from cplex.exceptions import CplexError
from inputdata import read_dat_file

#  Clase para los datos 
class ProbData:

    def __init__(self, filename):
        
        # lee los datos del archivo diet.dat
        self.foodCost, self.foodMin, self.foodMax, self.nutrMin, \
                       self.nutrMax, self.nutrPer = \
                       read_dat_file("../../data/diet.dat")

        # revisa si los datos no son inconsistentes
        if len(self.foodCost) != len(self.foodMin) or \
               len(self.foodCost) != len(self.foodMax) or \
               len(self.nutrMin)  != len(self.nutrMax) or \
               len(self.nutrMin)  != len(self.nutrPer):
            print "ERROR: Data file '%s' contains inconsistent data\n" % filename
            raise Exception("data file error")

        for np in self.nutrPer:
            if len(self.foodCost) != len(np):
                print "ERROR: Data file '%s' contains inconsistent data\n" % filename
                raise Exception("data file error")


def populatebyrow(prob, data):

    nFoods = len(data.foodCost)
    nNutrients = len(data.nutrMin)

    # we want to minimize costs
    prob.objective.set_sense(prob.objective.sense.minimize)

    # add variables to decide how much of each type of food to buy
    varnames = ["x"+str(j) for j in range(nFoods)]
    prob.variables.add(obj = data.foodCost,
                       lb = data.foodMin,
                       ub = data.foodMax,
                       names = varnames)

    # add constraints to specify limits for each of the nutrients
    for n in range(nNutrients):
        prob.linear_constraints.add(lin_expr = [[varnames,data.nutrPer[n]]],
                                    senses = ["R"],
                                    rhs = [data.nutrMin[n]],
                                    range_values = [data.nutrMax[n] - data.nutrMin[n]])
    
    
def populatebycolumn(prob, data):
    nFoods = len(data.foodCost)
    nNutrients = len(data.nutrMin)

    # we want to minimize costs
    prob.objective.set_sense(prob.objective.sense.minimize)

    # create empty constraints to be filled later
    rownames = ["r"+str(n) for n in range(nNutrients)]
    prob.linear_constraints.add(senses = ["R" * nNutrients],
                                rhs = data.nutrMin,
                                range_values = [data.nutrMax[n] - data.nutrMin[n]
                                                for n in range(nNutrients)],
                                names = rownames)

    # create columns
    for j in range(nFoods):
        prob.variables.add(obj = [data.foodCost[j]],
                           lb = [data.foodMin[j]],
                           ub = [data.foodMax[j]],
                           columns = [[rownames, [data.nutrPer[n][j]
                                                  for n in range(nNutrients)]]])

    
def diet(pop_method):
    try:
        # read the data in diet.dat
        data = ProbData("../../data/diet.dat")

        # create CPLEX object
        my_prob = cplex.Cplex()

        # populate problem
        if pop_method == "r":
            handle = populatebyrow(my_prob, data)
        if pop_method == "c":
            handle = populatebycolumn(my_prob, data)

        # solve problem
        my_prob.solve()
        
    except CplexError, exc:
        print exc
        return

    numrows = my_prob.linear_constraints.get_num()
    numcols = my_prob.variables.get_num()

    solution = my_prob.solution
    
    # solution.get_status() returns an integer code
    print "Solution status = " , solution.get_status(), ":",
    # the following line prints the corresponding string
    print solution.status[solution.get_status()]
    print "Objective value = " , solution.get_objective_value()
    
    x = solution.get_values(0, my_prob.variables.get_num()-1)
    for j in range(my_prob.variables.get_num()):
        print "Buy %d = %17.10g" % (j, x[j])

import sys

if __name__ == "__main__":
    if len(sys.argv) != 2 or sys.argv[1] not in  ["-r", "-c"]:
        print "Usage: diet.py -X"
        print "   where X is one of the following options:"
        print "      r          generate problem by row"
        print "      c          generate problem by column"
        print " Exiting..."
        sys.exit(-1)
    diet(sys.argv[1][1])
else:
    prompt = """Enter the letter indicating how the problem data should be populated:
    r : populate by rows
    c : populate by columns\n ? > """
    r = 'r'
    c = 'c'
    diet(input(prompt))
Referencias Introducción 
modelo_simple

No hay comentarios:

Publicar un comentario