IEEE14节点求解系统潮流matlab仿真( PQ分解法)
目录
程序运行结果与IEEE14节点标准测试系统数据比较
设置变量和相关参数
线路和变压器参数矩阵
节点参数矩阵
计算节点导纳矩阵
调整节点编号
计算节点导纳矩阵
PQ解偶迭代
发电机和负荷的注入功率
雅可比矩阵
迭代主程序
潮流计算
还原节点顺序
输出结果
程序运行结果与IEEE14节点标准测试系统数据比较
设置变量和相关参数
%节点并联对地之路参数
GND=[9,0,0.19;];
%对地支路数
ngnd=1;
%节点个数
nbb=14;
%线路支路数
ntr=20;
%最大迭代数
nitmax=1000;
%容许误差
error=1e-3;
%迭代次数初始化
nit=0;
%PQ、PV节点数
PQNodes=0;
PVNodes=0;
%节点注入功率
P0=zeros(1,nbb-1);
Q0=[];
%节点导纳矩阵
Y=zeros(nbb,nbb);
线路和变压器参数矩阵
% 入 出 电阻 电抗 电导 电纳 有无 变比
% 调压
% 变压器
%
tr=[1, 2, .01938, .05917, 0, .0264*2, 0, 1.000;2, 3, .04699, .19797, 0, .0219*2, 0, 1.000;2, 4, .05811, .17632, 0, .0187*2, 0, 1.000;1, 5, .05403, .22304, 0, .0246*2, 0, 1.000;2, 5, .05695, .17388, 0, .0170*2, 0, 1.000;3, 4, .06701, .17103, 0, .0173*2, 0, 1.000;4, 5, .01335, .04211, 0, .0064*2, 0, 1.000;7, 8, .00000, .17615, 0, .0000, 0, 1.000;7, 9, .00000, .11001, 0, .0000, 0, 1.000;9, 10, .03181, .08450, 0, .0000, 0, 1.000;6, 11, .09498, .19890, 0, .0000, 0, 1.000;6, 12, .12291, .15581, 0, .0000, 0, 1.000;6, 13, .06615, .13027, 0, .0000, 0, 1.000;9, 14, .12711, .27038, 0, .0000, 0, 1.000;10, 11, .08205, .19207, 0, .0000, 0, 1.000;12, 13, .22092, .19988, 0, .0000, 0, 1.000;13, 14, .17093, .34802, 0, .0000, 0, 1.000;4, 7, .00000, .20912, 0, .0000, 1, 0.978; 4, 9, .00000, .55618, 0, .0000, 1, 0.969; 5, 6, .00000, .25202, 0, .0000, 1, 0.932; ];
节点参数矩阵
%节点参数矩阵
%1是平衡节点
%2是pv节点
%3是pq节点node=[%编号 类型 发动机 负荷 电压幅值 电压相位 1, 1, 0, 0, 1.0600, 0; 2, 2, 0.40, 0.2170+0.1270i, 1.0450, 0; 3, 2, 0, 0.9420+0.1900i, 1.0100, 0; 4, 3, 0, 0.4780-0.0390i, 1.0000, 0; 5, 3, 0, 0.0760+0.0160i, 1.0000, 0; 6, 2, 0, 0.1120+0.0750i, 1.0700, 0; 7, 3, 0, 0.0000+0.0000i, 1.0000, 0; 8, 2, 0, 0.0000+0.0000i, 1.0900, 0; 9, 3, 0, 0.2950+0.1660i, 1.0000, 0; 10, 3, 0, 0.0900+0.0580i, 1.0000, 0;11, 3, 0, 0.0350+0.0180i, 1.0000, 0; 12, 3, 0, 0.0610+0.0160i, 1.0000, 0; 13, 3, 0, 0.1350+0.0580i, 1.0000, 0; 14, 3, 0, 0.1490+0.0500i, 1.0000, 0; ];
计算节点导纳矩阵
调整节点编号
%PQ
%PV
%平衡节点
node2=[];
node3=[];
for i=1:nbbif(node(i,2)==1)%平衡节点放最后node1=node(i,:);elseif(node(i,2)==2)%pv节点放中间node2=[node2;node(i,:)];PVNodes=PVNodes+1;elseif(node(i,2)==3)%PQ节点放开头node3=[node3;node(i,:)];PQNodes=PQNodes+1;end
end
node = [node3;node2;node1];%新节点
clearvars node1 node2 node3;
tr1=tr;
for k=1:nbb for i=1:ntrfor j=1:2if(tr1(i,j)==node(k,1))tr(i,j)=k;endendend
end
计算节点导纳矩阵
%2---------------------节点导纳矩阵:for i=1:ntrY(tr(i,1),tr(i,2))=-(1/tr(i,8))/(tr(i,3)+tr(i,4)*1i);Y(tr(i,2),tr(i,1))=-(1/tr(i,8))/(tr(i,3)+tr(i,4)*1i);Y(tr(i,2),tr(i,2))=Y(tr(i,2),tr(i,2))+1/(tr(i,3)+tr(i,4)*1i)+tr(i,6)/2*1i;Y(tr(i,1),tr(i,1))=Y(tr(i,1),tr(i,1))+(1/tr(i,8)^2)/(tr(i,3)+tr(i,4)*1i)+tr(i,6)/2*1i;
endfor cc=1:ngnd %对地并联支路的修正i=4;%GND(cc,1);Y(i,i)=Y(i,i)+GND(cc,2)-GND(cc,3)*1i;end
YR=real(Y);
YI=imag(Y);
PQ解偶迭代
PQ解偶后有:
发电机和负荷的注入功率
Q0=zeros(1,nbb-1-PVNodes);
for i=1:nbb-1P0(i)=node(i,3)-real(node(i,4));
end
for i=1:nbb-1-PVNodesQ0(i)=0-imag(node(i,4));
end
雅可比矩阵
B1=-YI;
B1(nbb,:)=[];
B1(:,nbb)=[];B2=-YI;
for i=0:nbb-1if(nbb-i<=nbb-1-PVNodes)break;endB2(nbb-i,:)=[];B2(:,nbb-i)=[];
end
迭代主程序
%迭代
while(nit<=nitmax)%DP=Delta P/U %DQ=Delta Q/UDP=zeros(1,nbb-1);DQ=zeros(1,nbb-1-PVNodes); DP1=zeros(1,nbb-1);DQ1=zeros(1,nbb-1-PVNodes);for i=1:nbb-1psum=0;for j=1:nbbpsum=psum+node(i,5)*node(j,5)*(YR(i,j)*cos(node(i,6)-node(j,6))+YI(i,j)*sin(node(i,6)-node(j,6))); %计算有功endDP1(i)=(P0(i)-psum);DP(i)=DP1(i)/node(i,5);endfor i=1:nbb-1-PVNodesqsum=0;for j=1:nbbqsum=qsum+node(i,5)*node(j,5)*YR(i,j)*sin(node(i,6)-node(j,6))-node(i,5)*node(j,5)*YI(i,j)*cos(node(i,6)-node(j,6)); %计算无功endDQ1(i)=(Q0(i)-qsum);DQ(i)=DQ1(i)/node(i,5);endMAX=max(max(abs(DP1)),max(abs(DP1)));if(MAX<error)break;endUDeltadelta=B1\DP';Deltadelta=UDeltadelta./node(1:nbb-1,5);DeltaU=B2\DQ';%更新deltafor i=1:nbb-1node(i,6)=node(i,6)+Deltadelta(i);end%更新Ufor i=1:nbb-1-PVNodesnode(i,5)=node(i,5)+DeltaU(i);endnit=nit+1;end
潮流计算
%计算发电机有功功率
psum=0;
for j=1:nbbpsum=psum+node(14,5)*node(j,5)*(YR(14,j)*cos(node(14,6)-node(j,6))+YI(14,j)*sin(node(14,6)-node(j,6))); %计算有功
end
node(14,3)=psum;
%计算发电机无功
for i=nbb-PVNodes:nbbqsum=0;for j=1:nbbqsum=qsum+node(i,5)*node(j,5)*YR(i,j)*sin(node(i,6)-node(j,6))-node(i,5)*node(j,5)*YI(i,j)*cos(node(i,6)-node(j,6)); %计算无功endnode(i,3)=node(i,3)+imag(node(i,4))*1i+qsum*1i;
end
还原节点顺序
node1=zeros(nbb,6);
for i=1:nbbfor j=1:nbbif(i==node(j,1))node1(i,:)=node(j,:);node1(i,6)=node1(i,6)*180/pi;endend
end
输出结果
fprintf("\n 节点号 节点类型 发动机功率 负载功率 电压幅值 相角\n")
format short
disp(node1);
程序下载地址:IEEE14节点PQ分解法求解系统潮流matlab仿真,程序结构简单,仿真结果几乎与原测试数据相同-Matlab文档类资源-CSDN下载
IEEE14节点求解系统潮流matlab仿真( PQ分解法)相关推荐
- matlab 电力系统分析 毕设,基于MATLAB的P-Q分解法电力系统潮流计算毕业设计
基于MATLAB的P-Q分解法电力系统潮流计算毕业设计 摘 要 电力系统潮流计算是研究电力系统稳定运行情况的一种重要的计算,在电力系统规划设计和现有电力系统运行方式的研究中,都需要利用它来定量地分析比 ...
- pq分解法matlab程序,基于MATLAB软件的PQ分解法潮流计算
基于MATLAB软件的PQ分解法潮流计算 基于MATLAB软件的P-Q分解法潮流计算 摘要 电力系统潮流计算是研究电力系统稳态运行情况的一种重要的分析计算,它根据给定的运行条件及系统接线情况确定整个电 ...
- matlab p q的确定,基于MATLAB软件的P-Q分解法潮流计算 (1)
基于MATLAB软件的P-Q分解法潮流计算 摘要 电力系统潮流计算是研究电力系统稳态运行情况的一种重要的分析计算,它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态:各母线的电压,各元 ...
- pq分解法潮流计算的matlab,第四节PQ分解法潮流计算
<第四节PQ分解法潮流计算>由会员分享,可在线阅读,更多相关<第四节PQ分解法潮流计算(17页珍藏版)>请在人人文库网上搜索. 1.第四节 PQ分解法潮流计算一 .PQ分解法的 ...
- 静电场的有限差分法与matlab 仿真课程设计,计算物理和MATLAB课程设计--自激振动系统的MATLAB仿真.doc...
东北石油大学课程设计任务书 课程 计算物理和MATLAB课程设计 题目 自激振动系统的MATLAB仿真 专业 姓名 学号 主要内容.基本要求.主要参考资料等 主要内容: 研究范?德?波耳(Van de ...
- matlab仿真生成信号程序,信号与系统的MATLAB仿真
用哪个MATLAB进行信号与系统仿真,有详细的实例程序及结果 1 信号与系统的MATLAB 仿真 一.信号生成与运算的实现 1.1 实现)3(sin )()(π±===t t t t S t f a ...
- dtmf信号系统的matlab仿真,dtmf信号系统的matlab仿真毕业设计
dtmf信号系统的matlab仿真毕业设计 DTMF 信号系统的 Matlab 仿真摘 要双音多频(Dual Tone Multi Frequency, DTMF)信号是音频电话中的拨号信号,由美国 ...
- matlab中电流调节器仿真模型,直流调速系统的matlab仿真(课程设计作业).docx
直流调速系统的matlab仿真(课程设计作业).docx 运动控制系统课程设计班级自动化0802姓名王有录学号0806050231摘要直流调速系统具有调速范围广.精度高.动态性能好和易于控制等优点,所 ...
- 对蓄电池matlab仿真,UPS系统的MATLAB仿真
科 学 论 坛 I■ China science and Technology Review UPS系统的 MATLAB仿真 刘春玉 (武汉 理工大学 4 300 70) [摘 要]UPS系统主要用 ...
最新文章
- ubuntu系统安装mysql(deb-bundle包)
- windows server 2012多用户用一个账号同时登陆
- Hibernate----Hibernate小配置
- Team Work(CF 932 E)[bzoj5093][Lydsy1711月赛]图的价值
- int size java_int size()
- 线性代数之——行列式及其性质
- 这个夏天,感动我的歌,感动我的你
- JavaScript(四)—— JavaScript 内置对象/JavaScript 简单数据类型与复杂类型
- Oracle10gR5GridControl InstallOnLinux
- php连接db2失败,PHP_浅谈DB2数据库故障处理及最佳实践, 当你在使用数据库时,可 - phpStudy...
- 牛客网-《剑指offer》-跳台阶
- JavaScript:三大家族
- java编程新手初学者入门学习,看什么书比较好
- 多个ai文件合并成pdf_AI打开多页PDF教程.pdf
- layui表单验证,表单提交的若干个方法
- python提取国家名称缩写_如何从python中的国家缩写中获取国家名称,并混合使用alpha_2和alpha_3 ccodes...
- 医院挂号系统代码_智慧医院中心是怎样做的?分诊叫号系统如何正确使用!
- 软件设计师——数据流图
- Android基于Facebook Rebound的动画效果框架Backboard demo (非常炫酷)
- rem适配布局 Less基础 插件 cssrem
热门文章
- 计算机导论知识点总结与试题
- 使用DCA1000和EVM板进行数据采集
- 论文翻译:Semantics-Guided Neural Networks for Efficient Skeleton-Based Human Action Recognition
- python合成图片_python图片合成的示例
- 高并发框架 Disruptor
- 专题05-python操作微信(wxpy)
- 超新星计算机网络技术就业方向,超新星网络影响力榜单:林彦俊第三,蔡徐坤跌至第二,第一力压登顶...
- 成兴光 | LED灯珠的封装形式
- 爱了爱了!java是不是比前端难
- IIS——实现局域网访问本地网站