没有合适的资源?快使用搜索试试~ 我知道了~
MATLAB潮流程序(IEEE14_直角坐标_牛拉法).doc
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
5星 · 超过95%的资源 2 下载量 185 浏览量
2022-07-06
01:23:01
上传
评论
收藏 41KB DOC 举报
温馨提示
试读
4页
MATLAB潮流程序(IEEE14_直角坐标_牛拉法)
资源推荐
资源详情
资源评论
MATLAB 潮流程序(IEEE14 直角坐标 牛拉法)
clear
baseMVA=100; %功率基值
%%读Data1中数据
load Data1.txt
Bus=Data1(:,1); %节点号
Vtype=Data1(:,5); %节点类型
Pload=Data1(:,8); %负载有功
Qload=Data1(:,9); %负载无功
Pgen=Data1(:,10); %发电机发出有功
Qgen=Data1(:,11); %发电机发出无功
Vset=Data1(:,13); %电压设定点
Qsh=Data1(:,17); %并联电容电纳标幺值
%%读Data2中数据
load Data2.txt
II=Data2(:,1);
JJ=Data2(:,2); %支路端点号
Ltype=Data2(:,5); %线路类型
R=Data2(:,6); %两点间电阻
X=Data2(:,7); %两点间电抗
B=Data2(:,8)/2; %线路对地电纳
K=Data2(:,14); %变压器非标准变压比
%%求导纳矩阵Y
y1=zeros(14);
y2=zeros(14);
y3=zeros(14);
lin=length(II); %支路数
for x=1:lin
switch Ltype(x)
case 1
y1(II(x),JJ(x))=1/(R(x)+i*X(x));
y1(JJ(x),II(x))=y1(II(x),JJ(x));
y3(II(x),JJ(x))=i*B(x);
y3(JJ(x),II(x))=i*B(x);
case 2
y1(II(x),JJ(x))=1/((R(x)+i*X(x))*K(x));
y1(JJ(x),II(x))=y1(II(x),JJ(x));
y2(II(x),JJ(x))=(1-K(x))/((R(x)+i*X(x))*K(x)^2);
y2(JJ(x),II(x))=(K(x)-1)/((R(x)+i*X(x))*K(x));
end
end
clear x
Y=zeros(14);
资源评论
- 阿J7172022-07-13感谢资源主分享的资源解决了我当下的问题,非常有用的资源。
- FugwAa2023-05-04感谢大佬,让我及时解决了当下的问题,解燃眉之急,必须支持!
老帽爬新坡
- 粉丝: 79
- 资源: 2万+
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功