发动机潮流值与出力值.docx

1、发动机潮流值与出力值探究发动机机组的线路潮流值与机组出力关系1.摘要我国电力系统市场化改革正在积极进行，店里市场输电管理显得尤为重要。而准确找出潮流值与出力之间的关系是解决输电效率的关键一步。对于问题一，根据表内所给数据，适当取出数据，利用控制变量法以及excel画图，做出散点图，发现各线路潮流值关于发电机组出力是近似线性关系，并且带有常数项。对于问题二，在解决问题一的基础上，设出线路潮流值关于机组出力的表达式，利用matlab软件中回归分析命令求解出各项机组系数以及常数项。并且带入原数据进行检验。关键词：控制变量 近似线性关系 回归分析2.问题重述某电网有8台发电机组，6条主要线路，表1和表

2、2中的方案0给出了各机组的当前出力和各线路上对应的有功潮流值，方案132给出了围绕方案0的一些实验数据，本文将研究下列问题：（1）从所给数据出发分析6条线路上的潮流值与各机组出力方案的可能关系；（2）由这些数据确定各线路上有功潮流关于各发电机组出力的近似表达式。 3.问题假设和变量说明3.1 问题假设：1. 假设在发电机组工作时，没有线路的损耗，并且所给数据正确没有较大误差。3.2变量说明：第i条线路的潮流值：第i条线路第j各发电机组的出力值：第i条线路第j各发电机组的出力值的系数：常数项4.模型建立及求解 . 4.1问题一：根据所给数据，取出一组出力值(除第i个发电机组的出力值不同其他相同)

3、和相应的线路潮流值。例如，165.81165.51167.93166.79133.02129.63158.77145.32即方案1,2,3,4的第一条线路的潮流值和第一条第一个发电机组的出力值将数据在excel中画出散点图，并画出趋势线和r值如下 值很接近于1，其他发电机出力与潮流值都类似，并且很接近于1.所以得出结论，6条线路上的潮流值与各机组出力方案的可能关系是近似线性关系。4.2问题二： 在问题一中得到了潮流值与各机组出力方案的可能关系是近似线性关系，所以设出表达式，使用matlab软件作回归分析，得出各个系数以及区间估计，残差，置信区间和检验回归数据的统计量。（代码详见附录6.2）线路

4、一区间估计为minmaxB1109.3757111.2173A110.081090.084593A120.0443180.05224A130.0516420.0543A140.1168410.123026A15-0.02737-0.02351A160.1193940.124628A170.1185450.124606A18-0.003350.0009相关系数=0.9995 线路二：区间估计为minmaxB2130.4149132.0429A21-0.05611-0.05301A220.1243490.131352A23-0.00120.001149A240.0305430.03601A250.

5、085140.088553A26-0.11475-0.11013A27-0.02161-0.01625A280.0968480.100605相关系数=0.9996 线路三：区间估计为minmaxB3-109.713-108.033A31-0.07114-0.06794A320.0580320.065258A33-0.15783-0.1554A34-0.01274-0.0071A350.1227340.126255A36-0.000270.004504A37-0.005280.000253A38-0.20333-0.19945相关系数=0.9998 线路四：区间估计为minmaxB476.673

6、3578.29002A41-0.036-0.03293A42-0.10589-0.09893A430.2039980.206331A44-0.02354-0.01811A45-0.01352-0.01013A460.0036560.00825A470.1422580.147579A480.074680.078411相关系数=0.9999线路五：minmaxB5132.1371133.8118A51-0.00110.0021A520.23970.2469A53-0.0658-0.0633A54-0.0439-0.0383A55-0.0670-0.0635A56 0.0680 0.0727A57-

7、0.0070-0.0015A58-0.0108-0.0070相关系数=0.9996线路六：minmaxB6119.7094121.6171A610.23600.2396A62-0.0643-0.0561A63-0.0793 -0.0765A640.08980.0962A650.04490.0489A66 -0.0026 0.0028A670.16280.1691A68-0.00150.0029相关系数=0.99985.检验分析 因为方案1-32是围绕方案0 给出的数据，所以将方案0的数据代入所求解的表达式中，求出误差。 线路一线路二线路三线路四线路五线路六计算值164.712140.8238-

8、144.205119.0412135.3803157.6206实际值164.78140.87-144.25119.09135.44157.69误差（%）0.0412670.0327960.0311270.0409770.0440790.04401误差均在0.1%以下。 6.模型推广 本体所涉及的线性回归模型，它是用试验值（样本值）对未知参数和作点估计和假设检验，从而建立y与之间的数量关系，它可以在处对y的值作预测与控制，即对y作区间估计，具有一定的预测估计作用。7.附录6.1题中所给数据表1 各机组出力方案 （单位：兆瓦，记作MW）方案机组123456780120731808012512581

9、.1901133.02731808012512581.1902129.63731808012512581.1903158.77731808012512581.1904145.32731808012512581.190512078.5961808012512581.190612075.451808012512581.190712090.4871808012512581.190812083.8481808012512581.190912073231.398012512581.1901012073198.488012512581.1901112073212.648012512581.19012120

10、73190.558012512581.190131207318075.85712512581.190141207318065.95812512581.190151207318087.25812512581.190161207318097.82412512581.190171207318080150.7112581.190181207318080141.5812581.190191207318080132.3712581.190201207318080156.9312581.190211207318080125138.8881.190221207318080125131.2181.1902312

11、07318080125141.7181.190241207318080125149.2981.19025120731808012512560.5829026120731808012512570.9629027120731808012512564.8549028120731808012512575.5299029120731808012512581.1104.8430120731808012512581.1111.2231120731808012512581.198.09232120731808012512581.1120.44表2 各线路的潮流值（各方案与表1相对应，单位：MW）方案线路123

12、4560164.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

13、.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.8413

14、5.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

15、.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附录6.2在matlab中新建变量x，y将表一数据复制到x中，分别把表二方案1-32各线路的潮流值复制到y中【a，b，c，d，e】=regress（y，x）