在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)
   ; T8 e8 x! h3 C2 A; ?: O9 A5 l
3. switch flag,
   ( N8 F2 X1 I+ N; e, u0 E
4. 
   
5.   case 0,
   
6.     [sys,x0,str,ts]=mdlInitializeSizes;
   
7.   case {1,4,9},
   2 q6 K- `" |# ~- b. m$ L" _
8.     sys=[];
   
9.   case 2,
   
10. %     sys=mdlUpdate(t,x,u);
    
11.     sys=mdlUpdate(t,x,u,yita);
    
12. 
    ! n1 t/ d6 ?) S: b) N) p
13.   case 3,
    
14.     sys=mdlOutputs(t,x,u,K);
    & O' P/ H( \3 `- f& X" l. V
15.   otherwise
    
16.     error(['Unhandled flag = ',num2str(flag)]);
    ; |4 z5 P; L) r/ q# k3 E
17. end
    6 _5 ?2 J3 R8 @3 Q# K
18. 
    9 ]4 u. ~2 C' e9 O% A% \; g3 R
19. function [sys,x0,str,ts]=mdlInitializeSizes()
    
20. 
    
21. sizes = simsizes;
    8 V  a& l% k' \/ B0 n& _
22. sizes.NumContStates  = 0;
    
23. sizes.NumDiscStates  = 3;
    . ~5 J0 F- a& G# m' f: @
24. sizes.NumOutputs     = 4;
    # p8 j: ^% b5 l: ]
25. sizes.NumInputs      = 4;
    
26. sizes.DirFeedthrough = 1;
    8 H0 Y) E2 T. a! a2 a1 ^# K! x
27. sizes.NumSampleTimes = 1;  
    
28. sys = simsizes(sizes);
    ; K: K* x1 l( h  P/ P
29. 
    $ z  Z5 h, Y' s- N) Z9 l
30. x0  = [9;0.01;0];
    - x1 q! ^$ \' Z' Q+ }; _) y: ?# ]
31. 
    
32. str = [];
    - a7 F1 H, V0 f) G* z! E
33. 
    1 @# F- ^9 ^2 v, s% E8 {3 N
34. ts  = [-1 0];
    3 N% |6 E  V8 P2 d& m7 e% \. {. y  h
35. 
    9 ]$ t' i7 z! ^, @
36. function sys=mdlUpdate(t,x,u,yita)
    ; }  \0 g0 ~3 F! B
37. % function sys=mdlUpdate(t,x,u)
    
38. sys=x+yita*u(1)*u(4)*(2*u(1)-u(2));
    
39. 
    
40. function sys=mdlOutputs(t,x,u,K)   
    * B0 D1 W  Q# P
41. xx=[u(1)-u(2) u(1) u(1)+u(3)-2*u(2)];   
    5 P/ w/ E$ I/ u; i2 r
42. sys =[u(4)+K*xx*x/sum(abs(x));x/sum(abs(x))];
    

复制代码

redplum

fortran语言怎么和MATLAB结合呢

QTB

学习了，谢谢楼主