潮流计算地MATLAB源程序

潮流计算地MATLAB源程序

ID:38401473

大小:33.50 KB

页数:9页

时间:2019-06-11

潮流计算地MATLAB源程序_第1页
潮流计算地MATLAB源程序_第2页
潮流计算地MATLAB源程序_第3页
潮流计算地MATLAB源程序_第4页
潮流计算地MATLAB源程序_第5页
资源描述:

《潮流计算地MATLAB源程序》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、实用标准文案%简单潮流计算的牛顿拉夫逊程序,相关的原始数据数据数据输入格式如下:%B1是支路参数矩阵,第一列和第二列是节点编号。节点编号由小到大编写%对于含有变压器的支路,第一列为低压侧节点编号,第二列为高压侧节点%编号,将变压器的串联阻抗置于低压侧处理。%第三列为支路的串列阻抗参数。%第四列为支路的对地导纳参数。%第五列为含变压器支路的变压器的变比%第六列为变压器是否含有变压器的参数,其中“1”为含有变压器,%“0”为不含有变压器。%B2为节点参数矩阵,其中第一列为节点注入发电功率参数;第二列为节点负荷功率

2、参数;第三列为节点电压参数;第六列为节点类型参数,其中“1”为平衡节点,“2”为PQ节点,“3”为PV节点参数。%X为节点号和对地参数矩阵。其中第一列为节点编号,第二列为节点对地%参数。n=input('请输入节点数:n=');n1=input('请输入支路数:n1=');isb=input('请输入平衡节点号:isb=');pr=input('请输入误差精度:pr=');B1=input('请输入支路参数:B1=');B2=input('请输入节点参数:B2=');X=input('节点号和对地参数:X='

3、);Y=zeros(n);Times=1;               %置迭代次数为初始值%创建节点导纳矩阵fori=1:n1   ifB1(i,6)==0      %不含变压器的支路       p=B1(i,1);       q=B1(i,2);       Y(p,q)=Y(p,q)-1/B1(i,3);       Y(q,p)=Y(p,q);       Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4);       Y(q,q)=Y(q,q)+1/B1(i,3)+0.5

4、*B1(i,4);   else               %含有变压器的支路       p=B1(i,1);       q=B1(i,2);       Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5));       Y(q,p)=Y(p,q);       Y(p,p)=Y(p,p)+1/B1(i,3);       Y(q,q)=Y(q,q)+1/(B1(i,5)^2*B1(i,3));   end文档大全实用标准文案endYOrgS=zeros(2*n-2,1);DetaS=z

5、eros(2*n-2,1);  %将OrgS、DetaS初始化%创建OrgS,用于存储初始功率参数h=0;j=0;fori=1:n           %对PQ节点的处理   ifi~=isb&B2(i,6)==2         h=h+1;       forj=1:n           OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3

6、))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));           OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));       end   endendfori

7、=1:n          %对PV节点的处理,注意这时不可再将h初始化为0   ifi~=isb&B2(i,6)==3       h=h+1;       forj=1:n           OrgS(2*h-1,1)=OrgS(2*h-1,1)+real(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))+imag(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(

8、B2(j,3)));           OrgS(2*h,1)=OrgS(2*h,1)+imag(B2(i,3))*(real(Y(i,j))*real(B2(j,3))-imag(Y(i,j))*imag(B2(j,3)))-real(B2(i,3))*(real(Y(i,j))*imag(B2(j,3))+imag(Y(i,j))*real(B2(j,3)));       end   end

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。