电力市场的输电阻塞管理—全国数学建模大赛

《电力市场的输电阻塞管理—全国数学建模大赛》由会员分享，可在线阅读，更多相关《电力市场的输电阻塞管理—全国数学建模大赛（42页珍藏版）》请在装配图网上搜索。

1、电力市场的输电阻塞管理摘要针对问题一，即对各线路上有功潮流关于各发电机组出力的研究，我们通过线性回归分析的思想，建立了的数学模型。然后利用MATLAB得到具体参数，从而确定各线路上有功潮流关于各发电机组出力的近似表达式。通过对置信区间的分析，检查回归变量对因变量的影响的显著程度。针对问题二，在输电阻塞发生时，阻塞费用计算应包括两部分：一、通过竞价取得发电权的发电容量（称序内容量）不能出力所引起的发电的成本提高和利润减少的部分；二、因一些在竞价中未取得发电权的发电容量（称序外容量）要在低于对应报价的清算价上出力所引起的额外出力的费用。为了公平起见，第一部分的阻塞费用不应由网方全部承担，具体折扣由

2、买卖双方协定，此处我们假定为50%，而第二部分的阻塞费用由网方全部承担。针对问题三，根据题目给出的电力市场交易原则，下一个时段预报的负荷需求和各机组的段容量、段价和爬坡速率的数据，可以得到下一个时段各机组的出力分配预案。针对问题四，按照题目给出的潮流限值，有问题一得到的各线路上有功潮流关于各发电机组出力的近似表达式可以检测问题三的出力分配预案是否会引起输电阻塞；在检测到发生输电阻塞时，根据各机组出力对各线路潮流值的影响比例大小，调整各机组出力分配方案，并给出与该方案相应的阻塞费用。针对问题五，根据对问题三和问题四的分析方法，可以得到相应结论。关键词：回归分析 出力分配预案 阻塞费用一、问题重述

3、设某电网有若干台发电机组和若干条主要线路，每条线路上的有功潮流（输电功率和方向）取决于电网结构和各发电机组的出力。电网每条线路上的有功潮流的绝对值有一安全限值，限值还具有一定的相对安全裕度（即在应急情况下潮流绝对值可以超过限值的百分比的上限）。如果各机组出力分配方案使某条线路上的有功潮流的绝对值超出限值，称为输电阻塞。当发生输电阻塞时，需要研究如何制订既安全又经济的调度计划。l 电力市场交易规则：1. 以15分钟为一个时段组织交易，每台机组在当前时段开始时刻前给出下一个时段的报价。各机组将可用出力由低到高分成至多10段报价，每个段的长度称为段容量，每个段容量报一个价（称为段价），段价按段序数单

4、调不减。在最低技术出力以下的报价一般为负值，表示愿意付费维持发电以避免停机带来更大的损失。2. 在当前时段内，市场交易-调度中心根据下一个时段的负荷预报,每台机组的报价、当前出力和出力改变速率，按段价从低到高选取各机组的段容量或其部分（见下面注释），直到它们之和等于预报的负荷，这时每个机组被选入的段容量或其部分之和形成该时段该机组的出力分配预案（初始交易结果）。最后一个被选入的段价（最高段价）称为该时段的清算价，该时段全部机组的所有出力均按清算价结算。注释：（a） 每个时段的负荷预报和机组出力分配计划的参照时刻均为该时段结束时刻。（b） 机组当前出力是对机组在当前时段结束时刻实际出力的预测值。

5、（c） 假设每台机组单位时间内能增加或减少的出力相同，该出力值称为该机组的爬坡速率。由于机组爬坡速率的约束，可能导致选取它的某个段容量的部分。（d） 为了使得各机组计划出力之和等于预报的负荷需求，清算价对应的段容量可能只选取部分。市场交易-调度中心在当前时段内要完成的具体操作过程如下：1、 监控当前时段各机组出力分配方案的执行，调度AGC辅助服务，在此基础上给出各机组的当前出力值。2、 作出下一个时段的负荷需求预报。3、 根据电力市场交易规则得到下一个时段各机组出力分配预案。4、 计算当执行各机组出力分配预案时电网各主要线路上的有功潮流，判断是否会出现输电阻塞。如果不出现，接受各机组出力分配预

