Skip to content

逻辑非

y=Not(x)={1,¬x0,x

离散

常量

m=max(x)

额外变量

y{0,1} :多项式非的逻辑值。

导出符号

y=y

数学模型

s.t.m(1y)x1yx

即:

y={1,x=00,x1

连续(非精确)

常量

m=max(x)

额外变量

b[0,1] :多项式的归一化值。

y{0,1} :多项式非的逻辑值。

导出符号

y=y

数学模型

s.t.x=mb1ybϵ(1y)b

即:

y={1,x=00,xϵ

连续(精确)

导出符号

y=Ulp(x)

Ulp(x) 可参考 一元分段线性函数,使用的采样点为 {(0,1),(pϵ,1),(p,0),(max(x),0)}

即:

y={1,0x<p0,xp

样例

离散

kotlin
import kotlinx.coroutines.*
import fuookami.ospf.kotlin.utils.math.*
import fuookami.ospf.kotlin.core.frontend.variable.*
import fuookami.ospf.kotlin.core.frontend.expression.polynomial.*
import fuookami.ospf.kotlin.core.frontend.expression.symbol.linear_function.*
import fuookami.ospf.kotlin.core.frontend.inequality.*
import fuookami.ospf.kotlin.core.frontend.model.mechanism.*
import fuookami.ospf.kotlin.core.backend.plugins.scip.*

val x = BinVar("x")
val not = NotFunction(x, name = "not")
val solver = ScipLinearSolver()

val model1 = LinearMetaModel()
model1.add(x)
model1.add(not)
model1.addConstraint(not)
model1.maximize(x)
val result1 = runBlocking { solver(model1) }
assert(result1.value!!.obj eq Flt64.zero)

val model2 = LinearMetaModel()
model2.add(x)
model2.add(not)
model2.addConstraint(!not)
model2.minimize(x)
val result2 = runBlocking { solver(model2) }
assert(result2.value!!.obj eq Flt64.one)

连续

kotlin
import kotlinx.coroutines.*
import fuookami.ospf.kotlin.utils.math.*
import fuookami.ospf.kotlin.core.frontend.variable.*
import fuookami.ospf.kotlin.core.frontend.expression.polynomial.*
import fuookami.ospf.kotlin.core.frontend.expression.symbol.linear_function.*
import fuookami.ospf.kotlin.core.frontend.inequality.*
import fuookami.ospf.kotlin.core.frontend.model.mechanism.*
import fuookami.ospf.kotlin.core.backend.plugins.scip.*

val x = URealVar("x")
x.range.leq(Flt64.one)
val not = NotFunction(x, name = "not")
val solver = ScipLinearSolver()

val model1 = LinearMetaModel()
model1.add(x)
model1.add(not)
model1.addConstraint(not)
model1.maximize(x)
val result1 = runBlocking { solver(model1) }
assert(result1.value!!.obj eq Flt64.zero)

val model2 = LinearMetaModel()
model2.add(x)
model2.add(not)
model2.addConstraint(!not)
model2.minimize(x)
val result2 = runBlocking { solver(model2) }
assert(result2.value!!.obj geq Flt64(1e-6))

完整实现请参考:

完整样例请参考: