谁能帮我指教一下改怎么编潮流计算

zhouyou529

clear;
9 m8 D+ g0 I5 X  @/ gYb=[6.25-18.75i  -5+15i        -1.25+3.75i    0             0
    -5+15i       10.834-32.5i  -1.667+5i     -1.667+5i    -2.5+7.5i
1 G& l0 `6 W' H    -1.25+3.75i  -1.667+5i     12.917-38.75i -10+30i      0
    0            -1.667+5i     -10+30i       12.917-38.75i  -1.25+3.75i
    0            -2.5+7.5i     0             -1.25+3.75i    3.75-11.25i];
Pknow=[0;0.2;-0.45;-0.4;-0.6];Qknow=[0;0.2;-0.15;-0.05;-0.1];e1=1.06;f1=0;%已知量
e=[1.06;1;1;1;1];f=[0;0;0;0;0];
( V8 g5 M9 b4 z8 j" Xk=0;%跌代次数
# v. d6 X6 ?3 z8 K2 M1 g- dPQX=[1;1;1;1;1;1;1;1];%修正量   
P=zeros(4,1);Q=zeros(4,1);I=zeros(4,1);YAKEBI=zeros(8,8);
1 ?8 z  p4 ?1 I( h) a7 H; mwhile (max(PQX)>=0.1|min(PQX)<=-0.1)
    k=k+1;
% Z5 m& Y* k: w    for n=1:4
        for m=1:5
7 r& K# E; f" y+ [1 A            P(n)=P(n)+e(n+1)*(real(Yb(n+1,m))*e(m)-imag(Yb(n+1,m))*f(m))+f(n+1)*(real(Yb(n+1,m))*f(m)+imag(Yb(n+1,m)*e(m)));
            Q(n)=Q(n)+f(n+1)*(real(Yb(n+1,m))*e(m)-imag(Yb(n+1,m))*f(m))-e(n+1)*(real(Yb(n+1,m))*f(m)+imag(Yb(n+1,m)*e(m)));
! c! C( J+ ^( Y2 k        end
        PQX(2*n-1)=Pknow(n+1)-P(n);
        PQX(2*n)=Qknow(n+1)-Q(n);
8 X3 g; `. u( c% l- S        I(n)=(P(n)-Q(n)*i)/(e(n+1)-f(n+1)*i);
! k; c6 A7 K! w! d" h. m    end
2 D& B. |/ a! ?6 |3 I: I$ b    for n=1:4
        for m=1:4
$ w' ?; u  O8 `* \% ^! y2 r+ t6 o            YAKEBI(2*n-1,2*m-1)=-imag(Yb(n+1,m+1))*e(n+1)+real(Yb(n+1,m+1))*f(n+1);%Hnm
: `7 B3 [* N' a- y6 `- E5 a            YAKEBI(2*n-1,2*m)=real(Yb(n+1,n+1))*e(n+1)+imag(Yb(n+1,n+1))*f(n+1);%Nnm
: I' d* m* \- ]& ]1 ?            YAKEBI(2*n,2*m-1)=-real(Yb(n+1,m+1))*e(n+1)-imag(Yb(n+1,m+1))*f(n+1);%Jnm
            YAKEBI(2*n,2*m)=-imag(Yb(n+1,m+1))*e(n+1)+real(Yb(n+1,m+1))*f(n+1);%Lnm
        end
% i: c; Q7 y! a& @, I6 e4 ?" ]        YAKEBI(2*n-1,2*n-1)=-imag(Yb(n+1,n+1))*e(n+1)+real(Yb(n+1,n+1))*f(n+1)+imag(I(n));%Hnn
0 ~4 x4 ]( y* m) a4 P        YAKEBI(2*n-1,2*n)=real(Yb(n+1,n+1))*e(n+1)+imag(Yb(n+1,n+1))*f(n+1)+real(I(n));%Nnn
        YAKEBI(2*n,2*n-1)=-real(Yb(n+1,n+1))*e(n+1)-imag(Yb(n+1,n+1))*f(n+1)+real(I(n));%Jnn
        YAKEBI(2*n,2*n)=-imag(Yb(n+1,n+1))*e(n+1)+real(Yb(n+1,n+1))*f(n+1)-imag(I(n));%Lnn
5 r9 T3 L2 ^* G7 `3 l; ?* G! s    end% build the yakebi matrix
/ Y8 F" V; x# t: p+ u' s    Dfe=inv(YAKEBI)*PQX; %修正fe
& x/ c' M' b5 q. K( v    for n=1:4
        f(n+1)=f(n+1)+Dfe(2*n-1);
* G2 y  W" Y2 m! z/ g        e(n+1)=e(n+1)+Dfe(2*n);
4 R) d  W; S1 h7 R/ c5 y" w    end
    f
    e
, X: \$ \8 Z# G# H/ Q4 t5 vend
" l% a! d% {( @8 ^+ o% a. ^4 p5 ?
& h" B9 W5 `' H- T6 m
9 k0 v! [) A' p/ q  T, M运行会出现
Warning: Matrix is singular to working precision.
(Type "warning off MATLAB:singularMatrix" to suppress this warning.)
> In f:\MATLAB6p5\work\yy.m at line 35

; `! ?9 g# |# J$ A我这个程序是按照陈珩那本稳态分析的例子编的，想试试。完全没考虑PV节点
! E/ J$ k7 G0 w这个程序要怎么改进啊？

wangmuyidong

我没有时间仔细给你检查，但是想告诉你，你的雅克比矩阵的形成太麻烦了，极有可能存在错误，并且你得程序也不通用，雅克比矩阵的形成要仔细分析其各个表达式的关系，可以用矩阵运算直接形成。而没有必要用循环。这样会非常耗时的。