求教，关于PSAT数据如何导入和其他一些问题

qaws2002

问题1：我按照PSAT中文说明书编写的程序，但是不能得到相应的结果，在额外附加载荷50%之后，运行后相应的系统P,Q数值是增加50%，但是电压voltages值没有相应的改变，计算结果的voltages和没有加载的电压值相同。
5 C/ q8 l, n8 F% h3 M3 J   程序如下：
initpsat
8 W; Y4 C. n$ v2 O+ a, R* Hclpsatreadfile = 0;
runpsat('d_006_mdl','data')
runpsat('pf')
voltages = DAE.y(1+Bus.n:2*Bus.n);
: J, p, R9 r$ n3 c+ h7 z! p; Kfor i = 1:10
! |, q# ~9 g0 {3 Z- \  PQ.store(:,[4,5]) = (1+i/20)*[0.9, 0.6; 1, 0.7; 0.9, 0.6];
  PV.store(:,4) = (1+i/20)*[0.9; 0.6];
' y8 l2 M1 K: K9 x- q6 Y9 D( \1 J1 D  runpsat('pf')
  voltages = [voltages, DAE.y(1+Bus.n:2*Bus.n)];
end
1 H* U9 k4 j% g2 M1 l  Sclosepsat
9 [! z0 @7 R; T. {disp(voltages)

1 E  U" F: R* I3 c- F0 d! D其中在循环语句中的voltages语句，着实看不懂。
9 k# o- J- P" q. e8 B; {" y
' z) J. |4 b  A) g% y# S1 k   问题二：  如何将大量的数据导入到psat中呢？
7 a7 r" _& I7 _0 O" q
) _% n9 A! _9 S& l/ x  请教各位高手， 在线求答复。     

  2014年3月20日

qaws2002

晚上我闲着没事，重新仔细的检查了一下，
1 c1 B- h  ?) @: l+ B* C问题1我解决了，主要是在读入数据的时候  clpsatreadfile = 0;  在clpsat后面需要加入一个“.”  这样程序读入时候的含义是   不读入原始数据，是我粗心了。
" U/ A; d- Z0 [8 u: R但是第二个问题我稍微有点思路了，郁闷了好几天。
但还是希望有人能详细回复我的第二个问题。

goodluckyang

回复 2# qaws2002

" C7 m, _6 h9 R8 L
: W0 U6 ^. ^( e: H" A" v$ ^第二个问题 怎么破？