电力市场的输电阻塞管理资料bubh.docx

电力市场场的输电电阻塞管管理摘要本文研究究了电力力市场的的输电阻阻塞管理理问题。电网公司在组织交易、调度和配送时，遵循的是电网“安全第一，兼顾经济”的原则，制订的电力市场交易规则，是按照购电费用最小的经济目标来运作的。问题一，我们利用给出的数据，采用多元线性回归来求取有功潮流与各机组出力之间的关系。而且用题中所给的数据进行了其进行相应的误差分析与灵敏度分析，残差很小，并且分析了所得线性函数中的常数项存在的物理背景。问题二，网方对发电商的供电取舍有着宏观调控的能力，在竞价中可能会造成发电商的一定损失。网方必须对这部分损失给与发电商一定补偿，即阻塞费用。阻塞补偿偿=+问题三，其其目标是是总费用用最小。又又因为各各机组均均受爬坡坡速率的的限制，因因此我们们可以建建立一线线性规划划模型，运运用Maatlaab工具具箱计算算得到各各机组的的出力。问题四，根根据给出出各线路路的潮流流值，检检验是否否会产生生输电阻阻塞。若若产生，则则按照输输电阻塞塞管理原原则，依依次采取取调整预预案、裕裕度输电电、拉闸闸限电来来保证输输电线路路安全运运行。计计算结果果如下：负荷需求求为9882.44MW时时的出力力分配方方案：线路123456预案潮流流值（MMW）173.74140.97-1500.833119.27136.81168.57调整后潮潮流值（MMW）164.99149.42154.99126.27132.00159.57清算价，阻阻塞费用用为47749元元。问题五，采用问题三相同的算法思想，计算各机组出力。对阻塞方案进行调整，得不到可行解。采用裕度输电，保证每条线路上的潮流绝对值超过限定值的百分比尽可能小，然后考虑经济最优，建立多目标规划模型，计算各机组出力。最后检验各线路的潮流值均未超过限值并计算得到阻塞费用。负荷需求求为10052.8MWW时的出出力分配配方案：机组 12345678预案出力力（mww）15081218.299.55135150102.1117调整后出出力（mmw）1538822899.5515215560.33117线路123456调整后潮潮流（mmw）173.41143.58155.21124.68135.3160.42清算价，阻阻塞费用用为8889.555元。一、问题题重述电力从生生产到使使用的四四大环节节发电、输输电、配配电和用用电是瞬瞬间完成成的。我我国电力市市场初期期是发电电侧电力力市场，采采取交易易与调度一一体化的的模式。电电网公司司在组织织交易、调调度和配配送时，必必须遵循循电网“安全第第一”的原则则，同时时要制订订一个电电力市场场交易规规则，按按照购电电费用最最小的经经济目标标来运作作。市场场交易-调度中中心根据据负荷预预报和交交易规则则制订满满足电网网安全运运行的调调度计划划，即各各发电机机组的出出力（发发电功率率）分配配方案。设设某电网网有若干干台发电电机组和和若干条条主要线线路，每每条线路路上的有有功潮流流（输电电功率和和方向）取取决于电电网结构构和各发发电机组组的出力力。电网网每条线线路上的的有功潮潮流的绝绝对值有有一安全全限值，限限值还具具有一定定的相对对安全裕裕度（即即在应急急情况下下潮流绝绝对值可可以超过过限值的的百分比比的上限限）。如如果各机机组出力力分配方方案使某某条线路路上的有有功潮流流的绝对对值超出出限值，称称为输电电阻塞。当当发生输输电阻塞塞时，需需要研究究如何制制订既安安全又经经济的调调度计划划，使各各线路的的潮流值值的绝对对值不超超过相对对安全裕裕度且阻阻塞费用用最小。二、模型型假设1、忽略略电流在在网络传传输中的的损失，即即线路有有功潮流流等于机机组出力力之和；2、机组组的爬坡坡速率是是一个总总体指标标，各个个段其变变化的速速率与之之不等；3、每条条线路上上的有功功潮流只只取决于于电网结结构和各各发电机机组的出出力，即即输电线线路与发发电机组组间满足足固定的的函数关关系；4、每个个时段的的负荷预预报和机机组出力力分配计计划的参参照时刻刻均为该该时段结结束时刻刻。三、问题题分析电网公司司在组织织交易、调调度和配配送时，遵遵循的是是电网“安全第第一，兼兼顾经济济”的原则则，制订订的电力力市场交交易规则则，是按按照购电电费用最最小的经经济目标标来运作作的。由由于输电电线路传传输容量量的限制制，有功功潮流的的绝对值值不能超超过其安安全限值值，否则则将引起起输电阻阻塞，危危及电网网安全，此此时必须须对各机机组的出出力分配配方案做做调整，以以消除阻阻塞。这这样就会会使部分分序内容容量不能能出力，只只能用序序外容量量在被迫迫低于其其报价的的清算价价上出力力。为了了解决利利益冲突突，网方方要为因因输电阻阻塞而不不能严格格执行原原出力分分配方案案付出代代价，对对发电商商潜在的的和实际际的利益益损失给给予补偿偿阻塞塞费用。一：有功功潮流与与各机组组出力之之间的关关系虽然有功功潮流与与各级组组之间没没有明确确的关系系，但是是我们可可以利用用给出的的大量数数据，基基于对数数据的统统计分析析去建立立模型。通通过对数数据的分分析，我我们采用用多元线线性回归归来求取取有功潮潮流与各各机组出出力之间间的关系系。同时时还对其其进行相相应的误误差分析析与灵敏敏度分析析。二：阻塞塞费用的的计算规规则在输电阻阻塞发生生时，要要调整预预案来避避免阻塞塞的现象象发生。但是调整之后，一部分序内容量不能出力，而一些在竞价中未取得发电权的发电容量（称序外容量）要在低于对应报价的清算价上出力，这样会导致网方与发电商之间的经济冲突。因此需要网方支付一定的阻塞费来作为赔偿，赔偿因涵盖以上两个方面。三：出力力分配方方案电力市场场应该以以安全第第一为原原则，并并按照一一下步骤骤实施方方案：1， 设计初步步的分配配方案，此此时因考考虑到段段容量，爬爬坡速率率的约束束。2， 计算检查查是否会会产生输输电阻塞塞。若不不发生阻阻塞，则则阻塞产产生的费费用为零零，如果果有阻塞塞产生，则则继续下下一步。3， 写出阻塞塞费用的的目标函函数，根根据已知知的约束束条件，建建立模型型求出最最优解。若若该模型型有解，则则即为最最合理的的分配方方案。若若无最优优解，则则继续下下一步。4， 因为会有有阻塞产产生，我我们只能能运用相相对安全全裕度输电电。由于于安全问问题，和和经济效效益，我我们在考虑使每每条线路路上潮流流的绝对对值超过过限值的的百分比比尽量小小的同时时，也应应该考虑虑考虑经经济最优优的方案案。5， 如果使用用安全裕裕度输电电仍不能能满足要要求，只只能采取取拉闸限限电的措措施来保保证安全全。四、模型型的建立立与求解问题一的的求解由表一、表表二给出出的各机机组出力力与相应应的各线线路有功功潮流的的大量实实验数据据，通过过对表一一中每四四组数据据（共八八项，每每项只有有一个变变量的变变化）进进行exxcell的简单单规划求求解，发发现其近近似为线线性关系系，又由由于与因因变量相相关的自自变量不不止一个个，故而而我们考考虑采用用最小二二乘准则则建立多多元线性性回归模模型对有有功潮流流值与各各机组的的关系进进行验证证。模型的建建立与求求解设为第条条线路上上的有功功潮流值值，为第第个机组组的出力力值，假假设它们们有如下下的线性性关系式式： 其中如果对变变量与自自变量同同时做次次观察得得组观察察值，采采用最小小二乘估估计求得得回归方方程： 其中中其中是待待估计的的回归系系数，为为随机误误差。为为了用矩矩阵表示示上式，令令，,于是建立立问题一一得线性性回归模模型如下:利用MAATLAAB统计计工具箱箱可得到到初步的的回归方方程,其其具体系系数矩阵阵如下:同时利用用统计工工具箱中中的reegreess命命令得出出各线路路的回归归系数估估计值及及其置信信区间（置置信水平平）、检检验统计计值的结结果统计计,其中中为回归归方程的的决定系系数（为为相关系系数），是统计量值，是与统计量对应的概率值。表一 对于线线路一的的回归模模型的结结果统计计回归系数数回归系数数的估计计值回归系数数的置信信区间110.477751099.54421 1111.4412990.088260.008088 0.084440.044780.004377 0.051180.05528 0.05114 00.055420.11199 0.11666 00.12231-0.002577-0.02777 -00.022370.122160.111900 0.124430.122200.111899 0.12551-0.000155-0.00337 00.00007表二 对于线线路二的的回归模模型的结结果统计计回归系数数回归系数数的估计计值回归系数数的置信信区间131.352211300.54461 1332.115800-0.005477-0.05663 -00.055320.122750.112400 0.13110-0.000011-0.00113 00.000100.033320.003044 0.036600.088670.008500 0.08884-0.111277-0.11550 -00.11104-0.001866 -00.02213 -0.0016000.09985 0.09666 00.10004表二 对于线线路二的的回归模模型的结结果统计计其余表格格在附录录中给出出。对表格中中结果的的分析：表一表二二显示，说明因变量基本可由模型确定，值远远超过检验的临界值，远远小于置信水平，因而我们可以认为线路一、二的回归模型是可用的，同理可知线路四六的回归模型也是可用的。模型的时时序残差差图分析析图一 第一条条线路的的时序残残差图 图二 第二条条线路的的时序残残差图其余线路路残差图图将在附附录中给给出。图形分析析：从图可以以看出残残差图中中大多数数点的残残差离零零点均较较近，且且残差的的置信区区间均包包含零点点，这说说明回归归得到的的结果能能较好的的符合原原始数据据，而一一些各别别点其置置信区间间不包含含零点，这这些点可可视为异异常点，发发现关系系式的线线性性显显著，说说明模型型的线性性拟合程程度较好好，模型型可用。从而得出出各线路路有功潮潮流值与与各发电电机组出出力的近近似关系系式如下下：对于所得得结果的的几点分分析：1、 各线路的的回归方方程中含含有常数数系数项项，可知当各各机组未未出力时时线路中中仍将具具有有功功潮流值值，所以以由基尔尔霍夫定定律可知知网络中中的线路路并不是是一个封封闭的环环路，可可能还有有若干条条外围输输电线路路，不影影响模型型的线性性关系。2、 方案中的的各线路路有功潮潮流值总总和并不不严格等等同于各各机组的的出力值值，这符符合实际际线路传传输具有有一定损损耗的事事实，有有利于模模型的推推广。3、 各个机组组的回归归系数相相差较大大，可能能造成额额外的阻阻塞费用用或资源源浪费，有有待协调调与改进进。问题二:阻塞费费用来源源：由题意可可知，网网方对发发电商的的供电取取舍有着着宏观调调控的能能力，即即发电商商有自己己制定段段容量、段段价的自自由，但但无法决决定是否否出售电电以及出出售多少少电给网网方；而而且网方方制定的的实际运运行方案案发电商商必须无无条件接接受。而电力市市场的付付费原则则是对一一个给定定的方案案，求出出其清算算价，清清算价即即为最后后付费标标准。在在网方进进行宏观观的出力力方案调调整后，由由于网方方具有宏宏观调控控的能力力可能造造成以下下情况：情况一：一些在在竞价中中由于价价钱太高高而未取取得发电电权的发发电容量量，由于于网络的的要求，而而被网方方强制购购买，而而这部分分发电容容量是在在低于所所期望的的清算价价下出卖卖电力，同同样也造造成发电电商的经经济损失失；情况二：一些通通过竞价价取得发发电权的的发电容容量，也也由于网网络自身身的限制制，而没没有全部部被网方方购买，造造成了发发电商的的一定损损失。因此制定定如下输输电阻塞塞费用计计算规则则：设序外容容量出力力部分最最低报价价为；序外出出力为：清算价价：第j组组,第ii段发电电机组序序内计划划出力，；：第j组组,第ii段发电电机组序序内实际际出力，；对于情况况一,网网方应支支付给发发电方的的赔偿为为:对于情况况二,网网方应支支付的阻阻塞赔偿偿为:故总的阻阻塞补偿偿=+问题三、四四、五:在我国电电力市场场是发电电侧市场场的国情情下，电电力市场场交易一一调度中中心作为为市场的的中间机机构，在在保证电电网经济济、安全全运行上上起着核核心作用用。它的的首要任任务就是是通过AAGC系系统得到到的各机机组的当当前出力力值、出出力上下下限、爬爬坡速率率以及下下一时段段的负荷荷需求预预报，按按照电力力市场交交易规则则，建立立一个竞竞价模型型，完成成下一时时段各机机组的最最优出力力分配预预案，使使得购电电费用最最小。我们在考考虑段容容量、爬爬坡速率率的约束束下，建建立目标标函数，要要求购电电费用最最小。在在以上分分析的基基础上建建立电力力市场分分配模型型：我们定义义符号如下下：:负荷需需求量；：清算价价；：为序外外容量出出力部分分最低报报价；：第j组组,第ii段发电电机组预预案出力力，；：第j组组,第ii段发电电机组实实际出力力，；：第j组组,第ii段的段段容量；：第i个个机组当当前出力力；：第i个个机组的的爬坡速速率；：第m条条线路上上的有功功潮流；：第m条条线路上上有功潮潮流的限限值；：第m条条线路的的相对安安全裕度度算法思想想：我们们先将所所有段的的段价按按照由低低到高进进行排序序，然后后按照排排好序的的序列选选取相应应的段容容量相加加，直到到所选取取的段容容量满足足要求的的负荷需需求量为为止。当负荷的的需求量量为9882.44MW时时，利用用mattlabb工具箱箱可以计计算得到到相应的的各机组组出力情情况机组12345678出力（MMW）1507918099.5512514095113.9此时的清清算价代入计算算输电线线路有功功潮流值值的表达达式得：线路123456潮流值173.74140.97-1500.833119.27136.81168.57超出百分分比（%）5.29950003.64424.0554由上表可可见，按按照预案案分配，会会导致线线路1、55、6发发生输电电阻塞，故需要按照步骤调整尽量避免阻塞发生，我们可以认为是阻塞费用为相应的目标函数。建立相应应的模型型：+我们采用用模拟退退火算法法，来求求解上述述问题。模模拟退火火算法的的基本思思想如下下：模拟退火算算法来源源于固体体退火原原理，将将固体加加温至充充分高，再再让其徐徐徐冷却却，加温温时，固固体内部部粒子随随温升变变为无序序状，内内能增大大，而徐徐徐冷却却时粒子子渐趋有有序，在在每个温温度都达达到平衡衡态，最最后在常常温时达达到基态态，内能能减为最最小。根根据Meetroopollis准准则，粒粒子在温温度T时时趋于平平衡的概概率为ee-E/(kT)，其中中E为温温度T时时的内能能，E为其其改变量量，k为为Bolltzmmannn常数。用用固体退退火模拟拟组合优优化问题题，将内内能E模模拟为目目标函数数值f，温温度T演演化成控控制参数数t，即即得到解解组合优优化问题题的模拟拟退火算算法：由由初始解解i和控控制参数数初值tt开始，对对当前解解重复“产生新新解计算目目标函数数差接受或或舍弃”的迭代代，并逐逐步衰减减t值，算算法终止止时的当当前解即即为所得得近似最最优解，这这是基于于蒙特卡卡罗迭代代求解法法的一种种启发式式随机搜搜索过程程。退火火过程由由冷却进进度表(Cooolinng SScheedulle)控控制，包包括控制制参数的的初值tt及其衰衰减因子子t、每每个t值值时的迭迭代次数数L和停停止条件件S。 模拟拟退火算算法的模模型模拟退退火算法法可以分分解为解解空间、目目标函数数和初始始解三部部分。模拟退退火的基基本思想想:(1) 初始始化：初初始温度度T(充充分大)，初始始解状态态S(是是算法迭迭代的起起点)， 每个TT值的迭迭代次数数L(2) 对kk=1，L做第(3)至第6步：(3) 产生新解S(4) 计算增量t=C(S)-C(S)，其中C(S)为评价函数(5) 若t<0则接受S作为新的当前解，否则以概率exp(-t/T)接受S作为新的当前解.(6) 如果满足终止条件则输出当前解作为最优解，结束程序。终止条件通常取为连续若干个新解都没有被接受时终止算法。(7) T逐渐减少，且T->0，然后转第2步。Matllab求求解，改改变调整整后的方方案机组12345678出力150.01228822880.224215297.114770117调整后线线路潮流流值：线路123456潮流值（MMW）164.99149.42154.99126.27132.00159.57超出百分分比000000此时的阻阻塞费为为47449元.同理计算算负荷需需求量为为10552.88MW时时各机组组出力情情况：机组 12345678出力(MMW)15081218.299.55135150102.1117此时清算算价。显然此时时会有线线路1，5，66产生输输电阻塞塞，按照照上述模模型对阻阻塞方案案进行调调整，但但是求不不出可行行的解，也也就是说说此时阻阻塞无法法避免，故故需要考考虑运用用裕度输输电。但但是调整整后的方方案，要要保证每每条线路路上的潮潮流绝对对值超过过限定值值的百分分比尽可可能小，也要考虑经济最优。因此，两者分别有一个权重，于是，我我们将问问题转化化为多目目标规划划问题：目标1：阻塞费费用最小小，计为为+。目标2：每条线线路上潮潮流的绝绝对值超超过限值值的百分分比尽量量小，设设F为此目目标的度度量函数数，F表示线路路上潮流流的绝对对值超过过限值的的百分比比的值。记为(MMax F)。由于两个个个目标标之间没没有必然然的联系系，求解解十分困困难，因因此我们们可以将将其统一一成为只只有一个个目标，我我们设阻阻塞费用用在该目目标中的的系数为为1，最最大潮流流超过限限制的比比例为，可可以得到到以下目目标：+在这里体体现出了了安全输输电与费费用最低低之间的的关系，越大，表明安全输电在目标函数中影响越大，反之则越小。我们在此认定，安全输电十分重要，优先安全输电，再考虑经济最优的情况，故取权值为100进行计算。此时，约约束条件件如下：当负荷需需求量时时，利用用mattlabb工具箱箱可以计计算得到到相应的的各机组组出力情情况：机组12345678出力1538822899.5515215560.33117计算此时时的的线线路的潮潮流值：线路123456潮流绝对对值173.41143.58155.21124.68135.3160.42潮流限值值165150160155132162潮流最大大值186.45177174.4172.05151.8184.68超出限值值的比例例5.1%0002.5%0经检验该该方案满满足所有有四个约约束条件件。故该该方案即即为问题题五所要要求调整整之后下下一个时时段各机机组的出出力分配配预案。根据问题题二得出出的结果果可计算算该调整整方案的的阻塞费费用，为为：8899.555元。模型评价价、 由题目所所给各线线路潮流流值只与与电网结结构和各各机组出出力有关关，通过过分析所所给大量量数据，找找出各线线路潮流流值与机机组出力力存在线线性关系系，从而而建立多多元线性性回归模模型；、 对所得结结果进行相应的的误差分分析与灵灵敏度分分析，验验证结果果的合理理性和准准确性；、 制定的阻阻塞费用用计算规规则简单单明了，并并公平地地对待序序内容量量不能出出力的部部分和报报价高于于清算价价的序外外容量出出力的部部分；、 以购电费费用最小小为目标标并考虑虑爬坡速速率约束束而建立立的线性性规划模模型，只只要给出出负荷需需求，即即可计算算得到下下一个时时段各机机组的出出力分配配预案，适适用性强强；、 对给出的的规划模模型，无无法精确确解出，我我们想到到使用模模拟退火火算法求求出近似似解；、 由于题中中只给出出了333组数据据，在利利用多元元线性回回归模型型求解时时会造成成一定的的误差，但但是通过过分析，误误差较小小，在接接受范围围以内。参考文献献：1姜姜启源 谢金星星 叶俊俊,数学学模型（第第三版），北京：高等教育出版社，2003；2赵赵东方 ，数数学模型型与计算算，北北京：科科学出版版社，220077；3李李志林 欧宜贵贵，数数学建模模及典型型案例分分析，北北京：化化学工业业出版社社，20006；附录：问题一mmatllab源源代码：MATLLAB程程序x=1120,73,1800,800,1225,1125,81.1,990133.02,73,1800,800,1225,1125,81.1,990129.63,73,1800,800,1225,1125,81.1,990158.77,73,1800,800,1225,1125,81.1,990145.32,73,1800,800,1225,1125,81.1,990120,78.5966,1880,880,1125,1255,811.1,90120,75.45,1800,800,1225,1125,81.1,990120,90.4877,1880,880,1125,1255,811.1,90120,83.8488,1880,880,1125,1255,811.1,90120,73,2311.399,800,1225,1125,81.1,990120,73,1988.488,800,1225,1125,81.1,990120,73,2122.644,800,1225,1125,81.1,990120,73,1900.555,800,1225,1125,81.1,990120,73,1800,755.8557,1125,1255,811.1,90120,73,1800,655.9558,1125,1255,811.1,90120,73,1800,877.2558,1125,1255,811.1,90120,73,1800,977.8224,1125,1255,811.1,90120,73,1800,800,1550.771,1125,81.1,990120,73,1800,800,1441.558,1125,81.1,990120,73,1800,800,1332.337,1125,81.1,990120,73,1800,800,1556.993,1125,81.1,990120,73,1800,800,1225,1138.88,81.1,990120,73,1800,800,1225,1131.21,81.1,990120,73,1800,800,1225,1141.71,81.1,990120,73,1800,800,1225,1149.29,81.1,990120,73,1800,800,1225,1125,60.5822,900120,73,1800,800,1225,1125,70.9622,900120,73,1800,800,1225,1125,64.8544,900120,73,1800,800,1225,1125,75.5299,900120,73,1800,800,1225,1125,81.1,1104.84120,73,1800,800,1225,1125,81.1,1111.22120,73,1800,800,1225,1125,81.1,998.0092120,73,1800,800,1225,1125,81.1,1120.44;>> ffor i=11:333x1(ii)=xx(i);end;>> ffor i=11:333x2(ii)=xx(i+33);end;>> ffor i=11:333x3(ii)=xx(i+66);end;>> ffor i=11:333x4(ii)=xx(i+99);end;>> ffor i=11:333x5(ii)=xx(i+1322);end;>> ffor i=11:333x6(ii)=xx(i+1655);end;>> ffor i=11:333x7(ii)=xx(i+1988);end;>> ffor i=11:333x8(ii)=xx(i+2311);end;y=1164.78,1400.877,-1144.25,1199.099,1335.444,1157.69165.81,1400.133,-1145.14,1188.633,1335.337,1160.76165.51,1400.255,-1144.92,1188.7,1355.333,1559.998167.93,1388.711,-1146.91,1177.722,1335.441,1166.81166.79,1399.455,-1145.92,1188.133,1335.441,1163.64164.94,1411.5,-1443.884,1118.43,1366.722,1557.222164.8,1141.13,-1444.007,1118.82,1366.022,1557.55165.59,1433.033,-1143.16,1177.244,1339.666,1156.59165.21,1422.288,-1143.49,1177.966,1337.998,1156.96167.43,1400.822,-1152.26,1299.588,1332.004,1153.6165.71,1400.822,-1147.08,1222.855,1334.221,1156.23166.45,1400.822,-1149.33,1255.755,1333.228,1155.09165.23,1400.855,-1145.82,1211.166,1334.775,1156.77164.23,1400.733,-1144.18,1199.122,1335.557,1157.2163.04,1400.344,-1144.03,1199.311,1335.997,1156.31165.54,1411.1,-1444.332,1118.84,1355.066,1558.226166.88,1411.4,-1444.334,1118.67,1344.677,1559.228164.07,1433.033,-1140.97,1188.755,1333.775,1158.83164.27,1422.299,-1142.15,1188.855,1334.227,1158.37164.57,1411.444,-1143.3,1119,1344.888,1558.001163.89,1433.611,-1140.25,1188.644,1333.228,1159.12166.35,1399.299,-1144.2,1119.1,1136.33,1577.599165.54,1400.144,-1144.19,1199.099,1335.881,1157.67166.75,1388.955,-1144.17,1199.155,1336.555,1157.59167.69,1388.077,-1144.14,1199.199,1337.111,1157.65162.21,1411.211,-1144.13,1166.033,1335.55,1554.226163.54,1411,-1144.16,1177.566,1335.444,1155.93162.7,1141.14,-1444.221,1116.74,1355.4,1544.888164.06,1400.944,-1144.18,1188.244,1335.44,1556.668164.66,1422.277,-1147.2,1120.21,1355.288,1557.665164.7,1142.94,-1448.445,1120.68,1355.166,1557.663164.67,1411.566,-1145.88,1199.688,1335.229,1157.61164.69,1433.844,-1150.34,1211.344,1335.112,1157.64;>> ffor i=11:333y1(ii)=yy(i);end;>> ffor i=11:333y2(ii)=yy(i+33);end;>> ffor i=11:333y3(ii)=yy(i+66);end;>> ffor i=11:333y4(ii)=yy(i+99);end;>> ffor i=11:333y5(ii)=yy(i+1322);end;>> ffor i=11:333y6(ii)=yy(i+1655);end;模型计算算及时序序残差图图的生成成n=333;m=8;>> XX=ooness(n,1),x1'',x22',xx3',x4'',x55',xx6',x7'',x88'>> b,bbintt,r,rinnt,ss=rregrresss(y11',XX,0.05);>> bb,biint,r,rrintt,srcopplott(r,rinnt)n=333;m=8;>> XX=ooness(n,1),x1'',x22',xx3',x4'',x55',xx6',x7'',x88'>> b,bbintt,r,rinnt,ss=rregrresss(y22',XX,0.05);>> bb,biint,r,rrintt,srcopplott(r,rinnt)n=333;m=8;>> XX=ooness(n,1),x1'',x22',xx3',x4'',x55',xx6',x7'',x88'>> b,bbintt,r,rinnt,ss=rregrresss(y33',XX,0.05);>> bb,biint,r,rrintt,srcopplott(r,rinnt)n=333;m=8;>> XX=ooness(n,1),x1'',x22',xx3',x4'',x55',xx6',x7'',x88'>> b,bbintt,r,rinnt,ss=rregrresss(y44',XX,0.05);>> bb,biint,r,rrintt,srcopplott(r,rinnt)n=333;m=8;>> XX=ooness(n,1),x1'',x22',xx3',x4'',x55',xx6',x7'',x88'>> b,bbintt,r,rinnt,ss=rregrresss(y55',XX,0.05);>> bb,biint,r,rrintt,srcopplott(r,rinnt)n=333;m=8;>> XX=ooness(n,1),x1'',x22',xx3',x4'',x55',xx6',x7'',x88'>> b,bbintt,r,rinnt,ss=rregrresss(y66',XX,0.05);>> bb,biint,r,rrintt,srcopplott(r,rinnt)残差图及及模型计计算结果果表格回归系数数回归系数数的估计计值回归系数数的置信信区间-1088.99928-1009.882077 -1108.16550-0.006944-0.07