pscad 自定义模块

zlw

pscad自带的例子中，和MATLAB互联仿真画磁滞回线的那个例子。matlab的程序能否解读一下，想把它做成c语言编的程序。
; ]. y- ^# ~4 {6 A  ?& @' xfunction [out_IF] = hyst_plot(I_F)
%
3 R* j7 |/ W' a! u7 C' I1 `( m% This function is generated by MATLAB-PSCAD component generator V1.0
; Y9 [5 y! G$ c4 c1 P%
3 {* q& B, G+ w, ^/ G; ~: {/ h% Input:
. N. W, J4 I3 r1 n/ I( f4 J8 s%        I_F of dimension 2 1
%
. @, h, f. j3 d9 c: d* @; Z% D%
/ X; I' x+ G  i) E% Output:
%        out_IF of dimension 1 1
5 {* F- }5 J" ~3 p( A& v% t+ g
global Ia flux my_nom

+ ]# L7 I0 |# d' N6 V. N: X: m9 ]9 `! dn2=15;                % points number
/ n& C9 M1 Z) b5 U0 U, ZFlim=2.5;                % Flux limit
) V! E$ Q" R2 }2 Y; Q! SIlim=0.1;                % Current limit
/ v% J6 n$ {* v
%if ~exist('my_nom')
4 Q+ y( Z# y  [  my_nom=1;
  Ia(n2)=0;
  a$ j% Q3 {5 n" V  flux(n2)=0;
%end
& J1 B- Z, i. @2 \figure(1);
for i=n2-1:-1:1
  Q# h+ t& M3 E0 `. o2 B0 [ Ia(i+1)=Ia(i);
flux(i+1)=flux(i);
& Q, _$ c& q% @8 N9 qend
, h6 J9 H* s5 N5 M6 h Ia(1)=I_F(1);
+ Z8 F) [) I4 j4 |9 c  z flux(1)=I_F(2);
% line(Ia,flux);
plot(Ia,flux);
  _+ c4 W1 b! q set (gca,'XLim',[-Ilim Ilim],'YLim',[-Flim Flim]);
" N7 I# }1 x0 P' Y pause(realmin);
8 _, ^' j3 S+ J  }
out_IF=0.0;