在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)
   # U) m# k- {; |& x/ X; i6 G! v
3. switch flag,
   
4. 
   
5.   case 0,
     w4 D  E' V. U4 ?9 |4 u/ K
6.     [sys,x0,str,ts]=mdlInitializeSizes;
   2 o( M/ |7 V6 J
7.   case {1,4,9},
   ' P( c7 v$ }/ g! [. w7 g. w
8.     sys=[];
   
9.   case 2,
   5 a0 u4 e; k# r- L/ m
10. %     sys=mdlUpdate(t,x,u);
    
11.     sys=mdlUpdate(t,x,u,yita);
    
12. 
    
13.   case 3,
    
14.     sys=mdlOutputs(t,x,u,K);
    
15.   otherwise
    
16.     error(['Unhandled flag = ',num2str(flag)]);
    
17. end
    
18. 
    " U+ v' {9 o& t+ |/ F  Y9 z# H2 U/ K& A
19. function [sys,x0,str,ts]=mdlInitializeSizes()
    8 s: \& T0 D" _1 R
20. 
    
21. sizes = simsizes;
    
22. sizes.NumContStates  = 0;
    
23. sizes.NumDiscStates  = 3;
    ( t0 J1 l  e" n" T1 \. J
24. sizes.NumOutputs     = 4;
    
25. sizes.NumInputs      = 4;
    
26. sizes.DirFeedthrough = 1;
    
27. sizes.NumSampleTimes = 1;  
      M7 d5 {4 A$ C" i9 f" d0 B
28. sys = simsizes(sizes);
    
29. 
    
30. x0  = [9;0.01;0];
    
31. 
    
32. str = [];
    - f# Z. B! l" ^, `% P: ^/ t
33. 
    8 g8 v+ p. G; _" E& u
34. ts  = [-1 0];
    2 C6 o3 r: q1 p1 I, u& l. v  u
35. 
    
36. function sys=mdlUpdate(t,x,u,yita)
    ; P% `' v4 r/ e7 c/ n
37. % function sys=mdlUpdate(t,x,u)
    
38. sys=x+yita*u(1)*u(4)*(2*u(1)-u(2));
    
39. 
    ) `8 r/ v8 L: q8 i
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))];
    

复制代码

redplum

fortran语言怎么和MATLAB结合呢

QTB

学习了，谢谢楼主