6、案；否则，按照如下原则实施阻塞管理：l 输电阻塞管理原则：（1） 调整各机组出力分配方案使得输电阻塞消除。（2） 如果（1）做不到，还可以使用线路的安全裕度输电，以避免拉闸限电（强制减少负荷需求），但要使每条线路上潮流的绝对值超过限值的百分比尽量小。（3） 如果无论怎样分配机组出力都无法使每条线路上的潮流绝对值超过限值的百分比小于相对安全裕度，则必须在用电侧拉闸限电。（4） 当改变根据电力市场交易规则得到的各机组出力分配预案时，一些通过竞价取得发电权的发电容量（称序内容量）不能出力；而一些在竞价中未取得发电权的发电容量（称序外容量）要在低于对应报价的清算价上出力。因此，发电商和网方将产生经济利

7、益冲突。网方应该为因输电阻塞而不能执行初始交易结果付出代价，网方在结算时应该适当地给发电商以经济补偿，由此引起的费用称之为阻塞费用。网方在电网安全运行的保证下应当同时考虑尽量减少阻塞费用。你需要做的工作如下：1. 某电网有8台发电机组，6条主要线路，表1和表2中的方案0给出了各机组的当前出力和各线路上对应的有功潮流值，方案132给出了围绕方案0的一些实验数据，试用这些数据确定各线路上有功潮流关于各发电机组出力的近似表达式。2. 设计一种简明、合理的阻塞费用计算规则，除考虑上述电力市场规则外，还需注意：在输电阻塞发生时公平地对待序内容量不能出力的部分和报价高于清算价的序外容量出力的部分。 3.

8、假设下一个时段预报的负荷需求是982.4MW，表3、表4和表5分别给出了各机组的段容量、段价和爬坡速率的数据，试按照电力市场规则给出下一个时段各机组的出力分配预案。4. 按照表6给出的潮流限值，检查得到的出力分配预案是否会引起输电阻塞，并在发生输电阻塞时，根据安全且经济的原则，调整各机组出力分配方案，并给出与该方案相应的阻塞费用。5. 假设下一个时段预报的负荷需求是1052.8MW，重复34的工作。二、问题分析我们进行了如下的分析：Step1: 首先我们明确问题中所需做的几个任务：1. 各线路上有功潮流关于各发电机组出力的近似表达式；2. 一种简明、合理的阻塞费用计算规则；3. 假设下一个时段

9、预报的负荷需求是982.4MW，给出下一个时段各机组的出力分配预案；4. 检测并调整分配预案，并计算相应的阻塞费用。Step2: 即对各线路上有功潮流关于各发电机组出力的研究，我们通过线性回归分析的思想，建立了公式的数学模型。然后利用MATLAB得到具体参数，从而确定各线路上有功潮流关于各发电机组出力的近似表达式。通过对置信区间的分析，检查回归变量对因变量的影响的显著程度。Step3: 在输电阻塞发生时，阻塞费用计算应包括两部分：一、通过竞价取得发电权的发电容量（称序内容量）不能出力所引起的发电的成本提高和利润减少的部分；二、因一些在竞价中未取得发电权的发电容量（称序外容量）要在低于对应报价的

10、清算价上出力所引起的额外出力的费用。为了公平起见，第一部分的阻塞费用不应由网方全部承担，具体折扣由买卖双方协定，此处我们假定为50%，而第二部分的阻塞费用由网方全部承担。Step4: 根据题目给出的电力市场交易原则，下一个时段预报的负荷需求和各机组的段容量、段价和爬坡速率的数据，可以得到下一个时段各机组的出力分配预案。Step5: 按照题目给出的潮流限值，有问题一得到的各线路上有功潮流关于各发电机组出力的近似表达式可以检测问题三的出力分配预案是否会引起输电阻塞；在检测到发生输电阻塞时，根据各机组出力对各线路潮流值的影响比例大小，调整各机组出力分配方案，并给出与该方案相应的阻塞费用。Step6:

