我正在尝试优化Excel中的三个参数,以便最小化实验值和理论值之间的误差。我在for循环中为每个参数使用Solver,一次一个。但是,我想迭代这个求解器for循环(循环内循环),直到实验值中的误差和理论值小于某个目标值。
我的实验值为$K25
我的理论值(根据我的模型方程计算)是$J$25
我需要优化的参数是$C$4
,$C$5
,$C$6
当我运行以下VBA代码时,$C$4
,$C$5
,$C$6
中的参数不会更改其初始值。但是,宏编译没有错误。任何人都可以帮助我吗?
以下是代码:
Sub Macro3()
Application.ScreenUpdating = False
SolverReset
Dim j As Integer
For j = 1 To 100 Step 1
If "$J$25" > "$K$25" Then
Dim i As Integer, s As String
For i = 4 To 6 Step 1
s = Format(i, "0")
SolverOk SetCell:="$J$25", MaxMinVal:=2, ValueOf:=0, ByChange:="$C$" & s, Engine:= _
1, EngineDesc:="GRG Nonlinear"
SolverOptions MaxTime:=0, Iterations:=1000000, Precision:=0.000001, Convergence _
:=0.00001, StepThru:=False, Scaling:=True, AssumeNonNeg:=True, Derivatives:=1
SolverOptions PopulationSize:=100, RandomSeed:=0, MutationRate:=0.075, Multistart _
:=False, RequireBounds:=True, MaxSubproblems:=0, MaxIntegerSols:=0, _
IntTolerance:=1, SolveWithout:=False, MaxTimeNoImp:=30
SolverOk SetCell:="$J$25", MaxMinVal:=2, ValueOf:=0, ByChange:="$C$" & s, Engine:= _
1, EngineDesc:="GRG Nonlinear"
SolverSolve (True)
SolverReset
Next i
End If
Next j
Application.ScreenUpdating = True
End Sub
答案 0 :(得分:1)
我不确定你是否需要在VBA中这样做,因为你正在寻找的正是解决者应该做的事情 - 修改一组参数,以便最大化/最小化其他东西!
因此,您需要做的就是在另一个单元格中插入公式=ABS(J25-K25)
。此单元格将显示实验值与理论值之间的差值。现在设置你的求解器,以便通过改变你的三个参数来最小化这个单元 - 你已经完成了! (请注意,您可以在“通过更改可变单元格”字段中提供多个单元格!)
如果你想坚持你的方法,这里是语法正确的代码。请注意,我还没有对它进行测试 - 但只是通过查看代码来纠正我可以发现的错误。它有望成为一个很好的起点。事实上,看看这种方法,我确信你最终会得到错误的结果,因为每次运行只优化一个变量 - 因此你永远不会考虑由两个或三个参数组合产生的任何效果!
无论如何,这是你的代码:
Sub RunSolver()
Dim j As Integer, i As Integer
Application.ScreenUpdating = False
SolverReset
For j = 1 To 100
Application.Statusbar = j & "/100"
If Range("$J$25") > Range("$K$25") Then
For i = 4 To 6
SolverOk SetCell:=Range("$J$25"), MaxMinVal:=2, ValueOf:=0, ByChange:=Range("$C$" & i), Engine:= _
1, EngineDesc:="GRG Nonlinear"
SolverOptions MaxTime:=0, Iterations:=1000000, Precision:=0.000001, Convergence _
:=0.00001, StepThru:=False, Scaling:=True, AssumeNonNeg:=True, Derivatives:=1
SolverOptions PopulationSize:=100, RandomSeed:=0, MutationRate:=0.075, Multistart _
:=False, RequireBounds:=True, MaxSubproblems:=0, MaxIntegerSols:=0, _
IntTolerance:=1, SolveWithout:=False, MaxTimeNoImp:=30
SolverSolve (True)
SolverReset
Next i
End If
Next j
Application.StatusBar = False
Application.ScreenUpdating = True
End Sub
答案 1 :(得分:1)
您可以仔细检查代码中的那行:
Engine:= 1, EngineDesc:="GRG Nonlinear"
可能您的目标函数是非线性,并且您认为您正在使用GRG非线性求解器,因为您在 EngineDesc 参数下提及它。这是不正确的。这只是一个描述参数。
您实际使用的解算器是 Simplex LP ,其值 1 。
更改为 2 以使用 GRG非线性求解器。