《基于matlab的电力系统潮流计算课程设计(毕业论文).docx》由会员分享，可在线阅读，更多相关《基于matlab的电力系统潮流计算课程设计(毕业论文).docx（43页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、课程设计论文基于MATLAB的电力系统潮流计算 目录：一、电力系统潮流计算概述3二、潮流计算方法概述321节点分类32.2牛顿拉夫逊法概要32.2.1牛顿拉夫逊法迭代原理322.1 牛顿法的框图及求解过程3三、课程设计任务33.1题目原始数据33.2课程设计要求3四、初步分析34.1节点设置及分类34.2 参数计算34.3等值电路绘制3五、潮流计算35.1 给定负荷下的潮流计算35.1.1 B1/B2矩阵的形成35.1.2 潮流调整及分析35.1.3 仿真比较35.2 变电所负荷变化时的潮流计算35.2.1 4个变电所负荷同时以2%的比例增大35.2.2 4个变电所负荷成比例2%下降35.2.
2、3 1、4号负荷下降，2、3号负荷上升35.2.4 仿真比较35.3断线潮流计算35.5.1 断开1、3节点间的一条支路35.3.2 断开1、5支路的一条线35.3.3 断开1、7支路的一条线35.3.4 断开2、9支路的一条线35.3.5 断开7、9支路的一条线35.3.6 仿真比较3六、潮流分析总结3七、心得体会3参考文献：3附录：3潮流计算课程设计摘要电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压。各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮
3、流计算来定量的分析比较供电方案或运行方式的合理性、可靠性和经济性。此外，电力系统的潮流计算也是计算机系统动态稳定和静态稳定的基础，所以潮流计算是研究电力系统的一种重要的计算。潮流计算在数学上是多元非线性方程组的求解问题，求解的方法有很多种，牛顿拉夫逊法是数学上解非线性方程组的有效方法，有较好的收敛性。本文应用了电力系统潮流计算仿真软件DDTR与利用程序计算的结果进行比较，使计算的结果更加准确。利用成形的程序对系统中出现的各种情况，例如负荷的变化以及线路上所发生的变化进行计算，并对母线上不满足范围的数据进行调控，使得系统处于一个较稳定的状态。关键词：牛顿拉夫逊法 MATLAB DDRTS 潮流计
4、算 一、电力系统潮流计算概述在电力系统的正常运行中，随着用电负荷的变化和系统运行方式的改变，网络中的损耗也将发生变化。要严格保证所有的用户在任何时刻都有额定的电压是不可能的，因此系统运行中个节点出现电压的偏移是不可避免的。为了保证电力系统的稳定运行，要进行潮流调节。随着电力系统及在线应用的发展，计算机网络已经形成，为电力系统的潮流计算提供了物质基础。它的发展主要围绕这样几个方面：计算方法的收敛性、可靠性；计算速度的快速性；对计算机存储容量的要求以及计算的方便、灵活等。牛顿-拉夫逊法是数学中解决非线性方程式的典型方法，有较好的收敛性。在解决电力系统潮流计算问题时，是以导纳矩阵为基础的，只要我们能
5、在迭代过程中尽可能保持方程式系数矩阵的稀疏性，就可以大大提高潮流计算的效率。二、潮流计算方法概述21节点分类常规的电力系统潮流计算中一般具有三种类型的节点：PQ、PV及平衡节点。一个节点有四个变量，即注入有功功率、注入无功功率，电压大小及相角。常规的潮流计算一般给定其中的二个变量：PQ节点（注入有功功率及无功功率），PV节点（注入有功功率及电压的大小），平衡节点（电压的大小及相角）。第一类称PQ节点：等值负荷功率PLi、QLi和等值电源功率PGi、QGi是给定的，从而注入功率P、Q是给定的，待求的则是节点电压的大小Ui和相位角。属于这类节点的有按给定有功、无功率发电的发电厂母线和没有其他电源的
6、变电所母线。第二类称PV节点：等值负荷和等值电源的有功功率PLi、PGi是给定的，从而注入有功功率Pi是给定的。等值负荷的无功功率QLi和节点电压的大小Ui也是给定的。待求的则是等值电源的无功功率QGi，从而注入无功功率Qi和节点电压的相位角。有一定无功功率储备的发电厂和有一定无功功率电源的变电所母线都可以作为PV节点。第三类平衡节点：潮流计算时一般只设一个平衡节点。等值负荷功率PLs、QLs是给定的，节点电压的大小和相位也是给定的。担负调整系统频率任务的发电厂母线往往被选作为平衡节点。2.2牛顿拉夫逊法概要2.2.1牛顿拉夫逊法迭代原理已知一个变量X函数为： ，由适当的近似值出发，根据：反复
7、进行计算，当满足适当的收敛条件就是上面方程的根。这样的方法就是所谓的牛顿拉夫逊法。这一方法还可以做下面的解释，设第次迭代得到的解语真值之差，即的误差为时，则：把在附近对用泰勒级数展开上式省略去以后部分的误差可以近似由上式计算出来。比较两式，可以看出牛顿拉夫逊法的休整量和的误差的一次项相等。用同样的方法考虑，给出个变量的个方程：对其近似解得修正量可以通过解下边的方程来确定：式中等号右边的矩阵都是对于的值。这一矩阵称为雅可比（JACOBI）矩阵。按上述得到的修正向量后，得到如下关系这比更接近真实值。这一步在收敛到希望的值以前重复进行，一般要反复计算满足为预先规定的小正数，是第n次迭代的近似值。22
8、.1 牛顿法的框图及求解过程（1）用牛顿法计算潮流时，有以下的步骤：给这各节点电压初始值；将以上电压初始值代入公式，求修正方程的常数项向量 ；将电压初始值在带入上述公式，求出修正方程中系数矩阵的各元素。解修正方程式；修正各节点电压，；将，在带入方程式，求出；检验是否收敛，即如果收敛，迭代到此结束，进一步计算各线路潮流和平衡节点功率，并打印输出结果。如果不收敛，转回进行下次迭代计算，直到收敛为止。（2） 牛顿拉夫逊法计算程序框图启动输入原始数据形成节点导纳矩阵分解各节点初始电压的实部和虚部迭代次数K=0求PQ节点的,，求PV节点的，置节点号i=0雅克比矩阵是否形成，in求得雅克比矩阵各元素增大节
9、点号i=i+1把雅克比矩阵单位化求解修正方程，得，求解最大修正量，是否收敛 回带各电压新值，K=K+1计算输出电压大小及相角，节电功率及支路损耗停止否是否是三、课程设计任务3.1题目原始数据1.系统图：两个发电厂分别通过变压器和输电线路与四个变电所相连。变电所1变电所2变电所3变电所435kV母线10kV母线35kV母线10kV母线一次侧电压220kV一次侧电压220kV母线1母线3母线2线路长为100km线路长为50km线路长为90km线路长为100km2*QFQ-50-22*QFS-50-2TQN-100-22*TQN-100-2。线路长为80km。电厂一电厂二2、发电厂资料：母线1和2为
10、发电厂高压母线，发电厂一总装机容量为（ 300MW ），母线3为机压母线，机压母线上装机容量（100MW），最大负荷和最小负荷为50MW和20MW；发电厂二总装机容量（200MW）。3、变电所资料：（一）变电所1、2、3、4低压母线的电压等级分别为：35KV 10KV 35KV 10KV（二）变电所的负荷分别为：60MW 50MW 40MW 60MW （三） 每个变电所的功率因数均为cos=0.85；（四） 变电所1和变电所3分别配有两台容量为75MVA的变压器，短路损耗414KW，短路电压（%）=16.7；变电所2和变电所4分别配有两台容量为63MVA的变压器，短路损耗为245KW，短路电压
11、（%）=10.5；4、输电线路资料：发电厂和变电所之间的输电线路的电压等级及长度标于图中，单位长度的电阻为，单位长度的电抗为，单位长度的电纳为。3.2课程设计要求 1. 对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图。2. 输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷情况下的潮流计算，并对计算结果进行分析。3. 跟随变电所负荷按一定比例发生变化，进行潮流计算分析。1) 4个变电所的负荷同时以2%的比例增大；2) 4个变电所的负荷同时以2%的比例下降3) 1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的负荷同时以2%的比例上升；4. 在不同的负荷情况下
12、，分析潮流计算的结果，如果各母线电压不满足要求，进行电压的调整。（变电所低压母线电压10KV要求调整范围在9.5-10.5之间；电压35KV要求调整范围在35-36之间）5. 轮流断开支路双回线中的一条，分析潮流的分布。（几条支路断几次）6. 利用DDRTS软件，进行绘制系统图进行上述各种情况潮流的分析，并进行结果的比较。 7. 最终形成课程设计成品说明书。四、初步分析4.1节点设置及分类根据系统图可知此系统为两端供电网路，为了潮流计算的完整性，将母线1，2设为节点1，2，将变电所1、2、3、4的高低压侧分别设为节点3、4、5、6、7、8、9、10。并且，将节点1设为平衡节点，将节点2设为PV
13、节点，其余节点设为PQ节点。则网络共具有10个节点和9条支路。4.2 参数计算根据题目原始资料，计算发电厂、变压器及线路的有名值参数。（1）运用下列公式计算变压器低压侧参数： （2）计算线路参数 将参数整理成如下表格：首端号末端号阻抗有名值电纳有名值138.5+j20.1j2.780e-4154.25+j10.05j1.390e-4177.65+j18.09j2.502e-4298.5+j20.1j2.780e-4340.05455+j1.650240563.375e-3+j0.100830780.05455+j1.650240796.8+j16.08j2.224e-49103.375e-3+
14、j0.100830（3）变电所负荷分别为：变电所1 =60+j37.185 变电所2 =50+j30.987变电所3 =40+j24.790 变电所4 =60+j37.185(4）计算变压器分接头变比变压器有5个抽头，电压调节范围为2*2.5%， 对应的分接头开始时设变压器高压侧接主接头,降压变压器5个分接头时变比K以备调压时选用。 4.3等值电路绘制根据题目给定的参数和要求，手绘形成了网络的等值电路图,并且在图中标注出了各参数和节点编号。60+37.18466j发电厂一1:k21:k31357468=0U=235P=200MW8.5+20.1j4.25+10.05j7.65+18.09j0.
15、05455+1.65024j0.05455+1.65024j0.00373+0.10084j1:k150+30.98722j40+24.78977j0.000139j0.00025 j0.00025j0.000139j0.01405+0.03322j0.00028j0.00028j0.00022j0.00022 j8.5+20.1j发电厂二U=23560+37.18466j图4.3 等值电路图4.4计算方法利用牛顿拉夫逊法进行求解，用MATLAB软件编程，可以求解系统潮流分布。根据题目的不同要求对参数进行调整，通过调节变压器变比和发电厂的电压，求解出合理的潮流分布，最后用DDTRS进行潮流分析
16、，将两者进行比较。五、潮流计算5.1 给定负荷下的潮流计算 5.1.1 B1/B2矩阵的形成根据所求参数，以及B1、B2矩阵的含义，列写B1、B2矩阵如下： 138.5+j20.1j5.560e-4100154.25+j10.05j2.780e-4100B1=177.65+j18.09j5.004e-4100298.5+j20.1 j5.560e-4100340.05455+j1.650240231/38.511563.375e-3+j0.100830231/1111780.05455+j1.650240231/38.511796.8+j16.08 j4.448e-41009103.375e-
17、3+j0.10083 0231/11110 0 235 0 1200 0 235 0 30 0 220 0 2B2=0 60+j37.18466 35 0 20 0 220 0 20 50+j30.98722 10 0 20 0 220 0 20 40+j24.78977 35 0 20 0 220 0 20 60+j37.18466 10 0 2 5.1.2 潮流调整及分析将编好的的程序(程序见附录)输入MATLAB中运行，可以得到如下运行结果:节点12345电压235.0000235.0000230.241736.5113232.9385节点678910电压10.7756232.16913
18、7.5076231.395110.6310由表观察到节点4,6,8,10的电压值都较正常范围偏高，因此调节发电厂电压和变压器分接头，经过多次调整，使节点电压均满足要求，调节结果如下：调整过程：调节方法发电厂一发电厂二分接头一分接头二分接头三分接头四未调整235235231231231231第一次230230231231231231第二次225225231231231231第三次225225225.5231231231调整结果：电压12345678910未调整235.00235.00230.2436.51232.9310.77232.1637.50231.3910.63第一次230.00230.
19、00225.0435.59227.8710.52226.9636.60226.1510.37第二次225.00225.00219.8334.66222.8010.27221.7535.70220.91 10.11第三次225.00225.00219.8635.62222.8010.27221.7535.70222.9110.11由上表观察到，先后进行发电机机端电压调节和变压器分接头的调节使得节点4,6,8,10的节点电压均在允许范围内。调节前后各支路损耗对比如下：功率损耗未调整调整后(1,3)支路0.712464-28.4049i0.798149-25.6253i(1,5)支路0.251484
20、-14.6238i0.279321-13.2764i(1,7)支路1.90622-22.796i2.13376-19.9238i(2,9)支路6.33378-15.2602i6.8637-11.41i(3,4)支路0.20389+6.1682i0.21416+6.4789i(5,6)支路0.1113+3.0048i0.12238+3.3037i(7,8)支路0.085869+2.5977i0.094746+2.8663i(7,9)支路2.77942-17.3235i3.02607-14.6339i(9,10)支路0.16467+4.4453i0.18207+4.915i总损耗12.5491-8
21、2.1924i13.7144-67.3056i由此可以看出，随着发电机端电压的降低，网络有功损耗逐渐增大，符合实际情况，而最终调节使得电压在规定范围内。调整后的各支路的功率如下表所示：各条支路的首端功率Si各条支路的末端功率SjS(1,3)=61.0123+18.0382iS(3,1)=-60.2142-43.6635iS(1,5)=50.4017+21.0145iS(5,1)=-50.1224-34.2909iS(1,7)=-87.69965+67.51637iS(7,1)=89.83341-87.4402iS(2,9)=200-43.72836iS(9,2)=-193.1363+32.31
22、837iS(3,4)=60.2142+43.6635iS(4,3)=-60-37.1847iS(5,6)=50.1224+34.2909iS(6,5)=-50-30.9872iS(7,8)=40.0947+27.656iS(8,7)=-40-24.7898iS(7,9)=-129.9282+59.78417iS(9,7)=132.9542-74.41806iS(9,10)=60.1821+42.0997iS(10,9)=-60-37.1847i支路的末端功率加上支路损耗等于该支路的首段功率，符合实际情况。且得到各节点的功率 (节点号从小到大排列)：0.2371+1.0657i 2.0000-0
23、.4373i -0.0000+0.0000i -0.6000-0.3718i -0.0000-0.0000i -0.5000-0.3099i -0.0000+0.0000i -0.4000-0.2479i 0.0000+0.0000i -0.6000-0.3718i5.1.3 仿真比较由MATLAB得到电压、功角和具体的潮流分布图，调整好以后的数据如图：在DDTRS中绘制出电路图，仿真并进行调节以后可以得到潮流分布图：在用DDTRS绘制的系统潮流分布图中，各个节点电压在调试之后均能满足要求，均达到系统电压等级标准。将PSAT仿真得到的潮流图中个点电压的有名值与MATLAB计算得到的各节点电压进
24、行比较：电压12345678910DDTRS224.99224.99219.8535.61222.7910.27221.7435.70221.9010.11MATLAB225.00225.00219.8635.62222.8010.27221.7535.70222.9110.11由上述比较可观察到: 应用DDTRS仿真和MATLAB计算的到的结果大致相同，可以认为此潮流计算是正确的。5.2 变电所负荷变化时的潮流计算 5.2.1 4个变电所负荷同时以2%的比例增大由于各个变电所的负荷同时以2%的比例增大，所以此时的各负荷变为： 变电所1 =61.2+j37.92835 变电所2 = 51+j3
25、1.60696变电所3 =40.8+j25.28557变电所4 =61.2+j37.92835相应的，B2矩阵变为： 0 0 235 0 1200 0 235 0 30 0 220 0 2B2=061.2+j37.92835 35 0 20 0 220 0 20 51+j31.60696 10 0 20 0 220 0 2040.8+j25.28557 35 0 20 0 220 0 20 61.2+j37.92835 10 0 2将原程序中的B2矩阵替换后，运行并进行调节，经过调节后达到要求，对应的调节方法和结果见下表：调节方法：调节方法发电厂一发电厂二分接头一分接头二分接头三分接头四未调整
26、235235231231231231调整后225225225.5231231231调整结果： 电压12345678910未调整235.00235.00230.1036.44232.8810.76232.0837.46231.2910.61调整后225.00225.00219.7135.55222.7410.26221.6635.66220.8010.09由上表观察到，调节后的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗，记录于下表：功率损耗未调整调整后(1,3)支路0.746602-28.3064i0.836494-25.5166i(1,5)支路0.26271
27、7-14.5939i0.291806-13.2435i(1,7)支路1.86419-22.8856i2.08903-20.0197i(2,9)支路6.32203-15.275i6.8524-11.4236i(3,4)支路0.2129+6.4407i0.22369+6.767i(5,6)支路0.116+3.1315i0.12756+3.4437i(7,8)支路0.089529+2.7084i0.098809+2.9891i(7,9)支路2.738-17.4024i2.98121-14.7207i(9,10)支路0.17175+4.6366i0.18995+5.128i总损耗12.5237-81.
28、5461i13.6909-66.5963i由电压和功率损耗可以观察到，最后一次的电压和有功损耗均符合题目要求，可认为是合理的。具体的潮流分布如下：各条支路的首端功率Si各条支路的末端功率SjS(1,3)=62.2602+19.1788iS(3,1)=-61.4237-44.6954iS(1,5)=51.4194+21.8072iS(5,1)=-51.1276-35.0507iS(1,7)=-85.78861+67.73745iS(7,1)=87.87764-87.75718iS(2,9)=200-42.57042iS(9,2)=-193.1476+31.14682iS(3,4)=61.4237
29、+44.6954S(4,3)=-61.2-37.9284iS(5,6)=51.1276+35.0507iS(6,5)=-51-31.607iS(7,8)=40.8988+28.2747iS(8,7)=-40.8-25.2856iS(7,9)=-128.7764+59.48247iS(9,7)=131.7577-74.20318iS(9,10)=61.39+43.0564iS(10,9)=-61.2-37.9284i各节点的功率S为(节点号从小到大排列)：0.278 +1.0872i 2.0000-0.4257i 0.0000-0.0000i -0.6120-0.3793i 0.0000+0.0
30、000i -0.5100 0.3161i -0.000-0.0000i -0.4080-0.2529i 0.0000+0.0000i -0.6120-0.3793i5.2.2 4个变电所负荷成比例2%下降由于各个变电所的负荷同时以2%的比例下降，所以此时的各负荷变为： 变电所1 =58.8+j36.44097 变电所2 = 49+j30.36747变电所3 =39.2+j24.29398变电所4 =58.8+j36.44097相应的，B2矩阵变为： 0 0 235 0 1200 0 235 0 30 0 220 0 2B2=058.8+j36.44097 35 0 20 0 220 0 20
31、49+j30.36747 10 0 20 0 220 0 2039.2+j24.29398 35 0 20 0 220 0 20 58.8+j36.44097 10 0 2将程序中的B2矩阵替换后，运行并进行调节，此系统经过调节后达到要求，对应的调节方法和结果见下表：调节方法：调节方法发电厂一发电厂二分接头一分接头二分接头三分接头四未调整235235231231231231调整后225225225.5231231231调整结果： 电压12345678910未调整235.00235.00230.3836.57232.9810.78232.2537.54231.4910.64调整后225.0022
32、5.00220.0135.69222.8510.28221.8335.74221.0110.12由上表观察到，调节后的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗，记录于下表：功率损耗未调整调整后(1,3)支路0.679291-28.5011i0.760896-25.7314i(1,5)支路0.24052-14.653i0.267137-13.3086i(1,7)支路1.94936-22.7038i2.17972-19.825i(2,9)支路6.34587-15.2445i6.87539-11.3954i(3,4)支路0.19512+5.9027i0.2048
33、9+6.1982i(5,6)支路0.10671+2.8808i0.11731+3.1668i(7,8)支路0.082294+2.4896i0.09078+2.7463i(7,9)支路2.82114-17.2439i3.07126-14.5462i(9,10)支路0.15775+4.2586i0.17437+4.7072i总损耗12.5781-82.8145i13.7417-67.988i由电压和功率损耗可以观察到，最后一次的电压和有功损耗均符合题目要求，可认为是合理的。具体的潮流分布如下：各条支路的首端功率Si各条支路的末端功率SjS(1,3)=59.7658+16.9078iS(3,1)=-
34、59.0049-42.6392iS(1,5)=49.3844+20.2257iS(5,1)=-49.1173-33.5343iS(1,7)=-89.60849+67.30186iS(7,1)=91.7882-87.12685iS(2,9)=200-44.88005iS(9,2)=-193.1246+33.48465iS(3,4)=59.0049+42.6392iS(4,3)=-58.8-36.441iS(5,6)=49.1173+33.5343iS(6,5)=-49-30.3675iS(7,8)=39.2908+27.0402iS(8,7)=-39.2-24.294iS(7,9)=-131.0
35、79+60.08661iS(9,7)=134.1502-74.63283iS(9,10)=58.9744+41.1482iS(10,9)=-58.8-36.441i各节点的功率S为(节点号从小到大排列)：0.1954+1.0444i 2.0000-0.4488i -0.0000+0.0000i -0.5880-0.3644i 0.0000+0.0000i -0.4900-0.3037i -0.0000+0.0000i -0.3920-0.2429i -0.0000-0.0000i - 0.5880-0.3644i5.2.3 1、4号变电所负荷成2%下降，2、3号变电负荷成2%上升此时的各负荷变
36、为：变电所1 =58.8+j36.44097 变电所2 = 51+j31.60696变电所3 =40.8+j25.28557变电所4 =58.8+j36.44097相应的，B2矩阵变为： 0 0 235 0 1200 0 235 0 30 0 220 0 2B2=058.8+j36.44097 35 0 20 0 220 0 20 51+j31.60696 10 0 20 0 220 0 2040.8+j25.28557 35 0 20 0 220 0 20 58.8+j36.44097 10 0 2将程序中的B2矩阵替换后，运行并进行调节，此系统经过调节后达到要求，对应的调节方法和结果见下表
37、：调节方法：调节方法发电厂一发电厂二分接头一分接头二分接头三分接头四未调整235235231231231231调整后225225225.5231231231调整结果： 电压12345678910未调整235.00235.00230.3836.57232.8810.76232.1537.48231.4410.64调整后225.00225.00220.0135.69222.7410.26221.7335.67220.9610.12由上表观察到，调节后的电压均在题目允许的范围内，对线路损耗进行分析，统计每次调整后各个支路的有功损耗，记录于下表：功率损耗未调整调整后(1,3)支路0.679291-28
38、.5011i0.760896-25.7314i(1,5)支路0.262717-14.5939i0.291806-13.2435i(1,7)支路1.92183-22.7577i2.15108-19.8814i(2,9)支路6.33938-15.253i6.8691-11.4033i(3,4)支路0.19512+5.9027i0.20489+6.1982i(5,6)支路0.116+3.1315i0.12756+3.4437i(7,8)支路0.08947+2.7066i0.098736+2.987i(7,9)支路2.81284-17.248i3.06221-14.552i(9,10)支路0.1578
39、3+4.2608i0.17446+4.7098i总损耗12.5745-82.352i13.7407-67.4728i由电压和功率损耗可以观察到，最后一次的电压和有功损耗均符合题目要求，可认为是合理的。具体的潮流分布如下：各条支路的首端功率Si各条支路的末端功率SjS(1,3)=59.7658+16.9078iS(3,1)=-59.0049-42.6392iS(1,5)=51.4194+21.8072iS(5,1)=-51.1276-35.0507iS(1,7)=-88.04441+67.85231iS(7,1)=90.19549-87.7337iS(2,9)=200-44.26567iS(9,2)=-193.1309+32.86236iS(3,4)=59.0049+42.6392iS(4,3)=-58.8-36.441iS(5,6)=51.1276+35.0507iS(6,5)=-51-31.607iS(7,8)=40.8987+28.2725iS(8,7)=-40.8-25.2856iS(7,9)=-131.0942+59.46118iS(9,7)=134.1564-74.01315iS(9,10