在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)
   
3. switch flag,
   8 j; F5 R2 h/ [9 q9 R% b/ w
4. 
   ( i) v# A* J! q" l
5.   case 0,
   ; ?+ V' k% w- ]/ p
6.     [sys,x0,str,ts]=mdlInitializeSizes;
   ) f, F$ C/ R+ j1 ^
7.   case {1,4,9},
   
8.     sys=[];
   
9.   case 2,
   
10. %     sys=mdlUpdate(t,x,u);
    9 F) w& J6 f0 V
11.     sys=mdlUpdate(t,x,u,yita);
    : z' E* E8 l, R" w( [1 |- L# Q
12. 
    2 ^" p/ C  L5 F
13.   case 3,
    
14.     sys=mdlOutputs(t,x,u,K);
    & W, h. r+ }. V' K( V
15.   otherwise
    * ~6 X! ]. s" U' N0 x9 G
16.     error(['Unhandled flag = ',num2str(flag)]);
    - s& |8 f. R# Q- ]/ i
17. end
    ' E% j' y2 N' I& r3 J
18. 
    
19. function [sys,x0,str,ts]=mdlInitializeSizes()
    ( D& s5 h- M( `5 e# e, j2 B
20. 
    
21. sizes = simsizes;
    
22. sizes.NumContStates  = 0;
    " E! j7 x  n+ s) P3 E' |9 @2 N6 a! g4 ]
23. sizes.NumDiscStates  = 3;
    
24. sizes.NumOutputs     = 4;
    " Z  Q0 a+ @1 j4 f
25. sizes.NumInputs      = 4;
    
26. sizes.DirFeedthrough = 1;
    . |5 o/ e( d8 K! ~
27. sizes.NumSampleTimes = 1;  
    
28. sys = simsizes(sizes);
    
29. 
    
30. x0  = [9;0.01;0];
    
31. 
    1 b# X# p! h9 X3 Q$ V
32. str = [];
    
33. 
    7 F( }3 o, v5 E4 h7 O/ V) L
34. ts  = [-1 0];
    
35. 
    9 o8 P4 u* K) q8 {
36. function sys=mdlUpdate(t,x,u,yita)
    * N7 i; i3 o$ O6 ]# I5 x" d& B0 }/ C
37. % function sys=mdlUpdate(t,x,u)
    9 G" N6 [* ~3 Z$ a9 d3 v
38. sys=x+yita*u(1)*u(4)*(2*u(1)-u(2));
    
39. 
    
40. function sys=mdlOutputs(t,x,u,K)   
    
41. xx=[u(1)-u(2) u(1) u(1)+u(3)-2*u(2)];   
    ) F& r1 e/ ~* o% B  J* t& ~# a
42. sys =[u(4)+K*xx*x/sum(abs(x));x/sum(abs(x))];
    3 j! J: b4 d6 U/ v$ f

复制代码

redplum

fortran语言怎么和MATLAB结合呢

QTB

学习了，谢谢楼主