没有合适的资源?快使用搜索试试~ 我知道了~
潮流计算matlab.doc
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 107 浏览量
2022-07-05
14:17:09
上传
评论
收藏 119KB DOC 举报
温馨提示
试读
20页
潮流计算matlab
资源推荐
资源详情
资源评论
电力系统大作业——用 matlab 潮流计算
08292007 电气 0809 班
%主程序:
[dfile,pathname]=uigetfile('ieee14.m','Select Data File');
if pathname==0
error('you must select a valid data file')
else
lfile=length(dfile);
% strip off .m
eval(dfile(1:lfile-2)); %打开数据文件
end
global n;
global m;
[nb,mb]=size(bus); %节点重新编号
[nl,ml]=size(line);
nSW=0; % 平衡节点数目
nPV=0; % PV 节点数目
nPQ=0; % PQ 节点数目
for i=1:nb, % nb 为总节点数
type=bus(i,6);
if type==3,
nSW=nSW+1; % 统计平衡节点数目
SW(nSW,:)=bus(i,:);
elseif type==2,
nPV=nPV+1; % 统计 PV 节点数目
PV(nPV,:)=bus(i,:);
else
nPQ=nPQ+1; % 统计 PQ 节点数目
PQ(nPQ,:)=bus(i,:);
end
end;
bus=[PQ;PV;SW];
newbus=[1:nb]';
f=bus(:,1);
nodenum=[newbus bus(:,1)];
bus(:,1)=newbus;
for i=1:nl
for j=1:2
for k=1:nb
if line(i,j)==nodenum(k,2)
line(i,j)=nodenum(k,1);
break
end
end
end
end
Y=y(bus,line); %形成节点导纳矩阵
K=0; %迭代次数初值
Kmax=10; %最大迭代次数
eps1=1.0e-10;
eps2=1.0e-10;
m=nPQ;n=nb;
Um=eye(m,m);
myf=fopen('output1.dat','w');
for K=1:Kmax
for i=1:m
for j=1:m
if i==j
Um(i,j)=bus(i,2);
end
end
end
b=dPQ(Y,bus);
C=jac(bus,Y);
dX=C\b';
dx=dX';
[nx,mx]=size(dx);
for i=1:n-1 %计算相角
bus(i,3)=bus(i,3)-dX(i,1);
end
B=dx(nx,n:mx)*Um; %计算电压差
bus(1:m,2)=bus(1:m,2)-B'; %计算电压值
dx(nx,n:mx)=B;
fprintf(myf,'--第%d 次迭代时雅可比矩阵--',K)
fprintf(myf, '\n');
for i=1:(n+m-1)
for j=1:(n+m-1)
fprintf(myf,'%8.6f', C(i,j));
fprintf(myf, ' ');
end
fprintf(myf, '\n');
end
fprintf(myf,'--第%d 次迭代时 dPQ 的误差--',K)
fprintf(myf, '\n');
for i=1:(n+m-1)
fprintf(myf,'%8.6e', b(1,i));
fprintf(myf, '\n');
end
fprintf(myf, '\n');
fprintf(myf,'--第%d 次迭代时 dx(误差)--',K)
fprintf(myf, '\n');
for i=1:(n+m-1)
fprintf(myf,'%8.6e', dX(i,1));
fprintf(myf, '\n');
end
fprintf(myf, '\n');
fprintf(myf,'第%d 次迭代后节点电压(仅 PQ 节点)',K)
fprintf(myf, '\n');
for i=1:m
fprintf(myf,'%8.6f', bus(i,2));
fprintf(myf, '\n');
end
fprintf(myf, '\n');
fprintf(myf,'第%d 次迭代后相角(角度)',K)
fprintf(myf, '\n');
for i=1:n
fprintf(myf,'%8.6f', bus(i,3)*180/pi);
fprintf(myf, '\n');
end
fprintf(myf, '\n');
if (max(abs(dx))<eps1)&(max(abs(b))<eps2) %判断是否达到计算精度
break;
end
end
%计算功率
P1=0;T=0; %计算平衡节点的有功
for j=1:n
T=bus(n,2)*bus(j,2)*(real(Y(n,j))*cos(bus(n,3)-bus(j,3))+imag(Y(n,j))*sin(bus(n,3)-bus(j,3)));
P1=P1+T;
end
bus(n,4)=P1;
for k=m+1:n %计算 PV 节点、平衡节点的无功
Q1=0;T=0;
for j=1:n
T=bus(k,2)*bus(j,2)*(real(Y(k,j))*sin(bus(k,3)-bus(j,3))-imag(Y(k,j))*cos(bus(k,3)-bus(j,3)));
Q1=Q1+T;
end
bus(k,5)=Q1;
end
bus(:,1)=f; %换回各节点、支路的初始编号
r=zeros(1,mb);
for t=1:n
for l=t+1:n
if bus(t,1)>bus(l,1)
r=bus(t,:);bus(t,:)=bus(l,:);bus(l,:)=r;
end
end
end
for i=1:nl
for j=1:2
for k=1:nb
if line(i,j)==nodenum(k,1)
line(i,j)=nodenum(k,2);
break
end
end
end
end
fclose(myf);
Pf=loss(bus,line); %计算支路潮流及损耗
%将节点导纳矩阵、节点潮流计算结果写入文件 output2
myf=fopen('output2.dat','w');
fprintf(myf, '--节点导纳矩阵--\n');
for k=1:n
for j=1:n
fprintf(myf,'%8.6f', real(Y(k,j)));
fprintf(myf, '+i*');
fprintf(myf,'%8.6f', imag(Y(k,j)));
fprintf(myf, ' ');
end
fprintf(myf, '\n');
end
fprintf(myf, '------------------牛顿-拉夫逊法潮流计算结果-------------\n');
fprintf(myf, '------------节点计算结果------------\n');
fprintf(myf, '-- 节点 节点电压 节点相角 注入有功功率(P) 注入无功功率(Q) 类型--\n');
for l=1:nb
for j=1:mb
if j==1|j==6
fprintf(myf, ' %8.1f ', bus(l,j));
elseif j==3
fprintf(myf, ' %8.6f ', bus(l,j)*180/pi);
else
fprintf(myf, ' %8.6f ', bus(l,j));
end
end
fprintf(myf, ' \n');
end
fprintf(myf, '--支路计算结果--\n');
fprintf(myf, '-- 节点(I) 节点(J) 线路功率 S(I,J) 线路功率 S(J,I) 线路损耗 dS(I,J)--\n');
for k=1:nl
for j=1:5
if j<=2
fprintf(myf,'%8.1f ', Pf(k,j));
fprintf(myf, ' ');
else
fprintf(myf,'%8.6f', real(Pf(k,j)));
fprintf(myf, '+i*');
fprintf(myf,'%8.6f', imag(Pf(k,j)));
fprintf(myf, ' ');
end
剩余19页未读,继续阅读
资源评论
oligaga
- 粉丝: 52
- 资源: 2万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功