电力负荷分配matlab程序

function [leastCostPt]=ecoDisp（coeffPminPmaxUPload）

% 本函数用于求解某个时段内各投入机组的最优发电功率和最小损耗=====================>>经济调度

% leastCost是该时段内所有机组的最小发电成本========>>经济调度
% Pt为该时段内所有机组的最优发电功率

% coeff为所有机组的成本系数矩阵，设为二次函数形式[  二次项系数   一次项系数    常数项系数    ]
% 在此假定发电机成本函数为二次函数，coeff=[c2c1c0]为发电机成本函数的系数
% PminPmax为所有机组的发电功率下限和上限向量（列向量）
% U 为由动态规划求解得到的所有机组在该时段内的运行状态（列向量）
% Pload为该时段内的负荷（标量）

n=size（coeff1）;% n 为机组数包括未投运的机组和已投运机组
c2=zeros（n1）;% 二次项系数
c1=zeros（n1）;% 一次项系数
c0=zeros（n1）;% 常数项系数
c2=coeff（:1）;% 二次项系数
c1=coeff（:2）;% 一次项系数
c0=coeff（:3）;% 常数项系数

leastCost=0;
Pt=zeros（n1）;%各机组在该时段内的发电功率

lambda=10; % 设定拉格朗日乘子的初值（可以改变）
epsilon=0.1; % 允许的误差
lostPower=0.1;% 负荷平衡约束条件的偏差
sumPower=0; % 各机组的发电功率之和

maxPower=0;% 投运机组的最大发电功率之和
for i=1:n
   maxPower = maxPower + U（i）*Pmax（i）; % 计算投运机组的最大发电功率
end  

if（ Pload<=maxPower）                  % 如果系统负荷大于投运机组的最大发电功率，则不必进行计算
   whil