复杂网络NR法潮流分析计算

复杂网络N-R法潮流解析计算复杂网络N-R法潮流解析计算30/30复杂网络N-R法潮流解析计算大连理工大学专业学位硕士学位论文编号课程设计2012级本科）题目：复杂网络N-R法潮流解析与计算的设计系（部）院：物理与机电工程学院专业：电气工程及其自动化作者姓名：指导教师：职称：副教授完成日期：2015年6月30日二○一五年六月常用软件课程设计大连理工大学专业学位硕士学位论文河西学院本科生课程设计任务书设计题目复杂网络N-R法潮流解析与计算的设计作者姓名学院、专业、年级物电学院电气工程专业12级指导教师姓名、职副教授任务下达日期2015年5月20日称一、设计资料1.系统图的确定选择六节点、环网、两电源和多引出的电力系统，简化电力系统图如图1所示，等值阻抗图如图2所示。运用以直角坐标表示的牛顿-拉夫逊计算如图1所示系统中的潮流分布。计算精度要求各节点电压的误差或修正量不大于105。1.8+j0.41.6+j0.8③④0.08+j0.3L4L2L50.06+j0.0250.1+j0.35②L3⑤0.04+j0.252+j13.7+j1.31.05:11.05:1U=1.05U=1.05①GG⑥图1电力系统图常用软件课程设计大连理工大学专业学位硕士学位论文1.8+j4.01.6+j0.8③④5j0.25j0.255230..0j0.250jj++16.00.0j0.25j0.25j0.25②⑤2+j10.04+j0.253.7+j1.31.05:11.05:1j0.03①U11.05⑥o0图2电力系统等值阻抗图各节点的初值及阻抗参数该系统中，节点①为平衡节点，保持U1=1.05+j0为定值，节点⑥为PV节点，其他四个节点都是PQ节点。给定的注入电压标幺值、线路阻抗标幺值、线路阻抗标幺值、输出功率标幺值和变压器变比标幺值如图2所示的说明。表1各节点电压标幺值参数U1U2U3U4U5U61.051.001.001.001.001.05表2线路、变压器阻抗标幺值线路TL2L3L4L5L6阻抗j0.030.06+j0.0250.04+j0.250.08+j0.300.1+j0.35j0.015表3节点输出功率节点②③④⑤⑥功率2+j11.8+j0.401.6+j0.83.7+j1.35注：各PQ节点的电压取1是为了方便计算和最后考据程序的正确性。常用软件课程设计大连理工大学专业学位硕士学位论文二、设计的基本要求3.1设计及计算说明书（1）说明书要求书写整齐，井井有条，表达正确、语言正确。（2）计算书内容：为各设计内容最后成就、确定供应依照进行的技术解析、论证和定量计算，如。（3）计算书要求：计算无误，解析论证过程简单了然，各设计内容列表汇总。3.2图纸（1）绘制解析所需的必要图纸（2）图纸要求：用标准符号绘制，部署均匀，设备符号大小合适，清楚雅观。三、论文(设计)进度安排阶段论文（设计）各阶段名称起止日期1熟悉设计任务书、设计题目及设计背景资料5.20～5.252查阅相关资料5.26～5.273阅读设计要求必读的参照资料5.28～5.294书写设计说明书5.30～6.155小组争辩思疑6.21～6.226上交设计成就6.30四、需收集和阅读的资料及参照文件（指导教师指定）[1]:陈珩.电力系统稳态解析（第三版）[M]，北京，中国电力初版社，2007：何仰赞.温增银.《电力系统解析》第三版[M]，武汉,华中科技大学初版社，2002[3]：陈悦.《电气工程毕业设计指南电力系统分册》[M]，北京，中国水利水电初版社，2008：：教研室意见负责人签字：年代日常用软件课程设计大连理工大学专业学位硕士学位论文目录1牛顿－拉夫逊法归纳11.1牛顿-拉夫逊法基根源理11.2牛顿--拉夫逊法潮流求解过程22手算潮流计算52.1用图中数据和等值网络形成节点导纳矩阵Y5B2.2设各节点电压初始值为：62.3用公式62.4求取雅可比矩阵72.5求△修正量矩阵72.6计算修正各节点电压83计算机算法潮流计算83.1牛顿—拉夫逊法的程序框图83.2结果显示9总结18附件19参照文件23常用软件课程设计大连理工大学专业学位硕士学位论文牛顿－拉夫逊法归纳电力系统潮流计算是电力系统解析中的一种最基本的计算，是对复杂电力系统正常和故障条件下稳态运行状态的计算。潮流计算的目标是求取电力系统在给定运行状态的计算。即节点电压和功率分布，用以检查系统各元件可否过负荷。各点电压可否满足要求，功率的分布和分配可否合理以及功率耗费等。对现有电力系统的运行和扩建，对新的电力系统进行规划设计以及对电力系统进行静态和暂态牢固解析都是以潮流计算为基础。潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算的模型和方法有直接影响。实质电力系统的潮流技术那主要采用牛顿-拉夫逊法。1.1牛顿-拉夫逊法基根源理牛顿--拉夫逊法(简称N—R法)在数学上是求解非线性代数方程式的有效方法。其要点是把非线性方程式的求解过程变成屡次地对相应的线性方程式进行求解的过程。即平时所称的逐次线性化过程。对于非线性代数方程组：=0即(1-1-1)在待求量x的某一个初始估计值x(0)周边，将上式张开成泰勒级数并略去二阶及以上的高阶项，得到以下的经线性化的方程组：(1-1-2)上式称之为牛顿法的修正方程式。由此可以求得第一次迭代的修正量(1-1-3)将x(0)和x(0)相加，获取变量的第一次改进值x(1)。接着就从x(1)出发，重复上述计算过程。所以从必然的初值x(0)出发，应用牛顿法求解的迭代格式为：(1-1-4)(1-1-5)上两式中：'f(x)是函数f(x)对于变量xJk为迭代次数。的一阶偏导数矩阵，即雅可比矩阵—有上式可见，牛顿法的核心即是屡次形式并求解修正方程式。牛顿法当初始估计值x(0)和方程的精确解足够凑近时，收敛速度特别快，拥有平方收敛特点。牛顿潮流算法突出的优点是收敛速度快，若选择到一个较好的初值，算法将拥有平方收敛特性，一般迭代4~5次便可以收敛到一个特别精确的解。而且其迭代次数与所计算网络的规模基本没关。牛顿法也拥有优异的收敛可靠性，对于对以节点导纳矩阵为基础的高斯法呈病态的系统，牛顿法也能可靠收敛。牛顿法所需的内存量及每次迭代所需时间均较高斯法多。牛顿法的可靠收敛取决于有一个优异的启动初值。若是初值选择不当，算法有可能根本不收常用软件课程设计大连理工大学专业学位硕士学位论文敛或收敛到一个无法运行的节点上。对于正常运行的系统，各节点电压一般均在额定值周边，偏移不会太大，而且各节点间的相位角差也不大，所以对各节点可以采用一致的电压初值(也称为平直电压)，如假定：或(1-1-6)这样一般能获取满意的结果。但若系统因无功紧张或其他原因以致电压质量很差或有重载线路而节点间角差很大时，仍用上述初始电压就有可能出现问题。解决这个问题的方法可以用高斯法迭代1~2次，以此迭代结果作为牛顿法的初值。也可以先用直流法潮流求解一次以求得一个较好的角度初值，尔后转入牛顿法迭代。1.2牛顿--拉夫逊法潮流求解过程以下谈论的是用直角坐标形式的牛顿—拉夫逊法潮流的求解过程。当采用直角坐标时，潮流问题的待求量为各节点电压的实部和虚部两个重量由于平衡节点的电压向量是给定的，所以待求两共需要2(n-1)个方程式。事实上，除了平衡节点的功率方程式在迭代过程中没有拘束作用以外，其他每个节点都可以列出两个方程式。对PQ节点来说和是给定的，所以可以写出（1-2-1）对PV节点来说，给定量是和，所以可以列出（1-2-2）求解过程大体可以分为以下步骤：1）形成节点导纳矩阵2）将各节点电压设初值U，3）将节点初值代入相关求式，求出修正方程式的常数项向量4）将节点电压初值代入求式，求出雅可比矩阵元素5）求解修正方程，求修正向量6）求取节点电压的新值（7）检查可否收敛，如不收敛，则以各节点电压的新值作为初值自第3步重新开始进行狭义次迭代，否则转入下一步（8）计算支路功率分布，PV节点无功功率和平衡节点注入功率。常用软件课程设计大连理工大学专业学位硕士学位论文以直角坐标系形式表示迭代计算式采用直角坐标时,节点电压相量及复数导纳可表示为:(1-2-3)将以上二关系式代入上式中,张开并分开实部和虚部;假定系统中的第1,2,,m号为P—Q节点,m+1,m+2,,n-1为P—V节点,依照节点性质的不同样,获取以下迭代计算式:对于PQ节点(1-2-4)i=1,2,,m对于PV节点(1-2-5)⑶对于平衡节点平衡节点只设一个,电压为已知,不拜会迭代,其电压为:(1-2-6)修正方程式(2-3-5)和(2-3-6)两组迭代式共包括2(n-1)个方程.选定电压初值及变量修正量符号此后代入式(2-3-5)和(2-3-6),并将其按泰勒级数张开,略去,二次方程及今后各项,获取修正方程如下:(1-2-7)=

