MATLAB代码："Stackelberg博弈下的光伏用户群内部电价优化定价模型
光伏用户群里最头疼的问题就是"电怎么分"。各家屋顶上的太阳能板白天发电用不完，想卖给邻居
还得走电网过一道手，一来二去成本上去了，谁都没积极性。这就像菜市场里卖菜的想直接卖给隔壁餐馆，
结果非得让中间商赚差价，搁谁都不乐意。
这时候运营商这个"中间商"突然开窍了：与其被政策卡脖子，不如自己搞个内部交易市场。我们来
看MATLAB里怎么用Stackelberg博弈玩转这个局。先看运营商定价的核心代码：
```matlab
function [price,profit] = operator_optimize(user_response)
cvx_begin
variable p(n_users)
maximize(sum(p.*user_response) - 0.5*sum_square(p))
subject to
p >= 0.35
p <= 0.6
cvx_end
price = p;
profit = cvx_optval;
end
```
这段代码用CVX工具包求解运营商的最优定价。目标函数里`sum(p.*user_response)`是售电收入，
后面那个二次项其实是防止电价波动过大的平滑处理。注意约束条件里的0.35-0.6元/度，这个价格区间
可不是随便定的——比电网收购价高但比零售价低，既让用户有动力内部交易，又能保证运营商有赚头。
用户端的需求响应更有意思。每家每户都跟精算师似的，拿着内部电价算小账：
```matlab
classdef UserModel
properties
base_load
pv_generation
end
methods
