如何解决R中有两种情况的SDE?

如何解决R中有两种情况的SDE?,r,differential-equations,stochastic,R,Differential Equations,Stochastic,我想用R解以下随机微分方程: \分形{dx}{dt}=f(x)+sigma*dW f(x)=a+bx+cx^2(对于x\leq 1)f(x)=a+bx(对于x>1) 及 西格玛=d^2 其中(a、b、c和d为常数) 我尝试使用: f = expression(a+bx+cx^2) s = expression(d^2) solution <- sde.sim(X0=0.6, t0=0, N=2000, delta=0.01, drift = f, sigma = s ) f=表达式(a

我想用R解以下随机微分方程:

\分形{dx}{dt}=f(x)+sigma*dW

f(x)=a+bx+cx^2(对于x\leq 1)f(x)=a+bx(对于x>1)

西格玛=d^2

其中(a、b、c和d为常数)

我尝试使用:

f = expression(a+bx+cx^2)
s = expression(d^2)

solution <- sde.sim(X0=0.6, t0=0, N=2000, delta=0.01, drift = f, sigma = s )
f=表达式(a+bx+cx^2)
s=表达式(d^2)
解决方案1)


很抱歉,数学表达式包含得不好。我不知道如何在这里写latex。

可能是这样的where
(x可能是这样的where
(x
f = expression(1+ 2 * x + (x <= 1) * 3*x^2)
s = expression(2^2)

solution <- sde.sim(X0=0.6, t0=0, N=2000, delta=0.01, drift = f, sigma = s)