网站大量收购闲置独家精品文档,联系QQ:2885784924

实验三源码.doc

  1. 1、本文档共12页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验三源码

源码如下: disp(注意:输入“var=1时选择三相短路电流计算;输入“var=2“时选用计算曲线计算三相短路电流;输入其他数字将选择不对称故障计算。); var=input(var=); switch var case 1 % 三相短路电流计算 NF=input(请输入短路点的数目:NF=); n=input(请输入节点数:n=); nI=input(请输入支路数:nI=); B=input(请输入支路参数形成的矩阵:B=); V0=input(请输入由各节点的初电压标幺值形成的列矩阵:V0=); D=input(请输入由短路号、短路点阻抗组成的矩阵:D=); m=0;Z=zeros(n);V=zeros(n);I=zeros(nI); for k1=1:nI; p=B(k1,1);q=B(k1,2); if B(k1,6)==0 k=1./B(k1,5); else k=B(k1,5); end if p==0 if qm Z(q,q)=B(k1,3);m=m+1; else for i=1:m, Z(i,m+1)=-Z(i,q);Z(m+1,i)=-Z(q,i); end Z(m+1,m+1)=Z(q,q)+B(k1,3); for i=1:m for j=1:m Z(i,j)=Z(i,j)-Z(i,m+1)*Z(m+1,j)./Z(m+1,m+1); end Z(i,m+1)=0; end for i=1:m+1 Z(m+1,i)=0; end end else if qm for i=1:m Z(i,q)=Z(i,p)*k;Z(q,i)=Z(p,i)*k; end Z(q,q)=k^2*Z(p,p)+k^2*B(k1,3); m=m+1; else for i=1:m Z(i,m+1)=k*Z(i,p)-Z(i,q); Z(m+1,i)=k*Z(p,i)-Z(q,i); end Z(m+1,m+1)=k^2*Z(p,p)+Z(q,q)-2*k*Z(p,q)+k^2*B(k1,3); for i=1:m for j=1:m Z(i,j)=Z(i,j)-Z(i,m+1)*Z(m+1,j)./Z(m+1,m+1); end Z(i,m+1)=0; end for i=1:m+1 Z(m+1,i)=0; end end end end for k=1:NF I(D(k,1),D(k,1))=V0(D(k,1),1)./(Z(D(k,1),D(k,1))+D(k,2)); ft=num2str(D(k,1)); ts1=(点短路时); ts2=(电流的标幺值 If=); dn=strcat(ft,ts1,ts2); disp(dn); disp(I(D(k,1),D(k,1))); for i=1:n V(i,i)=V0(i,1)-I(D(k,1),D(k,1))*Z(i,D(k,1)); end for i=1:nI if B(i,6)==0 k=B(i,5); else k=1./B(i,5); end p=B(i,1);q=B(i,2); if p==0 e=0;b=B(i,3); I(i,i)=(e-V(q,q)

文档评论(0)

haihang2017 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档