哈工大电力系统分析短路及潮流计算实验上机程序汇总.docx

1、哈工大电力系统分析短路及潮流计算实验上机程序汇总上机实验实验一 节点导纳阵生成和短路电流计算 实验二 简单系统的牛顿法潮流计算姓名：班级： 141班学号：实验说明：本程序的电路结构来自翁增银、何仰赞主编的电力系统分析的例题实验一 节点导纳阵生成和短路电流计算一、实验目的根据所给的电力系统，编制短路电流计算程序，通过计算机进行调试，最后完成一个切实可行的电力系统计算应用程序。通过自己设计电力系统计算程序加深对电力系统分析的理解，同时加强计算机实际应用能力。二、实验内容1、编写数据输入、形成导纳阵程序2、电力系统短路计算实用公式的计算程序及编制和调试三、实验过程1、程序代码及说明%清屏clear%

2、读数据fid=fopen(node5.txt,r);A=fscanf(fid,%d,2);B=fscanf(fid,%f,5,A(2,1);fclose(fid);B=B;%求节点导纳矩阵Y=zeros(A(1,1); % 形成55的0阵%= %导纳阵元素计算for a=1:1:A(2,1)m=B(a,1);n=B(a,2);if B(a,5)0Y(m,m)=Y(m,m)+1/(B(a,3)+j*B(a,4);Y(n,n)=Y(n,n)+1/(B(a,3)+j*B(a,4)/(B(a,5)2);Y(m,n)=-1/(B(a,3)+j*B(a,4)/B(a,5);Y(n,m)=Y(m,n);el

3、seif B(a,5)0I(m,n)=(B(a,5)*V(1,m)-V(1,n)/(B(a,3)+j*B(a,4);elseif B(a,5)0I(m,n)=(V(1,m)-V(1,n)/(B(a,3)+j*B(a,4);endend%=fid=fopen(ans.txt,w);fprintf(fid,Y_matrixn);%= %输出导纳阵for i=1:1:5for j=1:1:5k=Y(i,j);re=real(k);fprintf(fid,%8.4f,re);im=imag(k);if im=0fprintf(fid,+%8.4fjt,im);endendfprintf(fid,n);

4、end%=%= %输出阻抗阵，导纳阵求逆fprintf(fid,Z_matrixn);for i=1:1:5for j=1:1:5k=Z(i,j);re=real(k);fprintf(fid,%8.4f,re);im=imag(k);if im=0fprintf(fid,+%8.4fjt,im);endendfprintf(fid,n);end%=%= %打印输出短路电流fprintf(fid,If=);re=real(If);fprintf(fid,%8.4f,re);im=imag(If);if im=0fprintf(fid,+%8.4fjt,im);endfprintf(fid,nV

5、n);%=%= %打印输出节点电压for i=1:1:5fprintf(fid,Note%d:V%d=,i,i);k=V(1,i);re=real(k);fprintf(fid,%8.4f,re);im=imag(k);if im=0fprintf(fid,+%8.4fjt,im);endfprintf(fid,n);end%=%= %输出打印支路电流fprintf(fid,Ibrn);for i=1:1:5for j=1:1:4k=I(i,j);re=real(k);im=imag(k);if(re=0|im=0)fprintf(fid,Branch%d-%d:I%d%d=,i,j,i,j)

6、;fprintf(fid,%8.4f,re);if im=0fprintf(fid,+%8.4fjn,im);endendendend%=fclose(fid); %关闭文件附：node5.txt2、程序输出结果节点导纳阵：节点阻抗阵：短路电流：If= 0.0001+ 0.1082j节点电压：各支路电流：四、实验总结这是我的第一次上机实验，感觉稍微有点难，主要还是在工具软件C语言或者MATLAB的运用上，但是我相信，以后学习中，我会努力掌握的，这是我把理论应用于实际中的必要桥梁！实验二 简单系统的牛顿法潮流计算一、实验目的根据所给的电力系统，编制牛顿法潮流计算程序，通过计算机进行调试，最后完成

7、一个切实可行的电力系统计算应用程序。通过自己设计电力系统计算程序加深对电力系统分析的理解，同时加强计算机实际应用能力。二、实验内容电力系统潮流计算的计算程序设计及编制和调试。三、程序框图四、实验过程1、实验程序及说明clearfid=fopen(node4.txt,r); %打开输入数据A=fscanf(fid,%f,8); %读8个数B=fscanf(fid,%f,5,A(2,1);C=fscanf(fid,%f,3,(A(1,1)-1);fclose(fid);B=B;C=C;B(2,5)=1/B(2,5);Y=zeros(A(1,1); %得44的0阵%= %与前一实验同法求导纳阵for

8、 a=1:1:A(2,1)m=B(a,1);n=B(a,2);if B(a,5)0Y(m,m)=Y(m,m)+1/(B(a,3)+j*B(a,4);Y(n,n)=Y(n,n)+1/(B(a,3)+j*B(a,4)/(B(a,5)2);Y(m,n)=-1/(B(a,3)+j*B(a,4)/B(a,5);Y(n,m)=Y(m,n);else if B(a,5)0Q(2*i-1,1)=C(i,3);Q(2*i,1)=0;elseQ(2*i-1,1)=1;Q(2*i,1)=0;endQ(2*A(3,1)-1,1)=A(4,1);Q(2*A(3,1),1)=0;endfid=fopen(answer.t

9、xt,w);fprintf(fid,=节点电压V=n); fprintf(fid,迭代计数tt V1=e1+jf1tttt V2=e2+jf2ttttV3=e3+jf3ttn);%=求W阵= W=zeros(2*(A(1,1)-1),50);for x=1:1:50 %设置迭代次数为50次for i=1:1:A(7,1)k=C(i,1);p=0;q=0;m=0;n=0;for j=1:1:A(1,1)g=real(Y(k,j);b=imag(Y(k,j);e=Q(2*j-1,1);f=Q(2*j,1);p=p+g*e-b*f;q=q+g*f+b*e;m=m+g*e-b*f;n=n+g*f+b*

10、e;endW(2*k,x)=C(k,2)-Q(2*k-1,1)*p-Q(2*k)*q;W(2*k-1,x)=C(k,3)-Q(2*k)*m+Q(2*k-1,1)*n;End%=PQ节点=%=PV节点= for l=1:1:A(8,1)k=C(l+A(7,1),1);p=0;q=0;m=0;n=0;for j=1:1:A(1,1)g=real(Y(k,j);b=imag(Y(k,j);e=Q(2*j-1,1);f=Q(2*j,1);p=p+g*e-b*f;q=q+g*f+b*e;endW(2*k,x)=C(k,2)-Q(2*k-1,1)*p-Q(2*k)*q;W(2*k-1,x)=(C(k,3)

11、2-(Q(2*k-1,1)2+Q(2*k)2);End%=%= %比较是否符合条件Max=0;for i=1:1:2*(A(1,1)-1)Max=max(abs(W(i,x),Max);endif Max=0fprintf(fid,%8.4f+%8.4fjtt ,Q(2*k-1,1),Q(2*k,1);elsefprintf(fid,%8.4f%8.4fjtt,Q(2*k-1,1),Q(2*k,1);endendfprintf(fid,n);end%=%=平衡点功率=k=A(3,1);v=0;j=sqrt(-1);for b=1:A(1,1)m=conj(Y(k,b);p=Q(2*b-1,1)

12、+Q(2*b,1)*j;n=conj(p);v=v+m*n;endWp=(Q(2*k-1,1)+j*Q(2*k,1)*v;end%=节点电压V= %fid=fopen(answer.txt,w);fprintf(fid,节点电压Vn);for i=1:A(1,1)-1k=C(i,1);if Q(2*k,1)=0fprintf(fid,node%d:%8.4f+%8.4fjn,k,Q(2*k-1,1),Q(2*k,1);elsefprintf(fid,node%d:%8.4f%8.4fjn,k,Q(2*k-1,1),Q(2*k,1);endendfprintf(fid,平衡点功率P+jQ=%8.

13、4f+%8.4fj,real(Wp),imag(Wp);fclose(fid);附注：node4.txt %输入数据2、程序输出结果五、实验总结本次实验是我把应用理论知识的重要实践，经过实验，我有两点感想，首先，作为一名工科学生，应该能熟练运用C语言和MATLAB等工具，其次，理论如果不用于实践，就永远不知道理论是用来干什么的，学到头一直是满脑子的浆糊，所以，以我的切身经历建议，把这门实验放在跟课程平行的时间上进行，这样不仅有利于实验开展，也有利于学生更加深刻地学习！/p (关闭)2.在图标菜单点击全选和自动调节缩放 3. 在功能目录表(图9的X)选择 镜像单元 4. 在 形式 选择栏选择 复

14、制 5. 在 镜像平面 选择栏选择 z-x 平面 6. 确认 y : 0 后点击 键 7. 在Icon Menu点击 全选 8. 在 镜像平面 选择栏选择 y-z 平面 9. 确认 x : 0 点击 10. 点击 键 Toggle on X 图 3.9 完成模型的腹板部分开口部详细分析把临时输入的梁单元扩展为板单元来完成如图11的圆形开口部的翼缘、垂直和水平加劲板以及梁的翼缘。1. 在图标菜单点击2. 图标菜单点击 标准视图 设定全局坐标系3. 选择 工作 选项(图10的)4. 在 特性值截面 双击截面号999(曲管形补强 加劲板)5. 在单元工具条点击 将单元扩展6. 在 扩展类型 选择栏选

15、择 线单元 平面单元7. 在 厚度 选择栏选择 10 : 0.0100008. 在 生成形式 选择栏选择 复制和移动9. 在 等间距 的 dx, dy, dz 输入栏输入 0, - 0.1, 010. 在 复制次数 输入栏输入 3 输入11. 点击 键选择属性 12. 在图标菜单点击13. 在 选择属性 选择 截面14. 选择截面号 998(垂直、水平加劲板)15. 点击16. 点击 键 键17. 在 厚度 选择栏选择 15 : 0.01500018. 点击 键19. 与步骤1216相同的方法选择截面号 997(梁的翼缘)20. 在 厚度 选择栏选择 40 : 0.040000 21. 在 复制数量 输入栏输入 422. 点击23. 点击 键 键15例题16