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

gaoxunhui · 发表于 2009-3-10 11:53:31

支持一下！

momingrui · 发表于 2009-3-10 12:08:26

这是我写的潮流的一个函数
function jacobian=jacobian(V,Y1,dij,nodeNum,PVI,balanceNode)

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

H(balanceNode,:)=0;
N(balanceNode,:)=0;
J(balanceNode,:)=0;
L(balanceNode,:)=0;
H(:,balanceNode)=0;
N(:,balanceNode)=0;
J(:,balanceNode)=0;
L(:,balanceNode)=0;

H(balanceNode,balanceNode)=100;
L(balanceNode,balanceNode)=100;

L(PVI,:)=0;
J(:,PVI)=0;
L(:,PVI)=0;
N(PVI,:)=0;

PV_Diag=sparse(PVI,PVI,100,nodeNum,nodeNum);
L=L+PV_Diag;
jacobian=[H,J;N,L];

senny123 · 发表于 2009-3-10 13:47:20

也是，先谢啦

blueray2008 · 发表于 2009-3-15 15:54:26

THANK YOU

senny123 · 发表于 2009-3-18 09:30:55

正在编写而之中

new2009 · 发表于 2009-3-19 13:06:26

fortran很适合科学计算类的问题，特别是计算量比较大、比较复杂的问题
' V; [: T% l) S5 @很多年以前学过点皮毛
U; S" T, Y( r1 N3 p6 V/ [kingbat 发表于 2009-2-23 22:00

; i% s+ T+ x% c6 H. n \( ^
登录/注册后可看大图

PSCAD里面用的是fortran语言，我以前学C语言的，不知道里面能不能用C编程？

8willson · 发表于 2009-3-21 11:09:42

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

lhw606 · 发表于 2009-4-3 00:17:51

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

cfs112xc · 发表于 2009-4-15 19:59:24

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

zlyang · 发表于 2009-4-20 14:14:21

我自己也写过，没成功

账号		自动登录	找回密码
密码			立即加入

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