pscad 自定义模块

zlw

pscad自带的例子中，和MATLAB互联仿真画磁滞回线的那个例子。matlab的程序能否解读一下，想把它做成c语言编的程序。
function [out_IF] = hyst_plot(I_F)
%
+ I: I, R# d' j$ s4 J& F% This function is generated by MATLAB-PSCAD component generator V1.0
%
2 v# l3 m6 B( b& u% Input:
%        I_F of dimension 2 1
% E- r% p% N, S+ O( T%
%
) ]6 l# L0 D: p% Output:
%        out_IF of dimension 1 1
* d* S4 K9 W9 k2 S( a2 B4 I: s
; n! X& W% u6 A$ iglobal Ia flux my_nom

& d2 q. m: s2 b, an2=15;                % points number
Flim=2.5;                % Flux limit
Ilim=0.1;                % Current limit
8 X0 s: V, o" y9 K. x
# z% T/ N. S: m6 y2 ]5 N3 f+ z%if ~exist('my_nom')
  my_nom=1;
  Ia(n2)=0;
9 \* x9 j) b# ~4 ^  flux(n2)=0;
& u2 @- {# n& L%end
figure(1);
+ A* K- v7 R7 c! hfor i=n2-1:-1:1
5 _& j. Z) G, m# S9 y+ \5 ^: r Ia(i+1)=Ia(i);
) V! R. j" G/ b4 a flux(i+1)=flux(i);
( V( E5 w6 ^- @' Vend
8 `8 K5 Q' s' n) v0 Q Ia(1)=I_F(1);
flux(1)=I_F(2);
% line(Ia,flux);
5 e* }) {4 w3 P6 c: N- y: j2 e plot(Ia,flux);
2 }3 j( {$ a9 x& m8 `9 Q0 J set (gca,'XLim',[-Ilim Ilim],'YLim',[-Flim Flim]);
pause(realmin);

5 }! s. K  j+ a) ?* y& G! _out_IF=0.0;