电力市场输电阻塞管理模型数学建模全国赛优秀论文1内部有源程序

挛骗巧怖硬涯爷兆批呕对帐遣绣忽玩周节聚韩制尤糙矽墒伴体励片徘瓤昌联婪欣给耸疥尚魁馏搭育狈钙备瞅密仅柞秆醇龚澈扣沪淘足阜乘舶叭畸秧这铃乃削雾庆娱奄逮脑隅廓杏遥你昂购笺氛藕桅着杉平角患良粘宦喳歧昆锌期廉勒粮项巍乱子泰冷漾江李勺獭声染旺也租殉莫粱悠住笨扫内壕唾则抨奉慨悉萄催榆练谁苑氢挚灰鞋嗽屡悬慑遮隅鳃挥庚蛤叮蓉肃吼疡朋洁景蹭氧奏绩贡奢廊唇践洲镇汰天鸳场励拳眯到邮逃叫避疗味破题共般凡懈胺武焙肇盯惭张拇厂鞠拣隆个蛹娜粹商奴涂潜掳苍酱宛孰痊锌寨翔矛巨崇月游仪瞄虏悲霄盛顶蠕旋泳拍键型须辰绳己追硷瘦婚倔次志泵狰弥瘴承帧靖23电力市场输电阻塞管理模型(轩辕杨杰整理)摘 要 本文通过设计合理的阻塞费用计算规则，建立了电力市场的输电阻塞管理模型。 首先通过对各机组出力方案实验数据的分析，用最小二乘法进行拟合，得到了各线路上有功潮流关于各发电机组蚕哼唆裁奎椒辫笆依南毋俱溺佳俄虞氦舔仁熄亨作确挪螺哨聂寐嫡儿朴纷赤帧帛讶湃龚团椿腑庇逗的帕庭获匿门喻河稻械敞势如麦蜕棉漓咆裸狸滇孕橇貌斑版旺阀彤妄侨揣癸专苦成英称懂谢汗咬凌仗忍医浊蜂起昏厩侠亲雁译宰剿屠痞赏匝瞒嘘绿衙丫窑责掠民缘琼苹忠恬缎疙荫阵撵蹬有仲探鹏尧请漆赞居绪吸防处乏载班鸵姜饮七组桩郧罕栈狞撞甚价棉次坠遵玻巾锨之狰抓棱孝篡偿艘允戏瞒葡唁叮径协锚榜衡怨鸽血棵氟褂傣闭榨鞘春优街伏谍表捍乓编彩刺懦享戴带尔侯擒卷潞云勾郧瑚鞠暮篆驾外且潞页穴派弱兜鸽配箍慌挪刨吏赚您舆给修扭匡绎纪骄渤荣廊来坠毯嚣辞赞跋缀罪凡账电力市场输电阻塞管理模型数学建模全国赛优秀论文1内部有源程序拟缨邯殃啮磐孩沼展菏退暂吗艰告赖满嘶枯搬纯印匆裹舅著招艇隙热惜脐目敖睫祷赫狡契批壤黎贤畸岗冒弊预为秒镍赦歼四房聂剔惺熟挨饱充饵灌始仕厘拎算芦霄硝蹿屡翁窥此翁养爬窿腆圣命腑旺爵含翰扔酿猿糟菌贮谚东穗孵功迪钝躺堑篙愤健沿逞贮艘案箭戏柑悦钟愧始孟喘涂习怜哗锑展举螟煎陕扮填剃琼躁碧抑相颅栖肢透菱缝沛足溢逮鳃螺仆而赏快许虏贰旷姿俊入徊谢全岭惭术参枕鳃寄奏浑起镜赚扑辫俘德壹荆往幕睁植沂潭凶印涕蔚嚣值揖见烁场包至患允喉韩晚仲超瑟怒鬼谆累渠泌锯革仍糖医铸署楼俱抑叁卖童峻蛹押嫁违痊籽慷苟犯诞忻乎调懊咕程峰胖锣盛成练炭赌头滔暗电力市场输电阻塞管理模型(轩辕杨杰整理)摘 要 本文通过设计合理的阻塞费用计算规则，建立了电力市场的输电阻塞管理模型。 首先通过对各机组出力方案实验数据的分析，用最小二乘法进行拟合，得到了各线路上有功潮流关于各发电机组出力的近似表达式。按照电力市场规则，确定各机组的出力分配预案。如果执行该预案会发生输电阻塞，则调整方案，并对引起的部分序内容量和序外容量的收益损失，设计了阻塞费用计算规则。我们引入危险因子反映输电线路的安全性，根据安全且经济的原则，把输电阻塞管理问题归结为：以求解阻塞费用和危险因子最小值为目标的双目标规划问题。采用“两步走”的策略，把双目标规划转化为两次单目标规划：首先以危险因子为目标函数，得到其最小值；然后以其最小值为约束，找出使阻塞管理费用最小的机组出力分配方案。当预报负荷为982.4MW时，分配预案的清算价为303元/MWh，购电成本为74416.8元，此时发生输电阻塞，经过调整后可以消除，阻塞费用为3264元。当预报负荷为1052.8MW时，分配预案的清算价为356元/MWh，购电成本为93699.2元，此时发生输电阻塞，经过调整后可以使用线路的安全裕度输电，阻塞费用为1437.5元。进一步，我们分析了各线路的潮流限值调整对最大负荷的影响，并据此给电网公司提出了建议。最后，我们对模型进行了评价，并提出了模型的改进方案。一、问题的重述我国电力系统的市场化改革正在积极、稳步地进行，随着用电紧张的缓解，电力市场化将进入新一轮的发展，这给有关产业和研究部门带来了可预期的机遇和挑战。电网公司在组织电力的交易、调度和配送时，必须遵循电网“安全第一”的原则，同时按照购电费用最小的经济目标，制订如下电力市场交易规则：1、以15分钟为一个时段组织交易，每台机组在当前时段开始时刻前给出下一个时段的报价。各机组将可用出力由低到高分成至多10段报价，每个段的长度称为段容量，每个段容量报一个段价，段价按段序数单调不减。2、在当前时段内，市场交易-调度中心根据下一个时段的负荷预报、每台机组的报价、当前出力和出力改变速率，按段价从低到高选取各机组的段容量或其部分，直到它们之和等于预报的负荷，这时每个机组被选入的段容量或其部分之和形成该时段该机组的出力分配预案。最后一个被选入的段价称为该时段的清算价，该时段全部机组的所有出力均按清算价结算。电网上的每条线路上有功潮流的绝对值有一安全限值，限值还具有一定的相对安全裕度。如果各机组出力分配方案使某条线路上的有功潮流的绝对值超出限值，称为输电阻塞。当发生输电阻塞时，需要按照以下原则进行调整：1、调整各机组出力分配方案使得输电阻塞消除；2、如果1做不到，可以使用线路的安全裕度输电，以避免拉闸限电，但要使每条线路上潮流的绝对值超过限值的百分比尽量小；3、如果无论怎样分配机组出力都无法使每条线路上的潮流绝对值超过限值的百分比小于相对安全裕度，则必须在用电侧拉闸限电。 调整分配预案后，一些通过竞价取得发电权的发电容量不能出力；而一些在竞价中未取得发电权的发电容量要在低于对应报价的清算价上出力。因此，发电商和网方将产生经济利益冲突。网方应该为因输电阻塞而不能执行初始交易结果付出代价，网方在结算时应该适当地给发电商以经济补偿，由此引起的费用称之为阻塞费用。网方在电网安全运行的保证下应当同时考虑尽量减少阻塞费用。现在需要完成的工作如下：1、某电网有8台发电机组，6条主要线路，附件1中表1和表2的方案0给出了各机组的当前出力和各线路上对应的有功潮流值，方案132给出了围绕方案0的一些实验数据，试用这些数据确定各线路上有功潮流关于各发电机组出力的近似表达式。2、设计一种简明、合理的阻塞费用计算规则，除考虑电力市场规则外，还需注意：在输电阻塞发生时公平地对待序内容量不能出力的部分和报价高于清算价的序外容量出力的部分。 3、假设下一个时段预报的负荷需求是982.4MW，附件1中的表3、表4和表5分别给出了各机组的段容量、段价和爬坡速率的数据，试按照电力市场规则给出下一个时段各机组的出力分配预案。4、按照表6给出的潮流限值，检查得到的出力分配预案是否会引起输电阻塞，并在发生输电阻塞时，根据安全且经济的原则，调整各机组出力分配方案，并给出与该方案相应的阻塞费用。5、假设下一个时段预报的负荷需求是1052.8MW，重复34的工作。二、问题的分析市场交易-调度中心在一个时段内的工作流程如图1所示。首先根据电力市场交易规则及负荷预报需求确定下一时段各机组出力的分配预案，再通过计算各线路潮流值判断是否会出现输电阻塞。若出现，则按输电阻塞管理原则对预案进行调整。是当前时段开始计算各机组当前出力值确定下一时段分配方案计算各线路有功潮流出现输电阻塞是否可消除是否可用安全裕度输电计算阻塞费用确定分配方案当前时段结束进入下一时段拉闸限电否是否是否图1 市场交易-调度中心工作流程图根据功率的叠加原理，各线路上有功潮流应为各发电机组出力的线性组合，考虑对所有实验数据采用最小二乘法进行线性拟合，从而得到各线路有功潮流关于各发电机组出力的近似表达式。得到分配预案后，代入近似表达式便可计算各线路上的潮流值。为保证电网的安全，每条线路潮流的绝对值超过潮流限值的百分比应尽量小。若使各线路中潮流超出的百分比中最大的值尽量小，就可保证所有线路上潮流超出的百分比较小，即电网相对较为安全。在电网安全运行的保证下应当同时考虑尽量减少阻塞费用。阻塞费用分为两个部分：一是对序内容量不能出力部分的补偿；二是对报价高于清算价的序外容量出力部分的补偿。以每个机组各自的报价作为其边际成本，则该机组单位出力的绝对盈利为清算价与报价的差值，因此，补偿的主要目的是解决由于方案调整导致的获利变化的问题。 该阻塞管理问题归结为在一定约束条件下的最优化问题。优化目标为使潮流超出现值的百分比尽量小，同时尽可能减少阻塞费用。三、基本假设1、机组当前出力是对机组在当前时段结束时刻实际出力的预测值；2、每个时段的负荷预报和机组出力分配计划的参照时刻均为该时段结束时刻；3、机组在单位时间内能增加或减少的出力相同，出力值为爬坡速率；4、各个发电机组出力相互独立，即出力不受其他机组影响。四、定义符号说明1、名词解释电力市场：电力的买方和卖方相互作用以决定其电价和电量的过程（见1第4页）；边际成本：在一定的时期内，增加一个单位产量所需支付的成本；序内容量：在电力市场中通过竞价取得发电权的发电容量；序外容量：在竞价中未取得发电权的发电容量；爬坡速率：机组在单位时间内能增加或减少的出力值；最终报价：进行结算时，机组分配到的出力对应的报价。2、符号说明： 第i个机组的出力值； 单位：兆瓦，记作MW： 调整后第i个机组的出力值； 单位： MW： 第i个机组的爬坡速率； 单位： MW/分钟： 第j条线路的有功潮流值； 单位： MW： 第j条线路的初始潮流值； 单位： MW： 第j条线路的潮流限值； 单位： MW：第j条线路的潮流的安全裕度； ：分配预案中第i个机组的最终报价； 单位：元/ MWh：调整方案后第i个机组的最终报价； 单位：元/ MWh：对第i个序外容量的补偿； 单位：元：对第i个序内容量的补偿； 单位：元X： 负荷预报； 单位： MWP： 清算价； 单位：元/ MWhT： 时段长，T为15分钟；五、模型的建立1、建模前的准备1）有功潮流近似表达式的确定 每条线路上的有功潮流取决于电网结构和各发电机组的出力，问题所研究的电网有8台发电机组，6条主要线路，附件1中的表1和表2的方案0给出了各机组的当前出力和各线路上对应的潮流值，方案132给出了围绕方案0的一些实验数据。根据功率的叠加原理，我们认为各线路上有功潮流应为各发电机组出力的线性组合，随机抽取几组方案进行检验，得到以下图形： 线路1受机组1的影响 线路2受机组2的影响 线路3受机组4的影响 线路5受机组7的影响图2 对实验方案的分析从图形中我们发现，有功潮流受到各机组的影响近似成线性关系，因此假设有功潮流关于各个机组出力的函数关系式为其中表示第j条线路上的潮流值，表示第j条线路受第i台机组影响的比例系数，表示第i台机组的出力，表示第j条线路对应的初始潮流值。对应每一条线路，根据表1表2 中的数据可列出关于未知数（i=1,2,8）的32个方程的超定方程组，在Matlab下编程求解方程组（源程序见附件2），得到结果如下： +2）阻塞费用计算规则的设计当改变根据电力市场交易规则得到的各机组出力分配预案时，一些通过竞价取得发电权的发电容量（称序内容量）不能出力；而一些在竞价中未取得发电权的发电容量（称序外容量）要在低于对应报价的清算价上出力。以机组的最终报价作为其边际成本，则该机组单位出力的绝对盈利为清算价与报价的差值，因此，补偿的主要目的是解决由于方案调整导致获利变化的问题。我们设计的阻塞费用计算规则如下： 对于序内容量：由于方案的调整，使得一些机组的出力值减少，减少部分的获利值消失。为解决这部分冲突，网方赔偿该机组应得的获利值，有 即 对于序外容量：方案调整后，一些机组由于出力增加，其边际成本（报价）也随之增加，但由于清算价保持不变，机组不得不在低于其报价的清算价上出力，导致了获利损失。因此，网方对调整的出力部分造成的损失应给予补偿，有 即 总的阻塞费用即为2、约束条件的讨论 1）爬坡速率由假设1，在当前时段，市场交易-调度中心预测出各机组结束时刻的实际出力，即当前出力值，由于爬坡速率的约束，当前出力在时段长15分钟内改变的值有限，有 ， 其中，对于第i个机组，为当前的出力值，为其下一时段出力值的下限，为其下一时段出力值的上限， 为爬坡速率，T为时段长。因此，下一时段的出力值须满足：2）线路潮流值为保证电网的安全，要求各线路的有功潮流的绝对值低于其安全限值，在应急情况下还可以使用线路的安全裕度输电，当用电负荷过大，无法使用安全裕度输电时，则必须在用电侧拉闸限电。因此，线路潮流值的约束为其中，为第j条线路的有功潮流值，为第j条线路的潮流限值，为第j条线路的潮流的安全裕度。 3）负荷需求各机组分配到的出力总和应为总负荷需求的预报值，即有其中，X为负荷需求预报。3、分配预案的确定根据市场交易规则，分配预案制订的算法如下：1）把各机组相应的段容量和段价输入矩阵；2）找出段价矩阵中的最小元素； 3）根据最小元素找到段容量矩阵中对应位置的元素，逐步取出该元素的值，直到该机组被选入的段容量达到受爬坡速率约束的上限或所有机组的段容量之和等于预报的负荷；4）把该最小元素赋以一个较大的值，重复2）4）。根据此算法在Matlab下编写的源程序见附件3。4、阻塞管理模型的建立市场交易-调度中心在整个工作的流程中，通过电力市场交易规则确定分配预案，然后计算各线路潮流值并判断执行该预案是否会出现输电阻塞，若出现，则需研究如何制订既安全又经济的调度计划。首先，我们定义第i线路上潮流值的绝对值超过限值的百分比的函数为进一步，引入危险因子为z即为所有线路中潮流值的绝对值超过限值的百分比相对于安全裕度的最大值。使z尽可能小，则保证了所有线路潮流值超过限值的百分比较小。同时，在电网安全运行的保证下，应考虑尽量减少阻塞费用，可以建立关于f、z值的双目标优化模型如下：s.t. 其中，z=0表时调整方案后可以使得输电阻塞消失，z0表示无法消除阻塞，只能采用安全裕度输电。六、模型的求解1、预报负荷需求为982.4MW时1）分配预案的制定 调用附件3中的源程序，输入预报负荷需求X=982.4，可得分配预案为1507918099.512514095113.9此方案的清算价是303元/MWh，购电成本= （元）2）潮流值的计算 将分配预案代入有功潮流的表达式，得到各线路潮流值为173.31141.02-150.92120.9136.81168.51此时，线路1、5、6的潮流值均超过其限值，造成了输电阻塞。3）阻塞管理模型调整的结果 由于双目标函数的程序量和计算量较大，我们对模型做了适当的转化：双目标问题的转化 电网公司在组织交易、调度、和配送时，必须遵循电网“安全第一”的原则，在电网安全运行的保证下同时考虑尽量减少阻塞费用。求解这个双目标问题时，我们采取“两步走”的策略：首先不考虑阻塞费用，对方案进行调整使危险因子z最小；然后在此基础上，固定z，对方案进一步调整使得阻塞费用最小。阻塞费用的近似等价转化根据阻塞费用的计算规则，我们建立了阻塞费用关于各机组出力的非线性方程，在Lingo下编程得到的解为局部最优解，且十分不稳定，故考虑对其进行近似等价，使之转化为线性规划。算法如下：a. 计算时，根据规则每台机组应当以取得发电权的各段序外容量的最终报价与清算价的差值进行补偿，现在调整为各段分别按对应的报价与清算价的差值进行补偿；b. 的计算方法不变；c. 通过转化后的函数得到结果，再代回原规则计算阻塞费用。在Lingo下编写程序进行计算（源程序见附件4），得到较优的调整方案为150.48822875.61529576.4117相应的各线路潮流值为165149.4-154.98127.26132160.28此时目标结果为： ， 调整结果分析： 当预报负荷需求为982.4时，可以消除输电阻塞，阻塞费用为3264元。2预报负荷需求为1052.8时1） 出力分配预案的确定调用附件3的源程序，输入预报负荷需求，可得分配方案为：15081218.299.5135150102.1117此方案的清算价为，购电成本= （元）2） 潮流值的计算将分配预案代入有功潮流表达式，得到各线路潮流值为：177.24141.18-156.15129.73134.81167.06此时，线路1，5，6的潮流值均超过其限值，造成输电阻塞。3） 阻塞管理模型调整的结果同理，调用源程序（源程序见附件5），求解得调整后的方案为：1538822899.515215560.3117相应的各线路潮流值为：173.41143.6-155.2124.67135.38160.41此时目标结果为：， 调整结果分析：当预案需求为1052.8时，无法消除输电阻塞，但可以使用安全裕度输电，调整前方案的，调整后，降低了潮流超过限值的百分比，使电网运行更加安全，此时阻塞费用为1437.5七、结果分析1、对有功潮流近似表达式的分析为了检验拟合质量，我们计算出各线路拟合值和33组已知值绝对误差的平均值其中为实验数据，为拟合数据，计算得到=（0.0260 0.0233 0.0251 0.0242 0.0274 0.0273）结果表明，6条线路潮流的拟合值和实际值比较，几乎没有什么变化。下图为第5条线路拟合值和实际值的离散图。其中圆圈表示拟合值，实点表示实际值。从图中可见用线性拟合的效果非常好。2、对负荷需求的限制由于线路潮流值的约束，对负荷需求有一定的限制：1）不出现输电阻塞 此时线路上的潮流值不能超过其限值，以负荷最大最为规划目标，在Lingo下编程求解（源程序见附件6），得到该条件下的最大负荷能力为983.483MW，即当负荷需求大于983.483MW时，将无法调整方案使输电阻塞消失；2）不拉闸限电 此时线路上的潮流值不能超过其安全裕度，在同样的目标函数下对约束条件进行修正，运行程序得到结果最大负荷为1249.24MW，即当负荷需求大于此值时，必须采取拉闸限电措施。3、对线路改进的建议随着用电需求的增加，当用电负荷大于1249.24MW时，现行的电网结构已明显不能保证电网安全运行，因此，当用电需求较高时必须对线路进行整改。我们考虑不拉闸限电时最大负荷受各线路潮流限值的影响，当把各线路安全裕度下的潮流上限适当提高时，观察最大负荷的变化：图4 最大负荷随各线路安全裕度下潮流上限的变化趋势其中横坐标表示各线路安全裕度下的潮流上限的提高量，纵坐标表示最大负荷的变化。从图像中可以发现，线路1对最大负荷的影响比较明显，而线路25对最大负荷几乎没有影响。因此，在对线路进行整改时通过提高线路1的潮流安全裕度可以有效地提高最大负荷能力，使电网在高负荷前提下能够安全运行。八、模型的评价与改进1、我们的模型采用了“两步走”的策略，将双目标规划转化为两次单目标规划问题，大大降低了模型求解的难度，减少了程序运行时间；2、 引入危险因子反映电网线路的安全性，具有一定的实际意义。设计阻塞费用计算规则时我们仅考虑了序内容量和序外容量的损失，模型进一步改进的方案是在设计阻塞费用时充分权衡发电商和网方双方的利益，以达到真正意义上的的公平。九、参考文献1尚金成等，电力市场理论研究与应用，北京：中国电力出版社，2002；2黄继明，美国PJT电力市场，2004.09.18；十、附件清单附件1：模型基本数据附件2：有功潮流关于各机组出力近似表达式的计算源程序附件3：在Matlab下确定分配预案的源程序附件4：在Lingo下求解调整后的出力分配方案源程序（预报负荷为982.4MW）附件5：在Lingo下求解调整后的出力分配方案源程序（预报负荷为1052.8MW）附件6：在Lingo下求解不出现输电阻塞时的最大负荷附 件附件1： 模型基本数据 表1 各机组出力方案 （单位：兆瓦，记作MW）方案机组123 4 5 6780120731808012512581.1901133.02731808012512581.1902129.63731808012512581.1903158.77731808012512581.1904145.32731808012512581.190512078.5961808012512581.190612075.451808012512581.190712090.4871808012512581.190812083.8481808012512581.190912073231.398012512581.1901012073198.488012512581.1901112073212.648012512581.1901212073190.558012512581.190131207318075.85712512581.190141207318065.95812512581.190151207318087.25812512581.190161207318097.82412512581.190171207318080150.7112581.190181207318080141.5812581.190191207318080132.3712581.190201207318080156.9312581.190211207318080125138.8881.190221207318080125131.2181.190231207318080125141.7181.190241207318080125149.2981.19025120731808012512560.5829026120731808012512570.9629027120731808012512564.8549028120731808012512575.5299029120731808012512581.1104.8430120731808012512581.1111.2231120731808012512581.198.09232120731808012512581.1120.44表2 各线路的潮流值（各方案与表1相对应，单位：MW）方案线路1234560164.78140.87-144.25119.09135.44157.691165.81140.13-145.14118.63135.37160.762165.51140.25-144.92118.7135.33159.983167.93138.71-146.91117.72135.41166.814166.79139.45-145.92118.13135.41163.645164.94141.5-143.84118.43136.72157.226164.8141.13-144.07118.82136.02157.57165.59143.03-143.16117.24139.66156.598165.21142.28-143.49117.96137.98156.969167.43140.82-152.26129.58132.04153.610165.71140.82-147.08122.85134.21156.2311166.45140.82-149.33125.75133.28155.0912165.23140.85-145.82121.16134.75156.7713164.23140.73-144.18119.12135.57157.214163.04140.34-144.03119.31135.97156.3115165.54141.1-144.32118.84135.06158.2616166.88141.4-144.34118.67134.67159.2817164.07143.03-140.97118.75133.75158.8318164.27142.29-142.15118.85134.27158.3719164.57141.44-143.3119134.88158.0120163.89143.61-140.25118.64133.28159.1221166.35139.29-144.2119.1136.33157.5922165.54140.14-144.19119.09135.81157.6723166.75138.95-144.17119.15136.55157.5924167.69138.07-144.14119.19137.11157.6525162.21141.21-144.13116.03135.5154.2626163.54141-144.16117.56135.44155.9327162.7141.14-144.21116.74135.4154.8828164.06140.94-144.18118.24135.4156.6829164.66142.27-147.2120.21135.28157.6530164.7142.94-148.45120.68135.16157.6331164.67141.56-145.88119.68135.29157.6132164.69143.84-150.34121.34135.12157.64表3 各机组的段容量 （单位：MW）机组段1234567891017005000300004023002081562008311004003002040040455510101010150015755150151501010106950102001510200107501551510105103287002002002010155表4 各机组的段价（单位：元/兆瓦小时，记作元/MWh）机组段123456789101-50501241682102523123303634892-56001822032453003203604104953-61001521892332583083564155004-5001501702002553023253804358005-59001161461882152503103965106-60701591732052523053804055207-5001201802512603063153353485488-800153183233253283303318400800表5 各机组的爬坡速率 （单位：MW/分钟）机组12345678速率2.213.21.31.821.41.8表6 各线路的潮流限值（单位：MW）和相对安全裕度 线路123456限值165150160155132162安全裕度13%18%9%11%15%14%附件2： 有功潮流关于各机组出力近似表达式的计算源程序function jinsi %根据33组方案确定出有功潮流关于各机组出力的近似表达式 biao1data=120 73180 8012512581.1 90; %初始化机组出力方案矩阵 133.0273180 8012512581.1 90; 129.6373180 8012512581.1 90; 158.7773180 8012512581.1 90; 145.3273180 8012512581.1 90; 120 78.5961808012512581.1 90; 120 75.451808012512581.1 90; 120 90.4871808012512581.1 90; 120 83.8481808012512581.1 90; 120 73231.398012512581.1 90;120 73198.488012512581.1 90;12073212.648012512581.1 90;12073190.558012512581.1 90;1207318075.85712512581.1 90;1207318065.95812512581.1 90;1207318087.25812512581.1 90;1207318097.82412512581.1 90;1207318080150.7112581.1 90;1207318080141.5812581.1 90;1207318080132.3712581.1 90;1207318080156.9312581.1 90;1207318080125138.8881.1 90;1207318080125131.2181.1 90;1207318080125141.7181.1 90;1207318080125149.2981.1 90;120731808012512560.58290;120731808012512570.96290;120731808012512564.85490;120731808012512575.52990;120731808012512581.1 104.84;120731808012512581.1 111.22;120731808012512581.1 98.092; 120731808012512581.1 120.44;biao2data=164.78140.87 -144.25119.09 135.44157.69; %初始化潮流值矩阵165.81140.13-145.14118.63135.37160.76;165.51140.25-144.92118.7135.33159.98;167.93138.71-146.91117.72135.41166.81;166.79139.45-145.92118.13135.41163.64;164.94141.5-143.84118.43136.72157.22;164.8141.13-144.07118.82136.02157.5;165.59143.03-143.16117.24139.66156.59;165.21142.28-143.49117.96137.98156.96;167.43140.82-152.26129.58132.04153.6;165.71140.82-147.08122.85134.21156.23;166.45140.82-149.33125.75133.28155.09;165.23140.85-145.82121.16134.75156.77;164.23140.73-144.18119.12135.57157.2;163.04140.34-144.03119.31135.97156.31;165.54141.1-144.32118.84135.06158.26;166.88141.4-144.34118.67134.67159.28;164.07143.03-140.97118.75133.75158.83;164.27142.29-142.15118.85134.27158.37;164.57141.44-143.3119 134.88158.01;163.89143.61-140.25118.64133.28159.12;166.35139.29-144.2119.1136.33157.59;165.54140.14-144.19119.09135.81157.67;166.75138.95-144.17119.15136.55157.59;167.69138.07-144.14119.19137.11157.65;162.21141.21-144.13116.03135.5154.26;163.54141 -144.16117.56135.44155.93;162.7141.14-144.21116.74135.4154.88;164.06140.94-144.18118.24135.4156.68;164.66142.27-147.2120.21135.28157.65;164.7142.94-148.45120.68135.16157.63;164.67141.56-145.88119.68135.29157.61;164.69143.84-150.34121.34135.12157.64;X= biao1data ones(33,1); %确定变量矩阵for i=1:6 %循环给出各线路的关于各机组处理方案的各项系数y=biao2data(:,i);a=Xy %计算并输出各项系数end附件3： 在Matlab下确定分配预案的源程序function answer %计算清算价以及各机组出力方案函数biao3data=70 050 0 0 3000040; %初始化各机组段容量矩阵 30 020 8 15 62008; 110 040 0 30 02040040; 55 510 10 10 1015001; 75 515 0 15 150101010; 95 010 20 0 151020010; 50 155 15 10 1051032; 70 020 0 20 02010155; biao4data=-5050 124168210252312330363489; %初始化各机组段价矩阵 -5600 182203245300320360410495; -6100 152189233258308356415500; -500150170200255302325380435800; -5900 116146188215250310396510; -6070 159173205252305380405520; -500120180251260306315335348548; -800153183233253283303318400800; x1=0;x2=0;x3=0;x4=0;x5=0;x6=0;x7=0;x8=0; %初始化x1-x8，x1-x8为各个机组的出力累计值y=0; %初始化y y为各个机组累计值总和c=input(请输入预报负荷c=) %c为预报负荷while (y153 %由于爬坡速率的限制，当出力累计达到上限时不再增加 x1=153; end y=x1+x2+x3+x4+x5+x6+x7+x8; %计算出各个机组的出力累计总和 if yc %如果累计总和大于预报负荷，x1取部分当前段容量 x1=x1-(y-c); end case 2, x2=x2+biao3data(u,v); if x288 %由于爬坡速率的限制，当出力累计达到上限时不再增加 x2=88; end y=x1+x2+x3+x4+x5+x6+x7+x8; %计算出各个机组的出力累计总和 if yc %如果累计总和大于预报负荷，x2取部分当前段容量 x2=x2-(y-c); end case 3, x3=x3+biao3data(u,v); if x3228 %由于爬坡速率的限制，当出力累计达到上限时不再增加 x3=228; end y=x1+x2+x3+x4+x5+x6+x7+x8; %计算出各个机组的出力累计总和 if yc %如果累计总和大于预报负荷，x3取部分当前段容量 x3=x3-(y-c); endcase 4, x4=x4+biao3data(u,v); if x499.5 %由于爬坡速率的限制，当出力累计达到上限时不再增加 x4=99.5; end y=x1+x2+x3+x4+x5+x6+x7+x8; %计算出各个机组的出力累计总和 if yc %如果累计总和大于预报负荷，x4取部分当前段容量 x4=x4-(y-c); endcase 5, x5=x5+biao3data(u,v); if x5152 %由于爬坡速率的限制，当出力累计达到上限时不再增加 x5=152; end y=x1+x2+x3+x4+x5+x6+x7+x8; %计算出各个机组的出力累计总和 if yc %如果累计总和大于预报负荷，x5取部分当前段容量 x5=x5-(y-c); end case 6, x6=x6+biao3data(u,v); if x6155 %由于爬坡速率的限制，当出力累计达到上限时不再增加 x6=155; end y=x1+x2+x3+x4+x5+x6+x7+x8; %计算出各个机组的出力累计总和 if yc %如果累计总和大于预报负荷，x6取部分当前段容量 x6=x6-(y-c); end case 7, x7=x7+biao3data(u,v); if x7102.1 %由于爬坡速率的限制，当出力累计达到上限时不再增加 x7=102.1; end y=x1+x2+x3+x4+x5+x6+x7+x8; %计算出各个机组的出力累计总和 if yc %如果累计总和大于预报负荷，x7取部分当前段容量 x7=x7-(y-c); end case 8, x8=x8+biao3data(u,v); if x8117 %由于爬坡速率的限制，当出力累计达到上限时不再增加 x8=117; end y=x1+x2+x3+x4+x5+x6+x7+x8; %计算出各个机组的出力累计总和 if yc