Julia中的约束
在Julia中，约束（constraints）是一种用于限制变量取值范围的机制。它们在优化问题、线性规划、非线性规划等领域中起着重要作用。
约束可以分为等式约束和不等式约束。等式约束要求某个表达式的值等于一个常数，而不等式约束则要求某个表达式的值满足某种关系，如大于、小于等。
在Julia中，可以使用JuMP（Julia for Mathematical Programming）库来处理约束。JuMP是一个优化建模语言，它提供了一种简洁的方式来定义优化问题，并支持多种优化求解器。
下面是一些常见的约束类型及其应用场景：
- 线性约束：线性约束是指约束条件中的表达式是线性的。它们在线性规划和整数规划等问题中经常出现。例如，限制某个变量的取值范围为一个线性不等式。
- 非线性约束：非线性约束是指约束条件中的表达式是非线性的。它们在非线性规划和最优化问题中常见。例如，限制某个变量的取值范围为一个非线性不等式。
- 等式约束：等式约束要求某个表达式的值等于一个常数。它们在最小二乘法和参数估计等问题中常用。例如，限制某个变量的取值满足一个等式。
- 不等式约束：不等式约束要求某个表达式的值满足某种关系，如大于、小于等。它们在优化问题中广泛应用。例如，限制某个变量的取值范围为一个不等式。
在处理约束时，可以使用JuMP库中的@constraint宏来定义约束。具体使用方法可以参考JuMP的官方文档（https://jump.dev/JuMP.jl/stable/constraints/）。
腾讯云提供了一系列与云计算相关的产品，如云服务器、云数据库、云存储等。这些产品可以帮助用户快速搭建和管理云计算环境，提供稳定可靠的计算和存储能力。具体产品信息可以参考腾讯云官方网站（https://cloud.tencent.com/）。
请注意，本回答仅涵盖了Julia中的约束概念及其应用，不涉及其他云计算品牌商。
相关·内容
有没有人能帮我写这段代码：
struct GenericPoint{T<:Real}
x::T
y::T
end
居里昂{T<:Real}中的<:是什么意思？
浏览 1提问于2020-11-20得票数 3
在帖子的启发下，我想问同样的问题，但有一个陷阱。
考虑：
struct HappyImmutable
a::Int
end
我不能改变这个领域。当然，它是不可改变的：
julia> struct HappyImmutable
a::Int
end
julia> b = HappyImmutable(1)
julia> b.a
1
julia> b.a = 2
ERROR: setfield! immutable struct of type HappyImmutable cannot be changed
Stacktrace:
[1] setp
浏览 0提问于2019-08-20得票数 1
1回答
当使用类型别名时，Julia (v1.0)不会在方法定义之外传播空闲类型参数：
const RT{R<:Real} = Type{R}
rt(::RT{R}) where R = R
rt2(::Type{R}) where {R<:Real} = R
# there's more logic here for other subtypes
这对于R<:Real来说很好
julia> rt(Int), rt2(Int) # works for both
(Int64, Int64)
但是rt()可以自由地接受非<:Real输入。
julia> r
浏览 0提问于2019-02-12得票数 2
回答已采纳
1回答
当函数最小化时，我试图防止一个值超过使用向量中的另一个值。
例如，这段代码优化了init_x
using Optim
f(x) = (1.0-x[1])^2+100.0*(-x[2]^2)^2
init_x = [0.0,0.0]
res = Optim.minimizer(Optim.optimize(f, init_x))
返回：
2-element Vector{Float64}:
0.9999860561511018
-0.0006072923701139431
在优化过程中，如何防止res (0.9999860561511018)的第一个值超过res的第二个值(-0.0006
浏览 1提问于2022-01-01得票数 2
回答已采纳
1回答
有没有人能解释一下如何更新Julia1.0 function _encode_zigzag{T <: Integer}(n::T)
num_bits = sizeof(T) * 8
(n << 1) ⊻ (n >> (num_bits - 1))
end 以及与以下各项的区别： function _encode_zigzag(n::Integer)
num_bits = sizeof(T) * 8
(n << 1) ⊻ (n >> (num_bits - 1))
end
浏览 11提问于2019-05-11得票数 0
回答已采纳
1回答
我对编译知之甚少，但我惊讶地发现，Julia编译器并没有优化几个过程。
让我们来看看Julia (这是一个即时编译器)，让我们来考虑这两段代码，它们本质上都是一样的。
n=1
@time for i = 1:10^8 n=n+1 end
elapsed time: 3.535394087 seconds (0 bytes allocated)
n=1
@time n=n+10^8
elapsed time: 6.599e-6 seconds (112 bytes allocated)
为什么现代编译器无法理解这个长循环只会将10^8添加到n中。
下面的例子更让我吃惊
n=1
@time fo
浏览 1提问于2015-01-28得票数 2
回答已采纳
2回答
朱莉娅的数理逻辑
、、
Jump是否支持在模型中直接使用数理逻辑？
我知道它支持条件指示器，但是这些条件条件呢？
例如：
1- bi-conditional indicator是一个二元变量，与当前解决方案相关的约束的满意和不满意状态相关联。
δi = 1 ⟺ sum(j, a_ij. x_j) <= b_i
2- Either-or Constraints确保至少满足两个约束中的一个。
3- If-Then-Else Constraints描述了一种情况，在这种情况下，如果约束C1成立，则必须保持约束C2；否则必须保持约束C3
浏览 7提问于2021-05-31得票数 1
2回答
由于，.4 * .4 = 0.16000000000000003在朱莉娅。我想得到0.16的数学正确答案，以一种CPU-高效的方式。我知道圆形()是可行的，但这需要事先知道答案所占的小数位数，所以它不是一个通用的解决方案。
浏览 2提问于2015-05-15得票数 10
回答已采纳
1回答
我对参数化方法是如何在Julia中起作用感到困惑。我也希望参数是正确的词在这里使用。我已经阅读了有关方法的文档，但仍然不清楚为什么会发生以下错误。如果我按以下方式定义一个函数
function Bazinga{T<:Real}(mat ::Union{Array{T,2},SparseMatrixCSC})
mat^4
end
然后跑
Penny = sparse(randn(10,10))
Bazinga(Penny)
我得到了
ERROR: MethodError: `Bazinga` has no method matching Bazinga(::SparseMatrix
浏览 2提问于2016-08-22得票数 3
回答已采纳
2回答
我正在使用Julia/跳转来实现一个算法。在第一部分中，我定义了一个连续变量的模型，并对线性模型进行了求解。我做了一些其他的计算，在此基础上，我在模型中添加了一些约束，然后我想要解决同样的问题，但是要使用整数变量。我不能使用convert()函数，因为它不接受变量。
我试图再次将变量定义为整数，但模型似乎没有考虑它！我在这里提供了一个示例代码：
m = Model()
@defVar(m, 0 <= x <= 5)
@setObjective(m, Max, x)
@addConstraint(m, con, x <= 3.1)
solve(m)
println(getVa
浏览 3提问于2015-01-04得票数 2
回答已采纳
1回答
我有几个布尔矩阵，我想要一个结果矩阵，它表明这些矩阵的那个位置上的任何元素是否为真。在Julia语言中，有没有一个函数可以让我对任意数量的矩阵进行元素级OR运算？ # My data
a = Bool[1 0; 1 1]
b = Bool[0 0; 1 1]
c = Bool[0 0; 0 0]
d = Bool[0 0; 1 1]
# Arrays of Bool Arrays
z1 = [a]
z2 = [a, b]
z3 = [b, c, d]
z4 = [a, b, c, d]
z100 = [rand(Bool, 2, 2) for i in 1:100]
# Expected
浏览 42提问于2021-11-03得票数 6
1回答
约束数组Julia
、、、、
我试图得到一系列的约束，但我不断地得到不同类型的错误，我不知道我做错了什么。
数据数组为77x9 (整数值)。
foods是大小为77x1的列向量，是数组类型的变量。
下是一个9x1的列向量，是一个整数向量。
我应该有9个约束
这就是我所拥有的
model2 = Model()
@variable(model2, foods[i=1:77] >= 0) # Quantity of food
for i ∈ 1:9
for j ∈ 1:77
@constraint(model2, c2[i], sum(data[j][i]*foods[j])<=lower[i])
end
浏览 6提问于2022-02-08得票数 2
回答已采纳
你能帮帮我吗?为什么在我运行下面的代码时会出现“type VariableIndex has no field head”的错误呢？ using JuMP,CPUTime, Distributions, Ipopt,Juniper,Gurobi
#-----Model parameters--------------------------------------------------------
sig, C1, c0 = 2, 300, 10;
E, landa, T0, T1, T2, gam1, gam2, a1, a2, a3, ap = 0.05, 0.01, 0, 2, 2,
浏览 33提问于2021-11-17得票数 1
回答已采纳
命令"Pkg.add("ConditionalJuMP)"有效。 然而，当我尝试命令(using ConditionalJuMP)时，我得到了这个错误。有人能给我提个建议吗？ ?
浏览 8提问于2020-10-06得票数 0
我在为朱莉娅的优化而奋斗。
我过去经常使用Matlab，但我试图用朱莉娅代替。
下面是我编写的代码。
using Optim
V = fill(1.0, (18,14,5))
agrid = range(-2, stop=20, length=18)
dgrid = range(0.01, stop=24, length=14)
#zgrid = [0.5; 0.75; 1.0; 1.25; 1.5]
zgrid = [0.7739832502827438; 0.8797631785217791; 1.0; 1.1366695315439874; 1.2920176239404275
浏览 3提问于2021-10-25得票数 1
回答已采纳
我试图在木星笔记本上使用朱莉娅在VS代码中的多线程，但我发现的所有帮助都表明，VS代码扩展具有一个很好的GUI，它允许您在扩展的设置中更改"Julia: Num Thread“下的线程数，如下所示：
然而，这对我来说是不可用的(Windows 10、Julia 1.7.2和VS Code Julia Extensionv1.6.6)。相反，它说“编辑在settings.json”，这是我完全不熟悉的。对我来说是这样：
值得注意的是，打开VS代码Powershell并使用$env:JULIA_NUM_THREADS=4更改环境变量确实有效.但只有终端在VS代码和留下我的VS代码朱庇特
浏览 7提问于2022-03-24得票数 1
回答已采纳
2回答
循环遍历JuMP约束中的变量子集
、、、、
我正在尝试实现一个弧流问题，其中我在一个数组中有一组弧线。每条弧都是由from/to节点组成的自定义数据结构。我想添加一个约束，其中我只包含从特定节点(1)开始的弧，类似于： @constraint(m, sum(x[a] for a in arcs; a.from==1) == 1) 这不起作用。处理这个问题的最佳方法是什么？有没有一种方法可以在不预先计算每个节点的所有输出弧线的情况下做到这一点？如果是这样，有没有办法添加额外的条件？提前感谢 贝尔纳多
浏览 10提问于2020-10-20得票数 2
回答已采纳
1回答
当类型推断错误(::Any in @code_warntype printout)时，我的理解是函数调用是动态分派的。换句话说，在运行时，将检查参数的类型，以找到具体参数类型的专门化(MethodInstance)。需要在运行时而不是编译时这样做会带来性能成本。
(编辑:最初，我说的是类型检查和专业化查找之间的“多重分派找到合适的方法”，但实际上我不知道这个部分是否在运行时发生。似乎只有在不存在有效的专门化且需要编译的情况下才需要这样做。)
在只需要检查一个参数的具体类型的情况下，是否有可能进行更快的动态单次分派，例如在某种类型的专门化查找表中？我只是找不到一种方法来访问和调用MethodIn
浏览 3提问于2022-01-02得票数 5
我想从朱莉娅那里检索目前正在运行的朱莉娅翻译的路径。在Python中，这可以通过sys.executable来实现。
浏览 2提问于2017-04-24得票数 3
回答已采纳
2回答
这是为了简化我问的问题的一部分
我想编写一些代码，这些代码可以保证在满足特定条件的类型上工作。假设今天我编写了一些代码：
immutable Example
whatever::ASCIIString
end
function step_one(x::Example)
length(x.whatever)
end
function step_two(x::Int64)
(x * 2.5)::Float64
end
function combine_two_steps{X}(x::X)
middle = step_one(x)
result = step_
浏览 5提问于2016-02-25得票数 8
我已经在我的macOS上安装了Julia1.7，因为我不能直接从终端调用它，所以我想将它添加到路径中，为此，我使用以下命令
ln -s /Applications/Julia-1.7.app/Contents/Resources/julia/bin/julia /usr/local/bin/julia
然而，结果是拒绝了权限。
有谁有解决办法吗？
干杯
浏览 4提问于2022-05-05得票数 0
我有大约50个函数，应该只消耗正数。现在，我正在用"if“检查每次输入的数字是否为零：
function grof(x::Int)
if (x % 2) == 0
println("good")
else
throw("x is not an even number!!!!!!!!!!!!! Stupid programmer!")
end
end
理想情况下，我希望有一个自动生成此类型的数据类型，即
function grof(x::EvenInt)
println("good&#
浏览 2提问于2020-04-09得票数 4
回答已采纳
3回答
我正在朱莉娅中构建一个两个元素数组，其中每个子数组有一个不同的类型(一个是Int64的向量，另一个是Float32的数组)。
代码下面自动地将我想成为Int64的元素转换为Float32，这正是我不想要的：
my_multitype_array = [ collect(1:5), rand(Float32,3) ]
生成的数组自动将第一个数组中的Int64(通过collect(1:5)定义)转换为Float32，生成的my_multitype_array具有2-element Array{Array{Float32,1}}类型。如何强制它使第一个子数组保持Int64s？在用值填充my_mult
浏览 3提问于2020-07-20得票数 4
回答已采纳
1回答
在Julia中对这个约束进行建模时遇到了一些问题
现在我有
for k=1:v
for j = 2:nodes
@constraint(model,sum(x[i,j,k] for i=1:nodes) == sum(x[j,i,k] for i=1:n))
end
end
其中节点是客户集，n包括仓库和仓库克隆。K是车辆的数量，i是开始节点，j是结束节点。
浏览 0提问于2017-01-06得票数 2
在C++中，您可以编写：
double x = 0.0, y = 0.0, z = 0.0;
而不是：
double x = 0.0; double y = 0.0; double z = 0.0;
若要定义相同类型的三个变量，请在本例中使用double。
在“朱莉娅”一书中，相似之处是：
x::Float64 = 0.0, y::Float64 = 0.0, z::Float64 = 0.0
是否有一个语法糖，像上面的C++，允许我省略这三种类型的说明符，当它们都一样的时候用一个替换它们？
注：，我知道我能写
x = y = z = 0.0
让朱莉娅推断出这种类型。这不是我想要的。我想明确地
浏览 3提问于2016-05-14得票数 3
我正在使用Python从web上抓取数据。我想用这些数据在Julia中运行计算。
有没有可能在Julia中调用函数并返回其结果，或者我只是导出到CSV并以这种方式加载数据会更好？
浏览 2提问于2017-10-09得票数 3
我正在尝试比较Matlab和Julia的速度和性能。我正在看一个代码，它对给定负载下的连续体结构进行拓扑优化。我正在查看的代码是公共代码topopt88.m：
本质上，它是一种迭代算法，在每次迭代中求解一个Ax=b系统(x=A\b)，其中A取决于结构设计(它是有限元刚度矩阵)，并在每次迭代中更新。
在Julia中，相同的代码运行速度比Matlab慢。我在Julia中做了一些代码优化，在函数定义中声明类型，尽可能多地使用函数，避免全局变量，并实现我在互联网上找到的其他技巧。但是Julia仍然比相同的Matlab代码慢(在概念步骤上是相同的)。
我的问题是:既然Matlab系统求解"\“是
浏览 9提问于2020-09-28得票数 2
我正在尝试将函数f(x::Array{Float64, 1}) -> Array{Float64, 2}应用于Julia中的m x n数组中的每一行(如果相关的话，我使用的是v1.1 )。
现在，我希望可以简单地应用mapslices，如下所示：
# toy example of f
f = (x -> randn(length(x), length(x)))
A = randn(100, 50)
# intent: apply f to every row in A and collect the result into a 100 x 50 x 50 matrix.
resu
浏览 1提问于2021-05-18得票数 2
回答已采纳
我正在尝试在Julia中实现OLS回归作为一种学习练习。我想要的一个特性是将公式作为参数(例如，'formula = Y ~ x1 + x2'，其中Y、x1和x2是DataFrame中的列)。这是一个。
如何将公式/表达式“映射”到正确的DataFrame列？
浏览 2提问于2016-10-06得票数 2
1回答
在Julia中，可以创建的线程数量是否仍然受到物理可用内核数量的限制？或者，Julia lang自3月份以来是否发生了变化，从而将线程的数量限制为逻辑核心的数量？
浏览 22提问于2020-11-26得票数 3
1回答
我有不同的数学表达式，必须在运行时进行评估。目前，该任务是通过用等效值替换符号并对结果进行eval来完成的。(可以通过任何现有的符号包来完成)
现在，请参阅Julia-lang中模块的定义：
Julia中的模块是独立的可变工作空间，即它们引入了一个新的全局范围。模块允许您创建顶级定义(也称为全局变量)，而无需担心代码与其他人的代码一起使用时的名称冲突。
有了朱莉娅的力量，
我想知道是否有可能在运行时m=Module()中创建匿名模块，并将它们用作计算表达式eval(m, :(a+b))的作用域。
但我根本找不到将变量加载到运行时模块的方法。虽然我可以得到以下结果：
julia>
浏览 4提问于2016-10-12得票数 3
回答已采纳
1回答
我在我的Windows 10笔记本电脑上安装了Julia编程语言。我想将它添加到系统环境变量中的PATH变量中，但是我无法找到可执行文件。你能帮我解决这个问题吗？
浏览 2提问于2020-05-31得票数 4
回答已采纳
但是，在尝试了不同的选项之后，我尝试在Python中运行朱莉娅.jl文件，这些选项都没有工作。
我试过使用PyJulia。导入Julia并定义Julia对象。没有实现。
有没有人将Python和Julia的技术混合在一起，并取得了成功？(从Python中运行Julia )
浏览 0提问于2018-04-10得票数 12
回答已采纳
例如，我有一个基本类型的抽象类型，我希望实现一个类型工厂，它在这个抽象类型下动态地创建具体类型，并使用用户输入参数给出的名称(和其他类型特征)。
abstract type AbstractFoo end
function TypeFactory(typename::String, supertype::DataType)
...
return ConcreteSubtype
end
函数TypeFactory使用typename和supertype参数，并创建一个属于supertype的具体类型，其名称与typename相同。
我知道这种类工厂在Python (例如)中实现
浏览 1提问于2018-02-07得票数 3
回答已采纳
function main()
(k, a, b, a1, b1) = (BigInt(2), BigInt(4), BigInt(1), BigInt(12), BigInt(4))
while true
(p, q, k) = (k*k, BigInt(2)*k+BigInt(1), k+BigInt(1))
(a, b, a1, b1) = (a1, b1, p*a+q*a1, p*b+q*b1)
(d,d1) = ( div(a,b),div(a1,b1) )
while d == d1
浏览 1提问于2018-08-13得票数 2
在REPL中调试或运行julia代码时，我通常会看到显示... at ./REPL[161]:12 [inlined]...的错误消息。我猜，数字161是指161-th在REPL中的值。所以我的问题是，我们可以在julia的提示符中显示这个数字吗，即用julia [161]>代替julia>
浏览 3提问于2018-01-02得票数 5
回答已采纳
1回答
假设有一个向量是一个元组$a$，我想在julia中定义一个函数p(x)=x^a。
例如，如果a=(1,2,3)，则结果函数为x^1 *y^2 * z^3。
对于任何元组，我都希望有一个通用的方法，但是我不知道批准的符号。在我的代码中，我有一个元组数组，我想为数组中的每个元组定义一个单项。
非常感谢你们的合作。
浏览 3提问于2022-11-01得票数 3
回答已采纳
我知道有几篇文章介绍了如何在以前的julia版本中清除工作区和变量。 但是，我使用的是Julia v"1.0.3“和Atom的最新更新。我还有最新版本的Julia "1.1.0“ 不过，我已经试过了 x=49
workspace()
clear!(:x)
# and even
clearconsole() # working but just clear the concle
Control + l # working but just clear the concle 他们中没有一个人为我清除了变量。同样，workspace()和clear!(:x)返回错误消息 ERRO
浏览 125提问于2019-03-11得票数 1
1回答
我有一个混合整数编程问题。目标函数是向量中最大变量值的极小化。变量的上限为5，问题如下：
m = Model(solver = GLPKSolverMIP())
@objective(m, Min, max(x[i] for i=1:12))
@variable(m, 0 <= x[i] <= 5, Int)
@constraint(m, sum(x[i] for i=1:12) == 12)
status = solve(m)
max变量不是julia JuMP语法的一部分。所以我把问题修改为
t=1
while t<=5 && (status == :No
浏览 1提问于2017-10-02得票数 3
回答已采纳
3回答
我在Ipython中准备了一个小的演示，在这里我想展示在Julia中进行并行操作是多么的容易。
它基本上是蒙特卡洛Pi计算
问题是我不能让它在IPython (木星)笔记本中并行工作，它只使用一个。
我创办了朱莉娅：julia -p 4
如果我在REPL中定义了函数并在那里运行，它就能正常工作。
@everywhere function compute_pi(N::Int)
"""
Compute pi with a Monte Carlo simulation of N darts thrown in [-1,1]^2
Returns esti
浏览 0提问于2015-06-23得票数 8
回答已采纳
问题:如何使用检查变量的字节大小？
我尝试过的内容：在Matlab中，whos()函数提供了这个信息，而在朱莉娅中，它只提供变量名称和模块。在朱莉娅手册中浏览标准库时，sizeof()看起来很有希望，但它只提供了规范二进制表示的大小，而不是当前变量。
浏览 8提问于2015-02-09得票数 16
回答已采纳
1回答
是否有可能在朱莉娅0.6中获得unicode字符的名称？
在python中，可以通过
s = u'\N{SECTION SIGN}'
有相似之处吗？我知道我能做到
s = '\u00a3'
但我想用名字来做。
谢谢
浏览 0提问于2018-05-11得票数 3
回答已采纳
如何将向量中的字符或字符串串联在一个字符串中，从而使"a“、"b”、"c“成为"abc"？
我试过用vcat，hcat，但似乎什么都不管用.谢谢
回答
join(["a", "b", "c"])
浏览 0提问于2018-08-28得票数 8
回答已采纳
2回答
我是新来的朱莉娅。我开发了几行代码，以便从Python或R中找不到的包中获得所需的结果。现在，我正在尝试使该文件易于访问，并将代码包装在Python或R中。以前有人这样做过吗？我已经尝试了几种方法，但没有找到任何有用的方法。
要做到这一点，最简单的方法是只需几行代码调用.jl文件，运行它(然后将代码添加到julia的.txt文件中)，然后在代码完成时提醒您。
任何帮助都将不胜感激。R是更好的方法，在这一点上，Python也会受到欢迎。
浏览 0提问于2020-12-01得票数 2
1回答
我非常喜欢的想法，它是介于Tuple和完整的、用户定义的复合类型之间的中间地带。
我知道如何构建一个命名元组并访问它的一个字段。
julia> nt = (a=1, b=2.0)
(a = 1, b = 2.0)
julia> nt.a
1
然而，我不知道更多，甚至不知道是否可能做更多的事情。我正在考虑许多方法，我们可以操作普通元组(通常涉及飞溅)，并怀疑其中一些方法是否也适用于命名元组。例如，如何：
从字段和值列表中动态构建NamedTuple，增长NamedTuple，即向其添加新的字段-值对(在不可变的意义上)--在现有的命名为tuple<code>f 212&
浏览 3提问于2020-03-27得票数 8
我试着理解“朱莉娅”中copy()和deepcopy()之间的区别。根据我在中看到的内容，deepcopy()似乎复制了值，然后创建了一个与我复制的原始对象无关的新对象。那部分说得通。对于下列对象之间的关系，我更感到困惑：
julia> a = [1,2,3]
3-element Array{Int64,1}:
1
2
3
julia> b = copy(a)
3-element Array{Int64,1}:
1
2
3
julia> a == b
true
julia> isequal(a,b)
true
也许这只是我上面选择的一个糟糕的例子，但我
浏览 17提问于2020-02-01得票数 6
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
即时通信 IM
对象存储活动推荐
腾讯云开发者
