|
马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
[x]=xlsread('pqinput.xls','A2:A2');& i+ N# W+ K8 Q7 B0 m- b( z3 K
[y]=xlsread('pqinput.xls','B2:B2');) S4 T% c2 ?, Q) ]9 D# ]
e=xlsread('pqinput.xls','B4:B4');
: |( j- t+ V- c# J; |9 L[point]=xlsread('pqinput.xls','D3:H100');5 A$ j$ _4 _) @
[zhilu]=xlsread('pqinput.xls','J3:R100');: E% s B$ ~0 I- ]5 ~' S7 {1 X
TYPE=zeros(x,1);
& Q3 ]$ x. W9 N& Y% T; b% Y3 f$ sU=zeros(x,1);: Q6 n1 B- f3 M3 H$ O+ R% @
a=zeros(x,1);
- j3 y K' ~; t+ f* O# e0 \4 MP=zeros(x,1); J3 }$ {: `' y. o$ n1 M
Q=zeros(x,1);) Q; l# q( B0 [: e
I=zeros(y,1);* r7 z2 p1 H4 ~7 K' O
J=zeros(y,1);
, o2 G" X- Z& C, w0 x: j5 P$ CRij=zeros(y,1);9 R8 h5 h& Y" R0 k/ Y
Xij=zeros(y,1);' K6 v/ }2 T5 y+ }- @0 l$ W
Zij=Rij+j*Xij;
' }- X- p: g8 X" O. P+ yY=zeros(x); # m9 V8 m6 C+ p% S$ z
G=zeros(x);
$ I/ @. V: M9 ^9 SB=zeros(x); 1 [# r {' h( j' Y4 q `
B0=zeros(y,1);! L7 e ~ a% X$ J8 U# n* w- s
RT=zeros(y,1);
& c- q G S% c9 VXT=zeros(y,1);
/ G' c8 T% X" ^( C& DZT=RT+j*XT;
! Y, ]0 G( a5 i/ q1 w6 FKT=zeros(y,1);
5 P: e1 v3 w- @" K0 R%------------------------------矩阵赋初值:
~9 G* B% T$ [/ q! y/ _, \ e9 OTYPE=point(:,1);% M7 D& r7 W |
U=point(:,2);9 ]% k( g8 r! |) d0 a/ }9 s% x9 }2 E
a=point(:,3);
9 c- s2 u1 k, s0 P# k; qP=point(:,4);0 K* [/ Q+ T3 j, Y. `
Q=point(:,5);, B3 H5 S! f3 p$ I/ L
I=zhilu(:,1);: t! |" T0 U3 b, c
J=zhilu(:,2);0 T- b6 a, a* z
Rij=zhilu(:,3);- I+ F+ k' Z( m8 ]# ^
Xij=zhilu(:,4);6 \) x* E) A7 x ?. L" E
Zij=Rij+j*Xij;
4 T8 M `, ?- R# H& B9 u% ^B0=zhilu(:,5);
/ T& Z$ @3 C* _& ^. y* }RT=zhilu(:,6);
% I: F( T, G" U1 @& Q4 MXT=zhilu(:,7);5 t& W, a6 C% y( X" Q5 G) g
ZT=RT+j*XT;
, t# i" A* |: c& ~# c6 XKT=zhilu(:,8);4 V: Q) k* P- c+ w4 u$ `
W=zhilu(:,9);1 v$ v0 B/ S( }( n: O; \& }
%------------------------------求节点导纳矩阵Y) C4 @ {9 y6 W- l6 z( q
for m=1:y ! \ C2 T1 b: ]3 u8 v
if KT(m)==0) v9 e2 S% G7 Z6 F4 f. }
Y(I(m),J(m))=-1/Zij(m);2 t8 a8 u2 j8 j4 g
Y(J(m),I(m))=-1/Zij(m);
# w! c$ b0 T( S else
' ?0 V/ P/ X+ r Y(I(m),J(m))=-1/(KT(m)*ZT(m));
3 `4 d H8 g0 l9 k) i! a( n# G Y(J(m),I(m))=-1/(KT(m)*ZT(m));
* }, `) ?% q I: ]* q& A: V& } end# m& E" u. @& E' o" W
end- H0 M6 `3 E f( g& Q3 n
for m=1:x 9 O' y- S; s4 G- ?+ I8 ]
for n=1:y. q/ j; ?8 N) @# A B5 y x; ]
if KT(n)==0
0 `( R) l$ M/ \' A) e if(I(n)==m|J(n)==m)
* z# @, [& G0 ~- S% r Y(m,m)=Y(m,m)-Y(I(n),J(n))+j*B0(n)/2;
2 F0 E& i; W! J end* K+ k0 p k, R) D
else; o- V6 h) j9 \+ O# d0 }6 d: p
if W==0
. X# L& ~0 g, n/ E- R/ l if I(n)==m
- x7 p$ R7 A: z4 ?& T Y(m,m)=Y(m,m)-Y(I(n),J(n))+(KT(n)-1)/KT(n)*(1/ZT(n));
! I7 m9 a. Y4 K2 W elseif J(n)==m
! Y+ X9 _' J* l4 X) b/ |6 H% d Y(m,m)=Y(m,m)-Y(I(n),J(n))+(1-KT(n))/(KT(n)^2)*(1/ZT(n)); - L+ i. \. O) I
end
7 Q0 }" ]+ r) H4 P else- b* c u: J5 f/ ?! k
if I(n)==m/ B" ]: P+ \0 m3 W4 ]: X& N3 \
Y(m,m)=Y(m,m)-Y(I(n),J(n))+(1-KT(n))/(KT(n)^2)*(1/ZT(n)); ) C3 v% ]4 M8 s
elseif J(n)==m
" O: x8 ^" m2 i$ A Y(m,m)=Y(m,m)-Y(I(n),J(n))+(KT(n)-1)/KT(n)*(1/ZT(n));
% R8 u: \7 |7 }% y else Y(m,m)=Y(m,m);
! u" k4 R" m: Q# Y+ B end% { e1 R; S# L) K
end" N% }1 Y$ t! ^5 Y9 `: |9 @9 }
end' ~6 i2 c0 Z$ p ?7 O, \
end" [ }4 T! j" T* ]& z
end- I2 N# D9 _2 C7 ^) {
G=real(Y);
- J' Y+ P0 P, Q, x. I%-----------------------求B'矩阵及其逆矩阵B1) L0 E! t; f8 S* V# W
B=imag(Y);& u" X' B& }( n, E
ph=find(TYPE(:,1)==3);
* G9 U! N b+ q( m- z BB=B;' g1 f+ E. g% f' m' W
BB(:,ph)=[];; U' z% x# b: Z% P- L5 f! n
BB(ph,:)=[];4 c; H6 F! L2 P8 K+ u
B1=BB;
0 n @ L3 [3 z# k B1=inv(B1);& t4 v4 e0 S1 J2 ?
disp('B矩阵')- V9 x9 C/ e5 z9 i; ^3 g% n; d; x
disp(BB)
& {% A: s8 w0 \2 T8 J8 [disp('B逆矩阵')8 C; M6 d$ A7 ^9 J* x( J
disp(B1)9 g: e2 c7 [/ v" _$ M
%-----------------------%求B''及其逆矩阵B2% r/ W3 L) T. s- C
phpv=find(TYPE(:,1)>1);# h( w( K1 }3 j6 P- h
BB=B;
8 Y9 C! L, F \" O z, X BB(:,phpv)=[];# R2 ~* ?: l% D
BB(phpv,:)=[];2 K4 }. N, z* [* u4 c5 a/ t% G
B2=BB; 3 Z" H6 ?9 f! K" S% g+ q \
B2=inv(B2);
+ x1 ~' L0 v, V ldisp(BB); S( j# o# A2 L* A8 C; m
disp(B2)( `% s' B9 J. g
%-------------计算各节点有功功率不平衡量deltaPi
! h: }+ ? @9 J8 ~- i k=0; 6 {; y9 h6 c) y h4 k( k
kp=1;* O+ H! C: N5 {" I: j/ T$ l! c
kq=1;
_" B* o* P1 }; ]+ ?) J4 d6 fwhile(((kp~=0)||(kq~=0))&&k<=20)
* x8 P! P9 o9 C5 z kp=1;
5 |8 K# Y( p3 g( K) `4 L kq=1;8 I7 u; ]( D) b
notph=find(TYPE(:,1)<3);( _6 @# j% z1 l, C
deltaPi=zeros(x-1,1);
! x7 l+ b! x& \. i+ C* h' |) A, H pq=find(TYPE(:,1)==1);! Y a/ `- D3 i; b o
pqnum=size(B2);
" W$ u* S( m. q e. f3 X pqnum=pqnum(1);, b# O; Q& C( x
deltaQi=zeros(pqnum,1);7 C& J( z1 [6 d0 G" z; u- ^- ~) V; g
for m=1:(x-1)
$ U: R# z8 q" s' O sum1=0;' g. g6 E2 p5 Q$ P
for n=1:x! }) F: B1 Y& F' n, N
sum1=sum1+U(notph(m))*U(n)*(G(notph(m),n)*cos(a(notph(m))-a(n))+B(notph(m),n)*sin(a(notph(m))-a(n)));
. s1 L4 k1 ]" O) K end
; ^7 T, z6 `( c* y deltaPi(m)=P(notph(m))-sum1;
8 N' q8 ?" I1 U end. Z& S$ @3 l0 r1 L
disp('有功功率不平衡量')# d6 x2 x; U. M; b! {! Y/ N! y, o! J
disp(k)' V- [2 x9 G" J* F1 N+ M, H# v
disp(deltaPi)
# e S2 {5 d* w max1=max(abs(deltaPi));
' M" R9 I$ K* d) i if max1<=e 1 n$ r ]3 W8 d8 {
kp=0;
7 u( x2 ?8 N& K3 _" G+ P0 D. | if kq==01 p% _. M5 y6 ^0 s; ~: @0 K* C7 \3 j
break V$ i. f' l) w. ?9 r0 y
else( h: F2 g6 M- ?! t5 V2 E
for m=1:pqnum
, M$ B9 Y2 y1 @: s3 \ sum2=0;; G9 z. h$ N( L9 N
for n=1:x, p3 e3 X" Z! C: _: P d
sum2=sum2+U(pq(m))*U(n)*(G(pq(m),n)*sin(a(pq(m))-a(n))-B(pq(m),n)*cos(a(pq(m))-a(n)));
b, N/ a# y8 m J4 C) X end6 F2 p" i) ~# \: N% A0 S
) u/ z- u- c. n/ ]" q
deltaQi(m)=Q(pq(m))-sum2;/ Q" v* @2 |" R3 E+ i
end' Y ^* l" g4 _1 E: Z& A
max2=max(abs(deltaQi));( \& }8 i2 U0 r
if max2<=e . @; I7 P+ [; Q0 e; Q) m4 l; L
kq=0;5 D% H& @ g% g- }
if kp==0
8 p9 L! a1 S" d7 @' P( z break7 E* P) h# v$ [5 u: o
else
$ _. i5 e4 n8 z& H- O% | k=k+1;; u: j; ~" ~" K0 I/ c3 T+ Q; S/ m- Q
end* a* \3 R7 a H3 u/ F4 i
else2 y) W B% y3 r7 `, P8 N* q
Uq=U;
( v. c: U# l/ V( C Uq(phpv)=[];
9 v2 p# g: `+ v" o! C Upq=Uq;
3 d+ V |0 s+ a' } deltaU=-B2*(deltaQi./Upq); e( w! X2 r# F8 R' w
for m=1:pqnum , A$ d$ C6 y) Y( _+ D$ I
U(pq(m))=U(pq(m))+deltaU(m); P5 |: d( G0 d. C1 n
end
$ a0 ~% Z' a% Q: X2 S kp=1;
$ O' P8 m) f1 ? k=k+1;
7 X2 L T8 f/ D4 C5 }( r end
" c( _* t1 v) U! G7 x end
- A) w$ m0 k# F( t, C else
9 e# E+ k/ P2 X# R+ v Up=U; 3 z' a" \% N$ z' M
Up(ph)=[];
0 t& n$ Z# T% ~ Unotph=Up;
3 ~/ h7 @( |" L$ f" H1 O" z) T# A deltaa=((-B1*(deltaPi./Unotph))./Unotph);
" e% p! `2 f" W6 [" R; m1 x0 Q2 ? for m=1:(x-1)5 L' W* S$ l5 u* T. Q0 | \+ R
a(notph(m))=a(notph(m))+deltaa(m);6 J+ j' {8 S) f
end
& V3 b3 o1 H2 ^disp('电压相角修正')+ U- D4 t, n0 h+ U4 U F7 e8 X
disp(k)& E! M1 e: N7 D1 G% u. J, f
disp(deltaa)/ U9 T& g' V6 x- M) `- t
disp(a)% B# ]& p9 C/ O
kq=1;
* w. | v0 ~3 X for m=1:pqnum5 u3 H1 c0 G) J+ F7 A
sum2=0;2 k( E% P+ j1 j5 C) ]% v4 D; k& Q
for n=1:x) p! t5 o5 e7 q; @
sum2=sum2+U(pq(m))*U(n)*(G(pq(m),n)*sin(a(pq(m))-a(n))-B(pq(m),n)*cos(a(pq(m))-a(n))); 7 B, |. j& O7 [* c* p# U
end3 j A! J5 u) g, I, q$ e8 d; K2 q
deltaQi(m)=Q(pq(m))-sum2;
4 U$ X' y# `4 A ~& I2 B end 6 l( c* ^4 A$ ~) b0 p+ t, t
disp('wu功功率不平衡量')" a& ?* y' D+ U9 |9 Z
disp(k)* | c# a d: Z4 {
disp(deltaQi) 4 y( f4 l* d- e: `! D
Uq=U;
- S* h3 q/ u4 j Uq(phpv)=[];0 o1 C6 } S: P/ _1 [; c: o
Upq=Uq;
6 F, |, g- {: }4 h, w9 d deltaU=-B2*(deltaQi./Upq);
4 h, p k8 A2 Z for m=1:pqnum
: b: s; {" p: t2 } U(pq(m))=U(pq(m))+deltaU(m);
- o4 w* |2 P1 E8 | end
3 B0 E2 Z& ~8 k7 t# J" hdisp('电压修正')7 W' H4 Z* v/ S
disp(k)5 x! u7 T. f+ q2 u
disp(deltaU)+ E3 V- H* P3 D
disp(U)' C& \* ]. R9 L' I9 ^" m, p; ^
kp=1;
* z8 k. h/ E0 j& G3 a0 D k=k+1;
* {8 I, k) \+ T; G# ^& ^ end
- g( Q G' e* R$ U" e& oend* k, S1 k; P2 f$ v$ a
sum3=0+j*0;; G+ y; m+ ?% c2 b
for m=1:x9 c1 h, _. L; @
sum3=sum3+conj(Y(ph,m))*(U(m)*cos(a(m))-i*U(m)*sin(a(m)))
% j2 x/ ^& @1 @7 d; k6 o" N end
9 ]9 }. n: g" t. T$ I8 |+ F Sph=(U(ph)*cos(a(ph))+j*U(ph)*sin(a(ph)))*sum3;1 K- K0 q2 Y+ K+ X! K* m
%-----------------求线路功率Sij和Sji& a8 [; |8 e* Z
Sij=zeros(y,1);
7 |, f+ m0 m- c$ C. X Sji=zeros(y,1);
+ E% J% c; n7 j7 u for m=1:y6 n0 H( l \% s3 u0 q7 h% e
if KT(m)==0 Sij(m)=(U(I(m))*cos(a(I(m)))+i*U(I(m))*sin(a(I(m))))*((U(I(m))*cos(a(I(m)))-i*U(I(m))*sin(a(I(m))))*(-i*B0(m)/2)+(U(I(m))*cos(a(I(m)))-i*U(I(m))*sin(a(I(m)))-U(J(m))*cos(a(J(m)))+i*U(J(m))*sin(a(J(m))))*conj(-Y(I(m),J(m)))); 0 f; }" u, |' q- }# v/ E$ }
Sji(m)=(U(J(m))*cos(a(J(m)))+i*U(J(m))*sin(a(J(m))))*((U(J(m))*cos(a(J(m)))-i*U(J(m))*sin(a(J(m))))*(-i*B0(m)/2)+(U(J(m))*cos(a(J(m)))-i*U(J(m))*sin(a(J(m)))-U(I(m))*cos(a(I(m)))+i*U(I(m))*sin(a(I(m))))*conj(-Y(I(m),J(m))));
1 ?; Q+ d- L$ G5 O% j; I else
: X/ C3 |) K- ]1 I8 u& d/ | if W==1 Sij(m)=(U(I(m))*cos(a(I(m)))+i*U(I(m))*sin(a(I(m))))*((U(I(m))*cos(a(I(m)))-i*U(I(m))*sin(a(I(m))))*((1-KT(m))/(KT(m))^2*conj(1/ZT(m)))+(U(I(m))*cos(a(I(m)))-i*U(I(m))*sin(a(I(m)))-U(J(m))*cos(a(J(m)))+i*U(J(m))*sin(a(J(m))))*conj(-Y(I(m),J(m))));
+ H6 C$ x5 q. c% H7 E4 S& O/ f Sji(m)=(U(J(m))*cos(a(J(m)))+i*U(J(m))*sin(a(J(m))))*((U(J(m))*cos(a(J(m)))-i*U(J(m))*sin(a(J(m))))*((KT(m)-1)/KT(m)*conj(1/ZT(m)))+(U(J(m))*cos(a(J(m)))-i*U(J(m))*sin(a(J(m)))-U(I(m))*cos(a(I(m)))+i*U(I(m))*sin(a(I(m))))*conj(-Y(I(m),J(m))));' C: W X1 W; {& R$ ?
else Sij(m)=(U(I(m))*cos(a(I(m)))+i*U(I(m))*sin(a(I(m))))*((U(I(m))*cos(a(I(m)))-i*U(I(m))*sin(a(I(m))))*((KT(m)-1)/KT(m)*conj(1/ZT(m)))+(U(I(m))*cos(a(I(m)))-i*U(I(m))*sin(a(I(m)))-U(J(m))*cos(a(J(m)))+i*U(J(m))*sin(a(J(m))))*conj(-Y(I(m),J(m))));0 ~0 A% Y* I- b- t
Sji(m)=(U(J(m))*cos(a(J(m)))+i*U(J(m))*sin(a(J(m))))*((U(J(m))*cos(a(J(m)))-i*U(J(m))*sin(a(J(m))))*((1-KT(m))/(KT(m))^2*conj(1/ZT(m)))+(U(J(m))*cos(a(J(m)))-i*U(J(m))*sin(a(J(m)))-U(I(m))*cos(a(I(m)))+i*U(I(m))*sin(a(I(m))))*conj(-Y(I(m),J(m))));4 l6 j i3 D: p4 g9 K" g
end1 K2 A# H/ f' m o- B' D7 Q6 `
end
7 [" C2 ], x3 I+ v, c( y1 R end
8 _0 [2 D0 d2 R/ ndeltaSij=Sij+Sji;
! ~. N% X. f, g3 h+ m- eS=zeros(x,1)+i*zeros(x,1);
6 @: \+ Z8 f: d! d2 {! efor b=1:x
' D9 v7 K1 Z7 r: r6 ^0 F6 Y. r$ c for m=1:y
; A9 E* ^/ W* a$ O2 O$ e, n if I(m)==b ) \' L$ h1 ]' k8 g: O: F
S(b)=S(b)+Sij(m);
, q- I' C! F6 D4 ~) d `' T else4 o2 e4 v) L" V7 E
if J(m)==b* O" ?3 h. ~8 T3 T/ Y, {
S(b)=S(b)+Sji(m);" c. S- V7 d$ c! e1 Q
else
; a2 y2 {- x# v# e" y1 w- b3 r! z S(b)=S(b);
2 N7 h* U2 P0 V0 c end! ]" j' c1 ?' c+ u( c1 s2 K! [
end
c( U: h# C3 Q) ^4 p! _& {5 l5 Dend$ x9 G4 a8 _; _1 b- c& y
end
5 Y2 i; _0 ^; Z, B; k; P: J; BP=real(S); 3 C. q8 m8 _! C+ W& ^8 f
Q=imag(S);
2 z/ r5 Q# Z2 H# k" o+ dsumdeltaS=sum(S);
6 ?3 H6 Z& O1 a" \/ T. r# q; Cfid=fopen('shiyan1 d! r! e0 Q" q& _% l, ~
fprintf(fid,' ******************************电力系统06 白云霄 0603010213********************************\n');jieguo.txt','wt');/ v2 v0 z, [: O$ P, T
fprintf(fid,' *****************潮流上机输出结果*************\n'); G, T: l6 V/ r; ^# z# j& `; S
fprintf(fid,' 迭代次数k为: %d \n',k);4 s! [/ d4 h4 T: C
fprintf(fid,' ================================================\n');
+ b9 N9 x# [5 }" }fprintf(fid,' 平衡节点%d的复功率Sph为: %f+j*(%f) \n',ph,real(Sph),imag(Sph));
3 U7 q9 ?8 I( n* X7 L Dfprintf(fid,' ================================================\n');# J- l0 P9 ^- n
fprintf(fid,' 节点电压U为: \n');
0 x" g. W8 r/ E/ Hfor m=1:x) H* u* ~. Y5 ^( g5 x c
fprintf(fid,' 第%d个节点电压: %f\n',m,U(m));3 `: i5 S m( R6 w+ e' w
end
3 U) h9 s7 V+ Z6 J2 z/ W8 X L0 ]; v) Bfprintf(fid,' ================================================\n');
9 ]! I0 t/ M# \, k7 m; @! ofprintf(fid,' 节点相角a为: \n');' Q6 R2 H9 Q4 f( v' D( J
for m=1:x- b, r" V$ [& w) s
fprintf(fid,' 第%d个节点相角: %f\n',m,a(m));( `1 ?/ v& q( g; k P: A
end/ |4 P* |) r: c& F' x3 U& @
fprintf(fid,' ================================================\n');
, Q; q0 F7 q; k3 zfprintf(fid,' 节点复功率S为: \n');. x- n. N9 e- b3 _7 f0 G
for m=1:x4 f& G/ Y& Z6 y7 X: m5 ?
fprintf(fid,' 第%d个节点复功率: %f+i*(%f)\n',m,real(S(m)),imag(S(m)));
2 D3 w) {2 \" T! f7 F8 |2 ]5 tend0 z' P# J- [, T9 F1 S
fprintf(fid,' ================================================\n');6 W# D' O4 h5 `& |9 o
fprintf(fid,' 节点有功功率P为: \n');; K, M! r% y, _) b$ g$ L# M# c
for m=1:x* M+ {' G- t+ M, m7 s
fprintf(fid,' 第%d个节点有功功率: %f\n',m,P(m));, S; j" |5 r, f+ X- Z( q
end( s7 V+ [1 k* h( H5 f5 o
fprintf(fid,' ================================================\n');
5 O* w" n8 d. |$ o( Rfprintf(fid,' 节点无功功率Q为: \n');
/ M1 \! \4 }+ q( R8 Nfor m=1:x
( x q/ j7 i: L1 e fprintf(fid,' 第%d个节点无功功率Q: %f\n',m,Q(m));" Q3 X' O$ n! t
end( S; @% s1 N7 R
fprintf(fid,' ================================================\n');
: i. }7 m" D0 i7 {% L' Ofprintf(fid,' 线路功率Sij和Sji为: \n');2 Y% `$ l9 d3 P, n" W0 u7 h$ R
for m=1:y
" j6 A& Z% Z* W9 F) E fprintf(fid,' 节点%d到节点%d的功率为: %f+i*(%f) \n',I(m),J(m),real(Sij(m)),imag(Sij(m)));
1 P _6 Y K* j0 Y1 D! z3 A r& X) f fprintf(fid,' 节点%d到节点%d的功率为: %f+i*(%f) \n',J(m),I(m),real(Sji(m)),imag(Sji(m)));' J9 o: T* s& y) w6 G, Q
end
4 k: D$ k( i- D; {, g( G6 T8 k! ?3 Rfprintf(fid,' ================================================\n');
* s4 ^6 F0 H4 m8 k' i4 zfprintf(fid,' 网络总损耗sumdeltaS为: %f+i*(%f) \n',real(sumdeltaS),imag(sumdeltaS));
) Q1 ^( {' l/ h6 E% S( \fprintf(fid,' ================================================\n');! ?2 b! y, U" P) i# o+ j* Y/ k
fprintf(fid,' 线路功率损耗deltaSij为: \n');! i. v; R- r. G0 C; d8 Q4 W6 x" D. R `
for m=1:y
: B- O, T+ P6 V, [6 I fprintf(fid,' %d--%d线路的功率损耗为: %f+i*(%f) \n',I(m),J(m),real(deltaSij(m)),imag(deltaSij(m)));, z4 I' \) F' {9 u! N
end |
|