在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)
   $ g! I( ~2 c) J
3. switch flag,
   : V* D$ Y: D4 L3 O) }
4. 
   ) U( h2 h5 k) w
5.   case 0,
   
6.     [sys,x0,str,ts]=mdlInitializeSizes;
   
7.   case {1,4,9},
   3 p" v- m! A8 V8 x3 r/ ^; O' h" ]+ k
8.     sys=[];
   
9.   case 2,
   - T# e# J+ I+ k( I. T  A! `
10. %     sys=mdlUpdate(t,x,u);
    0 ^$ g6 R* v1 X8 `$ x" c' m
11.     sys=mdlUpdate(t,x,u,yita);
      e. i+ M6 e& V0 [- E/ F. K, \% w$ ]; F
12. 
    2 z1 j6 x) M* u1 U" S/ v) `$ h) i
13.   case 3,
    
14.     sys=mdlOutputs(t,x,u,K);
    + l6 ^# U8 l% \% I& t
15.   otherwise
    6 j! j( y+ b# j2 |" X
16.     error(['Unhandled flag = ',num2str(flag)]);
    
17. end
    3 r$ {, L5 U* o: p; L3 Q) W( W
18. 
    
19. function [sys,x0,str,ts]=mdlInitializeSizes()
    
20. 
    % B" c! \( d( S/ N  B" L
21. sizes = simsizes;
    
22. sizes.NumContStates  = 0;
    
23. sizes.NumDiscStates  = 3;
    
24. sizes.NumOutputs     = 4;
    8 c+ B" D* {+ l% _% ]
25. sizes.NumInputs      = 4;
    6 l2 B( _( x8 v" ^4 m  ]7 r
26. sizes.DirFeedthrough = 1;
    
27. sizes.NumSampleTimes = 1;  
    
28. sys = simsizes(sizes);
    : q: g. B: x. e9 A' a9 I7 `  e
29. 
    5 L1 T) X1 K% G7 g! T/ J
30. x0  = [9;0.01;0];
    
31. 
    2 l! V1 w+ A9 A0 T- A" j
32. str = [];
    # z0 {5 W- Y% K2 F
33. 
    9 w7 u* X- `  S; g
34. ts  = [-1 0];
    
35. 
    6 D6 ]2 H! Q0 Y& D4 S9 G
36. function sys=mdlUpdate(t,x,u,yita)
    4 N' ]; A/ s: ?. e. `
37. % function sys=mdlUpdate(t,x,u)
    
38. sys=x+yita*u(1)*u(4)*(2*u(1)-u(2));
    
39. 
    1 k9 x. G, `" N, b, J* S
40. function sys=mdlOutputs(t,x,u,K)   
    
41. xx=[u(1)-u(2) u(1) u(1)+u(3)-2*u(2)];   
    
42. sys =[u(4)+K*xx*x/sum(abs(x));x/sum(abs(x))];
    6 |7 B4 h* G: f; ~% E+ K  i

复制代码

redplum

fortran语言怎么和MATLAB结合呢

QTB

学习了，谢谢楼主