pscad 自定义模块

zlw

pscad自带的例子中，和matlab互联仿真画磁滞回线的那个例子。matlab的程序能否解读一下，想把它做成c语言编的程序。
function [out_IF] = hyst_plot(I_F)
$ |: K+ b& s" u& S%
* b. h2 g% D9 q' Z" R& S" H: g% This function is generated by MATLAB-PSCAD component generator V1.0
%
% Input:
" o! N  \; @6 B  \7 d# m%        I_F of dimension 2 1
) r6 Y7 G' X; Z. q. j%
%
& s5 P8 A: }2 b5 O* o; D9 J3 {% Output:
%        out_IF of dimension 1 1
/ N6 }/ O2 G- x" [5 Q, W, u9 ]
global Ia flux my_nom

1 o, H8 M5 y2 S. l4 m5 U  pn2=15;                % points number
Flim=2.5;                % Flux limit
Ilim=0.1;                % Current limit
) m# f  W2 C+ w. ?' B$ t' q' F/ l
0 H3 O# L- X9 O%if ~exist('my_nom')
  my_nom=1;
. o- D+ ]3 U: @  Ia(n2)=0;
0 k5 u9 `; R" W  m7 {  flux(n2)=0;
%end
figure(1);
- |0 D6 b! f: {! zfor i=n2-1:-1:1
Ia(i+1)=Ia(i);
' }! }3 z$ b( w) u  _" q flux(i+1)=flux(i);
end
  J% w0 ]: F8 i9 o& E( O Ia(1)=I_F(1);
& t/ v! Y6 I% C! f flux(1)=I_F(2);
3 c4 B! {' O% j# ]! V% line(Ia,flux);
# L3 @0 ]( a5 k  U, _4 k2 d plot(Ia,flux);
7 _  e: l+ W8 i) T( K% w set (gca,'XLim',[-Ilim Ilim],'YLim',[-Flim Flim]);
pause(realmin);

5 w) G; b$ }. m7 m" K" d, x" Q% Rout_IF=0.0;