状态方程中的特性矩阵A～D随时间变化，如何定义？

skyismine

请教各位个问题，
关于状态空间法中的四个特性矩阵随时间变化时，A～D在s-function中如何定义呢？
问题是这样的，状态空间法中求解时A～D一般是常数矩阵，遇到耦合的问题，导致A~D是时变的，即为时间t的函数，这时候可以用其他的办法求解出来，但是，我就想，能否在s-function中直接定义A～D四个时变矩阵。我写的s-function代码贴出来是这样的：

1. function [sys,x0,str,ts] = stest(t,x,u,flag)
   
2. switch flag,
   
3. case 0,
   
4. [sys,x0,str,ts]=mdlInitializeSizes();
   
5. case 1,
   
6. sys=mdlDerivatives(t,x,u);
   
7. case 3,
   
8. sys=mdlOutputs(t,x,u);
   
9. case {2,4,9},
   
10. sys=[];
    
11. otherwise
    
12. error(['Unhandled flag = ',num2str(flag)]);
    
13. end
    
14. 
    
15. function [sys,x0,str,ts]=mdlInitializeSizes()
    
16. A=zeros(4,4);
    
17. D=zeros(2,2);
    
18. sizes = simsizes;
    
19. sizes.NumContStates = size(A,1);
    
20. sizes.NumDiscStates = 0;
    
21. sizes.NumOutputs = size(A,1)+size(D,1);
    
22. sizes.NumInputs = size(D,2);
    
23. sizes.DirFeedthrough = 1;
    
24. sizes.NumSampleTimes = 1; % at least one sample time is needed
    
25. sys = simsizes(sizes);
    
26. x0 = zeros(size(A,1),1);
    
27. str = [];
    
28. ts = [-1 0];
    
29. 
    
30. function sys=mdlDerivatives(t,x,u)
    
31. A=[1,3,t,t^2;2,3,5,1;t-1,3,1,5;3,t,1,3];
    
32. B=[2,1;3,t;3,;5,t];
    
33. sys = A*x+B*u;
    
34. 
    
35. function sys=mdlOutputs(t,x,u)
    
36. C=[1,3,t,t^2;t-1,3,1,5];
    
37. D=[1,t^2;2,t];
    
38. sys = [C*x+D*u;x];

复制代码

运行就会提示undefined function or variable
麻烦高手帮我解答一下，十分感谢。