三相短路短路电流计算matlab程序文档.pdf
《三相短路短路电流计算matlab程序文档.pdf》由会员分享,可在线阅读,更多相关《三相短路短路电流计算matlab程序文档.pdf(25页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、 电力系统三相短路计算 a main.m clear tim%打开文件 dfile,pathname=uigetfile(*.m ,Select Data File);if pathname=0 error(you must select a valid data file)else lfile=length(dfile);eval(dfile(1:lfile-2);end%定义输出文件 output_file=fopen(output.dat,w);%开始计时 tic;%求解节点导纳矩阵,其中 Ymatrix1 是考虑了变比,且支路未近似的导纳矩阵;Ymatrix2 是近似变比为 1,但是支
2、路未近似计算的节点导纳矩阵;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-j999999 Ib=ReviseBcurrent(fault,Bcurrent)
3、;%结束计时 tim=toc;fprintf(程序运行结果 );fprintf(计算完成,共用时%4.4fs,相关结果已保存在 output.dat,tim);%输出结果 fprintf_result(output_file,Ib);fprintf_result1(Ib);b FixY.m function fixY=FixY(Y,bus,fault)%对形成的导纳矩阵进行故障点的修正 nb,mb=size(bus);nf,mf=size(fault);fixY=Y;%对发电机节点导纳修正 for k=1:nb busType=bus(k,7);if(busType=1)fixY(bus(k,
4、1),bus(k,1)=fixY(bus(k,1),bus(k,1)+1/1i/bus(k,8);end end%对节点短路和支路短路的导纳矩阵进行修正 for k=1:nf nodeI=fault(k,1);nodeJ=fault(k,2);dis=fault(k,3);if(nodeI=0)fixY(nodeJ,nodeJ)=999999+1i*999999;continue;end if(nodeJ=0)fixY(nodeI,nodeI)=999999+1i*999999;continue;end if(dis=0)&(nodeI*nodeJ=0)fixY(nodeI,nodeI)=99
5、9999+1i*999999;continue;end if(dis=1)&(nodeI*nodeJ=0)fixY(nodeJ,nodeJ)=999999+1i*999999;continue;end if(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 en
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 三相 短路 电流 计算 matlab 程序 文档
限制150内