PSCAD 自定义没有报错 但是输出不正常

vince9501

#include"math.h"
typedef double real;
6 x/ w/ v( O5 Q% ]7 s1 C0 f- U( rvoid mp_ptd__(real* arg1,real* arg2,real* arg3,real* res)
2 `& S3 {/ Q! S1 B/ e) z& x5 f: u{

1 I* @9 p! P# `5 R1 v. d' Ereal V0=(*arg3);
real c1=0.000003;
real c2=0.07872;
real Voc=21.1;
* ~9 m1 m3 n8 E: ~6 Q4 qreal Isc=3.864;
7 _( U6 I$ w# K4 R0 d: |; rreal DT=(*arg1)-25;
" B  |: }+ s$ z. s8 w' H  _real DI=0.00065*((*arg2)/1000)*DT+((*arg2)/1000-1)*3.846;
( r! _0 r- T, }" u% M) Mreal DV=-0.08*DT-(-3.4872)*DI;
real esp=0.001;
real V[99]={0};
% d( P2 r& Y8 D  Lreal cn=1;
/ u- \* b6 k$ b$ {) T* V3 t# gint i=0;
5 @7 T: N1 f. ]while(fabs(cn)>esp*V)
{
i++;
V=V[i-1]-(Isc*(1+c1-c1*exp((V[i-1]-DV)/(c2*Voc))-c1*V[i-1]*exp((V [i-1]-DV)/(c2*Voc)/(c2*Voc))))/(-(2*c1*exp((V[i-1]-DV)/(c2*Voc))/ (c2*Voc)+V[i-1]*c1*exp((V[i-1]-DV)/(c2*Voc))/((c2*Voc)*(c2*Voc))));
cn=V-V[i-1];
(*res)=V;
& j! j, `  s% C% G }
' h! P5 x" a* s! o( Q: m}
$ N8 P# y- U4 |( L* [( r调试没有报错，但是得不到输出，请指教。

sunek

这是要做什么？保护的判据么？

vince9501

回复 2# sunek
9 s& J+ D4 p/ d' ?0 ~( R
/ P" F/ K. s0 T/ I; z4 Y0 z$ _
   信号跟踪 cn和esp 作为精度。。。