clc
clear
%初始值设置
V = [0.9; 1.0; 1.0; 1.05; 1.05];
d = [0;0; 0; 0; -pi/2];
Ps=[1.6; 2.0; 3.7; 5];
Qs=[0.8; 1; 1.3];
Y=zeros(5);
%节点1
Y(1,1)=1/(0.1+j*0.35)+1/(0.04+j*0.25)+j*0.25;
Y(1,2)=1/(0.04+j*0.25);Y(2,1)=Y(1,2);
Y(1,3)=1/(0.1+j*0.35); Y(3,1)=Y(1,3);
Y(1,4:5)=0;Y(4:5,1)=0;
% 节点2
Y(2,2)=1/(0.04+j*0.25)+j*0.25+1/(0.08+j*0.30)+j*0.25+1/(j*0.015);
Y(2,3)=1/(0.08+j*0.30);Y(3,2)=Y(2,3);
Y(2,4)=1/(j*0.015);Y(4,2)=Y(2,4);
Y(2,5)=0;Y(5,2)=0;
% 节点3
Y(3,3)=1/(0.1+j*0.35)+1/(0.08+j*0.30)+j*0.25+1/(j*0.03);
Y(3,4)=0;Y(4,3)=0;
Y(3,5)=1/(j*0.03);Y(5,3)=Y(3,5);
% 节点4
Y(4,4)=1/(j*0.015);
Y(4,5)=0;Y(5,4)=0;
% 节点5
Y(5,5)=1/(j*0.03);
Y %导纳矩阵
y= abs(Y)
t = angle(Y)
iter=0; %计数器
pwracur = 0.5; % 潮流计算求解精度
DC = 10; %设置偏差系数
% while max(abs(DC)) > pwracur
% iter = iter +1 %迭代计数
% 求P1,P2,P3,P4,Q1,Q2,Q3初值
P=[V(1)^2*y(1,1)*cos(t(1,1))+V(1)*V(2)*y(1,2)*cos(t(1,2)-d(1)+d(2))+ ...
V(1)*V(3)*y(1,3)*cos(t(1,3)-d(1)+d(3));
V(2)*V(1)*y(2,1)*cos(t(2,1)-d(2)+d(1))+V(2)^2*y(2,2)*cos(t(2,2))+ ...
V(2)*V(3)*y(2,3)*cos(t(2,3)-d(2)+d(3))+V(2)*V(4)*y(2,4)*cos(t(2,4)-d(2)+d(4));
V(3)*V(1)*y(3,1)*cos(t(3,1)-d(3)+d(1))+V(3)*V(2)*y(3,2)*cos(t(3,2)-d(3)+d(2))+...
V(3)^2*y(3,3)*cos(t(3,3))+V(3)*V(5)*y(3,5)*cos(t(3,5)-d(3)+d(5));
V(4)*V(2)*y(4,2)*cos(t(4,2)-d(4)+d(2))+V(4)^2*y(4,4)*cos(t(4,4))]
Q=[-V(1)*V(2)*y(1,2)*sin(t(1,2)-d(1)+d(2))-V(1)^2*y(1,1)*sin(t(1,1))- ...
V(1)*V(3)*y(1,3)*sin(t(1,3)-d(1)+d(3));
-V(2)*V(1)*y(2,1)*sin(t(2,1)-d(2)+d(1))-V(2)^2*y(2,2)*sin(t(2,2))- ...
V(2)*V(3)*y(2,3)*sin(t(2,3)-d(2)+d(3))-V(2)*V(4)*y(2,4)*sin(t(2,4)-d(2)+d(4));
-V(3)*V(1)*y(3,1)*sin(t(3,1)-d(3)+d(1))-V(3)*V(2)*y(3,2)*sin(t(3,2)-d(3)+d(2))-...
V(3)^2*y(3,3)*sin(t(3,3))-V(3)*V(5)*y(3,5)*sin(t(3,5)-d(3)+d(5))]
% 计算雅可比矩阵
J(1,1)=V(2)*y(1,2)*cos(t(1,2)-d(1)+d(2))+2*V(1)*y(1,1)*cos(t(1,1))+...
V(3)*y(1,3)*cos(t(1,3)-d(1)+d(3));
J(1,2)=V(1)*y(1,2)*cos(t(1,2)-d(1)+d(2));
J(1,3)=V(1)*y(1,3)*cos(t(1,3)-d(1)+d(3));
J(1,4)=V(1)*V(2)*y(1,2)*sin(t(1,2)-d(1)+d(2))+V(1)*V(3)*y(1,3)*sin(t(1,3)-d(1)+d(3));
J(1,5)=-V(1)*V(2)*y(1,2)*sin(t(1,2)-d(1)+d(2));
J(1,6)=-V(1)*V(3)*y(1,3)*sin(t(1,3)-d(1)+d(3));
J(1,7)=0;
J(2,1)=V(2)*y(2,1)*cos(t(2,1)-d(2)+d(1));
J(2,2)=V(1)*y(2,1)*cos(t(2,1)-d(2)+d(1))+2*V(2)*y(2,2)*cos(t(2,2))+...
V(3)*y(2,3)*cos(t(2,3)-d(2)+d(3))+V(4)*y(2,4)*cos(t(2,4)-d(2)+d(4));
J(2,3)=V(2)*y(2,3)*cos(t(2,3)-d(2)+d(3));
J(2,4)=-V(2)*V(1)*y(2,1)*sin(t(2,1)-d(2)+d(1));
J(2,5)=V(2)*V(1)*y(2,1)*sin(t(2,1)-d(2)+d(1))+V(2)*V(3)*y(2,3)*sin(t(2,3)-d(2)+d(3))+...
V(2)*V(4)*y(2,4)*sin(t(2,4)-d(2)+d(4));
J(2,6)=-V(2)*V(3)*y(2,3)*sin(t(2,3)-d(2)+d(3));
J(2,7)=-V(2)*V(4)*y(2,4)*sin(t(2,4)-d(2)+d(4));
J(3,1)=V(3)*y(3,1)*cos(t(3,1)-d(3)+d(1));
J(3,2)=V(3)*y(3,2)*cos(t(3,2)-d(3)+d(2));
J(3,3)=V(1)*y(3,1)*cos(t(3,1)-d(3)+d(1))+V(2)*y(3,2)*cos(t(3,2)-d(3)+d(2))+...
2*V(3)*y(3,3)*cos(t(3,3))+V(5)*y(3,5)*cos(t(3,5)-d(3)+d(5));
J(3,4)=-V(3)*V(1)*y(3,1)*sin(t(3,1)-d(3)+d(1));
J(3,5)=-V(3)*V(3)*y(3,2)*sin(t(3,2)-d(3)+d(2));
J(3,6)=V(3)*V(1)*y(3,1)*sin(t(3,1)-d(3)+d(1))+V(3)*V(2)*y(3,2)*sin(t(3,2)-d(3)+d(2))+...
V(3)*V(5)*y(3,5)*sin(t(3,5)-d(3)+d(5));
J(3,7)=0;
J(4,1)=0;
J(4,2)=V(4)*y(4,2)*cos(t(4,2)-d(4)+d(2));
J(4,3)=0;
J(4,4)=0;
J(4,5)=-V(4)*V(2)*y(4,2)*sin(t(4,2)-d(4)+d(2));
J(4,6)=0;
J(4,7)=V(4)*V(2)*y(4,2)*sin(t(4,2)-d(4)+d(2));
J(5,1)=-V(2)*V(1)*y(2,1)*sin(t(2,1)-d(2)+d(1))-2*V(1)*y(1,1)*sin(t(1,1))-...
V(3)*y(1,3)*sin(t(1,3)-d(1)+d(3));
J(5,2)=-V(1)*y(1,2)*sin(t(1,2)-d(1)+d(2));
J(5,3)=-V(1)*y(1,3)*sin(t(1,3)-d(1)+d(3));
J(5,4)=V(1)*V(2)*y(1,2)*cos(t(1,2)-d(1)+d(2))+V(1)*V(3)*y(1,3)*cos(t(1,3)-d(1)+d(3));
J(5,5)=-V(1)*V(2)*y(1,2)*cos(t(1,2)-d(1)+d(2));
J(5,6)=-V(1)*V(3)*y(1,3)*cos(t(1,3)-d(1)+d(3));
J(5,7)=0;
J(6,1)=-V(2)*y(2,1)*sin(t(2,1)-d(2)+d(1));
J(6,2)=-V(1)*y(2,1)*sin(t(2,1)-d(2)+d(1))-2*V(2)*y(2,2)*sin(t(2,2))-...
-V(2)*y(2,3)*sin(t(2,3)-d(2)+d(3))-V(2)*V(4)*y(2,4)*sin(t(2,4)-d(2)+d(4));
J(6,3)=-V(2)*V(3)*y(2,3)*sin(t(2,3)-d(2)+d(3));
J(6,4)=-V(2)*V(1)*y(2,1)*cos(t(2,1)-d(2)+d(1));
J(6,5)=V(2)*V(1)*y(2,1)*cos(t(2,1)-d(2)+d(1))+V(2)*V(3)*y(2,3)*cos(t(2,3)-d(2)+d(3))+...
V(2)*V(4)*y(2,4)*cos(t(2,4)-d(2)+d(4));
J(6,6)=-V(2)*V(3)*y(2,3)*cos(t(2,3)-d(2)+d(3));
J(6,7)=-V(2)*V(4)*y(2,4)*cos(t(2,4)-d(2)+d(4));
J(7,1)=-V(3)*y(3,1)*sin(t(3,1)-d(3)+d(1));
J(7,2)=-V(3)*y(3,2)*sin(t(3,2)-d(3)+d(2));
J(7,3)=-V(1)*y(3,1)*sin(t(3,1)-d(3)+d(1))-V(2)*y(3,2)*sin(t(3,2)-d(3)+d(2))-...
2*V(3)*y(3,3)*sin(t(3,3))-V(5)*y(3,5)*sin(t(3,5)-d(3)+d(5));
J(7,4)=-V(3)*V(1)*y(3,1)*cos(t(3,1)-d(3)+d(1));
J(7,5)=-V(3)*V(2)*y(3,2)*cos(t(3,2)-d(3)+d(2));
J(7,6)=V(3)*V(1)*y(3,1)*cos(t(3,1)-d(3)+d(1))+V(3)*V(2)*y(3,2)*cos(t(3,2)-d(3)+d(2))+...
V(3)*V(5)*y(3,5)*cos(t(3,5)-d(3)+d(5));
J(7,7)=0;
% 利用已知量的误差求解未知量误差
DP = Ps - P;
DQ = Qs - Q;
DC = [DP; DQ] %偏差系数
J
DX = J\DC; %求解所求误差
% 修正误差
V(1)=V(1)+DX(1);
V(2)=V(2)+DX(3);
V(3)=V(3)+DX(3);
d(1)=d(1) +DX(4);
d(2)=d(2) +DX(5);
d(3)=d(3) +DX(6);
d(4)=d(4) +DX(7);
% end
% 得出结果
P5=V(5)*V(3)*y(5,3)*cos(t(5,3)-d(5)+d(3))+V(5)^2*y(5,5)*cos(t(5,5))
Q4=-V(4)*V(2)*y(4,2)*sin(t(4,2)-d(4)+d(2))-V(4)^2*y(4,4)*sin(t(4,4))
Q5=-V(5)*V(3)*y(5,3)*sin(t(5,3)-d(5)+d(3))-V(5)^2*y(5,5)*sin(t(5,5))
V
d
潮流计算matlab代码.zip_3节点_3节点潮流_congresszzb_潮流_潮流计算
版权申诉
5星 · 超过95%的资源 157 浏览量
2022-07-15
20:43:55
上传
评论
收藏 3KB ZIP 举报
我虽横行却不霸道
- 粉丝: 75
- 资源: 1万+
最新资源
- 中医临床诊疗术语-证候同义词典
- NovAtelConvert-Setup
- MySql主备数据库配置
- BGP路由基本配置(可以互通)
- 《软件方法2024版》公开内容202405更新-epub版
- 适用于tensorflow-2.11.0 CUDA版本11.2的cuDNN8.1版本
- 5Y study学习平台2016计算机基础-综合测试(8)_哔哩哔哩_bilibili_2580252704.mp4
- (大赛作品)STM32实现的F072RB NUCLEO智能家居控制.zip
- STM32实现的数字示波器源码+数字信号处理教程、配套实例.zip
- 【cookie续续】【cookie续续】【cookie续续】
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