我正在尝试实现中点公式y[n+1/2] = y[n] + h/2 f (x[n], y[n])和y[n+1] = y[n] + h *f (x[n] + h/2, y[n + 1/2]) ? 我的函数是 function [ x, y ] = Midpoint_ODE ( f, xRange, yInitial, numSteps )
% f = name of file with functionError in Midpoi
浏览 27提问于2020-01-13得票数 0
回答已采纳
我有一个模型，我需要更改当前在我的Simulink模型中的每个块的采样时间。
问题是我有太多的块，这使得为每个单独的块修改这个参数变得很麻烦。有没有办法对一组块进行更改？
浏览 0提问于2014-07-14得票数 0
回答已采纳
我正在用Matlab来解一个微分方程。我想强制ode45采取恒定的步长，所以在求解方程时，它总是在T轴上递增0.01。我该怎么做呢？ode45一直在采取优化的、随机的步长，而我似乎想不出如何让它采用一致的、小的0.01步长。0.001,'Stats','on');
%clf;tspan=[to tf];
%tspan = t0:0.01:t
浏览 0提问于2010-03-31得票数 6
回答已采纳
然而，我使用的数学在收敛时足够挑剔，所以我希望对使用哪种算法有更多的控制。我尝试过使用ode (特别是vode)，但由于所有的移动部件，我很难对其进行调试。我真的对分解我正在使用的ode解算器的所有细节不感兴趣:我只想能够告诉odeint使用一个特定的算法并结束它。这是可能的吗，还是我需要调试ode？
浏览 0提问于2019-03-05得票数 2
6回答
下面是我使用方法对进行数值求解时获得的/代码。然而，结果与我的教科书结果不一致，有时甚至不一致得可笑。代码有什么问题？
浏览 0提问于2010-05-30得票数 0
回答已采纳
1回答
Matlab解算器中显示的弹跳球示例提供了一种在触发even (https://www.mathworks.com/help/matlab/math/ode-event-location.html)时终止积分的方法在此示例中，当height为0时，它将终止。也可以在坡度从正值变为负值(球飞行的顶点)或从负值变为正值(当球落地时)时终止。有没有办法在Math.Net RungeKutta.Fou
浏览 23提问于2020-09-21得票数 0
2回答
获取输出相对于输入的梯度
、、、
目前，我正试图用Py火炬实现一个ODE解算器，我的解决方案需要计算每个输出wtr到其输入的梯度。range(len(y)): #compute output grad wrt input 我想知道是否有一种更优雅的方法来计算批处理中的每个输出的梯度我试着用：
torch.autograd.backward(x,y,reta
浏览 12提问于2020-11-24得票数 0
回答已采纳
4回答
同样的问题由MATLAB的刚性求解器(ode15s和ode23s)很高兴地解决了，但是我不能使用它(即使是来自Python，因为MATLAB的Python绑定没有一个实现回调，而且我需要将一个函数传递给ODE求解器)。我正在尝试PyGSL，但它非常复杂。如有任何建议，将不胜感激。
编辑:我在PyGSL中遇到的具体问题是选择正确的step函数。其中有几个，但没有直接类似于ode15s或ode23s (
浏览 2提问于2010-01-18得票数 21
回答已采纳
1回答
在了解如何在八度/matlab中实现事件时，在微分方程的解算方面，我遇到了一些困难。dy = odefun(t,y)endfunction
[T,Y] = ode45(@odefun,[0 12],[1],options);
现在我想介绍一个随机事件。例如，在某个固定的时间步长上，我</
浏览 2提问于2014-02-01得票数 2
回答已采纳
我在拉格朗日力学中遇到了下面的微分方程系统。你能推荐一种数值方法吗，有相关的链接和参考资料，关于我如何解决它。另外，在Matlab或Mathematica上有没有更短的实现？
浏览 0提问于2012-09-07得票数 0
我在Matlab中有一个ODE模型，对此我感兴趣的是执行一些参数扫描。我正在尝试将下面的代码从Matlab移植到R initial_conditions(6)=sweep1(i); parameters(3)=sweep2(j);
[t,y] = ode23s(@(timespan, initial_conditions) MODEL(timespan循环相同，求解
浏览 0提问于2013-12-19得票数 0
回答已采纳
我正在尝试解决液体流出储罐的问题；利用连续性和能量守恒方程。最初，储罐有一些压力和一些焓，它们在2D查询表中用于获得其他流体属性。当液体离开储罐时，压力和热焓发生变化；因此，导致了我的问题。我在这里查看了很多类似的东西，但我看到的很多都是随时间变化的变量。所以如果时间< 2，那么x=0 ..如果time >=为2，那么x= 2，我试图找出一种重新排列方程的方法，使其成为时间的函数，但它不是时间的直接函数，比如y(t) = m*t + b，我
浏览 0提问于2019-05-07得票数 0
1回答
我正在尝试通过用Matlab编写代码来解决一个微分方程系统。我在这个论坛上发帖，希望有人能以某种方式帮助我。我有一个包含10个耦合微分方程的系统。我使用ode45来求解这两个系统，使用相同的初始条件。从t0到tfinal，两个系统的时间向量也是相同的。当我通过绘制差值来比较矩阵ye0和yeL中time=tfinal的第10个状态变量的值时，我发现它在一段时间内变为负值。(使用命令：plot(te0,ye0(:,10)-yeL(:,10)))。对于从t0
浏览 3提问于2012-06-24得票数 0
我试图为anaconda安装一个包，结果显示了这条消息：在那之后，它搜索冲突，几个小时后我取消了这些冲突。例如，推荐的修复方法是使用conda install python=3.7.5降级python，但是这样做也会出现同样的问题。我<
浏览 126提问于2021-06-05得票数 2
1回答
我正试图为一个学校项目解决一个ODE问题，我遇到了一个变量的问题。想知道我能不能得到帮助。我用的是ODE15s解算器。options = odeset('RelTol',2.22045e-14, 'AbsTol', 1e-6);所讨论的变量是条件语句块的一部分，调试<e
浏览 0提问于2018-04-09得票数 0
回答已采纳
1回答
我是PyMC ()的初学者，正在尝试构建一个具有动态组件的模型，本质上是在每次调用模型时求解一个小的常微分方程(ODE)系统。我尽我所能搜索了谷歌和(废弃的) PyMC邮件列表，得出了模型结构(pk_model.py)和调用者(pk_fit.py)：
我将非常感谢在这
浏览 3提问于2013-07-28得票数 0
我正在为安卓开发一个2D，水下，动作角色扮演游戏，使用Box2D作为物理引擎，主要用于碰撞检测，碰撞响应和游戏中角色在由墙壁，岩石和其他生物组成的环境中的移动。
我正在尝试实现角色动画(根据)。问题是，当发生冲突时，我想抑制冲突响应。有人告诉我应该在碰撞事件中设置'NoCollision‘，我假设你在Pre-Solve事件中通过b2Contact::SetEnabled(false)来做这件事。但是，libgdx JNI实现
浏览 1提问于2011-04-11得票数 1
1回答
我想集成一个控制问题，即形式为dx/dt = A.f(t)的常微分方程，其中x(t)是R^3中的函数，f(t)是R^4中的函数，A是矩阵3x4。在我的特殊情况下，f(t) = F'(t)，即函数F的时间导数。此外，F是1周期的。因此，在区间0，1上积分ODE应该再次得到起始位置。然而，像scipy.integrate的solve_ivp这样的方法根本不考虑这种周期性(我已经尝试了所有可能的方法，比如RK45, Radau, DOP853, LSODA
浏览 19提问于2021-01-22得票数 0
回答已采纳
1回答
我正在尝试用CNN制作数独解算器。我想知道是否有一种方法可以在不使用完全连接层的情况下计算输出拼图和之间的损失？如果是，是如何实现的？
浏览 1提问于2020-08-12得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
