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

gaoxunhui

支持一下！

momingrui

这是我写的潮流的一个函数
function jacobian=jacobian(V,Y1,dij,nodeNum,PVI,balanceNode)
( c" ?9 R0 n' X" \2 H

H=diag(V)*(Y1.*sin(dij))*diag(V)-diag(diag(V)*(Y1.*sin(dij))*(V));
N=(-diag(V)*(Y1.*cos(dij)))+diag((Y1.*sin(dij))*(V));
8 T- r6 c( q4 S' ]J=diag(V)*(Y1.*cos(dij))+diag((Y1.*cos(dij))*(V))*diag(V);
, x( R. p. H3 O2 UL=diag(V)*(Y1.*sin(dij))+diag((Y1.*sin(dij))*(V))*diag(V);
/ y  z3 S! C" @7 R% ^
4 ]' Q, E/ Q/ d7 i" }H(balanceNode,:)=0;
" n) y4 ^* f) |' W* wN(balanceNode,:)=0;
J(balanceNode,:)=0;
# Z9 J7 y! p% }2 k  ?' IL(balanceNode,:)=0;
H(:,balanceNode)=0;
N(:,balanceNode)=0;
J(:,balanceNode)=0;
L(:,balanceNode)=0;
! |' U5 W3 T( [7 {! ?8 Z" P
- O, c" a3 {. L7 i
2 |$ y  ?" c! k! DH(balanceNode,balanceNode)=100;
/ a2 u& u, o/ A! ?6 }+ j" }) kL(balanceNode,balanceNode)=100;
  ]+ K1 y. \% r
L(PVI,:)=0;
! D" ^+ o4 q/ KJ(:,PVI)=0;
* @, d! H" }6 v2 B, ^3 @# k% UL(:,PVI)=0;
( ^' M# m' w0 @$ v* g6 _; ^  H" GN(PVI,:)=0;


PV_Diag=sparse(PVI,PVI,100,nodeNum,nodeNum);
, A9 d' W, p3 B0 k% Q7 B! c2 Y) ?L=L+PV_Diag;
' R  S) v- @* R9 h: ~3 Xjacobian=[H,J;N,L];

senny123

也是，先谢啦

blueray2008

THANK YOU

senny123

正在编写而之中

new2009

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

9 @- D  h# `# f4 G' H' i5 p7 H/ O                               
登录/注册后可看大图

9 a" q) f/ f, C) U9 A6 @

$ e' l2 l/ M' sPSCAD里面用的是fortran语言，我以前学C语言的，不知道里面能不能用C编程？

8willson

怎么会是笑话呢~~~
4 Y  `$ S* z2 S' a* x7 s$ q好东西什么时候都受欢迎的~~

lhw606

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

cfs112xc

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

zlyang

我自己也写过，没成功