求解决啊！急!

mvpmqc

clc;
( x: i8 b4 o: K" ]- a1 d- H. d# jclear;
tic;
. ^+ X$ b% N2 }6 ?& S& {/ Mfilename=uigetfile('*.*');
5 \' Y( }* m, S a=textread(filename,'','delimiter',',','emptyvalue',0);%读取TXT文件
c=a(1,:);       %平衡节点行数
$ `* z; m# I. S7 G5 K* e n=a(1,1);    %节点数
3 \# C7 R: S9 ]3 k b=find(a(:,1)==0);
line=b(2)-b(1)-1;             %线路共有行数
knum=b(1)+line;               %线路结束行数
7 t+ t- U$ G$ P* \2 r  Z& z lineblock=a(b(1)+1:knum,:);   %线路的参数
- R0 Z3 y: Q. e- I lineNo=lineblock(1:line,1);
linei=lineblock(1:line,2);
linej=lineblock(1:line,3);
liner=lineblock(1:line,4);
linex=lineblock(1:line,5);
6 a9 D$ w+ m7 q8 w; F, Q lineb=lineblock(1:line,6);
; Q  D0 [; G/ g: J* `( M trans=b(3)-b(2)-1;          %变压器共有行数
8 b( ~5 ^1 _0 c k1=b(2)+1;                  %变压器开始行数
k2=b(2)+trans;              %结束行数
transblock=a(k1:k2,:);      %变压器的参数
# R9 h3 y- {8 O, m8 z4 D transNo=transblock(1:trans,1);  
transi=transblock(1:trans,2);
6 J  r1 _& K$ Z+ g8 }6 [3 U, k transj=transblock(1:trans,3);
% g! @& J. X7 P4 o7 P transr=transblock(1:trans,4);
transx=transblock(1:trans,5);
transk=transblock(1:trans,6);
1 Z# D) u$ l1 t! o, |( c& ?" s6 m branch=b(4)-b(3)-1;          %接地支路共有行数
' ^$ a6 {. f8 E; [ k3=b(3)+1;                   %接地支路开始行数
  f# D2 x: L" A- \# Q k4=b(3)+1+branch;            %结束行数  
( B; _; w+ g3 o0 j3 @9 e5 X branchblock=a(k3:k4,:);      %接地支路参数
branchNo=branchblock(1:branch,1);
branchi=branchblock(1:branch,2);
0 O4 @% X' D3 w3 C8 u# r8 h branchg=branchblock(1:branch,3);
' |; `5 s' z9 z$ P/ @3 [( ^- E9 M branchb=branchblock(1:branch,4);
power=b(5)-b(4)-1;            %节点功率共有行数
k5=b(4)+1;                    %节点功率开始行数
6 f6 t- z* C+ S- s9 Z. I k6=b(4)+power;                %结束行数
  S7 K# a1 }# [9 G- @6 ?" H1 E powerblock=a(k5:k6,:);        %节点功率参数
* F  V; s$ Z6 N0 K  S powerNo=powerblock(1:power,1);
6 }) U, H0 Z: }- h% |: a$ f* { PQi=powerblock(1:power,2);       %powi
! |% N. X1 U( z' c2 {$ [; w/ ? PG=powerblock(1:power,3);        %powpgi
QG=powerblock(1:power,4);        %powqgi
4 K0 m0 f5 u+ g8 b2 D- A- ?+ d0 w2 I) k PD=powerblock(1:power,5);        %powpdi
QD=powerblock(1:power,6);        %powqdi
PV=b(6)-b(5)-1;         %PV节点共有行数
5 ?7 C6 N  p7 ]) @. M8 d k7=b(5)+1;              %PV节点开始行数
k8=b(5)+PV;             %结束行数
PVblock=a(k7:k8,:);     %PV节点参数
) ]8 F" q! s' i- s" v PVNo=PVblock(1:PV,1);
PVi=PVblock(1:PV,2);
PVV=PVblock(1:PV,3);
8 t* I% l( a2 _2 Z7 W* `( j PVqmin=PVblock(1:PV,4);
8 f: p, L( b* \ PVqmax=PVblock(1:PV,5);
clear a;
; d4 D* [+ F( K& S# ^) k G=zeros(n,n);
  O4 D0 l) Q5 R% a B=zeros(n,n);
for k=(1:line);           %线路导纳矩阵
     i=linei(k);
4 H" u& Y) z6 a% i     j=linej(k);
+ Q4 e& |: `+ N  S0 [     r=liner(k);
     x=linex(k);
% A: u% H2 }! A     b=lineb(k);
$ X5 v/ ^/ t/ C5 Q/ @     GIJ=r/(r*r+x*x);
     BIJ=-x/(r*r+x*x);
; j. D, Y" m& d0 y     G(j,i)=-GIJ;
     G(i,j)=-GIJ;
) ^( U/ `) |/ F% Z5 n% V' h     B(j,i)=-BIJ;
2 F. `/ l" Y' I: ]/ @! ^     B(i,j)=-BIJ;
     b=lineb(k);
) Z' M) h; w+ F. I/ }% }( q     G(i,i)=G(i,i)+GIJ;
3 E. F1 v# F$ o4 q2 ]- C! S. G$ _& c     B(i,i)=B(i,i)+BIJ+b;
2 w0 g4 `, T( z2 [6 |) @3 m     G(j,j)=G(j,j)+GIJ;
     B(j,j)=B(j,j)+BIJ+b;
* `- X: X; E. C/ V8 k  m, h end
for k=(1:trans);     %变压器导纳矩阵
# Z+ Q+ a% `0 Z# l# S! ^     i=transi(k);
& k) x& R! X* l+ ], R/ T, B9 h     j=transj(k);
9 v+ n( e, ~% Z: W6 I0 y* a0 T     r=transr(k);
5 B0 C1 Z9 {4 s: }  U$ p     x=transx(k);
     j=transj(k);
- @4 P# @7 c# D. j4 R' V     k0=transk(k);
     GIJ=r/(r*r+x*x);
     BIJ=-x/(r*r+x*x);
0 U. _, G9 v7 w5 m' z7 R( \6 v     G(i,j)=-GIJ/k0;
  C) h" R% L, ?     G(j,i)=G(i,j);
     B(i,j)=-BIJ/k0;
     B(j,i)=B(i,j);