11、 根据对问题三和问题四的分析方法，可以得到相应结论。三、符号设定 第i条线路的有功潮流 第i个机组的出力 各部分阻塞费用的总和发电厂每度电的变动成本发电厂发电的固定成本给定的各机组的段容量出力预案中具体的各机组中的段容量给定的各机组的清算价出力预案中具体的各机组中的清算价四、模型建立与求解(一)、问题一的求解5.1 记第i条线路的有功潮流为，第i个机组的出力为，根据题目给出的32种方案的各机组的出力值及相对应的各线路的潮流值，基于上面的分析，我们首先利用来建立的近似表达式。模型准备为了大致分析与的关系，首先利用表1和表2的数据可以做出对的散点图，依此类推可分别得到对的散点图。图1 对的散点图模

12、型建立从图1可以发现，随着的增加，的值有比较明显的线性增长趋势，图中的直线是用线性模型 （5.1.1）拟合的，（其中是误差）同样经过分析，发现随着的增加，的值有比较明显的线性增长趋势，图中的直线是用线性模型 （5.1.2）综合上面的分析，结合上述模型简历如下的回归模型（5.1.3）（5.1.3）右端的称为回归变量（自变量），是给定机组出力，第一条线路的潮流值的平均值，其中的参数称为回归系数，由表1、表2的数据估计，影响的其他因素作用都包含在随机误差中。如果模型选择得合适，应大致服从均值为零的正态分布。依此类推，我们可以用同样的方法得出分别对的基本模型 （5.1.4）模型求解直接利用MATLAB

