Cleber, 

Olhando de "soslaio" as suas condições de contorno vê-se que não se trata de um problema de LP. . .

Prá complicar suas condições de contorno têm um grau de liberdade (a variável x é livre[?])...


2016-05-14 9:27 GMT-03:00 Cleber N.Borges <klebyn@yahoo.com.br>:
Prezados,
Bom dia!
Como na lista tem pessoas com conhecimento e experiência, peço algumas orientações... :-)

Imagino que o método de  "Programação Linear (LP)"
deve servir para resolver meu problema de achar os parâmetros de uma dada função objetivo.

Sei que no R tem a opção nativa "constrOptim" (poderia ser outra, posso instalar pacotes)

Entretanto eu não sei como elaborar minhas condições de contorno em formato matricial
que é o modo como o comando constrOptim trabalha (até onde consegui entender)

Agradeceria muito se alguém puder me ajudar,
respondendo se eu realmente estou certo em pensar que a LP é a ferramenta certa para meu caso
e /(ou que) me indique um material básico (introdutório).

Abaixo segue um exemplo simplificado do tipos de função que pretendo obter parâmetros.

Agradeço antecipadamente pela ajuda e por qualquer dica

Cleber Borges

#####################
#
fobjetivo <- function( x = 0.1 ){
k1 <- 0.007585776
k2 <- 6.16595e-08
k3 <- 2.137962e-13
k4 <- 1e-14
k <- k1 * k2 * k3 * k4
# Determinar as 6 variáveis: a, b, c, d, e, f
###################
# 6 Condições de contorno: ( Igualdades )
###################
k1 == a * d / c            # 1
k2 == a * e / d            # 2
k3 == a * f / e            # 3
k4 == a * b                # 4
4 * a == b + d + e + f     # 5
x     == c + d + e + f     # 6
return( c( a, b, c, d, e, f ) )
}











Livre de vírus. www.avast.com.

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.