1回答
我的代码附在图片中。如何将解y( x )绘制为x的函数？此BVP y(x)的解由命令dsolve({de，ics}) 给出 w := 20 L := 200
Base_TD
浏览 23提问于2021-03-10得票数 1
在宇宙膨胀的研究中，我遇到了二阶和二阶微分方程。完整的方程是非常混乱的，所以让我给出一个原型：在这种情况下，我通常做的是让数学来解决这个问题但是，有时mathematica会出现一些错误，例如“x=d的奇异性或刚度”。现在，我试图用一些基本语言(阅读python )来解决这些问题。但是，我们求解二阶
浏览 4提问于2017-08-01得票数 0
最近，我一直在对一些动力学系统进行模拟，其中所有的动力学变量都是相互依赖的。因此，为了模拟动态，我在小的时间步长dt<<1上执行了循环，并在每次迭代中改变了数量。我得到了很好的结果，但由于迭代过程缓慢，模拟可能需要相当长的时间。一般来说，我听说应该避免像我使用过的for循环，因为它们会大大减慢模拟速度。然而，另一方面，我不知道如何在没有迭代的情况下以小的时间步长进行模拟。因此，我问你:对于一个动态系统，其中每个量都必须在超小的时间步长内变化
浏览 1提问于2014-06-14得票数 1
计算图是使用操作符重载和工厂函数生成的，例如sum()，exp()等。我已经使用反向累加实现了梯度的自动微分。然而，我发现实现二阶导数(黑森导数)的自动微分要困难得多。我知道如何进行单独的第二部分梯度计算，但我想不出一种智能的方法来遍历图形并进行累加。有没有人知道好的文章，给出了二阶导数的自动微分算法，或者开源库实现了相同的算法，我可以试着学习一
浏览 2提问于2010-07-04得票数 9
回答已采纳
1回答
我正在求解上的一个ODE函数，但它返回一个错误。我是Stan的新手，并且已经阅读了手册(版本2.23)，我构建了一个变量函数，并在后面的参数部分尝试使用它。
浏览 9提问于2021-03-06得票数 0
1回答
有一个类似于CalculateProfit(decimal a, decimal b, float c, TimeSpan d)的函数，它的每个输入参数都具有minimum、maximum和initial它的输出是平滑的，但不是线性的，它有多个峰值和下降。我想用暴力强迫它的输入，找到最大可能的输出。如何在不尝试每一种可能的组合的情况下对其进行优化？也许是二进制搜索？我认为算法应该在开始时使用大的增量步骤来找到大多数<
浏览 0提问于2018-12-02得票数 0
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这样的方法根本不考虑这种周期性(我已经尝试了所有可能的
浏览 19提问于2021-01-22得票数 0
回答已采纳
我已经在网上找了一段时间了，为一个物理引擎寻找集成的方法，我正试图编写一个有趣的代码(我喜欢那里的书呆子:P)。我找到了欧拉法、RK4和Verlet (以及时间校正版本)。我也一直在努力想出一些我自己的方法。我想知道你是否知道其他你觉得直观或有帮助的东西。谢谢。
编辑:感谢您到目前为止的所有帮助。至于澄清:也许我指的是数值积分。显然，我不能在C++中使用无穷小的数字，所以我必须使用瞬时积分(这是我在一些阅读中听到的术
浏览 0提问于2012-02-24得票数 1
回答已采纳
1回答
我在思考数值分析中的一个基本问题。
众所周知，在离散常微分方程时，二阶方法比一阶方法更精确，因为二阶方法的截断误差为O(dx^2)，一阶方法的截断误差为O(dx^2)。在这种情况下，由于dx^2 = 100和dx = 10，二阶方法是否不像一阶方法那样精确？我们在处理大规模问题时会遇到这种情况，比如气候建模(云的大小可能是几公里)。
浏览 0提问于2016-11-18得票数 1
我在MATLAB中用有限元方法求解一个二阶微分方程，其中函数f的二阶导数写为：现在f上的运算可以转化为一个矩阵，然后我可以找到特征向量，然后这些特征向量就是给定微分方程的解。所有这些都适用于x值的均匀网
浏览 1提问于2016-03-10得票数 0
1回答
我正试着完成一项任务，但我真的不知道如何做问题所要求的。我不是在寻找一个完整的答案，而是对我需要使用/做什么来解决这个问题的理解。以下是一个问题：
我们被要求为一类零阶J0(x)的贝塞尔函数提供一个插值。(a)为插值点x1 = 1.0、x2 = 1.3、x3 = 1.6、x4 = 1.9、x5 = 2.2创建一个7位小数点的数据表。通过点x1，x2，x3拟合二次多项式.使用此插值估计J0(1.5)。BesselJ到底是做什么的？如何通过这三个点来拟合二次多项
浏览 4提问于2013-02-24得票数 0
回答已采纳
2回答
我并不是在任何时候都明确地知道u(t)，但在离散的时间步骤中，我知道它是从先前的计算中得到的。我在网上找到的Python求解器的每一个例子(例如，表示scipy.integrate.odeint和scipy.integrate.ode)都知道作为时间函数解析的导数的表达式。有没有一种方法叫这些(或其他微分方程解)，而不知道导数的解析表达式？
现在，我已经编写了自己的Runge解决程
浏览 3提问于2017-11-15得票数 2
回答已采纳
我的密码import sympy as spsolution at n=0 (pre-subs): Eq(f(t), C1 + C2*t)
我的问题通用n的解决方案形式似乎没有准确地描述n=0的特定解决方案形式。具体来
浏览 0提问于2018-11-14得票数 1
回答已采纳
1回答
我正在做一个多自由度动力学问题，使用二阶拉格朗日方程。我用渐近得到了运动方程。现在这些方程在计算导数后变得相当长，尽管符号简化并不能进一步简化它。我的问题实际上是如何从这里解决这个三个二阶常微分方程的系统。我不知道如何转换这些公式，以便将它们与scipy.odeint()一起使用。我想到了替换，但有很多符号。所以我在搜索phi0，phi1和phi2，以及它们的一阶和二阶导数。初始条件都是phi=0和dphi=0。我希望有一种不用从头开始就能解决这个问题的</em
浏览 37提问于2020-04-05得票数 1
我正在尝试绘制一个带电粒子在Reissner-Nordström黑洞(带电黑洞)周围的轨道。 我有三个二阶微分方程和三个一阶微分方程。由于问题的性质，每个导数都是关于适当的时间，而不是时间t。我用四阶龙格库塔方法来积分轨道。我的困惑，以及我最有可能犯错的地方，是因为当你有一个二阶耦合微分方程时，你通常会把它归结为两个一阶微分方程。然而，在我的问题中，我得到了三个一
浏览 43提问于2021-03-17得票数 2
基本上，我编写了一个微分方程解算器类，它将从“方程”类中提取方程，并使用rK4方法求解它。我遇到的主要问题是，如果不通过继承扩展并获得访问权限，或者在ODE类中创建公式方法的特定实例，我无法找到将方法发送到另一个类的方法。public class Equations {
public double pres
浏览 1提问于2014-10-25得票数 0
1回答
MATLAB金融数据算法
、、、
因此，我有一个庞大的excel电子表格的历史期权数据，标准普尔100在不同的日期，2010年和现在的日期。我试图在每一个日期找到股票的概率密度函数。找出这个函数的方法是取看涨价格的二阶导数(列在第5栏中)。我希望我的MATLAB脚本工作如下：
使用曲线工具箱样条拟合函数，在每个日期创建一条相对于交易价格的看涨价格曲线。利用微分函数求出二阶导数的</em
浏览 2提问于2014-11-26得票数 0
1回答
Julia是否支持边值微分代数方程？我有一个隐式的常微分方程，它的可变质量矩阵有时是奇异的，所以我必须使用DAEProblem。我的问题是x1(t)和x2(t)的两个耦合的二阶常微分方程，通过设置x1'(t) = y1(t)和x2'(t)=y2(t)，我已经将它们转换为四个一阶方程。我在域的开头和结尾有x1和x2的值，但没有y1或y2 anywhere的
浏览 30提问于2021-07-08得票数 2
如果我有一个用一阶微分方程表示的二阶微分方程，它描述了一个质量-弹簧-阻尼器系统，当我不知道一阶微分方程时，我怎么能用欧拉法来画这个方程呢？我想用MatLab来做这件事。我只想要一个简单的大纲，告诉你该怎么做。有什么想法吗？
浏览 1提问于2011-01-28得票数 1
回答已采纳
云服务器
ICP备案
云直播
对象存储
实时音视频
腾讯云开发者
