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

gaoxunhui

支持一下！

momingrui

这是我写的潮流的一个函数
function jacobian=jacobian(V,Y1,dij,nodeNum,PVI,balanceNode)
# i! q8 V# o3 z. F+ ^% B$ N

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));
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);

' D: O5 L; a. H8 NH(balanceNode,:)=0;
: {6 }4 O) X: A2 GN(balanceNode,:)=0;
8 Q: ^" e: o' z9 Q: n/ y0 @6 @J(balanceNode,:)=0;
: F  F  N4 P0 E7 oL(balanceNode,:)=0;
* E% t0 O+ T& Z3 R: n1 lH(:,balanceNode)=0;
$ B9 t1 O4 p! ?6 g, nN(:,balanceNode)=0;
J(:,balanceNode)=0;
" Q/ o: N2 H( a0 e) Z9 ^9 sL(:,balanceNode)=0;
% c0 C( ~! n1 |1 w1 M& e2 J/ \( S

& N3 R9 `, [* dH(balanceNode,balanceNode)=100;
L(balanceNode,balanceNode)=100;

0 W* [$ A. y1 L) R' }L(PVI,:)=0;
0 a- w4 N/ J. C/ c9 I+ i% GJ(:,PVI)=0;
; e0 F$ c; t) W2 |L(:,PVI)=0;
# C: T. A* w6 {( W1 pN(PVI,:)=0;
; x7 i" s" q. v/ t; m, N

PV_Diag=sparse(PVI,PVI,100,nodeNum,nodeNum);
L=L+PV_Diag;
, Q1 ?$ v  }" q/ Tjacobian=[H,J;N,L];

senny123

也是，先谢啦

blueray2008

THANK YOU

senny123

正在编写而之中

new2009

fortran很适合科学计算类的问题，特别是计算量比较大、比较复杂的问题
1 A# `3 r6 Y; d8 q! j很多年以前学过点皮毛
kingbat 发表于 2009-2-23 22:00

                               
登录/注册后可看大图

" w/ _# A! S0 `
8 b1 `2 M" |& k9 x, ~- z
+ g% N; e2 n" m# l' @( xPSCAD里面用的是fortran语言，我以前学C语言的，不知道里面能不能用C编程？

8willson

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

lhw606

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

cfs112xc

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

zlyang

我自己也写过，没成功