版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领
文档简介
end
实验一
clear
%n=input('请输入节点数:n=');
%nl=input('请输入支路数:nl=');
%B1=input('请输入参数矩阵:B1=');
%B2=input('请输入变压器参数矩阵:B2=');
%B3=input('请输入节点参数矩阵:B3=');n=5;
n1=5;
B1=[00.105i000;0.105i00.024+0.065i0.03+0.08i0;00.024+0.065i00.018+0.05i0;0
0.03+0.08i0.018+0.05i00.184i;0000.184i0];
B2=[00000;1.050000;00000;00000.96;00000];
B3=[00000;000.016i0.02i0;00.016i00.013i0;00.02i0.013i00;00000];
Y=zeros(n);fori=1:n
forj=1:n
ifi==j
fork=1:n
Y(i,j)=Y(i,j)+B3(i,k);
ifB2(i,k)~=0
Y(i,j)=Y(i,j)+1./(B2(i,k)^2*B1(k,i));
else
ifB2(k,i)~=0
Y(i,j)=Y(i,j)+1./B1(i,k);
else
ifB1(i,k)~=0
Y(i,j)=Y(i,j)+1./B1(i,k);
else
end
end
end
end
ifB2(i,j)~=0
Y(i,j)=-1./(B1(i,j)*B2(i,j));
else
ifB2(j,i)~=0
Y(i,j)=-1./(B1(i,j)*B2(j,i));
else
ifB1(i,j)~=0
Y(i,j)=-1./B1(i,j);
end
end
end
end
end
disp(Y);0
+5.6612i
0-9.5238i 0+9.0703i 0 0
0+9.0703i 9.1085-33.1001i -4.9990+13.5388i -4.1096+10.9589i 0
0 -4.9990+13.5388i 11.3729-31.2152i -6.3739+17.7054i 0
0 -4.1096+10.9589i -6.3739+17.7054i 10.4835-34.5284i 0
0 0 0 0+5.6612i
0-5.4348i
实验二-1clear
%n=input('请输入节点数:n=');
%nl=input('请输入支路数:nl=');
%B1=input('请输入参数矩阵:B1=');
%B2=input('请输入变压器参数矩阵:B2=');
%B3=input('请输入节点参数矩阵:B3=');n=4;
n1=3;
B1=[0.2i00.51i0;
04i0.59i0;
0.51i0.59i01.43i;
001.43i0];
B2=[0000;0000;0000;0000];
B3=[0000;0000;0000;0000];f=4;
Y=zeros(n);fori=1:n
forj=1:n
ifi==j
fork=1:n
Y(i,j)=Y(i,j)+B3(i,k);
ifB2(i,k)~=0
Y(i,j)=Y(i,j)+1./(B2(i,k)^2*B1(k,i));
else
ifB2(k,i)~=0
Y(i,j)=Y(i,j)+1./B1(i,k);
else
ifB1(i,k)~=0
Y(i,j)=Y(i,j)+1./B1(i,k);
end
else
end
end
ifB2(i,j)~=0
Y(i,j)=-1./(B1(i,j)*B2(i,j));
else
ifB2(j,i)~=0
Y(i,j)=-1./(B1(i,j)*B2(j,i));
else
ifB1(i,j)~=0
Y(i,j)=-1./B1(i,j);
end
end
end
end
end
end
disp(Y);
Z=inv(Y);
If=1/Z(f,f);If
0-6.9608i
0
0+1.9608i
0
0
0-1.9449i
0+1.6949i
0
0+1.9608i
0+1.6949i
0-4.3550i
0+0.6993i
0
0
0+0.6993i
0-0.6993i
If=
0-0.4890i
>>
实验二-2clear
z=[0.2iinf0.51iinf;inf4i0.59iinf;0.51i0.59iinf1.43i;infinf1.43iinf];
y=[0000;
0000;
0000;
0000];f=4;
Y=zeros(4,4);
fori=1:4
forj=1:4ifi==j
Y(i,j)=Y(i,j)else
Y(i,j)=-1.0/z(i,j)end
end
end
fori=1:4
forj=1:4
Y(i,i)=Y(i,i)+y(i,j)+1.0/z(i,j)
end
end
Z=inv(Y);
If=1/Z(f,f);If
Y=
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
Y=
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
Y=
0
0
0+1.9608i
0
0
0
0
0
0
0
0
0
0
0
0
0
Y=
0
0
0+1.9608i
0
0
0
0
0
0
0
0
0
0
0
0
0
Y=
0
0
0+1.9608i
0
0
0
0
0
0
0
0
0
0
0
0
0
Y=
0
0
0+1.9608i
0
0
0
0
0
0
0
0
0
0
0
0
0
Y=
0
0
0+1.9608i
0
0
0
0+1.6949i
0
0
0
0
0
0
0
0
0
Y=
0
0
0+1.9608i
0
0
0
0+1.6949i
0
0
0
0
0
0
0
0
0
Y=
0
0
0+1.9608i
0
0
0
0+1.6949i
0
0+1.9608i
0
0
0
0
0
0
0
Y=
0
0
0+1.9608i
0
0
0
0+1.6949i
0
0+1.9608i
0
0+1.6949i
0
0
0
0
0
Y=
0
0
0+1.9608i
0
0
0
0+1.6949i
0
0+1.9608i
0
0+1.6949i
0
0
0
0
0
Y=
0
0
0+1.9608i
0
0
0
0+1.6949i
0
0+1.9608i
0
0+1.6949i
0
0
0
0+0.6993i
0
Y=
0
0
0+1.9608i
0
0
0
0+1.6949i
0
0+1.9608i
0
0+1.6949i
0
0
0
0+0.6993i
0
Y=
0
0
0+1.9608i
0
0
0
0+1.6949i
0
0+1.9608i
0
0+1.6949i
0
0
0
0+0.6993i
0
Y=
0
0
0+1.9608i
0
0
0
0+1.6949i
0
0+1.9608i
0
0+1.6949i
0
0
0+0.6993i
0+0.6993i
0
Y=
0
0
0+1.9608i
0
0
0
0+1.6949i
0
0+1.9608i
0
0+1.6949i
0
0
0+0.6993i
0+0.6993i
0
Y=
0-5.0000i
0
0+1.9608i
0
0
0
0+1.6949i
0
0+1.9608i
0
0+1.6949i
0
0
0+0.6993i
0+0.6993i
0
Y=
0-5.0000i
0
0+1.9608i
0
0
0
0+1.6949i
0
0+1.9608i
0
0+1.6949i
0
0
0+0.6993i
0+0.6993i
0
Y=
0-6.9608i
0
0+1.9608i
0
0
0
0+1.6949i
0
0+1.9608i
0
0+1.6949i
0
0
0+0.6993i
0+0.6993i
0
Y=
0-6.9608i
0
0+1.9608i
0
0
0
0+1.6949i
0
0+1.9608i
0
0+1.6949i
0
0
0+0.6993i
0+0.6993i
0
Y=
0-6.9608i
0
0
0
0+1.9608i
0+1.6949i
0
0
0+1.9608i
0
0
+1.6949i0
0
0+0.6993i
0
+0.6993i0
Y=
0-6.9608i
0
0+1.9608i
0
0
0-0.2500i
0+1.6949i
0
0+1.9608i
0
0
+1.6949i0
0
0+0.6993i
0
+0.6993i0
Y=
0-6.9608i
0
0+1.9608i
0
0
0-1.9449i
0+1.6949i
0
0+1.9608i
0
0
+1.6949i0
0
0+0.6993i
0
+0.6993i0
Y=
0-6.9608i
0
0+1.9608i
0
0
0-1.9449i
0+1.6949i
0
0+1.9608i
0
0
+1.6949i0
0
0+0.6993i
0
+0.6993i0
Y=
0-6.9608i
0
0+1.9608i
0
0
0-1.9449i
0+1.6949i
0
0+1.9608i
0+1.6949i
0-1.9608i
0+0.6993i
0
0
0+0.6993i
0
Y=
0-6.9608i
0
0+1.9608i
0
0
0+1.9608i
0
0-1.9449i
+1.6949i
0+1.6949i
0-3.6557i
0
0+0.6993i
0
0
0+0.6993i
0
Y=
0-6.9608i
0
0+1.9608i
0
0
0+1.9608i
0
0-1.9449i
+1.6949i
0+1.6949i
0-3.6557i
0
0+0.6993i
0
0
0+0.6993i
0
Y=
0-6.9608i
0
0+1.9608i
0
0
0+1.9608i
0
0-1.9449i
+1.6949i
0+1.6949i
0-4.3550i
0
0+0.6993i
0
0
0+0.6993i
0
Y=
0-6.9608i
0
0+1.9608i
0
0
0+1.9608i
0
0-1.9449i
+1.6949i
0+1.6949i
0-4.3550i
0
0+0.6993i
0
0
0+0.6993i
0
Y=
0-6.9608i
0
0+1.9608i
0
0
0+1.9608i
0
0-1.9449i
+1.6949i
0+1.6949i
0-4.3550i
0
0+0.6993i
0
0
0+0.6993i
0
Y=
温馨提示
- 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 肇庆固化地面施工方案
- 融雪剂施工方案
- 江夏阳光房制作施工方案
- 2024-2030年中国制造服务业市场现状调查及发展前景研判报告
- 中国贮氢材料行业市场深度分析及投资战略研究报告
- 2024-2025年中国SAAS行业市场调研分析及投资战略咨询报告
- 量角尺行业深度研究报告
- 2025年运动器材项目可行性研究报告-20250102-070309
- 2025年中国电动物流车行业发展监测及发展趋势预测报告
- 航次运输贵金属合同样本
- 2024年四川省德阳市中考道德与法治试卷（含答案逐题解析）
- 施工现场水电费协议
- SH/T 3046-2024 石油化工立式圆筒形钢制焊接储罐设计规范（正式版）
- 六年级数学质量分析及改进措施
- 一年级下册数学口算题卡打印
- 真人cs基于信号发射的激光武器设计
- 【阅读提升】部编版语文五年级下册第三单元阅读要素解析 类文阅读课外阅读过关（含答案）
- 四年级上册递等式计算练习200题及答案
- 法院后勤部门述职报告
- 2024年国信证券招聘笔试参考题库附带答案详解
- 道医馆可行性报告
评论
0/150
提交评论