MATLAB编程1.docVIP

电力系统基础课程设计报告 电力系统潮流计算的MATLAB编程 【程序流程图】： 【源程序】： disp( -----------《电力系统潮流计算的MATLAB编程》--------------------); disp(请输入电力系统参数！); n=input(请输入节点数:n=); nl=input(请输入支路数:nl=); isb=input(请输入平衡节点号:isb=); pr=input(请输入误差精度:pr=); B1=input(请输入由各支路参数形成的矩阵:B1=); B2=input(请输入各节点参数形成的矩阵:B2=); Y=zeros(n);U=zeros(1,n);cta=zeros(1,n);P=zeros(1,n);Q=zeros(1,n);C=zeros(1,13);D=zeros(1,13); DP=zeros(1,13);DQ=zeros(1,9);UU=zeros(1,n);c=zeros(1,n);Si=zeros(nl);Sj=zeros(nl);DS=zeros(nl); %节点导纳矩阵Y： for i=1:nl p=B1(i,1);q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)+1j*B1(i,4)); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1/(B1(i,3)+1j*B1(i,4))+1j*B1(i,6)/2; Y(p,p)=Y(p,p)+1/(B1(i,3)+1j*B1(i,4))+1j*B1(i,6)/2; end G=real(Y);B=imag(Y); %为U、cta、P、Q赋初值： for i=1:n cta(i)=B2(i,6); %电压相角 U(i)=B2(i,5); %电压模值 P(i)=B2(i,3); Q(i)=B2(i,4); %节点PQ end %迭代开始： ICT1=0;IT2=1; while IT2~=0 IT2=0;t1=1;t2=1; for i=1:n if i~=isb C(i)=0; D(i)=0; for j1=1:n C(i)=C(i)+U(i)*U(j1)*(G(i,j1)*cos(cta(i)-cta(j1))+B(i,j1)*sin(cta(i)-cta(j1))); %求Pi D(i)=D(i)+U(i)*U(j1)*(G(i,j1)*sin(cta(i)-cta(j1))-B(i,j1)*cos(cta(i)-cta(j1))); %求Qi end DP(t1)=P(i)-C(i); %求1~13点的ΔP t1=t1+1; if B2(i,2)==1 DQ(t2)=Q(i)-D(i); %求1~9点的ΔQ t2=t2+1; end end end t1=t1-1;t2=t2-1; %此时t1=n-1=13,t2=PQ节点数=9 DPQ=[DP;DQ]; %误差矩阵 %判断误差精度，决定是否继续循环运算: for i=1:t1+t2 if max(abs(DPQ(i,:)))pr IT2=IT2+1; end end %雅克比矩阵程序： H=zeros(t1,t1);N=zeros(t1,t2);K=zeros(t2,t1);L=zeros(t2,t2); %此时t1=13,t2=9 %求雅克比矩阵中的H : for i=1:t1 for j1=1:t1 if j1~=isbj1~=i H(i,j1)=0-U(i)*U(j1)*(G(i,j1