
Stochastic differential equations with Simulink?
Posted:
Aug 6, 2013 10:03 AM


Hi,
is it possible to simulate stochastic differential equations with Simulink? The one I'm working with is an inverted pendulum with two PDcontrollers. The other controller has time delay in input and works only when certain criteria is met. This is the equation:
I*A''(t)=m*g*h*A(t)[P1*A(t)+D1*A'(t)+f1(A(td))+f2(A'(td))]+s*n(t),
f1(A(td))=P2*A(td) & f2(A'(td))=D2*A'(td), when A(td)*[A'(td)a*A(td)]>0 and f1(A(td))=0 & f2(A'(td))=0 otherwise.
A(t) is the pendulum angle, I,m,g,h,s,P1,D1,P2,D2,a are constants, n(t) is Gaussian noise (mean=0, variance=1), t is time and d is time delay.
After some reformulation I believe this is what should be simulated:
x(n+1)=x(n)+f(x(n),x(nk))dt+s*n(t)*sqrt(dt)
k=d/dt. The initial state is x0=[A(0),A'(0)]=[0.01,0].
Apparently the sqrt(dt) instead of just dt is the key problem here. If this is not doable with Simulink, any hint about how to make this happen is appreciated.
Thank you, Aino

