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

gaoxunhui

支持一下！

momingrui

这是我写的潮流的一个函数
function jacobian=jacobian(V,Y1,dij,nodeNum,PVI,balanceNode)
, A+ S9 j% D/ |/ E

H=diag(V)*(Y1.*sin(dij))*diag(V)-diag(diag(V)*(Y1.*sin(dij))*(V));
- V2 l* }# }% V. [& O1 |) ~N=(-diag(V)*(Y1.*cos(dij)))+diag((Y1.*sin(dij))*(V));
( Q/ a( Y/ w' I' M# |$ pJ=diag(V)*(Y1.*cos(dij))+diag((Y1.*cos(dij))*(V))*diag(V);
3 m5 l( [- s  C# ?; d7 IL=diag(V)*(Y1.*sin(dij))+diag((Y1.*sin(dij))*(V))*diag(V);

5 F4 h" w. V! D7 X- P' S& U3 k5 K: CH(balanceNode,:)=0;
N(balanceNode,:)=0;
' d5 _( z2 a/ ^* J% C9 RJ(balanceNode,:)=0;
) U- N- @, X6 U' w% b* A* S3 TL(balanceNode,:)=0;
H(:,balanceNode)=0;
0 I7 r1 n2 O3 e+ h4 I  LN(:,balanceNode)=0;
J(:,balanceNode)=0;
, w; Q. A; K9 M6 S8 u" bL(:,balanceNode)=0;


0 `9 G, A( K( t+ S3 ]2 N6 jH(balanceNode,balanceNode)=100;
L(balanceNode,balanceNode)=100;
+ s' x* x  w- K( R
L(PVI,:)=0;
J(:,PVI)=0;
L(:,PVI)=0;
2 m8 t2 a# l' E6 ?N(PVI,:)=0;


PV_Diag=sparse(PVI,PVI,100,nodeNum,nodeNum);
L=L+PV_Diag;
- B; ^6 l9 h# t2 Z5 wjacobian=[H,J;N,L];

senny123

也是，先谢啦

blueray2008

THANK YOU

senny123

正在编写而之中

new2009

fortran很适合科学计算类的问题，特别是计算量比较大、比较复杂的问题
+ a0 m" _: E+ q& l/ x很多年以前学过点皮毛
+ B% t% T5 I' Qkingbat 发表于 2009-2-23 22:00

$ P* E4 @  e2 u) D                               
登录/注册后可看大图

9 f( D" R+ d8 Y3 v
PSCAD里面用的是fortran语言，我以前学C语言的，不知道里面能不能用C编程？

8willson

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

lhw606

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

cfs112xc

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

zlyang

我自己也写过，没成功