潮流计算报告.docx
《潮流计算报告.docx》由会员分享,可在线阅读,更多相关《潮流计算报告.docx(37页珍藏版)》请在咨信网上搜索。
1、一、系统结构图:二、网络参数:1、支路参数:网络类型支路编号电压等级(kV)装机容量(MW)导线技术参数()()()环网1-222010013.6125.567.851-38.321130.552.243-510.2128.874.992-38.5105.428.361-47.579129.651.454-513.84125.312.78辐射网1-2242-3363-4482-5125-6442、节点参数:节点类型节点编号发电功率(MW)负荷视在功率环网10未知(平衡节点)21000(PV节点)3015+9.4i4027+6i5035.5+25.5i辐射网10未知(平衡节点)204+2i306
2、+3.2i403+1.44i504+3.2i602+1.1i三、时尚计算步骤图:四、matlab程序:clear;Un=input(请输入Un:); %输入所需额定电压PQ=%节点电压 有功 无功Un 0 0Un 4 2 Un 6 3.2Un 3 1.44Un 4 3.2Un 2 1.1 ;FT=%首端 末端 4 33 26 55 22 1;RX=% R X 4 83 64 41 2 2 4; NN=size(PQ,1); %节点数NB=size(FT,1); %支路数数V=PQ(:,1); %V初始电压相量maxd=1k=1while maxd0.0001 k=k+1; PQ2=PQ; %每
3、一次迭代各节点注入有功和无功相同 PL=0.0;for i=1:NB kf=FT(i,1); %前推始节点号 kt=FT(i,2); %前推终节点号 x=(PQ2(kf,2)2+PQ2(kf,3)2)/V(kf)/V(kf);%计算沿线电流 /平方A losss(i,1)=RX(i,1)*x; %计算线路有功损耗 /MW losss(i,2)=RX(i,2)*x; %计算线路无功损耗/MW PQ1(i,1)=PQ2(kf,2)+RX(i,1)*x; %计算支路首端有功/MW RX(i,1)*R PQ1(i,2)=PQ2(kf,3)+RX(i,2)*x; %计算沿支路无功/MW RX(i,2)*
4、X PQ2(kt,2)= PQ2(kt,2)+PQ1(i,1); %用PQ1去修正支路末端节点有功P 单位MW PQ2(kt,3)= PQ2(kt,3)+PQ1(i,2); %用PQ1去修正支路末端节点有功Q 单位Mvarend angle(1)=0.0;for i=NB:-1:1 kf=FT(i,2); %回代始节点号 kt=FT(i,1); %回代终节点号 dv1=(PQ1(i,1)*RX(i,1)+PQ1(i,2)*RX(i,2)/V(kf); %计算支路电压损耗纵分量dv1 dv2=(PQ1(i,1)*RX(i,2)-PQ1(i,2)*RX(i,1)/V(kf); %计算支路电压损耗横
5、分量dv2 V2(kt)=sqrt(V(kf)-dv1)2+dv22); %计算支路末端电压/kV angle(kt)=angle(kf)+atand(dv2/(V(kf)-dv1); %计算支路endmaxd=abs(V2(2)-V(2);V2(1)=V(1);for i=3:1:NN if abs(V2(i)-V(i)maxd; maxd=abs(V2(i)-V(i); end end fullloss(1,1)=0;%计算线路总损耗 fullloss(1,2)=0; finalPQ=max(PQ1); for i=1:NB fullloss(1,1)=fullloss(1,1)+loss
6、s(i,1); fullloss(1,2)=fullloss(1,2)+losss(i,2); enddisp(辐射网迭代次数:)kdisp(辐射网系统电压差精度:) maxddisp(辐射网系统末端节点有功和无功:)finalPQ %时尚分布 即支路首端时尚MVAdisp(辐射网系统总功率损耗:)fullloss %线路总损耗MVAdisp(辐射网系统各支路功率损耗:)losss %各支路损耗MVAdisp(辐射网系统各节点电压幅值:)V=V2 %节点电压模计算结果kVdisp(辐射网系统各节点电压相角:)angle %节点电压角度计算结果单位度endclcdisp(辐射网迭代次数:)kdi
7、sp(辐射网系统电压差精度:) maxddisp(辐射网系统末端节点有功和无功/MVA:)FinPQ=finalPQ(1,1)+finalPQ(1,2)*j %时尚分布 即支路首端时尚MVAdisp(辐射网系统总功率损耗/MVA:)Fulloss=fullloss (1,1)+fullloss(1,2)*j %线路总损耗MVAdisp(辐射网系统各支路功率损耗/MVA:)for(a=1:5)LOSS=losss (a,1)+losss(a,2)*j %各支路损耗MVAenddisp(辐射网系统各节点电压幅值/KV:)V=V2 %节点电压模计算结果kVdisp(辐射网系统各节点电压相角:)ang
8、le %节点电压角度计算结果单位度n=5; %input(节点数);nl=6; %input(支路数);isb=1; %input(平衡母线节点号);pr=0.000001; %input(误差精度:pr=);B1=1,2,13.6+125.5i,0.00006785i,1,0;1,3,8.321+130.5i,0.00005224i,1,0;3,5,10.2+128.8i,0.00007499i,1,0;2,3,8.5+105.4i,0.00002836i,1,0;1,4,7.579+129.6i,0.00005145i,1,0;4,5,13.84+125.31i,0.0000278i,1,
9、0; %input(由支路参数形成矩阵);B2=-FinPQ,0,Un,0,0,1;100,0,Un,Un,0,3;0,15+9.4i,Un,0,0,2;0,27+6i,Un,0,0,2;0,35.5+25.5i,Un,0,0,2; %input(各节点参数形成矩阵);Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl);%对各矩阵置零%修改部分ym=1;SB=100;UB=Un; %定义视在功率和电压基值if ym=0 %若不是标幺值 YB=SB./UB./UB; %定义导纳标幺值 BB1
10、=B1; BB2=B2; for i=1:nl B1(i,3)=B1(i,3)*YB; %切换为阻抗标幺值 B1(i,4)=B1(i,4)./YB; %切换为导纳标幺值 end disp(支路矩阵B1=); sparseB1=sparse(B1); disp(sparseB1) %输出标幺值稀疏矩阵B1 disp(-); for i=1:n B2(i,1)=B2(i,1)./SB; %切换为视在功率标幺值 B2(i,2)=B2(i,2)./SB; %切换为视在功率标幺值 B2(i,3)=B2(i,3)./UB; %切换为电压标幺值 B2(i,4)=B2(i,4)./UB; %切换为电压标幺值
11、B2(i,5)=B2(i,5)./SB; %切换为视在功率标幺值 end disp(节点矩阵B2=); sparseB2=sparse(B2); disp(sparseB2) %输出标幺值稀疏矩阵B2enddisp(-);% % %-for i=1:nl %支路数 if B1(i,6)=0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); %使左节点处于低压侧 end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5); %求解非对角元导纳 Y(q,p)=Y(p,q); %对角元两侧对称 Y(q,q)=Y(q,q)+
12、1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2; %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元1侧 end%求导纳矩阵disp(导纳矩阵 Y=);sparseY=sparse(Y);disp(sparseY) %输出导纳稀疏矩阵disp(-);%-G=real(Y);B=imag(Y); %分解出导纳阵实部和虚部 for i=1:n e(i)=real(B2(i,3); %给定i节点初始电压实部 f(i)=imag(B2(i,3); %给定i节点初始电压虚部 V(i)=B2(i,4); %PV节点电压给定模值 endfor
13、i=1:n %给定各节点注入功率 S(i)=B2(i,1)-B2(i,2); %i节点注入功率SG-SL B(i,i)=B(i,i)+B2(i,5); %i节点无功赔偿量 end%=P=real(S);Q=imag(S); %定义有功功率和无功功率ICT1=0;IT2=1;N0=2*n;N=N0+1;a=0; %定义迭代次数ICT1和不满足精度要求节点个数IT2while IT2=0 %仍有不满足精度要求节点 IT2=0;a=a+1; %IT2置零 for i=1:n if i=isb %非平衡节点 C(i)=0;D(i)=0; for j1=1:n C(i)=C(i)+G(i,j1)*e(j
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 潮流 计算 报告
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【精****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【精****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。