1

我有这个简短的样本..

sysclk.value = 320e6

ddsdiv = cp.Variable(integer = True,                    # DDS divisor
         name = 'ddsdiv')
ddsdivpos = cp.Variable(pos = True)                     # Constrain to be positive

ddsclk = cp.Variable(pos = True, name = 'ddsclk')       # Resulting clock

constraints = [
        ddsdiv == ddsdivpos,                            # Constrain DDS divisor to be positive
        sysclk == ddsclk * ddsdiv,                      # Compute DDS clock
        ]
objective = cp.Minimize(cp.abs(ddsclk - 10e6))
prob = cp.Problem(objective, constraints);

prob.solve()
print(ddsclk.value)

然而它说..

DCPError: Problem does not follow DCP rules. Specifically:
The following constraints are not DCP:
sysclk == ddsclk * ddsdiv , because the following subexpressions are not:
|--  ddsclk * ddsdiv

同样，如果我尝试（更天真），它也会失败..

        ddsclk == sysclk / ddsdiv,                      # Compute DDS clock

它还会发出 DCPError。

我不明白一个正数除以一个正整数怎么不能是我的凸（但我的数学不是很好:)

我在 MacOSX 10.14.6 上使用 CVXPY 1.0.25、Python 3.7.5。

谢谢。

4

1 回答 1

2

假设您的约束是x * y == 9wherex和yare (连续) 变量。该方程的解集必须是凸集。我们可以通过采用两个解决方案并检查这​​两个解决方案之间的所有点是否也在集合中来测试这是否是凸的。

例如，(1,9) 和 (9,1) 都是有效解。然而，这些解决方案之间的中点 (5,5)不是解决方案。我们找到了一个反例，因此这不是一个凸集。

最终，由于这不是凸的，因此不可能使其满足 DCP。您需要重新表述您的问题以使其凸出或其他更容易处理的问题。

于 2019-12-05T07:37:53.727 回答