pscad 自定义模块

zlw

pscad自带的例子中，和MATLAB互联仿真画磁滞回线的那个例子。matlab的程序能否解读一下，想把它做成c语言编的程序。
function [out_IF] = hyst_plot(I_F)
: L( {; V' i/ I/ p' h%
+ h% {; ~; F' @$ I% This function is generated by MATLAB-PSCAD component generator V1.0
3 d* p1 v2 M% v%
% Input:
+ j# n4 h4 s" D%        I_F of dimension 2 1
%
9 z5 M8 }- g' c3 w4 F%
3 _; X! d8 V. o2 I) s+ U% Output:
%        out_IF of dimension 1 1

global Ia flux my_nom
9 n+ x2 y. \6 K6 R0 c
n2=15;                % points number
Flim=2.5;                % Flux limit
Ilim=0.1;                % Current limit
- d% G. E! i  e7 ~) B  u
%if ~exist('my_nom')
/ P% D! [* {+ z- A, c# J" S  my_nom=1;
' {( _  |& j. {8 m  Ia(n2)=0;
; A) L. i0 k" Y0 a' `+ r  flux(n2)=0;
( c+ U' E1 |' x& E%end
3 P6 k8 Y( q) y8 P/ }8 x, ffigure(1);
# G9 d7 i. `% Q& i9 nfor i=n2-1:-1:1
Ia(i+1)=Ia(i);
- K# M6 L' ^+ K flux(i+1)=flux(i);
" ~6 Z4 b5 T# k8 [1 Z, t" Iend
0 R1 T  `7 M# B# M6 ~0 j( A Ia(1)=I_F(1);
* x; _0 a" L' k/ J8 t/ [) ^ flux(1)=I_F(2);
9 L( h' w) u$ }5 T1 c( F$ [% line(Ia,flux);
( N7 i6 N3 M  x plot(Ia,flux);
set (gca,'XLim',[-Ilim Ilim],'YLim',[-Flim Flim]);
pause(realmin);

out_IF=0.0;