- G) o* L7 ?- Z+ S7 f; Z     G(i,i)=G(i,i)+GIJ/k0/k0;
" v5 H, L4 d+ M- @$ b' i     B(i,i)=B(i,i)+BIJ/k0/k0;
     G(j,j)=G(j,j)+GIJ;
) a# ]" D9 D) n+ I     B(j,j)=B(j,j)+BIJ;
( k$ O1 V5 w. }- D: ]2 [2 j end
for k=(1:branch);      %接地支路导纳矩阵
7 ^! N7 C# F  h9 p" X4 _$ j     i=branchi(k);
     G(i,i)=G(i,i)+branchg(k);
     B(i,i)=B(i,i)+branchb(k);
end
6 ?# [; r* R6 v0 h% Y* J; C Y=sparse(zeros(n,n));
- s' u* `' ]# F4 O: V9 u A=sparse(zeros(n,n));
clear j;
# L; H0 V6 f' N5 g' K5 UX=G+j*B;  
X=sparse(X);
3 _5 @, y" h$ {. d5 Z0 VY=(abs(X));      %导纳矩阵完成
A=angle(X);     %相角
dp=sparse(PQi,1,(PG-PD)./100,n,n);   %节点注入有功功率
# x: G3 r" r$ Jdq=sparse(PQi,1,(QG-QD)./100,n,1);   %节点注入无功功率
V0=sparse(ones(1,n)');   %节点电压大小初值
4 k/ T6 O  @6 S+ QVA=sparse(zeros(n,1));     %相角初值
V0(c(2))=c(3);    %平衡节点电压
' v8 \/ U% j3 m1 a- }0 nV0(PVi)=PVV;     %pv节点电压
ee=1;        %精度
k=1;         %迭代次数
  while(ee>c(4)&&k<20)
ij=repmat(VA,1,n)-repmat(VA',n,1)-angle(X);  %角
V=V0.*(cos(VA)+j*sin(VA));  %节点电压复数形式
1 ?/ m: [0 k# K& I( C/ mdelt0=conj(X*V);
: n! b& k; {: S& Q% j1 a1 `' Idelt1=V.*delt0;   %代入节点电压求出的功率
deltp=dp-real(delt1);  %有功修正量
deltq=dq-imag(delt1);  %无功修正量
deltp(c(2))=0;          %有功修正量(处理平衡节点)
deltq(c(2))=0;          %无功修正量（处理平衡节点）
- r. A7 d* |3 E( v. S/ I. adeltq(PVi)=0;           %无功修正量（处理pv节点）
+ [/ T8 m# p6 L7 U( Sdelt=[deltp;deltq];     %p,q修正量数组
1 `7 u# }8 G4 w: a+ ?0 }h=diag(diag(V0)*(Y.*sin(ij))*V0)-diag(V0)*(Y.*sin(ij))*diag(V0);
1 v: Q) L# O/ znn=-diag(V0)*(Y.*cos(ij))-diag((Y.*cos(ij))*V0);         %n矩阵
, {. f+ E, A! f" z) t" ]jj=-diag(diag(V0)*(Y.*cos(ij))*V0)+diag(V0)*(Y.*cos(ij))*diag(V0);
l=-diag(V0)*(Y.*sin(ij))-diag((Y.*sin(ij))*V0);
%雅克比矩阵完成
; ^/ k0 h! V' W9 ]% fnn(:,PVi)=0;              %处理pv节点
9 j9 P0 E( \, h0 B" }  o& q* M2 Mjj(PVi,:)=0;              %处理pv节点
l(PVi,:)=0;               %处理pv节点
l(:,PVi)=0;               %处理pv节点
' R% K7 J' S1 Kl=l+sparse(PVi,PVi,1,n,n);  %形成l
+ R( O. J3 p- K5 }ykb=[h nn;jj l];          %初步形成雅可比矩阵
ykb(c(2),:)=0;          %处理平衡节点a
; ~4 F: T+ U7 f( E- ^/ Oykb(:,c(2))=0;          %处理平衡节点a
; \! Y! F' |& uykb(c(2),c(2))=1;       %处理平衡节点a
ykb(c(2)+n,:)=0;          %处理平衡节点v
ykb(:,c(2)+n)=0;          %处理平衡节点v
8 N0 T: L9 j8 {ykb(c(2)+n,c(2)+n)=1;  %处理平衡节点v
jie=ykb\delt;
4 l1 V: H2 ^) ]/ Z8 gdelta=jie(1:n);         %相角修正量
; e& g3 E5 _' f- @' ~2 vdeltv=jie((n+1):2*n);   %电压大小修正量
* _' h' [" m' p8 f+ l3 j1 Z' X- iVA=VA-delta;
) N8 ?$ y8 u9 vV0=V0-deltv;
ee=max(abs(deltv))
k=k+1;
' S- |: t" e) W, U  end
VA=VA.*180./pi            %相角
V0                      %电压幅值
; D7 j3 w8 |& p: @: f. c  utoc;
求解释，求解答，改了很多次都没用
/ A* J1 m& b" L; m% L1 J/ U( k' r??? Error using ==> minus
Matrix dimensions must agree.
) \. B# _. Z- M$ v" `9 `$ f. y7 {' UError in ==> liweiming2 at 122
  V5 j5 L4 x% G9 Q" ndeltp=dp-real(delt1);  %有功修正量

mingyu

编程高手解释一下哦！

sdad

回复 1# mvpmqc

- b  R5 l0 h) a3 u* O0 K4 z2 n0 t
# |5 ?9 Q  j2 g" V& S您的问题，有可能是读取的数据格式，通常会错误，多半原因是格式问题，

# `5 Q, l9 U! Z/ o2 P% ]除非你上传你的文本文件才知道错误在哪里。