电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

电力系统三相短路计算的MATLAB代码

电力系统三相短路计算的MATLAB代码_第1页
1/9
电力系统三相短路计算的MATLAB代码_第2页
2/9
电力系统三相短路计算的MATLAB代码_第3页
3/9
电力系统三相短路计算amain.mclear 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,但是支路未近似计算的节点导纳矩阵;Ymatrix3是近似变比为1,实行近似支路参数1的导纳矩阵;Ymatrix4是近似变比为1,实行近似支路参数2的导纳矩阵.Y = Ymatrix2(bus,line);%对故障点进行导纳修正fixY = FixY(Y,bus,fault);%求注入电流Iinj = Inode(bus,calcSettings);U = fixY\Iinj;%得到故障支路与其他支路电流Bcurrent = Ibranch( line,U,fault,Y );%假如发生支路三相短路,那么对应该支路的电流修正为-999999-j999999Ib = ReviseBcurrent( fault,Bcurrent );%结束计时tim=toc;fprintf(’\n程序运行结果');fprintf(’\n计算完成,共用时%4.4fs,相关结果已保存在output.dat\n’,tim);%输出结果fprintf_result(output_file, Ib);fprintf_result1(Ib);bFixY.mfunction 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,1),bus(k,1)) = fixY(bus(k,1),bus(k,1)) + 1/1i/bus(k,8); endend%对节点短路和支路短路的导纳矩阵进行修正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) = 999999+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、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

电力系统三相短路计算的MATLAB代码

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部