最优潮流简化梯度法计算程序-推荐下载

2最优潮流简化梯度法计算程序程序所用例子：节点参数：线路参数IJRX130.100.4140.150.6150.050.2230.050.2250.050.2340.100.4程序：clearall;clf;y=[2.1569-8.6274i0-0.5882+2.3529i-0.3922+1.5686i-1.1765+4.7059i;02.3530-9.4118i-1.1765+4.7059i0-1.1765+4.7059i;-0.5882+2.3529i-1.1765+4.7059i2.3529-9.4117i-0.5882+2.3529i0;-0.3922+1.5686i0-0.5882+2.3529i0.9804-3.9215i0;-1.1765+4.7059i-1.1765+4.7059i002.3530-9.4118i];%导纳矩阵g=real(y);b=imag(y);u(:,1)=[1;1;1;1;1];%各节点初始电压值15234T3、管路敷设技术通1过管线0不仅可以解决吊顶层配置不规范高中资料试卷问题22，而且可保障各类管路习题到位。在管路敷设过程1中，要加强看2222222222护1关于管路高中资料试卷连接管口处理高中资料试卷弯扁度固定盒位置.保护层防腐跨接地线弯曲半径标高等，要求技术交底。管线敷设技术中包含线槽、管架等多项方式，为解决高中语文电气课件中管壁薄、接口不严等问题，合理利用管线敷设技术。线缆敷设原则：在分线盒处，当不同电压回路交叉时，应采用金属隔板进行隔开处理；同一线槽内，强电回路须同时切断习题电源，线缆敷设完毕，要进行检查和检测处理。4、电气课件中调试技术对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料试卷相互作用与相互关系，根据生产工艺高中资料试卷要求，对电气设备进行空载与带负荷下高中资料试卷调控试验；对设备进行调整使其在正常工况下与过度工作下都可以正常工作；对于继电保护进行整核对定值，审核与校对图纸，编写复杂设备与装置高中资料试卷调试方案，编写重要设备高中资料试卷试验方案以及系统启动方案；对整套启动过程中高中资料试卷电气设备进行调试工作并且进行过关运行高中资料试卷技术指导。对于调试过程中高中资料试卷技术问题，作为调试人员，需要在事前掌握图纸资料、设备制造厂家出具高中资料试卷试验报告与相关技术资料，并且了解现场设备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。5、电气设备调试高中资料试卷技术电力保护装置调试技术，电力保护高中资料试卷配置技术是指机组在进行继电保护高中资料试卷总体配置时，需要在最大限度内来确保机组高中资料试卷安全，并且尽可能地缩小故障高中资料试卷破坏范围，或者对某些异常高中资料试卷工况进行自动处理，尤其要避免错误高中资料试卷保护装置动作，并且拒绝动作，来避免不必要高中资料试卷突然停机。因此，电力高中资料试卷保护装置调试技术，要求电力保护装置做到准确灵活。对于差动保护装置高中资料试卷调试技术是指发电机一变压器组在发生内部故障时，需要进行外部电源高中资料试卷切除从而采用高中资料试卷主要保护装置。3ct(:,1)=[0;0;0;0;0];%初始角度p0=[1;0;-0.6;-0.4;-0.6];%各节点的输入有功功率q0=[0.5;0;-0.3;-0.1;-0.2];%各节点的输入无功功率x(:,1)=[0;1;0;1;0;1;0];%控制变量的初始值e=0.0001;%迭代精度zu(:,1)=[1;1;1;1];%最优控制变量解对应为P1、U1、U2、Tfork=1:1000fori=1:5a=0;forj=1:5a=u(j,k)*(g(i,j)*cos(ct(i,k)-ct(j,k))+b(i,j)*sin(ct(i,k)-ct(j,k)))+a;enddp(i,k)=p0(i)-u(i,k)*a;endfori=3:5z=0;forj=1:5z=u(j,k)*(g(i,j)*sin(ct(i,k)-ct(j,k))-b(i,j)*cos(ct(i,k)-ct(j,k)))+z;enddq(i,k)=q0(i)-u(i,k)*z;endfori=1:5forj=1:5ifi~=jhh(i,j)=u(i,k)*u(j,k)*(g(i,j)*sin(ct(i,k)-ct(j,k))-b(i,j)*cos(ct(i,k)-ct(j,k)));jj(i,j)=-u(i,k)*u(j,k)*(g(i,j)*cos(ct(i,k)-ct(j,k))+b(i,j)*sin(ct(i,k)-ct(j,k)));nn(i,j)=u(i,k)*u(j,k)*(g(i,j)*cos(ct(i,k)-ct(j,k))+b(i,j)*sin(ct(i,k)-ct(j,k)));ll(i,j)=u(i,k)*u(j,k)*(g(i,j)*sin(ct(i,k)-ct(j,k))-b(i,j)*cos(ct(i,k)-ct(j,k)));elsec=0;d=0;form=1:5ifm~=ic=u(m,k)*(g(i,m)*sin(ct(i,k)-ct(m,k))-b(i,m)*cos(ct(i,k)-ct(m,k)))+c;d=u(m,k)*(g(i,m)*cos(ct(i,k)-ct(m,k))+b(i,m)*sin(ct(i,k)-ct(m,k)))+d;endendhh(i,i)=-u(i,k)*c;jj(i,i)=u(i,k)*d;nn(i,i)=u(i,k)*d+2*u(i,k)^2*g(i,i);ll(i,i)=u(i,k)*c-2*u(i,k)^2*b(i,i);endendendjkb=[hh(3,3)nn(3,3)hh(3,4)nn(3,4)hh(3,5)nn(3,5)hh(3,1);3、管路敷设技术通1过管线0不仅可以解决吊顶层配置不规范高中资料试卷问题22，而且可保障各类管路习题到位。在管路敷设过程1中，要加强看2222222222护1关于管路高中资料试卷连接管口处理高中资料试卷弯扁度固定盒位置.保护层防腐跨接地线弯曲半径标高等，要求技术交底。管线敷设技术中包含线槽、管架等多项方式，为解决高中语文电气课件中管壁薄、接口不严等问题，合理利用管线敷设技术。线缆敷设原则：在分线盒处，当不同电压回路交叉时，应采用金属隔板进行隔开处理；同一线槽内，强电回路须同时切断习题电源，线缆敷设完毕，要进行检查和检测处理。4、电气课件中调试技术对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料试卷相互作用与相互关系，根据生产工艺高中资料试卷要求，对电气设备进行空载与带负荷下高中资料试卷调控试验；对设备进行调整使其在正常工况下与过度工作下都可以正常工作；对于继电保护进行整核对定值，审核与校对图纸，编写复杂设备与装置高中资料试卷调试方案，编写重要设备高中资料试卷试验方案以及系统启动方案；对整套启动过程中高中资料试卷电气设备进行调试工作并且进行过关运行高中资料试卷技术指导。对于调试过程中高中资料试卷技术问题，作为调试人员，需要在事前掌握图纸资料、设备制造厂家出具高中资料试卷试验报告与相关技术资料，并且了解现场设备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。5、电气设备调试高中资料试卷技术电力保护装置调试技术，电力保护高中资料试卷配置技术是指机组在进行继电保护高中资料试卷总体配置时，需要在最大限度内来确保机组高中资料试卷安全，并且尽可能地缩小故障高中资料试卷破坏范围，或者对某些异常高中资料试卷工况进行自动处理，尤其要避免错误高中资料试卷保护装置动作，并且拒绝动作，来避免不必要高中资料试卷突然停机。因此，电力高中资料试卷保护装置调试技术，要求电力保护装置做到准确灵活。对于差动保护装置高中资料试卷调试技术是指发电机一变压器组在发生内部故障时，需要进行外部电源高中资料试卷切除从而采用高中资料试卷主要保护装置。4jj(3,3)ll(3,3)jj(3,4)ll(3,4)jj(3,5)ll(3,5)jj(3,1);hh(4,3)nn(4,3)hh(4,4)nn(4,4)hh(4,5)nn(4,5)hh(4,1);jj(4,3)ll(4,3)jj(4,4)ll(4,4)jj(4,5)ll(4,5)jj(4,1);hh(5,3)nn(5,3)hh(5,4)nn(5,4)hh(5,5)nn(5,5)hh(5,1);jj(5,3)ll(5,3)jj(5,4)ll(5,4)jj(5,5)ll(5,5)jj(5,1);hh(1,3)nn(1,3)hh(1,4)nn(1,4)hh(1,5)nn(1,5)hh(1,1)];%牛拉法的雅各比矩阵dd=[dp(3,k);dq(3,k);dp(4,k);dq(4,k);dp(5,k);dq(5,k);dp(1,k)];dx(:,k)=inv(jkb)*dd;xzh=[1000000;0u(3,k)00000;0010000;000u(4,k)000;0000100;00000u(5,k)0;0000001];x(:,k+1)=x(:,k)+xzh*dx(:,k);%潮流计算解u(:,k+1)=[u(1,k);u(2,k);x(2,k+1);x(4,k+1);x(6,k+1)];ct(:,k+1)=[x(7,k+1);0;x(1,k+1);x(3,k+1);x(5,k+1)];f1=0;fori=1:5f1=u(j,k+1)*(g(1,j)*cos(ct(1,k+1)-ct(j,k+1))+b(1,j)*sin(ct(1,k+1)-ct(j,k+1)))+f1;endpg1=u(1,k+1)*f1;f2=0;fori=1:5f2=u(j,k+1)*(g(2,j)*cos(ct(2,k+1)-ct(j,k+1))+b(2,j)*sin(ct(2,k+1)-ct(j,k+1)))+f2;endpg2=u(2,k+1)*f2;fdx=[(702*pg1+100)*u(1,k+1)*u(3,k+1)*(g(1,3)*sin(ct(1,k+1)-ct(3,k+1))-b(1,3)*cos(ct(1,k+1)-ct(3,k+1)))+(778*pg2+100)*u(2,k+1)*u(3,k+1)*(g(2,3)*sin(ct(2,k+1)-ct(3,k+1))-b(2,3)*cos(ct(2,k+1)-ct(3,k+1)));(702*pg1+100)*u(1,k+1)*(g(1,3)*cos(ct(1,k+1)-ct(3,k+1))+b(1,3)*sin(ct(1,k+1)-ct(3,k+1)))+(778*pg2+100)*u(2,k+1)*(g(2,3)*cos(ct(2,k+1)-ct(3,k+1))+b(2,3)*sin(ct(2,k+1)-ct(3,k+1)));(702*pg1+100)*u(1,k+1)*u(4,k+1)*(g(1,4)*sin(ct(1,k+1)-ct(4,k+1))-b(1,4)*cos(ct(1,k+1)-ct(4,k+1)))+(778*pg2+100)*u(2,k+1)*u(4,k+1)*(g(2,4)*sin(ct(2,k+1)-ct(4,k+1))-b(2,4)*cos(ct(2,k+1)-ct(4,k+1)));(702*pg1+100)*u(1,k+1)*(g(1,4)*cos(ct(1,k+1)-ct(4,k+1))+b(1,4)*sin(ct(1,k+1)-ct(4,k+1)))+(778*pg2+100)*u(2,k+1)*(g(2,4)*cos(ct(2,k+1)-ct(4,k+1))+b(2,4)*sin(ct(2,k+1)-ct(4,k+1)));(702*pg1+100)*u(1,k+1)*u(