Matlab - Expresar como condición que un número sea entero

Hola,

Estoy resolviendo un problema de optimización, cuyo resultado ha de ser un vector con elementos de la forma
v = x/N, donde x ha de ser necesariamente un valor de tipo integer positivo que pueda tomar como valores desde el 0 hasta N.

El problema es que no sé cómo expresar esta condición de forma matemática mediante Matlab. Estoy utilizando la función fmincon para resolver el problema de optimización:

https://es.mathworks.com/help/optim/ug/fmincon.html

Y no sé en qué parámetro 'encajar' esta condición adicional. Había pensado en hacerlo dentro de ceq, de la siguiente forma:


Pero no parece estar haciendo nada. ¿Alguien conoce alguna otra idea, por favor?

Si pudiera seros de utilidad, podéis encontrar el resto de mi código en StackOverflow en español:
http://es.stackoverflow.com/questions/56013/encontrar-la-distribuci%C3%B3n-de-1-en-filas-y-columnas-de-una-matriz-de-comprobaci

Muchas gracias de antemano por vuestra ayuda, y un cordial saludo. ¡Que paséis todos un buen fin de semana!

Cuando de optimización se trata con conjuntos discretos como lo son lo números enteros puedes usar algoritmos genéticos , con ellos es más fácil trabajar con números enteros, y es la solución más recomendada, si se te hace complicado trata de redondear tu respuesta al final

¡Hola de nuevo!

Sí, tienes mucha razón en tu observación. Creo que al final no me queda otro remedio y tendré que utilizar algoritmos genéticos... aunque el problema es que la última vez nunca llegó a terminar de converger (utilicé la función 'ga') y ahora he vuelto a intentarlo y parece ir por el mismo camino... ¿tienes alguna sugerencia?

La verdad es que mi código incluye muchas restricciones no lineales de desigualdad pero son completamente necesarias.

Mucahs gracias de nuevo por tu ayuda, y buen fin de semana a ti y al resto!


Aquí se incluyen aparte las restricciones no lineales:

te sugiero que plantees tu propio algoritmo genético donde puedas tener mayor control de tu función de evaluación y esos detalles, así estás haciendo un algoritmo genético muy general para tu problema , si me planteas el problema real en palabras ya sea por aquí o por interno puedo recomendarte un método de optimización más eficiente

Hola de nuevo!

Pues no sé si es posible, pero te puedo dejar la misma duda posteada en StackOverflow por aquí:

http://es.stackoverflow.com/questions/56013/encontrar-la-distribución-de-1-en-filas-y-columnas-de-una-matriz-de-comprobaci

Básicamente necesitaría encontrar un par de vectores h y v de longitudes N y k = N(1-r, 0<r<1) respectivamente que satisfagan las condiciones que aparecen en la imagen colgada. El plan es, más concretamente, hallar v mediante programación lineal como solución de un problema de optimización dado un cierto valor de h.

El plan es construir una matriz de comprobación de paridad H (básicamente una matriz de unos y ceros) con la siguiente información: sus dimensiones (k filas y N>k columnas) y la distribución de los pesos de sus filas (h) y sus columnas, (v). Esto es, si por ejemplo tenemos un valor h_2 = 3/20, quiere decir que 3/20 de todas las filas de H contienen dos '1' y el resto de sus elementos son '0'. Asimismo, si tuviéramos v_5 = 1/6 entonces 1/6 de las columas de H tienen cinco '1' y el resto de sus elementos son 0.

Cuando obtenga los valores de v y de h quiero construir la correspondiente matriz mediante un algoritmo que ya he finalizado. Para que funcione, es necesario que tanto los elementos de v como los de h sean de la forma x/N o y/k, siendo x e y numeros naturales (incluyendo el 0).

Creo que no me dejo nada... Pero si tienes alguna duda estoy a tu disposición, porque estoy desesperado por acabar con esto de una vez...

P.S: Este problema de optimización necesito resolverlo mediante programación lineal pero también necesitaré más adelante hacer lo propio mediante evolución diferencial bajo las mismas restricciones y minizando otra función que dependa tanto de v como de h (que pasarían a ocupar un único vector).

Muchas gracias por tu atención!! Eres muy amable!