如何在m文件中实现simulink中的du/dt实时求导功能？

zearo

不想用simlink环境仿真。在设计控制器的中间步骤中用到中间量的导数，在simulink中可以用du/dt模块轻松实现，但到m文件中如何实现呢？请指教！谢谢先。

eight

原帖由 zearo 于 2007-6-8 17:55 发表

                               
登录/注册后可看大图

不想用simlink环境仿真。在设计控制器的中间步骤中用到中间量的导数，在simulink中可以用du/dt模块轻松实现，但到m文件中如何实现呢？请指教！谢谢先。

help diff

zearo

例如
在如下系统方程中，假设控制 u 为 Y(2) 的导数，简单写成 u=diff(Y(2))，那么运算时 u=[]，不可用，如何解决呢？
主程序为
function DuffingM
P1 = 0.4;
P2 = -1.1;
P3 = 1.0;
P4 = 1.8;
w = 1.8;
tspan=[0:0.01:100];
Y0=[0 0];
[T,YY]=ode45('Duffing',tspan,Y0,[],P1,P2,P3,P4,w);

系统方程为：
function DY=Duffing(t,Y,flag,P1,P2,P3,P4,w)
switch flag
    case ''
       u = diff(Y(2))
     DY=[Y(2); u- P1*Y(2) - P2*Y(1) - P3*Y(1)^3 + P4*cos(w*t)];
end

eight

原帖由 zearo 于 2007-6-8 18:39 发表

                               
登录/注册后可看大图

例如
在如下系统方程中，假设控制 u 为 Y(2) 的导数，简单写成 u=diff(Y(2))，那么运算时 u=[]，不可用，如何解决呢？
主程序为
function DuffingM
P1 = 0.4;
P2 = -1.1;
P3 = 1.0;
P4 = 1.8;
w = 1.8; ...

请仔细看 diff 的例子。另外，从数学上说，要求某一点的导数，能脱离它的周围点吗？