ode

    0热度

    1回答

    对于微分方程mx'' + kx = 0(其中x''是x相对于t的双重导数),如何解决这个问题x(t)?我的意思是如何得到这个公式: x(t) = c1*cos(sqrt(k/m)*t) + c2*sin(sqrt(k/m)*t) 我试了一下: t, g, k, m, w0, a_0, b_0, c1, c2 = symbols('t g k m w0 a_0 b_0 c1 c2') x =

    1热度

    1回答

    我想根据位置条件改变移动机器人的输入轮速度? if d > 0 & d < 0.4 p.WR = 0.51; p.WL = 0.5; elseif d > 0.4 & d < 0.8 p.WR = 0.5; p.WL = 0.51; elseif d == 0.4 p.WR = 0.5; p.WL = 0.

    0热度

    1回答

    我是相当新的Matlab和这个ODE求解器,下面是我的代码: 的main.m format short; tspan=[0 5]; y0=[0.30;-0.30;0;-0;0]; [t,y]=ode23s(@(t,y) pend(t,y),tspan,y0); figure(1) %subplot(2,1,1); plot(t,y(:,1),t,y(:,2),'k--') s

    -1热度

    1回答

    二进制开关我在试图解决在R A ODE跑起来对抗一个问题。我有一个参数Q,当它比另一个参数h更大时停止流入它。直到它到达的地方发生切换,然后停止运行,并给我的消息的时间点的ODE正常工作: DLSODA- At current T (=R1), MXSTEP (=I1) steps taken on this call before reaching TOUT In abov

    0热度

    2回答

    嗨,大家好,我在我的ODE系列中遇到了一个我的方程式问题。我有两个完全相同的方程,然而他们出现了非常不同的答案。我的方程式给出不同的答案有什么原因吗? 这两个方程是管理e和r的方程。 library(sigmoid) parameters <- c( a = 0.032, b = (9/140), c = (5/1400), d =

    2热度

    1回答

    如何在scipy.integrate.ode中找到默认参数?具体而言,atol和rtol在积分器dopri5?可以使用set_integrator方法设置参数,但我看不到任何查询方式,如果我不设置它们,请查找它使用的是什么。

    2热度

    1回答

    我很清楚如何以最简单的形式使用scipy.integrate.ode.integrate(t)函数,但API读取它也需要两个可选参数,即step和relax。目前的documentation没有关于这些参数的信息,也没有在example中使用。我想知道,他们做了什么以及哪些情况下他们有用?

    1热度

    2回答

    我想解决使用bvp4c(边界条件)的微分方程系统。我使用两个for循环来解决bvp4c的功能,但问题是循环先完成,然后只有最后一个值用于第二循环的迭代,是否有办法让它们同时工作?这意味着第一个for循环的第一次迭代使用第二个循环的第一次迭代值(而不是最后一次)。感谢 function RTrajfoll(X,Y) clf; for i = 1:length(X)-1

    0热度

    1回答

    我在写一些代码来求解二阶微分方程,但是它给出了一个完全错误的结果。我认为问题在于以正确的方式表达欧拉方法。我也尝试了另一个二阶ODE,但我也未能近似y(x)。你能指出哪里可能出错吗?请参阅图形和代码: 求解ODE: y"(x)=-1001y'(x)-1000y(x) y(0)=1, y'(0)=0 解析解: y(x)=(1000exp(-x)-exp(-1000x))/999 重写为2

    2热度

    1回答

    我试图解决这个问题: 其中U是 这里: s=c*e(t)+e_dot(t) and e(t)=theta(t)-thetad(t) 和 e_dot(t)=theta_dot(t)-thetad_dot(t) 其中thetad(希塔需要的话)= SIN(T) - 待跟踪即信号!我试着第一次使用odeint,它在t = 0.4后给出的误差是θ(上面的微分方程的解)平稳地下降到了0, 0并留