pscad 自定义模块

zlw

pscad自带的例子中，和matlab互联仿真画磁滞回线的那个例子。matlab的程序能否解读一下，想把它做成c语言编的程序。
" ~8 q+ R' ]0 [. t/ ?. X. Q4 d; Z& F9 Kfunction [out_IF] = hyst_plot(I_F)
2 ?$ E- j# ~) j& H%
% This function is generated by MATLAB-PSCAD component generator V1.0
%
5 d% M" c2 z5 P' E, t- y% Input:
%        I_F of dimension 2 1
4 v8 v8 s8 H" M9 X%
%
6 I; U% r; A# b$ i, v% Output:
+ l; d( p2 s' K) j%        out_IF of dimension 1 1
  T, z0 c' P+ k% v$ a
) l$ l  ?% L6 n. H9 Cglobal Ia flux my_nom
. }& l$ A8 w$ z# Y, R2 H7 J- K
1 Y* L" z8 ]( ?- k+ o5 N0 f/ U/ An2=15;                % points number
Flim=2.5;                % Flux limit
$ U; W1 V9 d) p. `7 m5 VIlim=0.1;                % Current limit

%if ~exist('my_nom')
6 f, s0 \' G% Q7 _  my_nom=1;
8 b2 [8 Y0 g# F* q) q/ b6 @: c  Ia(n2)=0;
1 D1 h# Y6 e, i/ o! V7 {  M% |. Y2 k  flux(n2)=0;
%end
" ~/ v2 E+ ], i$ \figure(1);
  Y" y% M/ C1 o7 D& N5 x# hfor i=n2-1:-1:1
9 d& J" [, E3 v4 [8 V. M& Z' s Ia(i+1)=Ia(i);
" t5 b4 h1 D; j flux(i+1)=flux(i);
1 a" P9 V/ o  Zend
8 \1 e  F7 M. A- W Ia(1)=I_F(1);
/ J& V' T! m" a- D9 L flux(1)=I_F(2);
1 F4 {3 K$ q- t* q, ?' A% line(Ia,flux);
1 J% g/ X* l8 N$ N7 y# H plot(Ia,flux);
set (gca,'XLim',[-Ilim Ilim],'YLim',[-Flim Flim]);
# l/ g( J" x& c, |# d+ o pause(realmin);
" T) Q( H0 z$ F! D  v2 n: a, `
out_IF=0.0;