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
相应的初值问题明显不同。因此产生了不同的结果