P1Q1PmQmPm1Um21Pn1Un21常用软件课程设计大连理工大学专业学位硕士学位论文PPPPPPPP11111111e1f1emfmem1fm1en1fn1Q1Q1Q1Q1Q1Q1Q1Q1e1f1emfmem1fm1en1fn1PmPmPmPmPmPmPmPme1f1emfmem1fm1en1fn1QmQmQmQmQmQmQmQmJe1f1emfmem1fm1en1fn1Pm1Pm1Pm1Pm1PPm1PPm1m1m1e1f1emfmem1fm1en1fn1U2m1U2m1U2m1U2m1U2m1U2m1U2m1U2m1e1f1emfmem1fm1en1fn1Pn1Pn1Pn1Pn1Pn1Pn1Pn1Pn1e1f1emfmem1fm1en1fn1(1-2-8)U2n1U2n1U2n1U2n1U2n1U2n1U2n1U2n1e1f1emfmem1fm1en1fn1③雅可比矩阵各元素的算式式(1-2-8)中,雅可比矩阵中的各元素可经过对式(1-2-4)和(1-2-5)进行偏导而求得.当ji时,雅比矩阵中非对角元素为(1-2-9)当ji时,雅可比矩阵中对角元素为:Pin(GijejBijfj)GiieiBiifieij1Pin(GijfjBijej)GiifiBiieifij1Qin(GijfjBijej)GiifiBiieieij1(1-2-10)Qin(GijejBijfj)GiieiBiififjj1Ui2eej2iUi2ffi2i由式(3-2-9和(3-2-10)看出,雅可比矩阵的特点:阵中各元素是节点电压的函数,在迭代过程中,这些元素随着节点电压的变化而变化;常用软件课程设计大连理工大学专业学位硕士学位论文⒉导纳矩阵中的某些非对角元素为零时,雅可比矩阵中对应的元素也是为零.若Yij0,则必有;⒊雅可比矩阵不是对称矩阵;(i=q1,2,,n,i≠s)雅可比矩阵各元素的表示以下:=手算潮流计算2.1用图中数据和等值网络形成节点导纳矩阵YB节点导纳矩阵对角线上的元素为：非对角线上的元素为：常用软件课程设计大连理工大学专业学位硕士学位论文所以导纳矩阵为Z=j33.3333j31.74600000j31.74601.53174j37.41660.9077j3.782200.6240j3.9002000.9077j3.78221.7376j6.39420.8299j3.112000000.8299j3.11201.5846j5.25350.7547j2.6415000.6240j3.900200.7547j2.64151.3787j66.5103j63.49210000j63.4921j66.66672.2设各节点电压初始值为：U=e+fe1=1.05f1=0e2=1f2=0e3=1f3=0e4=1f4=0e5=1f5=0e6=1.05f6=02.3用公式对PQ和PV节点求取,得常用软件课程设计大连理工大学专业学位硕士学位论文25995816373956502.4求取雅可比矩阵1.531841.01560.90773.7822000.62403.90020033.81761.53173.78220.9077003.90020.6240000.90773.78221.73766.89420.82993.112000003.78220.90775.89421.73763.11200.82990000000.82993.11201.58465.75350.75472.641500J03.11200.82994.75351.58462.64150.75470000.62403.9002000.75472.64151.378773.2083063.49213.90020.62400002.64150.754759.81231.378763.49210000000066.6667063.4921000000002.102.5求△修正量矩阵常用软件课程设计大连理工大学专业学位硕士学位论文0.094700.123000.061180.604410.13159x0.637770.096250.158570.000000.087752.6计算修正各节点电压=1.05000=1.09470=0.93882=0.86841=1.09625=1.05000=0=-0.12300=-0.60441=-0.63777=-0.15857=-0.08775计算机算法潮流计算3.1牛顿—拉夫逊法的程序框图常用软件课程设计大连理工大学专业学位硕士学位论文3.2结果显示：Y=0-33.3333i0+31.7460i00000+31.7460i1.5317-37.4166i-0.9077+3.7822i0-0.6240+3.9002i00-0.9077+3.7822i1.7376-6.3942i-0.8299+3.1120i0000-0.8299+3.1120i1.5846-5.2535i-0.7547+2.6415i00-0.6240+3.9002i0-0.7547+2.6415i1.3787-66.5103i0+63.4921i00000+63.4921i0-66.6667i次数time=1雅可比矩阵JJ=0.031229.5393-1.5317-37.41660.90773.7822000.62403.900233.95280.0312-37.41661.53173.7822-0.9077003.9002-0.6240000.90773.2822-1.7376-6.39420.82993.112000004.2821-0.9077-6.39421.73763.1120-0.82990000000.86762.4800-1.5846-5.25350.75472.641500003.7441-0.7921-5.25351.58462.6415-0.7547000.62403.9002000.6858-0.7310-1.3787-66.5103003.9002-0.6240006.0140-0.8237-66.5103常用软件课程设计大连理工大学专业学位硕士学位论文1.378700000000070.000000000000-2.10000E=0.81370.83660.67130.66361.05001.0500F=-0.1896-0.2825-0.7053-0.6191-0.07140U=0.8137-0.1896i0.8366-0.2825i0.6713-0.7053i0.6636-0.6191i1.0500-0.0714i1.0500dU=0.18630.18960.16340.28250.32870.70530.33640.619100.0714PQ=-1.96881.2067-1.80000.1000-1.5623-0.1679-3.76892.07255.00000precision=5次数time=2雅可比矩阵JJ=7.606224.4159-8.3403-30.15691.45572.9056001.24723.055427.2498-4.4313-30.15698.34032.9056-1.4557003.0554-1.2472常用软件课程设计大连理工大学专业学位硕士学位论文003.48652.4260-3.2599-4.85841.57342.369000003.3894-0.1690-4.85843.25992.3690-1.57340000003.91870.4693-4.7692-2.40892.36971.240900002.5381-1.5855-2.40894.76921.2409-2.3697002.82862.2020003.7251-1.1059-42.0899-43.2856002.2020-2.8286003.6775-0.5472-43.285642.089900000000070.000000000000-2.1000-0.1429E=1.12611.13330.86980.73241.04761.0500F=-0.1055-0.2062-0.5336-0.4376-0.07140U=1.1261-0.1055i1.1333-0.2062i0.8698-0.5336i0.7324-0.4376i1.0476-0.0714i1.0500dU=-0.3124-0.0841-0.2967-0.0762-0.1985-0.1717-0.0688-0.18150.00240常用软件课程设计大连理工大学专业学位硕士学位论文PQ=-0.4396-0.1479-0.2762-0.4656-0.0873-0.9285-1.1649-0.69650-0.0051precision=1.1649次数time=3雅可比矩阵JJ=5.127434.8149-5.6738-41.97401.42144.1634001.11434.326236.6848-1.5735-41.97405.67384.1634-1.4214004.3262-1.1143003.20413.4717-3.2879-6.88801.58233.355600004.7263-0.4133-6.88803.28793.3556-1.58230000003.18370.8345-4.1815-3.72382.06591.894800003.6933-1.5811-3.72384.18151.8948-2.0659002.16362.5835004.3148-2.1031-30.1120-48.1109002.5835-2.1636005.31200.8976-48.110930.112000000000070.000000000000-2.0951-0.1429E=1.12211.13360.87160.71591.04761.0500F=-0.0117-0.1248-0.5167-0.4513-0.07140U=1.1221-0.0117i1.1336-0.1248i0.8716-0.5167i0.7159-0.4513i1.0476-0.0714i1.0500常用软件课程设计大连理工大学专业学位硕士学位论文dU=0.0040-0.0939-0.0004-0.0814-0.0019-0.01690.01650.01370.0000-0.0000PQ=0.0997-0.1346-0.08920.0231-0.14030.0157-0.16890.2752-0.0000-0.0000precision=0.2752次数time=4雅可比矩阵JJ=2.145834.7123-2.1557-41.96691.06274.2333000.74584.369036.53151.4043-41.96692.15574.2333-1.0627004.3690-0.7458003.03723.6423-2.7679-7.03181.32923.424300004.70640.0350-7.03182.76793.4243-1.32920000003.28520.7975-4.0959-3.76042.02281.912500003.7700-1.3777-3.76044.09591.9125-2.0228002.20692.5105004.6076-2.0746-31.0036-46.9932002.5105-2.2069005.17561.1428-46.9932常用软件课程设计大连理工大学专业学位硕士学位论文31.003600000000070.000000000000-2.0951-0.1429E=1.12011.13210.87310.71741.04761.0500F=-0.0113-0.1244-0.5167-0.4519-0.07140U=1.1201-0.0113i1.1321-0.1244i0.8731-0.5167i0.7174-0.4519i1.0476-0.0714i1.0500dU=0.0020-0.00040.0016-0.0004-0.0014-0.0001-0.00150.00060.00000.0000PQ=0.0024-0.00010.00780.0114-0.00070.0027-0.0056-0.00230.0000-0.0000precision=0.0114次数time=5雅可比矩阵JJ=2.135434.6468-2.1387-41.89161.05954.2260000.74304.361436.46831.4174-41.89162.13874.2260-1.0595004.3614-0.7430常用软件课程设计大连理工大学专业学位硕士学位论文003.03093.6468-2.7627-7.02231.32673.419700004.69050.0345-7.02232.76273.4197-1.32670000003.28810.8064-4.0978-3.76802.02371.916300003.7701-1.3767-3.76804.09781.9163-2.0237002.21002.5161004.6104-2.0690-31.0428-47.0929002.5161-2.2100005.17711.1403-47.092931.042800000000070.000000000000-2.0951-0.1429E=1.12011.13200.87310.71741.04761.0500F=-0.0113-0.1244-0.5167-0.4519-0.07140U=1.1201-0.0113i1.1320-0.1244i0.8731-0.5167i0.7174-0.4519i1.0476-0.0714i1.0500dU=1.0e-005*0.1966-0.04810.1552-0.0453-0.1877-0.1040-0.1272-0.1305-0.0000常用软件课程设计大连理工大学专业学位硕士学位论文0PQ=1.0e-004*-0.0319-0.01800.09390.14950.0312-0.0620-0.0088-0.031700.0000precision=1.4954e-005次数time=6雅可比矩阵JJ=2.135434.6467-2.1387-41.89151.05954.2260000.74304.361436.46821.4175-41.89152.13874.2260-1.0595004.3614-0.7430003.03093.6468-2.7627-7.02231.32673.419700004.69040.0345-7.02232.76273.4197-1.32670000003.28810.8064-4.0978-3.76812.02371.916300003.7701-1.3767-3.76814.09781.9163-2.0237002.21002.5161004.6104-2.0690-31.0427-47.0930002.5161-2.2100005.17711.1403-47.093031.042700000000070.000000000000-2.0951-0.1429E=1.12011.13200.87310.71741.04761.0500F=-0.0113-0.1244-0.5167-0.4519-0.07140常用软件课程设计大连理工大学专业学位硕士学位论文U=1.1201-0.0113i1.1320-0.1244i0.8731-0.5167i0.7174-0.4519i1.0476-0.0714i1.0500dU=1.0e-011*0.17480.30540.16910.2967-0.09460.13580.05680.0980-0.00000PQ=1.0e-010*-0.1365-0.01870.13130.1907-0.0803-0.03560.11800.017600.0000precision=1.9074e-011S=-19.5876+27.9344i线路流动功率S12=-0.0000+2.5240iS25=-0.0000+2.5240iS23=4.7133+1.7798iS34=3.4698+1.3835iS45=0.3937+1.3666iS56=常用软件课程设计大连理工大学专业学位硕士学位论文-28.1644+18.3564i功率耗费S=-19.5876+27.9344i总结近两周的课程设计时间是紧迫的，在此过程第一让我对所学知识进行重新回顾，理解了使电力系统牢固的运行，必定经过精美的设计和运算。在进行此次课程设计的过程中，加深了我对潮流计算的认识，特别是对牛顿拉夫逊潮流计算的求解思路有了比较透彻的理解。经过自己着手查资料，解析原理，绘制电路以致撰写文档感觉了设计的难度。综合检验和牢固自己学到的知识。而且在此次课程设计中，我发现了自己的基础知识不足。在今后的学习中应注意基础知识的掌握。此次设计中采用Matlab进行潮流计算，提高了正确性，使用方便，已获取成功应用。在此，衷心的感谢老师对我们的悉心指导，感谢老师!常用软件课程设计大连理工大学专业学位硕士学位论文附件计算机算法程序：y11=-33.33333j;y12=31.74603j;y13=0;y14=0;y15=0;y16=0;y21=31.74603j;y22=1.53174-37.41662j;y23=-0.90772+3.78215j;y24=0;y25=-0.62402+3.90016j;y26=0;y31=0;y32=-0.90772+3.78215j;y33=1.7376-6.39418j;y34=-0.8299+3.11203j;y35=0;y36=0;y41=0;y42=0;y43=-0.82988+3.11203j;y44=1.5846-5.25354j;y45=-0.75472+2.64151j;y46=0;y51=0;y52=-0.62402+3.90016j;y53=0;y54=-0.75472+2.64151j;y55=1.37874-66.5103j;y56=63.49206j;常用软件课程设计大连理工大学专业学位硕士学位论文y61=0;y62=0;y63=0;y64=0;y65=63.4921j;y66=-66.66667j;Y=[y11y12y13y14y15y16;y21y22y23y24y25y26;y31y32y33y34y35y36;y41y42y43y44y45y46;y51y52y53y54y55y56;y61y62y63y64y65y66]W=zeros(6);fori=1:6%导纳矩阵第一行放到最后一行，方便下面程序执行。W(1,i)=Y(1,i);endfori=1:6Y(1,i)=Y(2,i);Y(2,i)=Y(3,i);Y(3,i)=Y(4,i);Y(4,i)=Y(5,i);Y(5,i)=Y(6,i);Y(6,i)=W(1,i);endE(1)=1.00;E(2)=1.00;E(3)=1.00;E(4)=1.00;E(5)=1.05;%为方便程序计算，以下程序节点1放到最后当作6点，F(1)=0;F(2)=0;F(3)=0;F(4)=0;F(5)=0;G=real(Y);B=imag(Y);S(1)=-2-1j;S(2)=-1.8-0.4j;S(3)=-1.6-0.8j;S(4)=-3.7-1.3j;S(5)=5;P=real(S);Q=imag(S);k=0;precision=1;N1=5;whileprecision>0.00001E(6)=1.05;F(6)=0;form=1:N1forn=1:N1+1Pt(n)=(E(m)*(G(m,n)*E(n)-B(m,n)*F(n))+F(m)*(G(m,n)*F(n)+B(m,n)*E(n)));常用软件课程设计大连理工大学专业学位硕士学位论文Qt(n)=(F(m)*(G(m,n)*E(n)-B(m,n)*F(n))-E(m)*(G(m,n)*F(n)+B(m,n)*E(n)));enddP(m)=P(m)-sum(Pt);dQ(m)=Q(m)-sum(Qt);endform=1:N1forn=1:N1+1Bi(n)=G(m,n)*F(n)+B(m,n)*E(n);Ai(n)=G(m,n)*E(n)-B(m,n)*F(n);endH(m,m)=-sum(Bi)+B(m,m)*E(m)-G(m,m)*F(m);N(m,m)=-sum(Ai)-G(m,m)*E(m)-B(m,m)*F(m);J(m,m)=-sum(Ai)+G(m,m)*E(m)+B(m,m)*F(m);L(m,m)=sum(Bi)+B(m,m)*E(m)-G(m,m)*F(m);endform=1:N1JJ(2*m-1,2*m-1)=N(m,m);JJ(2*m-1,2*m)=H(m,m);JJ(2*m,2*m-1)=L(m,m);JJ(2*m,2*m)=J(m,m);endform=1:N1forn=1:N1ifm==nelseH(m,n)=B(m,n)*E(m)-G(m,n)*F(m);N(m,n)=-B(m,n)*F(m)-G(m,n)*E(m);J(m,n)=B(m,n)*F(m)+G(m,n)*E(m);L(m,n)=B(m,n)*E(m)-G(m,n)*F(m);JJ(2*m-1,2*n)=H(m,n);JJ(2*m-1,2*n-1)=N(m,n);常用软件课程设计大连理工大