我正在尝试用python中的四阶龙格-库塔法求解椭圆型微分方程。在执行之后，我得到了应该获得的实际绘图的一小部分，并伴随着一个错误，说明： "RuntimeWarning: invalid value encountered in double_scalars" import]=x[i]+dx#print(y)plt.plot(x,y)plt.ylabel('Theta')
plt.grid
浏览 16提问于2021-01-21得票数 2
回答已采纳
2回答
为了说明具有无穷解族的一阶常微分方程的柯西问题，我将绘制参数化的解，并有可能通过滑块控制参数值。完整的柯西问题是y'=sqrt(|y|)，y(0)=0，参数解是y_c(x):={0，if c=>x；(x-c)^2/4，if x=>c}。所以我会得到一个y=y_c(x)的曲线图，用一个滑块来控制c的值。
谢谢。
浏览 1提问于2012-04-15得票数 4
回答已采纳
我有一个非线性的一阶常微分方程，可以简化为其中i= sqrt(-1)，g( t )是t的某个给定函数。y(t)是我要找的函数，但我不知道如何处理这个方程，因为有conj(y)，有人知道如何在MATLAB中求解这个问题吗？编辑:t是真实的。
浏览 4提问于2015-03-20得票数 0
我花了一周多的时间来解决问题：我必须使用牛顿-拉夫森方法来得到u的近似解。我有这样做的脚本，但我需要“线性化”这个非线性常微分方程。k1-k4不是常量。在每个网格点(x=1-100)上，它们得到一个不同的值，该值被计算出来。初始条件是u(0)=0。
浏览 1提问于2013-03-21得票数 0
我想解一个形式为的非线性常微分方程 Theta2 = (C + j(Theta2))**-1 * (f(t) - g(Theta1) -h(Theta0)) 其中f()、g()、h()和j()是已经定义的函数Theta2和Theta1是Theta0对时间t的二阶和一阶导数。我一直在使用SciPy.odeint函数求解不带j(Theta2)项的方程，代码如下： from scipy.integrate import odeint
def g
浏览 17提问于2019-09-12得票数 0
我有一个关于绘制x(t)的问题，知道dx/dt等于下面的表达式，下面的微分方程的解。在t=0时，x的值是0。syms x
dxdt = -(1.0*(6.84e+45*x^2 + 5.24e+32*x - 2.49e+42))/(2.47e+39*x + 7.12e+37) 我想画出这个一阶非线性微分方程的解解析解涉及复数，因此这并不重要，因为这个方程模拟了现实生活中的过程，但Matlab可以使用数值方法求解该<
浏览 26提问于2021-11-06得票数 0
回答已采纳
我想解决一个由4个耦合的微分方程组成的系统，其中包含python (sympy)：solution = dsolve_system(eqs, ics={cP1(0): 0, cE1(0): cE1_0, cE2(0):cE2_0, cE3(0):cE3_0})非常感谢!
浏览 7提问于2020-12-21得票数 0
我正在尝试用数值方法解决一个边值问题，我的朋友问我，求解器是否适用于这些条件。求解器的页面没有给出我们拥有的条件，即bc(y(a),y(b), p) = 0，但是我们问题的形式是y(0) = some constant value和y'(b) = 0，给出了我们的Neumann条件，你需要重写函数来有一个像打靶法一样的一阶降阶吗
浏览 29提问于2019-12-09得票数 1
回答已采纳
我应该写一个MATLAB函数，它取一阶常微分方程，形式为y‘(T)= a*y(t) +b，初始点y(t0)=y0作为输入，并计算解的前15个点。还会绘制前15个点的解决方案曲线。我们要解的方程是:y‘(T)= 4*y(t)+1，初始点y(0)=0。
对于这个函数，我写了下面的代码，但这给了我一个关于y的错误。我应该如何正确地实现euler函数？
浏览 0提问于2012-12-17得票数 0
回答已采纳
如果我有系统动力学，为什么我不能自己手动集成它并得到结果。ode45函数的确切用途是什么？
浏览 6提问于2020-02-16得票数 0
我有一个项目，我的微分方程和线性代数课程，我必须使用计算机代数系统来解决线性，常系数常微分方程和非线性常系数的常微分方程。我必须在分析和数值上证明这一点。我的计划是在分析部分使用dsolve函数的两个函数，以及我通过学到的一个名为ODE1的函数。该函数的定义是用欧拉法逼近微分方程的解。这些函数都使用相同的输入参数，因此每个输入都可以定义一次，并且所有函数都将了解要使用的输入(可能将函数嵌套在一个调用
浏览 1提问于2018-06-18得票数 0
1回答
我正努力实施这些方案：
这就是我尝试过的：这个有什么问题吗？如何使用稀疏操作计算此值？
我试过这样做：如何实现基于新x计算现有x的后向部分？用除法可以吗？
full(L) =
0.3333 -1.0000 0.333
浏览 8提问于2019-11-24得票数 0
1回答
这个问题问的是一个系统，它有两个质量，每个质量连接两个弹簧，看起来像这样：外部弹簧具有弹簧常数kb，而内弹簧具有常数k。我写了一些代码来寻找系统的正常模式，我得到的频率是0.160，我非常确定这是正确的。然而，当我试图用线性代数解决同样的问题时，我得到了-1.2和-1.0的特征值。我也不太明白如何编写代码，因为到目前为止，我的运动方程不适合任何python代码示例，我也不知道如何将矩阵d^2/dt^2X传递给eigh()函数。运动方程是：
d^2/dt^
浏览 2提问于2020-11-22得票数 0
2回答
我有一个由两个一阶常微分方程组成的系统，它们是非线性的，因此很难在封闭形式下解析求解。我想把这个常微分方程组的数值解拟合到一个数据集上。我的数据集只针对属于ODE系统的两个变量中的一个。我该怎么做呢？没有帮上忙，因为只有一个变量。
浏览 0提问于2016-10-23得票数 1
我有一个联立方程计算器，它包含两个方程的字符串，它应该输出到X和Y的解。这是它的代码。transformations=transformations))sol = solve(eqs_sympy)
final_lst的一个例子：，sol只是输出一个空白列表，为什么会这样呢？
浏览 2提问于2021-08-08得票数 1
回答已采纳
3回答
我是提亚诺的新手，我试图实现一个反应扩散系统的数值积分器--这个版本的：
目前，我的表达方式是：import theano.tensor as
浏览 7提问于2016-02-21得票数 8
回答已采纳
1回答
我有一个关于解一阶复微分方程的问题。我使用Runge-Kutta，但答案似乎不正确。这是我的方程式：常微分方程的结果是可以的，但对于复杂的方程，它似乎不是正确的。我用gfortran编译了我的代码，并用gnuplot绘制了没有i的常微分方程的图表。这张图对于方程的解是正确的。但我在方程式中加入了i，图中只有一条直线。另
浏览 1提问于2016-01-20得票数 2
我目前正在学习一些拉格朗日力学，这在很大程度上进行得很好，但我决定尝试一个计算机练习，我不知道如何让Python求解这些方程。我试图解决的二阶非线性常微分方程是：Y的初始值是y0 = 0。
求解一个普通的一阶方程很容易，因为你只需创建一个函数，将它设置为某个值，然后使用odeint。没问题。但是现在我得到了一个二阶方程，我不知道如何告诉odeint或函数，我要找的是两次
浏览 26提问于2019-08-05得票数 2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
云直播
对象存储
实时音视频活动推荐
腾讯云开发者
