在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,
   : a* c( x3 F- f3 t4 K
4. 
   
5.   case 0,
   
6.     [sys,x0,str,ts]=mdlInitializeSizes;
   
7.   case {1,4,9},
   
8.     sys=[];
   2 S# A& y8 G, J( p4 j
9.   case 2,
   
10. %     sys=mdlUpdate(t,x,u);
    + g$ F$ r& g; S- U) t3 U/ A0 G
11.     sys=mdlUpdate(t,x,u,yita);
    . H) Y) m; W1 v2 q+ m
12. 
    
13.   case 3,
    
14.     sys=mdlOutputs(t,x,u,K);
    + d- m  ~* h( j" ?2 p8 ^* e, a; r  q
15.   otherwise
    
16.     error(['Unhandled flag = ',num2str(flag)]);
    
17. end
    
18. 
    * ~% o  g' \, p/ Q8 [9 t
19. function [sys,x0,str,ts]=mdlInitializeSizes()
    8 K- _! X2 R+ p
20. 
    
21. sizes = simsizes;
    7 O+ c! R0 a. I; a5 d
22. sizes.NumContStates  = 0;
    " e  o; ^: h6 ]! _
23. sizes.NumDiscStates  = 3;
    
24. sizes.NumOutputs     = 4;
    * t. w7 N# c4 D+ |) u
25. sizes.NumInputs      = 4;
    - T' l  ^! D. h' E9 l
26. sizes.DirFeedthrough = 1;
    
27. sizes.NumSampleTimes = 1;  
    
28. sys = simsizes(sizes);
    : b( ]/ J0 ]  S7 J1 L
29. 
    ; p* p& p; v! T0 x
30. x0  = [9;0.01;0];
    # t. w9 {8 T( A. C& t# ?1 ^
31. 
    ; C& F4 V# R3 F" {& N
32. str = [];
    $ e4 o5 W% U( J( m. |. m2 _9 f
33. 
    
34. ts  = [-1 0];
    * f# X' C$ {5 y+ n4 C7 v  y
35. 
    . _' R# A* I* _: [0 V
36. function sys=mdlUpdate(t,x,u,yita)
    " ?; _+ v, [8 ~4 K2 W
37. % function sys=mdlUpdate(t,x,u)
    
38. sys=x+yita*u(1)*u(4)*(2*u(1)-u(2));
    
39. 
    9 S4 u, X6 K2 ^# Q: N
40. function sys=mdlOutputs(t,x,u,K)   
    
41. xx=[u(1)-u(2) u(1) u(1)+u(3)-2*u(2)];   
    " l* `8 Z4 Q. @+ ^& U" r  S2 E
42. sys =[u(4)+K*xx*x/sum(abs(x));x/sum(abs(x))];
    # Y7 Q2 N" \0 O2 W

复制代码

redplum

fortran语言怎么和MATLAB结合呢

QTB

学习了，谢谢楼主