我正在尝试用python中的odeint来解决ODE。我正在做一个物理项目。我的目标是在球面坐标中计算轨迹的方向和角度，当我把物体的速度，我的位置，着陆点。例如，当我设置参数- 800m/s，(S50°E135°)(我的位置)，向南10公里(着陆点)。然后结果是这样的-拍摄180°(0°是北)，海拔37°。当我知道初始值时，我可以计算轨迹，我可以计算结果，它的轨迹，但我不知道怎么做才能完成我的</e
浏览 36提问于2019-11-25得票数 0
我想在COMSOL中求解具有特定边界条件的Schrödinger方程。作为一个常微分方程，schrödinger方程是这样写的( 1D)：其中E是将由边界条件确定的未知常数。到目前为止，所有用于求解微分方程的模板都包含一些通用形式，其中您必须在每个项之前指定常量的值。这不适用于上面的特征
浏览 40提问于2016-04-04得票数 2
1回答
Solve_ivp是一个初值问题的求解函数。，method=‘RK45 45’，t_eval=None，dense_output=False，events=None，vectorized=False，args=None，**选项)参数有趣，系统右侧可调用.呼叫签名很有趣。或
浏览 2提问于2020-02-18得票数 3
1回答
我想用Julia语言实现一个通用的Runge-Kutta步长函数。在Python中，我可以只传递一个函数作为这个RK4函数作为输入获取的参数之一。如果我在Julia中这样做，会有性能损失吗？我的函数如下所示： Vs_k1 = rhs(Vs,Ps) Vs_k2 = rhs((Vs,Ps) Vs_next = Vs+(1/
浏览 2提问于2016-06-27得票数 2
1回答
是否有一个包含Gauss求解器的Python2.7软件包，用于包含3个以上未知数的线性代数方程组？下面给出一个简单的例子，说明我想要解决的问题。如果没有可用的模板或包，是否有可能在python中解决这个问题？如果是的话，请你给我建议一下最好的办法。谢谢。具有三个未知数(x，y，z)的三个线性代数方程的例子：
X
浏览 1提问于2014-04-06得票数 1
回答已采纳
1回答
有谁能给我一些建议，如何在Python中解决在其中实现了时延的ODE？我似乎不知道如何使用scipy.integrate.odeint来做到这一点。我要寻找的内容应该是：b = 1/50a = 0.8tau = 0.5
u = [b其中N(t)，N(t - tau)等表示函数的时间自变量。有没有一个很好的库来解决这些类型的方程
浏览 2提问于2017-03-07得票数 7
我在Jupyter笔记本上运行Python 3。键入：并不像您期望的那样只返回goodbye。相同的代码在Python控制台和终端的IPython笔记本中都可以正常工作，那么为什么它在Jupyter中不能工作呢？
浏览 5提问于2019-03-31得票数 1
下面是用于常微分方程参数的最小二乘拟合的代码。使用了Python的“最小化”和“最小二乘”函数。已经尝试了不同的方法和ODE解算器/步骤(scipy ode/odeint)。这是一个在MATLAB中很容易解决的问题，但是Python总是返回最初的猜测。我希望你找出一个编码错误，否则我会对Python优化函数感到失望。Obj表示目标(剩余平方和)，ode函数(一阶
浏览 23提问于2019-02-08得票数 0
我试图将线性回归方法应用于9个样本的数据集，该数据集具有大约50个特性，使用python。但我可以得到HR和NNLS具有不同解的直觉，而LR和闭形OLS具有相同的目标函数，即最小化给定样本中观测值与由一组特征的线性函数预测的差值的平方之和。train_features) 对于
浏览 7提问于2017-09-27得票数 1
1回答
我写了一些代码来解决二阶微分方程，但它给出了一个完全错误的结果。我认为问题在于以正确的方式表达Euler方法。我也尝试了另一个二阶常微分方程，但在近似y(x)时也失败了。你能指出错误在哪里吗？解算ODE：解析解决方案：重写为2个ODE的系统；
v'
浏览 1提问于2017-04-25得票数 0
1回答
这是我尝试过的代码：import matplotlib as mpl plt.show()
orbit()这显然是错误的，因为我认为太阳的轨道是椭圆形的。
浏览 9提问于2022-11-03得票数 -1
2回答
从已知物理中模拟不确定性
、、、、
我有一个方程式：其中，a是一个输入常数，\lambda是一个依赖于s的非线性项.我知道\lambda\left(s\right)的真正解决方案是\sin \left( s\right) \cos\left(s\right)
我用真实解和方程生成了S的数据，不包括一系列初始条件的非线性项和a's的范围。solve_ivp(foo, (t[0], t[-1])
浏览 0提问于2023-05-25得票数 4
我试图演示如何使用系统传递函数的定义和python-control模块来“求解”微分方程初值问题的“解”(模拟解)。事实上，我真的是个新手。我有一个简单的微分作为例子：y'' - 4y' + 13y = 0，具有这些初始条件：y(0) = 1和y'(0) = 0。所以，在Python中，我正在编写这段代码(请注意，y_ans是这个微分IVP作为的答案。)
plt.plot(T,
浏览 14提问于2020-10-29得票数 1
回答已采纳
所以我有四个节点具有已知值，我想通过插值找到一个位于这四个节点之间的值，但这四个节点不是一个矩形。*我在这里发现，，我应该解下面的一组方程：，其中Ni是：。最后，这将导致求解以下形式的一组方程：
a*x+b
浏览 14提问于2018-03-02得票数 2
回答已采纳
(TLDR:当您将方程作为纯Python列表传递时，SymPy的linsolve函数无法求解通过将有限差分方法应用于ODE BVP而生成的线性方程组，但可以将该方程列表放入SymPy的Matrix函数中这可能是一个需要修复的错误，特别是考虑到在SymPy文档中，它们提供的the example让它们将列表作为参数传递给linsolve。) 我有一个边值问题的常微分
浏览 33提问于2021-08-02得票数 0
1回答
我的自变量是时间。我正在做一个物理项目，在这个项目中，我应该通过求解一个常微分方程来得到一个矩阵函数(矩阵的每个元素都是时间的函数)。我使用了numpy的ODE求解器，它给出了矩阵函数在不同时间求值的数值结果。现在有了这些矩阵和时间，我想为每个元素矩阵找到一个与时间相关的表达式，以获得与时间相关的矩阵。我听说我应该做的是找到一个hat矩阵，我猜预测值或拟合值将是我的7x7矩阵，响应值将是时间数组。那
浏览 2提问于2015-06-21得票数 4
7回答
从字符串求解方程到C++
、、、、
我想知道是否有人有关于如何做一些听起来简单但看起来不像编程的事情的信息或经验。其想法是:给出一个包含等式的字符串，例如："2*x = 10“(这很简单，但它可能会变得非常复杂，如sqrt(54)*35=x^2；以此类推……)程序将返回x=5，并可能给出他是如何到达那里的日志。对于信息，有这个网站()，它用PHP做同样的事情，但不是开源的。
谢谢你的帮助！
浏览 0提问于2010-05-17得票数 3
回答已采纳
2回答
扬起力量
、、
= 3.299804.5 ^-4.5 = 0.00114959详细信息
您可以编写用于控制台的函数或类似的构造如果不能使用控制台输入，则可以假设两个数字都通过标准输出或写入文件保存在变量和ouptut中。输出必须正确到至少4个有意义的数字。您可以假设a和b都是非零的。明显超过1分钟的运行时间是不可接受的。最少的字节数将获胜。
浏览 0提问于2014-08-04得票数 12
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
云直播活动推荐
腾讯云开发者
