解广润-电力系统过电压-集中参数等值计算例题编程.docVIP

下载本文档

23
0
约3.14千字
约 4页
2017-05-20 发布于河南
举报
版权申诉

解广润-电力系统过电压-集中参数等值计算例题编程.doc

1、本文档共4页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

解广润-电力系统过电压-集中参数等值计算例题编程

本程序是使用matlab编制的一个函数，用于计算解广润—《电力系统过电压》P66页电路例题，计算方法是集中参数电容等值计算电路法，在matlab环境下运行该函数，输入观察时间的步长（步长选0.2及以下的数值最为合理）和观察的持续时间，即可得到各节点电压随时间的变化的图像和数据。经过与解广润—《电力系统过电压》P80-82页的计算结果比对，计算结果完全一致。（等效电路和计算公式本文已给出）节点电压方程 i12(t)=U1(t)/300+I12(t-2) i21(t)=U2(t)/300+I21(t-2) i24(t)=U2(t)/300+I24(t-2) i42(t)=U4(t)/300+I42(t-2) i14(t)=U1(t)/300+I14(t-2) i41(t)=U4(t)/300+I41(t-2) i23(t)=U2(t)/300+I23(t-4) i32(t)=U3(t)/300+I32(t-4) I12(t-2)=-U2(t-2)/300-i21(t-2) I21(t-2)=-U1(t-2)/300-i12(t-2) I24(t-2)=-U4(t-2)/300-i42(t-2) I42(t-2)=-U2(t-2)/300-i24(t-2) I14(t-2)=-U4(t-2)/300-i41(t-2) I41(t-2)=-U1(t-2)/300-i14(t-2) I23(t-2)=-U3(t-4)/300-i32(t-4) I32(t-2)=-U2(t-4)/300-i23(t-4) function U=gdy s=input(请输入步长（μs），按回车键确认); T=input(请输入持续时间的上限（μs），按回车键确认); t=0;%t为时间 k=0;%k为迭代次数 for i=1:1:8 I(i,1)=0; end %I(1,1)至I(8,1)值依次为等效电路中的I12，I21,I24,I42,I14,I41,I23,I32 for i=1:1:4 for j=1:1:4 y(i,j)=0; end end y(1,1)=1/150; y(2,2)=1/100; y(3,3)=10000+1/300; y(4,4)=1/150; %y为导纳矩阵 for t=0:s:T k=k+1; if tT break end a(1)=t-2; a(2)=t-2; a(3)=t-2; a(4)=t-2; a(5)=t-2; a(6)=t-2; a(7)=t-4; a(8)=t-4; t=t+s; for i=1:1:6 if a(i)0 I(i,1)=0; else b=k-ceil(2/s); I(i,1)=II(i,b); end end for i=7:1:8 if a(i)0 I(i,1)=0; else b=k-ceil(4/s); I(i,1)=II(i,b); end end Ie(1,1)=-I(5,1)-I(1,1); Ie(2,1)=-I(7,1)-I(3,1)-I(2,1); Ie(3,1)=10000-I(8,1); Ie(4,1)=-I(4,1)-I(6,1); u=pinv(y)*Ie; u(3,1)=1; i12=u(1,1)/300+I(1,1); i21=u(2,1)/300+I(2,1); i24=u(2,1)/300+I(3,1); i42=u(4,1)/300+I(4,1); i14=u(1,1)/300+I(5,1); i41=u(4,1)/300+I(6,1); i23=u(2,1)/300+I(7,1); i32=u(3,1)/300+I(8,1); II(1,k)=(-1/300)*u(2,1)-i21; II(2,k)=(-1/300)*u(1,1)-i12; II(3,k)=(-1/300)*u(4,1)-i42; II(4,k)=(-1/300)*u(2,1)-i24; II(5,k)=(-1/300)*u(4,1)-i41; II(6,k)=(-1/300)*u(1,1)-i14; II(7,k)=(-1/300)*u(3,1)-i32; II(8,k)=(-1/300)*u(2,1)-i23; %II的时间变量为t，而I的时间变量为t-τ U(2,k)=u(1,1); U(3,k)=u(2,1); U(4,k)=u(3,1); U(5,k)=u(4,1); end for i=1:1:k U(1,i)=(i-1)*s; end for i=1:1:k ti(i)