实验结果 文档
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
n=input('请输入节点数:n=');
nl=input('请输入支路数:nl=');
%isb=input('请输入平衡母线节点号:isb=');
%pr=input('请输入误差精度:pr=');
B1=input('请输入由各支路参数形成的矩阵:B1=');
%B2=input('请输入各节点参数形成的矩阵:B2=');
X=input('请输入由节点号和接地支路参数形成的矩阵:X=');
Y=zeros(n); e=zeros(1,n);f=zeros(1,n);V=zeros(1,n); O=zeros(1,n);S1=zeros(nl);
% zeros(m,n)产生m×n的零矩阵,zeros(n)产生n×n的全0方阵。
for i=1:n %节点数
if X(i,2)~=0; %matlab关系符号:>(大于),>=(大于等于),<(小于),<=(小于等于),==(等于)~=(不等于)
p=X(i,1);
Y(p,p)=1./X(i,2); %接地支路, “./”点除代表矩阵对应元素相除
end
end
for i=1:nl %支路数
if B1(i,6)==0 %折算到哪一侧的标志,0非标准变比在q侧,1非标准变比在在p侧p=B1(i,1);q=B1(i,2); % B1(i,1), B1(i,2)为支路编号
else
p=B1(i,2);q=B1(i,1);
end
Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5));
%计算节点p和q间的互导纳,B1(i,3)为支路阻抗,B1(i,5)为变压器支路变比k, B1(i,5)=1时不是变压器支路
Y(q,p)=Y(p,q);
Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2;
%对角元, 节点q的自导纳, B1(i,4)为支路对地容抗
Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; %对角元,节点p的自导纳
end
%输出导纳矩阵
disp('导纳矩阵Y=');
disp(Y);
请输入节点数:n=5
请输入支路数:nl=5
请输入由各支路参数形成的矩阵:B1=[1 2 0.03i 0 1.05 0;2 3 0.08+0.3i 0.5i 1 0;2 4 0.1+0.35i 0 1 0;3 4 0.04+0.25i 0.5i 1 0;3 5 0.015i 0 1.05 1];
请输入由节点号和接地支路参数形成的矩阵:X=[1 0;2 0;3 0;4 0;5 0];
实验结果:
导纳矩阵Y=
Columns 1 through 4
0 -33.3333i 0 +31.7460i 0 0
0 +31.7460i 1.5846 -35.7379i -0.8299 + 3.1120i -0.7547 + 2.6415i
0 -0.8299 + 3.1120i 1.4539 -66.9808i -0.6240 + 3.9002i
0 -0.7547 + 2.6415i -0.6240 + 3.9002i 1.3787 - 6.2917i
0 0 0 +63.4921i 0 Column 5
0 +63.4921i
0 -66.6667i