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

gaoxunhui

支持一下！

momingrui

这是我写的潮流的一个函数
function jacobian=jacobian(V,Y1,dij,nodeNum,PVI,balanceNode)
8 y$ p$ V( m( t" y
8 D9 u" Z- u3 T! y" m
& u6 ^2 R6 G. a" F8 m: PH=diag(V)*(Y1.*sin(dij))*diag(V)-diag(diag(V)*(Y1.*sin(dij))*(V));
) a9 Z9 Z/ P/ {( S( Z; G: kN=(-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);
0 T% ^" Q3 _! O. A+ w6 L  v3 F  @
H(balanceNode,:)=0;
N(balanceNode,:)=0;
J(balanceNode,:)=0;
L(balanceNode,:)=0;
H(:,balanceNode)=0;
N(:,balanceNode)=0;
  P# k* [% R2 gJ(:,balanceNode)=0;
6 m" m& k* |7 b4 ?3 b1 HL(:,balanceNode)=0;

- b1 e" G3 S% y' S" m) L
8 B) X- _. U5 c% H+ tH(balanceNode,balanceNode)=100;
L(balanceNode,balanceNode)=100;

+ Y* y- N+ k& H- P0 t' M4 oL(PVI,:)=0;
J(:,PVI)=0;
' Z6 i" o1 W: i2 `L(:,PVI)=0;
( p: s8 a. P+ C2 ?  p% X' \N(PVI,:)=0;


PV_Diag=sparse(PVI,PVI,100,nodeNum,nodeNum);
" }  M6 `0 h4 Q' rL=L+PV_Diag;
3 Z- S6 [: L' x! T. O- Z) r; bjacobian=[H,J;N,L];

senny123

也是，先谢啦

blueray2008

THANK YOU

senny123

正在编写而之中

new2009

fortran很适合科学计算类的问题，特别是计算量比较大、比较复杂的问题
很多年以前学过点皮毛
kingbat 发表于 2009-2-23 22:00

# ~! Q4 T1 @6 S  _" g. P5 F3 ?                               
登录/注册后可看大图

. i: b; A2 U5 x( @: q* d' n

* E8 w0 O& N" @- ], S/ l4 uPSCAD里面用的是fortran语言，我以前学C语言的，不知道里面能不能用C编程？

8willson

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

lhw606

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

cfs112xc

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

zlyang

我自己也写过，没成功