13、统计工具箱中的命令regress求解，使用格式为： b,bint,r,rint,stats=regress(y1(:),x0(:) x1(:) x2(:) x3(:) x4(:) x5(:) x6(:) x7(:) x8(:)其中输入y为模型（5.1.3）中y的数据，x为对应于回归系数的数据矩阵，置信水平。得到模型（5.1.3）的回归系数估计值及其置信区间（置信水平）的结果见下表（5.1）。参数参数值估计参数置信区间0110.4775 109.5421 111.412910.0826 0.0808 0.084420.0478 0.0437 0.051830.0528 0.0514 0.0542

14、40.1199 0.1166 0.12315-0.0257 -0.0277 -0.02376 0.12160.1190 0.12437 0.12200.1189 0.12518 -0.0015-0.0037 0.0007表（5.1） 模型（5.1.3）的计算结果同样，以此类推我们可以得出其他各线路上有功潮流关于各发电机组出力的近似表达式： （5.1.5）yi: =0 1 2 3 4 5 6 7 8其中，y1: =110.4775 0.0826 0.0478 0.0528 0.1199 -0.0257 0.1216 0.1220 -0.0015y2: =131.3521 0.0547 0.127

15、5 -0.0001 0.0332 0.0867 -0.1127 -0.0186 0.0985y3: =-108.99281 -0.0694 0.0620 -0.1565 -0.0099 0.01247 0.0024 -0.0028 -0.2012y4: =77.6116 -0.0346 -0.1028 0.2050 -0.0209 -0.0120 0.0057 0.1452 0.0763y5: =133.1334 -0.0003 0.2428 -0.0647 -0.0412 -0.0655 0.0700 -0.0039 -0.0092y6: =120.8481 0.2376 -0.0607 -

16、0.0781 -0.0929 0.0466 -0.0003 0.1664 0.0004结果分析表（5.1）的回归系数给出了模型（5.1.3）中的估计值，检查他们的置信区间发现，只有8置信区间包含零点，表示回归变量8（对因变量的影响）不是太显著，但我们仍将8保留在模型中。(二)、问题二的求解5.2 在输电阻塞发生时，阻塞费用（记为W）计算应包括两部分：一、通过竞价取得发电权的发电容量（称序内容量）不能出力所引起的发电的成本提高（记为W1）和利润减少的部分（记为W2）；二、因一些在竞价中未取得发电权的发电容量（称序外容量）要在低于对应报价的清算价上出力所引起的出力的额外费用（记为W3）。记给定的各

17、机组的段容量为，记出力预案中具体的各机组中的段容量为，记给定的各机组的清算价为，记出力预案中具体的各机组中的清算价为。模型准备为了公平起见，第一部分的阻塞费用不应由网方全部承担，具体折扣由买卖双方协定，此处我们假定为50%，而第二部分的阻塞费用由网方全部承担。即 （5.2）模型建立5.2.1 计算发电的成本提高W1发电厂的发电成本分为两部分：固定成本：包括厂房设备、保险、材料存贮、设备运行等，不随发电 度数的变化而变化；平均变动成本：包括原材料、燃料、工资等，随发电度数的变化而变化。所以由给定的各机组的段容量为，出力预案中具体的各机组中的段容量为，发电的成本提高W1的值确定为 （5.2.1）5

18、.2.2 计算发电的利润减少W2由给定的各机组的段容量为，出力预案中实际的机组中的段容量为，故少发电为则发电的利润减少W2的值确定为 （5.2.2）5.2.3 计算出力的额外费用W3由给定的各机组的清算价为，出力预案中具体的各机组中的清算价为，故多出的费用为则发电的出力的额外费用W3的值确定为 （5.2.3）5.2.4 计算阻塞费用W由（5.2）可得阻塞费用的值为 （5.2.4）(三)、问题三的求解5.3 已知下一阶段的负荷量，我们根据电力市场交易规则给出相应的出力分配预案。按段价从低到高选取各机组的段容量或其部分，直到它们之和等于预报的负荷，这时每个机组被选入的段容量或其部分之和形成该时段该

19、机组的出力分配预案。以下表格中的红色部分就是按照电力市场交易规则所取的各机组的段容量对应的段价。机组段123456789101-50501241682102523123303634892-56001822032453003203604104953-61001521892332583083564155004-5001501702002553023253804358005-59001161461882152503103965106-60701591732052523053804055207-5001201802512603063153353485488-8001531832332532833033

20、18400800因考虑到各机组的最后出力总和应该在最低出力和到最高出力这个范围内，以及最后的总负荷量的大小，所以可能有些段容量只能取其中的一部分。以下就是按照上面的两个约束条件作出调整后的结果。因为第四机组的最高出力只能是99.5，所以第四机组的第六段只能出力9.5，这样才能保证不超过最高出力。当清算价为302（第四机组第六段）时，段容量之和为978.5，离目标负荷只差3.9，所以在下一段段价为303（第八机组第七段）只要取3.9就能达到目标值了。机组12345678总出力当前出力120731808012512581.190最高出力1538822899.5152155102.1117最低出力8

21、75813260.5989560.163第一出力70301105575955070555第二出力/段5/215/320第三出力/段15/215第四出力/段50/340/35/295第五出力/段20/310/315/510/35/320/311010/420/4第七出力/段30/68/415/530/510/415/615/615/420/5158出力和一150731809012514085110953第九出力/段6/69.5/63.9/710/529.4出力和1507918099.5125143995110982.4所以下一阶段出力分配预案为：机组12345678出力1507918099.51

22、25143995110清算价为：303(4) 、问题四的求解5.4.1判断是否有输电阻塞现象 利用问题一中求出的各线路上有功潮流关于各发电机组出力的近似表达式，将各机组出力的具体数据代入计算，得出各线路上的潮流值。然后与题目所给的各线路潮流限值相比较，观察是否有输电阻塞现象。计算线路一上的有功潮流：y1=110.4775+0.0826*150+0.0478*79+0.0528*180+0.1199*99.5-0.0257*125+0.1216*143.9+0.1220*95-0.0015*110y1 =173.7885 165则线路一上有阻塞现象。依次我们可以得出y2=131.3521-0.0

23、547*150+0.1275*79-0.0001*180+0.0332*99.5+0.0867*125-0.1127*143.9-0.0186*95+0.0985*110y2=140.1930150所以线路二上无阻塞现象y3=-108.99281-0.0694*150+0.0620*79-0.1565*180-0.0099*99.5+0.01247*125+0.0024*143.9-0.0028*95-0.2012*110y3 =-164.1538 160所以线路三上无阻塞现象y4=77.6116-0.0346*150-0.1028*79+0.2050*180-0.0209*99.5-0.01

24、20*125+0.0057*143.9+0.1452*95+0.0763*110y4 =120.6281132所以线路五上有阻塞现象y6=120.8481+0.2376*150-0.0607*79-0.0781*180-0.0929*99.5+0.0466*125-0.0003*143.9+0.1664*9+0.0004*110y6 =135.7147162所以线路六上无阻塞现象综上所述，线路一和线路五会出现输电阻塞现象5.4.2.1调整出力分配方案 在发生输电阻塞时，根据安全且经济的原则，调整各机组出力分配方案。首先我们根据各线路上有功潮流关于各发电机组出力的近似表达式得出各机组出力对线路潮

25、流值的影响大小，如下表（表中1代表影响最大，8代表影响最小）机组12345678线路146537218线路251864273线路345263871线路453167824线路581453276线路615436827因为我们调整出力分配方案的原则是做尽量少的调整而能使得输电阻塞消除，所以我们要进行调整的首要对象是对线路一和五影响较大而对其他线路影响较小的机组。其中线路一超出范围较大，所以优先考虑对线路一的调整。由表可看出机组6、7的出力大小对线路一的潮流值的影响最大。所以可以使机组6、7的出力减小来达到线路一的潮流值变小。而要保持最后的出力总和不变，则应有相应机组的出力增加。机组5、8的变化对线路

26、一潮流值的影响较小，并且是负相关的。所以我们可以相应的增加机组5,、8的出力来保持出力总和不变，同时还可以减小线路一的潮流值。在调整的过程中，应时刻注意各机组的出力总和应该在最小出力和最大出力范围内，并且相对考虑一下段价的大小，尽量避免段价过高，使得阻塞费用变大。在调整的过程中，我们发现调整后使得线路一的潮流值小于潮流限值，相应的线路五的潮流值也会小于潮流限值，反而线路二很容易就超出潮流限值的范围。所以接下来还需要做微小的调整，使得所有线路的潮流值都在潮流限值的范围内，即无输电阻塞。所以调整后的出力方案为：机组12345678调整前1507918099.512514095113.9调整后153

27、86.522890.81529560.1117最高出力1538822899.5152155102.1117最低出力875813260.5989560.1635.4.2.2 检验是否存在输电阻塞现象y1=110.4775+0.0826*153+0.0478*86.5+0.0528*228+0.1199*90.8-0.0257*152+0.1216*95+0.1220*60.1-0.0015*117y1 =164.9776165 y2=131.3521-0.0547*153+0.1275*86.5-0.0001*228+0.0332*90.8+0.0867*152-0.1127*95-0.0186

28、*60.1+0.0985*117y2 =149.8820150y3=-108.99281-0.0694*153+0.0620*86.5-0.1565*228-0.0099*90.8+0.01247*152+0.0024*95-0.0028*60.1-0.2012*117y3 = -172.4142160y4=77.6116-0.0346*153-0.1028*86.5+0.2050*228-0.0209*90.8-0.0120*152+0.0057*95+0.1452*60.1+0.0763*117y4 = 124.6390155y5=133.1334-0.0003*153+0.2428*86

29、.5-0.0647*228-0.0412*90.8-0.0655*152+0.0700*95-0.0039*60.1-0.0092*117y5 =130.9803132y6=120.8481+0.2376*153-0.0607*86.5-0.0781*228-0.0929*90.8+0.0466*152-0.0003*95+0.1664*60.1+0.0004*117y6 =142.8104165则线路一上阻塞现象依次我们可以得出y2=131.3521-0.0547*150+0.1275*81-0.0001*218.2+0.0332*99.5+0.0867*135-0.1127*150-0.0

30、186*102.1+0.0985*117y2 =141.1811 150所以线路二上无阻塞现象y3=-108.99281-0.0694*150+0.0620*81-0.1565*218.2-0.0099*99.5+0.01247*135+0.0024*150-0.0028*102.1-0.2012*117y3 =-171.2970 160所以线路三上无阻塞现象y4=77.6116-0.0346*150-0.1028*81+0.2050*218.2-0.0209*99.5-0.0120*135+0.0057*150+0.1452*102.1+0.0763*117y4 =129.7333 132所

31、以线路五上有阻塞现象y6=120.8481+0.2376*150-0.0607*81-0.0781*218.7-0.0929*99.5+0.0466*135-0.0003*150+0.1664*102.1+0.0004*117y6 =148.5687 165 y2=131.3521-0.0547*153+0.1275*88-0.0001*228+0.0332*99.5+0.0867*152-0.1127*113.2-0.0186*102.1+0.0985*117y2 =147.5298 150 y3=-108.99281-0.0694*153+0.0620*88-0.1565*228-0.00

32、99*99.5+0.01247*152+0.0024*113.2-0.0028*102.1-0.2012*117y3 =-172.4812 160 y4=77.6116-0.0346*153-0.1028*88+0.2050*228-0.0209*99.5-0.0120*152+0.0057*113.2+0.1452*102.13+0.0763*117y4 =130.5095 132y6=120.8481+0.2376*153-0.0607*88-0.0781*228-0.0929*99.5+0.0466*152-0.0003*113.2+0.1664*102.1+0.0004*117y6 =

33、148.8944 x0=1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1; x1=120 133.02 129.63 158.77 145.32 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120 120;x2=73 73 73 73 73 78.596 75.45 90.487 83.848 73 73 73 73 73 73 73 73 73 73

34、 73 73 73 73 73 73 73 73 73 73 73 73 73 73;x3=180 180 180 180 180 180 180 180 180 231.39 198.48 212.64 190.55 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180 180;x4=80 80 80 80 80 80 80 80 80 80 80 80 80 75.857 65.958 87.258 97.824 80 80 80 80 80 80 80 80 80 80 80 80 80 8

35、0 80 80;x5=125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 150.71 141.58 132.37 156.93 125 125 125 125 125 125 125 125 125 125 125 125;x6=125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 125 138.88 131.21 141.71 149.29 125 125 125 125 125 125 125 12

36、5;x7=81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 81.1 60.582 70.962 64.854 75.529 81.1 81.1 81.1 81.1;x8=90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 90 104.84 111.22 98.092 120.44; y1=

37、164.78 165.81 165.51 167.93 166.79 164.94 164.8 165.59 165.21 167.43 165.71 166.45 165.23 164.23 163.04 165.54 166.88 164.07 164.27 164.57 163.89 166.35 165.54 166.75 167.69 162.21 163.54 162.7 164.06 164.66 164.7 164.67 164.69; y2=140.87 140.13 140.25 138.71 139.45 141.5 141.13 143.03 142.28 140.82

38、 140.82 140.82 140.85 140.73 140.34 141.1 141.4 143.03 142.29 141.44 143.61 139.29 140.14 138.95 138.07 141.21 141 141.14 140.94 142.27 142.94 141.56 143.84; y3=-144.25 -145.14 -144.92 -146.91 -145.92 -143.84 -144.07 -143.16 -143.49 -152.26 -147.08 -149.33 -145.82 -144.18 -144.03 -144.32 -144.34 -14

39、0.97 -142.15 -143.3 -140.25 -144.2 -144.19 -144.17 -144.14 -144.13 -144.16 -144.21 -144.18 -147.2 -148.45 -145.88 -150.34; y4=119.09 118.63 118.7 117.72 118.13 118.43 118.82 117.24 117.96 129.58 122.85 125.75 121.16 119.12 119.31 118.84 118.67 118.75 118.85 119 118.64 119.1 119.09 119.15 119.19 116.

40、03 117.56 116.74 118.24 120.21 120.68 119.68 121.34; y5=135.44 135.37 135.33 135.41 135.41 136.72 136.02 139.66 137.98 132.04 134.21 133.28 134.75 135.57 135.97 135.06 134.67 133.75 134.27 134.88 133.28 136.33 135.81 136.55 137.11 135.5 135.44 135.4 135.4 135.28 135.16 135.29 135.12; y6=157.69 160.7

41、6 159.98 166.81 163.64 157.22 157.5 156.59 156.96 153.6 156.23 155.09 156.77 157.2 156.31 158.26 159.28 158.83 158.37 158.01 159.12 157.59 157.67 157.59 157.65 154.26 155.93 154.88 156.68 157.65 157.63 157.61 157.64; b,bint,r,rint,stats=regress(y1(:),x0(:) x1(:) x2(:) x3(:) x4(:) x5(:) x6(:) x7(:) x

42、8(:)b = 110.4775 0.0826 0.0478 0.0528 0.1199 -0.0257 0.1216 0.1220 -0.0015bint = 109.5421 111.4129 0.0808 0.0844 0.0437 0.0518 0.0514 0.0542 0.1166 0.1231 -0.0277 -0.0237 0.1190 0.1243 0.1189 0.1251 -0.0037 0.0007r = 0.0615 0.0159 -0.0040 0.0088 -0.0201 -0.0458 -0.0356 0.0362 -0.0267 -0.0016 0.0158 0.0083 -0.0