自己写的matlab潮流分析程序，高斯，牛拉，PQ法

gaoxunhui

支持一下！

momingrui

这是我写的潮流的一个函数
! N6 b) s2 p5 Y% s, tfunction jacobian=jacobian(V,Y1,dij,nodeNum,PVI,balanceNode)

" \1 k; j7 k% B8 H3 ~5 [$ c* G* c% O
H=diag(V)*(Y1.*sin(dij))*diag(V)-diag(diag(V)*(Y1.*sin(dij))*(V));
+ I+ g' m4 k- _: ?& E& ^( o& MN=(-diag(V)*(Y1.*cos(dij)))+diag((Y1.*sin(dij))*(V));
J=diag(V)*(Y1.*cos(dij))+diag((Y1.*cos(dij))*(V))*diag(V);
L=diag(V)*(Y1.*sin(dij))+diag((Y1.*sin(dij))*(V))*diag(V);
5 P, A3 _1 e4 a1 ?, G  w8 P9 ~6 r
5 K6 c6 D8 l" |0 j" ^; y5 @) J. tH(balanceNode,:)=0;
) q1 ^) M1 v# o; q/ [  |N(balanceNode,:)=0;
J(balanceNode,:)=0;
L(balanceNode,:)=0;
H(:,balanceNode)=0;
- p3 l! d) l0 k" iN(:,balanceNode)=0;
J(:,balanceNode)=0;
L(:,balanceNode)=0;
7 |. c' d: X: i/ z0 i* r& Y
- y/ P! H  f4 ]8 i) ?
7 y* h% z$ g& Y! ?$ F: g( k7 A8 rH(balanceNode,balanceNode)=100;
* d+ }! [: a( ~L(balanceNode,balanceNode)=100;
3 h4 S8 t: G6 G. O9 z! s3 M+ b
L(PVI,:)=0;
J(:,PVI)=0;
L(:,PVI)=0;
N(PVI,:)=0;
9 G! ]2 G# k  L( @
& a8 S% j6 n. D8 a
PV_Diag=sparse(PVI,PVI,100,nodeNum,nodeNum);
L=L+PV_Diag;
jacobian=[H,J;N,L];

senny123

也是，先谢啦

blueray2008

THANK YOU

senny123

正在编写而之中

new2009

fortran很适合科学计算类的问题，特别是计算量比较大、比较复杂的问题
很多年以前学过点皮毛
5 U( O7 e# Y" V$ s2 F* Skingbat 发表于 2009-2-23 22:00

  Y, Q/ v; ^, q, q9 o                               
登录/注册后可看大图

* L7 M7 P$ C* l
$ n4 M& e. h! X" k* k: n( d( A
( @3 F7 W+ v5 o6 Y$ A# kPSCAD里面用的是fortran语言，我以前学C语言的，不知道里面能不能用C编程？

8willson

怎么会是笑话呢~~~
好东西什么时候都受欢迎的~~

lhw606

支持原创~~~~~~~~~~~~

cfs112xc

我的论文是用fortran编的潮流计算 花了一两个月 晕 用matlab还是比较快

zlyang

我自己也写过，没成功