%节点导纳矩阵的求解
%B1为5列,依次为i,j,r,x,k(节点号1,节点号2,电阻,电抗,变压器变比>1)
clear;
a=input('请输入节点数:'); %a是网络节点数
b=input('请输入支路数:'); %b是网络的支路数
fid=fopen('B1.txt','r');
A=fscanf(fid,'%f %f',[5 inf]); %支路数矩阵,有5列
B=A'; %支路数据矩阵
dn=zeros(a,a);
for i=1:b
if(B(i,5)==1) %线路支路
for n=1:2
for k=1:a
if(B(i,n)==k)
dn(k,k)=dn(k,k)+B(i,3)+j*B(i,4);
end
end
end
for l=1:b
for m=1:5
dn(B(l,1),B(l,2))=B(l,3)+j*B(l,4);
dn(B(l,2),B(l,1))=B(l,3)+j*B(l,4);
end
end
else %变压器支路
dn(B(i,2),B(i,2))=dn(B(i,2),B(i,2))+B(i,3)+j*B(i,4);
dn(B(i,1),B(i,1))=dn(B(i,1),B(i,1))+(1/(B(i,5)^2))*(B(i,3)+j*B(i,4));
dn(B(i,1),B(i,2))=-(1/B(i,5))*(B(i,3)+j*B(i,4));
dn(B(i,2),B(i,1))=-(1/B(i,5))*(B(i,3)+j*B(i,4));
end
end
%disp('节点数:');
%a
%disp('支路数:');
%b
%disp('导纳矩阵为:');
%dn
jddn.rar_导纳_节点导纳
版权申诉
192 浏览量
2022-09-24
02:12:00
上传
评论
收藏 573B RAR 举报
刘良运
- 粉丝: 68
- 资源: 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续续】
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