在PSCAD中用fortran编程怎样实现单神经元PID控制器

我家雪

大家好，最近一直在忙仿真，用到单神经元PID控制，在MATLAB中可以用S函数搭建其模型（见附件），但在pscad中用fortran语言不太会，请求高手指点，很急！各位朋友帮帮忙！下面是在MATLAB中的单神经元PID的S函数编写的文件，结合。

1. function [sys,x0,str,ts] =Sfunction(t,x,u,flag,yita,K)
   
2. % function [sys,x0,str,ts] =Sfunction(t,x,u,flag,K)
   $ J/ N/ |; T2 u  K. |
3. switch flag,
   
4. 
   
5.   case 0,
   
6.     [sys,x0,str,ts]=mdlInitializeSizes;
   
7.   case {1,4,9},
   
8.     sys=[];
   5 p7 o+ M+ T6 S# s
9.   case 2,
   
10. %     sys=mdlUpdate(t,x,u);
    / o1 @5 ^2 r: B' _8 R8 y
11.     sys=mdlUpdate(t,x,u,yita);
    
12. 
    $ z/ X  [* w! q- S: a
13.   case 3,
    ) v% D) |' g- o, n: o" N
14.     sys=mdlOutputs(t,x,u,K);
    # e! F) R* r! d# c* G' _; l. P
15.   otherwise
    
16.     error(['Unhandled flag = ',num2str(flag)]);
    1 t' H# `  \) ~
17. end
    ! u7 A/ A* s7 ~7 f
18. 
    6 K2 n/ B- d+ t( N2 s, T
19. function [sys,x0,str,ts]=mdlInitializeSizes()
    ! q4 Q( a5 p+ p2 J
20. 
    4 r% l0 `. b% k" ^
21. sizes = simsizes;
    
22. sizes.NumContStates  = 0;
    + @3 P1 V: N+ {' I* `1 K/ L+ {
23. sizes.NumDiscStates  = 3;
    & w6 H  r0 c$ e! s; w# T7 X+ w
24. sizes.NumOutputs     = 4;
    
25. sizes.NumInputs      = 4;
    6 e# O2 F4 ^6 }3 I- m/ ?$ j
26. sizes.DirFeedthrough = 1;
    
27. sizes.NumSampleTimes = 1;  
    
28. sys = simsizes(sizes);
    
29. 
    
30. x0  = [9;0.01;0];
    
31. 
    
32. str = [];
    
33. 
    
34. ts  = [-1 0];
    
35. 
    , V9 _5 C5 o. b0 X
36. function sys=mdlUpdate(t,x,u,yita)
    " |3 b/ Q2 o  F$ r
37. % function sys=mdlUpdate(t,x,u)
    ( W' {( r5 A+ H* O6 j( U
38. sys=x+yita*u(1)*u(4)*(2*u(1)-u(2));
    5 y5 A' |1 O" T. X& P& t2 l) [4 U
39. 
    % V6 @' H" x) V$ _4 g
40. function sys=mdlOutputs(t,x,u,K)   
    : ]7 ?$ D! T+ V/ Q' e. T; o
41. xx=[u(1)-u(2) u(1) u(1)+u(3)-2*u(2)];   
    ( G4 x' S9 \' P* x1 ^+ A
42. sys =[u(4)+K*xx*x/sum(abs(x));x/sum(abs(x))];
    

复制代码

redplum

fortran语言怎么和MATLAB结合呢

QTB

学习了，谢谢楼主