资源描述:
《matlab牛顿拉夫逊法算潮流分析(matlab newton raphson method for power flow analysis)》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、matlab牛顿拉夫逊法算潮流分析(MatlabNewtonRaphsonmethodforpowerflowanalysis)%ofthisprocedureisthefunctionoftheNewtonRaphsonmethodforpowerflowcalculation%B1matrix:1,branchheadnumber,2,terminalnumber;3,branchimpedance;4,linetogroundelectricity(ortransformeradmittance);%5,branchratio;6,branchheadendKsidein
2、1,the1sideis0;Thetransformerparametersof%7andline/transformeridentification(0/1)arenormalizedtotheterminalsidewhenthefirstbranchofthebranchisattheKside,andthe0iscalculatedtothefirstsideofthe1side%B2matrix:1,thenode,thenodepowergenerator;2loadpower;%3PQ,theinitialvalueofthegivennodevoltage,or
3、balancenodeandPVnodevoltagevalue%4,thenodereactivepowercompensationshuntcapacitor(sense)ofthesusceptance%5,nodeclassificationlabel:1isabalancednode(shouldbenode1);2isPQnode;3isPVnode;Clear;Isb=1;%input("pleaseinputbalancebusnumber:isb=');Pr=1e-5;%input("pleaseinputerror:pr=');%--------------
4、-------------------------------------N=10;%input("enternodenumber:n=');Nl=10;%input("inputbranchnumber:nl=');B1=[123.4+12.8i1.4e-4i100;145.1+19.2i2.1e-4i100;234.25+16i1.75e-4i100;454.25+16i7e-4i100;135.1+19.2i2.1e-4i100;645.95+22.4i9.8e-4i100;271.78+53.89i038.5/23101;381.49+48.02i011/23101
5、;491.49+48.02i011/23101;5102.46+70.17i038.5/23101]%input("pleaseinputbythebranchparametermatrixformed:B1=');B2=[0225.501;0022002;0022002;0022002;0022002;120023103;061.11+37.87i3502;047.53+29.46i1002;054.32+33.66i1002;040.74+25.25i3502]%input("pleaseinputparametersofeachnodematrixformed:B2=')
6、;%-------------------------------------------------------------%n=4;%input("enternodenumber:n=');%nl=4;%input("inputbranchnumber:nl=');%B1=[124+16i0100;%134+16i0100;%232+8i0100;%241.49+48.02i011/11001]%input("pleaseinputbythebranchparametermatrixformed:B1=');%B2=[0011501;%0011002;%020+4i11
7、002;%010+6i1002]%input("pleaseentereachnodeparametermatrixformed:B2=');%-------------------------------------------------------------Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(NL);%%%-----------admittancematrix------