![三相短路短路电流计算matlab程序文档_第1页](http://file4.renrendoc.com/view/3b8f6a1b60c1553a010464bf41c551d9/3b8f6a1b60c1553a010464bf41c551d91.gif)
![三相短路短路电流计算matlab程序文档_第2页](http://file4.renrendoc.com/view/3b8f6a1b60c1553a010464bf41c551d9/3b8f6a1b60c1553a010464bf41c551d92.gif)
![三相短路短路电流计算matlab程序文档_第3页](http://file4.renrendoc.com/view/3b8f6a1b60c1553a010464bf41c551d9/3b8f6a1b60c1553a010464bf41c551d93.gif)
![三相短路短路电流计算matlab程序文档_第4页](http://file4.renrendoc.com/view/3b8f6a1b60c1553a010464bf41c551d9/3b8f6a1b60c1553a010464bf41c551d94.gif)
![三相短路短路电流计算matlab程序文档_第5页](http://file4.renrendoc.com/view/3b8f6a1b60c1553a010464bf41c551d9/3b8f6a1b60c1553a010464bf41c551d95.gif)
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
电力系统三相短路计算main.macleartim%打开文件[dfile,pathname]=uigetfile(ifpathname==0'*.m','SelectDataFile');error('youmustselectavaliddatafile')elselfile=length(dfile);eval(dfile(1:lfile-2));end%定义输出文件output_file=fopen('output.dat','w');%开始计时tic;
%求解节点导纳矩阵,其中Ymatrix1是考虑了变比,且支路未近似的导纳矩阵;Ymatrix2是近似变比为1,但是支路未近似计算的节点导纳矩阵;Ymatrix3是近似变比为1,采取近似支路参数1的导纳矩阵;Ymatrix4是近似变比为1,采取近似支路参数2的导纳矩阵。Y=Ymatrix2(bus,line);%对故障点进行导纳修正fixY=FixY(Y,bus,fault);%求注入电流Iinj=Inode(bus,calcSettings);%得到故障支路与其他支路电流Bcurrent=Ibranch(line,U,fault,Y);%如果发生支路三相短路,那么对应该支路的电流修正为-999999-j999999Ib=ReviseBcurrent(fault,Bcurrent);%结束计时tim=toc;
fprintf('
程序运行结果');fprintf('
计算完成,共用时%4.4fs,相关结果已保存在output.dat
',tim);%输出结果fprintf_result(output_file,Ib);fprintf_result1(Ib);FixY.mbfunctionfixY=FixY(Y,bus,fault)%对形成的导纳矩阵进行故障点的修正[nb,mb]=size(bus);[nf,mf]=size(fault);fixY=Y;%对发电机节点导纳修正fork=1:nbbusType=bus(k,7);if(busType==1)
fixY(bus(k,1),bus(k,1))=fixY(bus(k,1),bus(k,1))+1/1i/bus(k,8);endend%对节点短路和支路短路的导纳矩阵进行修正fork=1:nfnodeI=fault(k,1);nodeJ=fault(k,2);dis=fault(k,3);if(nodeI==0)fixY(nodeJ,nodeJ)=999999+1i*999999;continue;endif(nodeJ==0)fixY(nodeI,nodeI)=999999+1i*999999;continue;
endif(dis==0)&&(nodeI*nodeJ~=0)fixY(nodeI,nodeI)=999999+1i*999999;continue;endif(dis==1)&&(nodeI*nodeJ~=0)fixY(nodeJ,nodeJ)=999999+1i*999999;continue;endif(dis~=1)&&(dis~=0)&&(nodeI*nodeJ~=0)fixY(nodeI,nodeI)=fixY(nodeI,nodeI)-fixY(nodeI,nodeJ)/dis;fixY(nodeJ,nodeJ)=fixY(nodeJ,nodeJ)-fixY(nodeI,nodeJ)/(1-dis);fixY(nodeI,nodeJ)=0;fixY(nodeJ,nodeI)=0;end
endendfprintf_result.mcfunction[output_args]=fprintf_result(output_file,Ib)%将得到的短路电流输入到输出文件中[n,m]=size(Ib);fprintf(output_file,'No.No.vectorofIvalueofI
');fork=1:nI=Ib(k,1);J=Ib(k,2);I01=real(Ib(k,3));I02=imag(Ib(k,3));I1=Ib(k,4);
if(I02>=0)fprintf(output_file,'%3d%3d%10.6f+j%10.6f%10.6f',I,J,I01,I02,I1);endif(I02<0)I02=abs(I02);fprintf(output_file,'%3d%3d%10.6f-j%10.6f%10.6f',I,J,I01,I02,I1);endfprintf(output_file,'
');endendfprintf_result1.mdfunction[output_args]=fprintf_result1(Ib)%UNTITLED???ú′?′|ê?è?oˉêy??òa
[n,m]=size(Ib);fprintf('No.No.vectorofIvalueofI
');fork=1:nI=Ib(k,1);J=Ib(k,2);I01=real(Ib(k,3));I02=imag(Ib(k,3));I1=Ib(k,4);if(I02>=0)fprintf('%3d%3d%10.6f+j%10.6f%10.6f',I,J,I01,I02,I1);endif(I02<0)I02=abs(I02);fprintf('%3d%3d%10.6f-j%10.6f%10.6f',I,J,I01,I02,I1);end
fprintf('
');endendIbranch.mefunctionBcurrent=Ibranch(line,U,fault,Y)%计算短路电流%记录短路故障参数,如短路节点,如为支路短路,记录距离节点的距离%此段计算采用的支路参数未近似,如果计算近似的时候需要修改[nl,ml]=size(line);Bcurrent=zeros(nl+1,4);faultI=fault(1,1);faultJ=fault(1,2);dis=fault(1,3);faultNode=0;if(faultI==0)
faultNode=faultJ;endif(faultJ==0)faultNode=faultI;endif(dis==1)&&(faultI*faultJ~=0)faultNode=faultJ;endif(dis==0)&&(faultI*faultJ~=0)faultNode=faultI;endif(faultNode~=0)Bcurrent(nl+1,1)=faultNode;Bcurrent(nl+1,2)=faultNode;Iij=0;
Iij1=0;end%计算非故障支路的短路电流fork=1:nli=line(k,1);j=line(k,2);Ui=U(i);ifj~=0Uj=U(j);elseUj=0;endifline(k,2)==0Ym=line(k,5)+1i*line(k,6);Iij=Ui*Ym;
Iij1=abs(Iij);endifline(k,2)~=0Zt=line(k,3)+1i*line(k,4);Yt=1/Zt;Ym=line(k,5)+1i*line(k,6);Iij=(Ui-Uj)*Yt+Ui*Ym;Iij1=abs(Iij);endBcurrent(k,1)=i;Bcurrent(k,2)=j;Bcurrent(k,3)=Iij;Bcurrent(k,4)=Iij1;end%如果为节点短路,修正短路点的电流大小
if(faultNode~=0)Bcurrent(nl+1,1)=faultNode;Bcurrent(nl+1,2)=faultNode;Ifault=0;branchCurrent=0;fork=1:nlI=line(k,1);J=line(k,2);if(I*J==0)continue;endbranchCurrent=(U(I)-U(J))/(line(k,3)+1i*line(k,4));if(I==faultNode)Ifault=Ifault-branchCurrent;elseif(J==faultNode)
Ifault=Ifault+branchCurrent;endendendBcurrent(nl+1,3)=Ifault;Bcurrent(nl+1,4)=abs(Bcurrent(nl+1,3));end%如果为支路短路,修正短路支路的短路电流大小if(dis~=0)&&(dis~=1)&&(faultI*faultJ~=0)Bcurrent(nl+1,1)=faultI;Bcurrent(nl+1,2)=faultJ;Bcurrent(nl+1,3)=U(faultI)*Y(faultI,faultJ)/dis+U(faultJ)*Y(faultI,faultJ)/(1-dis);Bcurrent(nl+1,4)=abs(Bcurrent(nl+1,3));endend
Inode.mffunctionIinj=Inode(bus,calcSettings)%计算节点注入电流[nb,mb]=size(bus);Iinj=zeros(nb,1);fork=1:nbbusType=bus(k,7);if(calcSettings(1)==1)v=1;elsev=bus(k,2);end%对发电机节点电流进行修正if(busType==1)Iinj(bus(k,1),1)=Iinj(bus(k,1),1)+v/1i/bus(k,8);
endendendReviseBcurrent.mgfunctionIb=ReviseBcurrent(fault,Bcurrent)%如果发生支路短路,对原来的计算电流进行修正,使该支路短路电流输出为-999999-j999999clearfaultIfaultJdis[nt,mt]=size(Bcurrent);Ib=zeros(nt,mt);faultI=fault(1,1);faultJ=fault(1,2);dis=fault(1,3);fork=1:nt-1i=Bcurrent(k,1);
j=Bcurrent(k,2);Ib(k,:)=Bcurrent(k,:);if(faultI*faultJ~=0)&&(dis~=1)&&(dis~=0)&&(i==faultI)&&(j==faultJ)Ib(k,1)=i;Ib(k,2)=j;Ib(k,3)=-999999-1i*999999;Ib(k,4)=-999999;endif(faultI*faultJ~=0)&&(dis~=1)&&(dis~=0)&&(i==faultJ)&&(j==faultI)Ib(k,1)=i;Ib(k,2)=j;Ib(k,3)=-999999-1i*999999;Ib(k,4)=-999999;end
Ib(nt,:)=Bcurrent(nt,:);endYmatrix1.mhfunctionY=Ymatrix1(bus,line)%考虑变压器,并且支路参数不近似的节点导纳矩阵[nb,mb]=size(bus);[nl,ml]=size(line);Y=zeros(nb,nb);fork=1:nlI=line(k,1);J=line(k,2);Zt=line(k,3)+1i*line(k,4);Yt=1/Zt;Ym=line(k,5)+1i*line(k,6);K=line(k,7);
if(K==0)&&(J~=0)Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt+Ym;Y(I,J)=Y(I,J)-Yt;Y(J,I)=Y(I,J);endif(K==0)&&(J==0)Y(I,I)=Y(I,I)+Ym;endifK>0Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt/(K*K);Y(I,J)=Y(I,J)-Yt/K;Y(J,I)=Y(I,J);end
ifK<0Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+K*K*Yt;Y(I,J)=Y(I,J)+K*Yt;Y(J,I)=Y(I,J);endendendYmatrix2.mifunctionY=Ymatrix2(bus,line)%考虑变压器变比近似为1,支路参数不等效[nb,mb]=size(bus);[nl,ml]=size(line);Y=zeros(nb,nb);fork=1:nl
I=line(k,1);J=line(k,2);Zt=line(k,3)+1i*line(k,4);Yt=1/Zt;Ym=line(k,5)+1i*line(k,6);ifJ~=0Y(I,I)=Y(I,I)+Yt+Ym;Y(J,J)=Y(J,J)+Yt+Ym;Y(I,J)=Y(I,J)-Yt;Y(J,I)=Y(I,J);endifJ==0Y(I,I)=Y(I,I)+Ym;endend
endYmatrix3.mjf
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024版公司的限制性股权激励协议
- 全国较好的城市绿化经验做法
- 2024版购销油罐合同
- 2024年江苏客运驾驶员考试题及答案解析
- 2024年武汉客运上岗证考试多少道题
- 2024年山西客运运输从业资格证模拟考试试题
- 2024年西安客运资格考试考题
- 2024年阜新道路旅客运输驾驶员从业资格模拟试题
- 2023版高考地理二轮总复习第2部分技能提升篇技能2区域认知
- 房屋自装用工合同模板
- 丰收宝使用说明
- 招商银行智慧营销体系规划方案(2022年-2023年)
- T-AHMS 0003-2023 气象文物价值分类指南
- 给孩子玩手机的协议书
- 上海通用汽车有限公司
- 河北单招考试五类职业适应性测试试题+答案
- 线上商城微信小程序开发说明书
- von frey丝K值表完整版
- 汽车维修保养服务单
- GB/T 24690-2018袋泡茶
- GB/T 19867.1-2005电弧焊焊接工艺规程
评论
0/150
提交评论