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

gaoxunhui

支持一下！

momingrui

这是我写的潮流的一个函数
7 H9 K$ ]$ N" ffunction jacobian=jacobian(V,Y1,dij,nodeNum,PVI,balanceNode)
' y- @# z6 W" j% S" I* i
2 q/ S1 r" k* r" M
% [) x9 B  }8 J% \: VH=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));
J=diag(V)*(Y1.*cos(dij))+diag((Y1.*cos(dij))*(V))*diag(V);
' R. ^2 e2 F/ i- `L=diag(V)*(Y1.*sin(dij))+diag((Y1.*sin(dij))*(V))*diag(V);
3 K: R; O- h5 `
H(balanceNode,:)=0;
4 u0 b) a* [2 ^* Z; h4 C$ o" w* @N(balanceNode,:)=0;
J(balanceNode,:)=0;
L(balanceNode,:)=0;
H(:,balanceNode)=0;
( }: N) `" e( k# ]7 p) ^$ X( U2 Y3 ^N(:,balanceNode)=0;
J(:,balanceNode)=0;
L(:,balanceNode)=0;
$ J5 H$ u: G* N  \6 ?+ g# E* C$ y0 m
3 A( H2 B* g/ p) \6 v  e" S
H(balanceNode,balanceNode)=100;
/ H- r# `7 F4 H0 {! H9 ~: B5 ML(balanceNode,balanceNode)=100;

L(PVI,:)=0;
! A5 r, |3 x+ j& o8 dJ(:,PVI)=0;
9 L# P; l7 H( H) H/ a8 YL(:,PVI)=0;
: P3 d" z+ M' i# \/ x( s% M$ J$ YN(PVI,:)=0;

8 ^5 S1 a0 i1 ?) j9 u
PV_Diag=sparse(PVI,PVI,100,nodeNum,nodeNum);
L=L+PV_Diag;
9 p. @9 c% d/ S% u/ e" |  Z- \jacobian=[H,J;N,L];

senny123

也是，先谢啦

blueray2008

THANK YOU

senny123

正在编写而之中

new2009

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

                               
登录/注册后可看大图

. o- U4 S# K3 {6 ?; F$ m

$ A; ]: y. [5 O5 S  \9 yPSCAD里面用的是fortran语言，我以前学C语言的，不知道里面能不能用C编程？

8willson

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

lhw606

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

cfs112xc

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

zlyang

我自己也写过，没成功