极坐标法潮流计算教材Word格式文档下载.docx
- 文档编号:14883860
- 上传时间:2022-10-25
- 格式:DOCX
- 页数:12
- 大小:53.13KB
极坐标法潮流计算教材Word格式文档下载.docx
《极坐标法潮流计算教材Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《极坐标法潮流计算教材Word格式文档下载.docx(12页珍藏版)》请在冰豆网上搜索。
试求:
采用极坐标下的牛顿-拉夫逊计算图1网络的潮流分布。
三、任务和要求任务:
熟练掌握计算机语言,并采用计算机编程进行下列计算:
根据电力系统网络推导电力网络数学模型,写出节点导纳矩阵;
掌握潮流计算的数值求解方法(节点导纳矩阵,修正方程),开发系统潮流计算的计算程序。
要求:
1.手工计算,手写,采用A4纸,得出计算结果。
2.编写程序:
它包括程序源代码;
程序说明;
部分程序的流程图;
程序运行结果,电子版。
注:
1.此表由指导教师填写,经系、教研室审批,指导教师、学生签字后生效;
2.此表1式3份,学生、指导教师、教研室各1份
四、参考资料和现有基础条件(包括实验室、主要仪器设备等)
[1]何仰赞等.电力系统分析[M].武汉:
华中理工大学出版社,2002.3
[2]西安交通大学等.电力系统计算[M].北京:
水利电力出版社,1993.12
五、进度安排
2010年12月20日:
下达课程设计的计划书,任务书,设计题目及分组情况。
2010年12月21日-23日:
学生完成潮流计算的手工计算。
2010年12月24日:
讲述课程设计编程的思路、要求。
举例:
用MATLAB软件编写的部分程序。
2010年12月25日-30日:
学生编写程序。
2011年1月1日-3日:
上机调试程序,得出正确结果。
2011年1月4日-5日:
整理课程设计报告。
2011年1月6日:
学生答辩。
六、教研室审批意见
教研室主任(签字):
年月日
七|、主管教学主任意见
主管主任(签字):
八、备注
学生(签字):
指导教师(签字):
设计主题
题目一:
在下图所示的简单电力系统中,系统中节点1、2为PQ节点,节点3为PV节点,节点4为平衡节点,已给定s10.4j0.3,s20.3j0.2,P30.4,V31.02,V41.05,40,网络各元件参数的标幺值如表1所示,给定电压的初始值如表2所示,收敛系数0.00001。
采用极坐标下的牛顿-拉夫逊计算图示网络的潮流计算。
表1网络各元件参数的标幺值
支路
电阻
电抗
1输电线路yc
2c
变压器变比k
1—2
0.02
0.06
0.01
—
1—3
0.03
2—3
0.07
2—4
0.0
0.05
0.9625
3—4
表2各节点电压(初值)标幺值参数
节点i
1
2
3
4
Ui(0)ei(0)jfi(0)
1.00+j0.0
1.0+j0.0
1.05+j0.0
3潮流计算流程图
本次课程设计采用极坐标下的牛顿-拉夫逊计算网络的潮流计算。
其牛顿-拉夫逊潮流计算程序框图如下所示。
图3.1极坐标下的牛顿-拉夫逊潮流计算程序框图
4手工计算插入手写的潮流计算过程
5MATLAB程序设计
5.1程序%电力系统极坐标下的牛顿-拉夫逊法潮流计算disp('
电力系统极坐标下的牛顿-拉夫逊法潮流计算:
'
);
clear
n=input('
请输入结点数:
n='
n1=input('
请输入PV结点数:
n1='
n2=input('
请输入PQ结点数:
n2='
isb=input('
请输入平衡结点:
isb='
pr=input('
请输入精确度:
pr='
K=input('
请输入变比矩阵:
K='
C=input('
请输入支路阻抗矩阵:
C='
y=input('
请输入支路导纳矩阵:
y='
U=input('
请输入结点电压矩阵:
U='
S=input('
请输入各结点的功率:
S='
Z=zeros(1,n);
N=zeros(n1+n2,n2);
L=zeros(n2,n2);
QT1=zeros(1,n1+n2);
form=1:
n
forR=1:
C(m,m)=C(m,m)+y(m,R);
ifK(m,R)~=0
C(m,m)=C(m,m)+1/(C(m,R)/(K(m,R)*(K(m,R)-1)));
C(R,R)=C(R,R)+1/(C(m,R)/(1-K(m,R)));
C(m,R)=C(m,R)/K(m,R);
C(R,m)=C(m,R);
end
form=1:
ifm~=R
Z(m)=Z(m)+1/C(m,R);
ifm==R
Y(m,m)=C(m,m)+Z(m);
else
Y(m,R)=-1/C(m,R);
enddisp('
结点导纳矩阵:
disp(Y);
disp('
迭代中的雅克比矩阵:
G=real(Y);
B=imag(Y);
O=angle(U);
U1=abs(U);
k=0;
PR=1;
P=real(S);
Q=imag(S);
whilePR>
pr
n2
UD(m)=U1(m);
n1+n2
PT(R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)-O(R))+B(m,R)*sin(O(m)-O(R)));
PT1(m)=sum(PT);
PP(m)=P(m)-PT1(m);
PP1(k+1,m)=PP(m);
QT(R)=U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R))-B(m,R)*cos(O(m)-O(R)));
QT1(m)=sum(QT);
QQ(m)=Q(m)-QT1(m);
QQ1(k+1,m)=QQ(m);
PR1=max(abs(PP));
PR2=max(abs(QQ));
PR=max(PR1,PR2);
H(m,m)=U1(m)^2*B(m,m)+QT1(m);
H(m,R)=-U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R))-B(m,R)*cos(O(m)-O(R)));
N(m,m)=-U1(m)^2*G(m,m)-PT1(m);
N(m,R)=-U1(m)*U1(R)*(G(m,R)*cos(O(m)-O(R))+B(m,R)*sin(O(m)-O(R)));
J(m,m)=U1(m)^2*G(m,m)-PT1(m);
J(m,R)=U1(m)*U1(R)*(G(m,R)*cos(O(m)-O(R))+B(m,R)*sin(O(m)-O(R)));
L(m,m)=U1(m)^2*B(m,m)-QT1(m);
L(m,R)=-U1(m)*U1(R)*(G(m,R)*sin(O(m)-O(R))-B(m,R)*cos(O(m)-O(
R)));
JJ=[HN;
JL];
disp(JJ);
PQ=[PP'
;
QQ'
];
DA=-inv(JJ)*PQ;
DA1=DA'
OO(m)=DA1(m);
form=n:
n1+n2+n2
UU1(m-n1-n2)=DA1(m);
UD2=diag(UD);
UU=UU1*UD2;
O(m)=O(m)+OO(m);
endform=1:
U1(m)=U1(m)+UU(m);
o(k+1,m)=180/pi*O(m);
u(k+1,m)=U1(m);
k=k+1;
b(m)=U1(m)*cos(O(m));
c(m)=U1(m)*sin(O(m));
U=b+i*c;
PH1(R)=U(isb)*conj(Y(isb,R))*conj(U(R));
PH=sum(PH1);
C1(m,R)=1/C(m,R);
C1(m,m)=C(m,m);
if(C(m,R)~=inf)&
(m~=R)
SS(m,R)=U1(m)^2*conj(C1(m,m))+U(m)*(conj(U(m))-conj(U(R)))*conj(C1(m,R));
迭代中的△P:
'
disp(PP1);
迭代中的△Q:
disp(QQ1);
disp('
迭代中相角:
disp(o);
迭代中电压的模:
disp(u);
平衡结点的功率:
disp(PH);
全部线路功率分布:
disp(SS);
matlab默认输出结果保留4位小数,可在显示屏上输入>
>
fomatlong此时小数点后面保留14位小数
5.2程序结果请输入结点数:
n=4请输入PV结点数:
n1=1请输入PQ结点数:
n2=2请输入平衡结点:
isb=4请输入精确度:
pr=0.00001请输入变比矩阵:
K=[0000;
0000
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 坐标 潮流 计算 教材