Graph wolfram alpha和倍频程下不同的常微分方程图

Graph wolfram alpha和倍频程下不同的常微分方程图,graph,octave,Graph,Octave,我用这个倍频程代码来解微分方程 # Define the right-hand side of the equation: xvall= -11 ;#xvall xvalu= 10 ;#xvalu range=5000; function ret=f(x,t);ret= t ;end; # ywill be the values of the function at these moments of time. t=linspace(xvall,xvalu,range); y=lsode (

我用这个倍频程代码来解微分方程

 # Define the right-hand side of the equation:
xvall= -11 ;#xvall
xvalu= 10 ;#xvalu
range=5000;
function ret=f(x,t);ret= t ;end;

# ywill be the values of the function at these moments of time.
t=linspace(xvall,xvalu,range);
y=lsode ('f', 2, linspace(xvall,xvalu,range));
y
plot(t,y);
我得到了这样的图表

但是当同样的条件传递给wolfram alpha时 我得到了y值从60到0的图表

图形是 为什么图形在两种情况下表现不同。
要指定常微分方程的初值问题,需要定义初始条件。对于八度音程,您指定了
x(-11)=2
,因为
xvall=-11
,对于Wolfram Alpha,您指定了
y(0)=2
。这就是为什么有两种不同的解决方案

倍频程

八度音阶的
lsode(f,x_0,ts)
解决了以下初值问题

 dx/dt = t
 x(t_0) = x_0
 t in ts
此处
ts
被指定为间隔
[t\u 0,t\u 1]
中的一组点。您已经指定了
t\u 0=-11
t\u 1=10

在封闭形式下,这个问题的解决方案是
x=(t^2-117)/2

Wolfram

对于Wolfram,您使用了半形式语法:

Runge-Kutta method, dy/dx = x, y(0) = 2, from -11 to 10, h = 0.25
在封闭形式下,这个问题的解决方案是
y=(x^2+4)/2


相应的初值问题明显不同。因此结果不同。

要指定常微分方程的初值问题,需要定义初始条件。对于八度音程,您指定了
x(-11)=2
,因为
xvall=-11
,对于Wolfram Alpha,您指定了
y(0)=2
。这就是为什么有两种不同的解决方案

倍频程

八度音阶的
lsode(f,x_0,ts)
解决了以下初值问题

 dx/dt = t
 x(t_0) = x_0
 t in ts
此处
ts
被指定为间隔
[t\u 0,t\u 1]
中的一组点。您已经指定了
t\u 0=-11
t\u 1=10

在封闭形式下,这个问题的解决方案是
x=(t^2-117)/2

Wolfram

对于Wolfram,您使用了半形式语法:

Runge-Kutta method, dy/dx = x, y(0) = 2, from -11 to 10, h = 0.25
在封闭形式下,这个问题的解决方案是
y=(x^2+4)/2


相应的初值问题明显不同。因此结果不同。

要指定常微分方程的初值问题,需要定义初始条件。对于八度音程,您指定了
x(-11)=2
,因为
xvall=-11
,对于Wolfram Alpha,您指定了
y(0)=2
。这就是为什么有两种不同的解决方案

倍频程

八度音阶的
lsode(f,x_0,ts)
解决了以下初值问题

 dx/dt = t
 x(t_0) = x_0
 t in ts
此处
ts
被指定为间隔
[t\u 0,t\u 1]
中的一组点。您已经指定了
t\u 0=-11
t\u 1=10

在封闭形式下,这个问题的解决方案是
x=(t^2-117)/2

Wolfram

对于Wolfram,您使用了半形式语法:

Runge-Kutta method, dy/dx = x, y(0) = 2, from -11 to 10, h = 0.25
在封闭形式下,这个问题的解决方案是
y=(x^2+4)/2


相应的初值问题明显不同。因此结果不同。

要指定常微分方程的初值问题,需要定义初始条件。对于八度音程,您指定了
x(-11)=2
,因为
xvall=-11
,对于Wolfram Alpha,您指定了
y(0)=2
。这就是为什么有两种不同的解决方案

倍频程

八度音阶的
lsode(f,x_0,ts)
解决了以下初值问题

 dx/dt = t
 x(t_0) = x_0
 t in ts
此处
ts
被指定为间隔
[t\u 0,t\u 1]
中的一组点。您已经指定了
t\u 0=-11
t\u 1=10

在封闭形式下,这个问题的解决方案是
x=(t^2-117)/2

Wolfram

对于Wolfram,您使用了半形式语法:

Runge-Kutta method, dy/dx = x, y(0) = 2, from -11 to 10, h = 0.25
在封闭形式下,这个问题的解决方案是
y=(x^2+4)/2

相应的初值问题明显不同。因此产生了不同的结果