Punishing Term for Smoothing Optimization Results

Question

Mathias Dirksmeier 2018-8-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/416068-punishing-term-for-smoothing-optimization-results

评论： Mathias Dirksmeier 2018-8-24

采纳的回答： Matt J

image.png

在 MATLAB Online 中打开

Hi all,

this might be an interesting brainteaser: I am optimizing electric storage dispatch (fmincon solver) with a function to be optimized as follows:

fun = @(x)sum(x.*p)
,p is given as a price vector, x is the storage dispatch in kWh.

There are some constraints to this so that optimization results in a vector x. So far so good.

Now, plotting x reveals its jumping behavior (cp. attachment). However, I would rather have it smoothed just as or at least alike the price vector (Clearing Price, also cp. attachment). The jumping seems arbitrary (so far, there is no punishment for jumping).

How do I punish for jumping in my optimization function? I could not think of a way to include this into my constraints as I do not want to define strong constraints rather a little punishment (just strong enough to prevent jumping when it does not serve to optimality but still allow big leaps if necessary).

My idea so far: Add a punishment term to the optimization function like this:

fun = @(x)sum(x.*p+(x(n+1)-x(n))/1000)

That way, there should be a little punishment for all changes in x. The problem, that would require a loop in my optimization function and I have no idea if that's possible at all. I am really looking forward to your ideas and the following discussion.

Thanks in advance, Mathias

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Matt J 2018-8-24

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/416068-punishing-term-for-smoothing-optimization-results#answer_333949

在 MATLAB Online 中打开

fun = @(x) x(:).'*p(:) + beta*norm(diff(x)).^2 ;

3 个评论
显示 1更早的评论隐藏 1更早的评论

Matt J 2018-8-24

编辑：Matt J 2018-8-24

The penalty term norm(diff(x)).^2 is the same as sum(x(n+1)-x(n))^2. It is a very traditional roughness penalty.

beta is a coefficient you must define to control the weight of the penalty term. The error message is because you have not defined it, so Matlab thinks you mean to invoke the command beta().

Mathias Dirksmeier 2018-8-24