问题标签 [nonlinear-equation]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
267 浏览

r - 如何求解 R 和 uniroot 函数中的非线性方程

我怎样才能找到 x 的值使用 R 的方程

a(x^b)+k*log(1+(x^c))+log(1-(u^(1/d)))=0 对于值

u=0.1,c=0.8,k=1.2,d=1.5,a=0.9,b=1 他们得到的答案是 0.0539。但我没有得到这个答案。

0 投票
1 回答
117 浏览

python - 如何使用 python 从复杂的非线性系统中获取所有根(复杂和实数)?

我正在尝试解决复杂的非线性系统。问题是有很多根,其中一些是复杂的根。从这些根中,我只需要选择实部在 [0,1] 之间并且没有复杂部分的那个(例如:0.23+0i)。前任:

根1:1.02 + 2i

根2:0.23 + 1.23i

根 3:0.23 + 0i

...

这是我的系统:tau1 和 tau2 是我需要找到的变量。方程是 t1 和 t2,它们取决于 tau1 和 tau2

我尝试使用 fsolve,但使用 fsolve 我无法获得复杂的部分。

有没有办法在python中做到这一点?

非常感谢你的帮助!

0 投票
0 回答
127 浏览

julia - 似乎无法让 NLsolve 在 Julia 中收敛。你能提出一些建议吗?

我正在尝试使用 Julia 解决经济学中的生命周期问题,但我在使用 NLsolve 时遇到了麻烦。该模型归结为尝试解决两个二等式系统,以找到每个工作期间的最佳休闲时间和资本存量。退休后的经济代理人设置休闲= 1,我只需要求解一个资本非线性方程。这部分工作正常。它正在解决似乎崩溃的两个方程系统。

因为我对 Julia / 编程很陌生,所以任何建议都会非常有帮助。此外,我们将不胜感激有关代码各个方面的建议/要点/建议。该模型从最后一个时间段向后求解。

我的尝试

代码似乎在第 35 期中断,并出现错误“在非线性系统的解析过程中,对以下方程的评估导致了一个非有限数:[1,2]”但是在第 37 期,解决方案似乎变得很奇怪。

0 投票
2 回答
187 浏览

matlab - 求解一个有 16 个变量的线性方程

我正在尝试求解一个 16 变量线性方程。我在 Matlab 中尝试过,但它说找不到解决方案。有没有一种有效的方法来解决方程。

在 MATLAB 中,我按照下面提到的方式进行了尝试,

警告:无法找到明确的解决方案。有关选项,请参阅帮助。

在 sym/solve 中(第 317 行)

0 投票
1 回答
61 浏览

algorithm - 如何为依赖于 x 的非线性薛定谔方程中的初始条件编写子程序

我正在通过分步傅里叶法求解非线性薛定谔方程: i df/dz+1/2* d^2f/dX^2+|f|^2*f=0

使用初始条件:f=q*exp(-(X/X0)^24)。

但是对于|x|<1,我必须使用q=1的条件,否则q=0。因此,我为横向变量 x 编写了以下子程序(涉及横向变量的代码摘录):

但似乎“f”出错了,它是一个 1 x 2048 向量,所有条目都为零。我没有得到预期的结果。如果我的初始条件只有 f=q*exp(-(X/X0).^24), q=1, 就很简单了,但是在上面的条件下(q=1 for |x|<1, 否则, q =0) 我需要做什么?任何帮助都感激不尽。提前致谢。

0 投票
1 回答
48 浏览

plot - 如何绘制具有非线性微分方程的振动系统的频率响应?

我有一个用于 3 自由度振动系统的非线性微分方程系统。 微分方程组

首先,我想绘制 y、y_L 和 y_R 与时间的关系(对于 Omega 的给定值),然后我想绘制不同数量的 Omega 的域(y、y_L 和 y_R 的最大值)。不幸的是,我不擅长 Octave。我在 Octave 中编写了以下代码(基于其中一位用户提供的示例),但它以这个错误结束:“匿名函数体必须是单个表达式”。

如果有人可以帮助我,我将不胜感激。

这是代码:

0 投票
0 回答
563 浏览

python - scipy.optimize.fsolve 和 sympy 的求解函数的区别

我想弄清楚 Optimize.fsolve 和 SymPys 求解函数之间的区别?**计算是否存在差异,或者对使用哪种方法有限制?** Scipy fsolve 是否比 Numpy 求解函数更准确?

我发现的明显区别是 Scipy.Optimize.Fsolve 可以求解给定猜测值的非线性方程,而 SymPys 求解函数可用于求解包含符号数学变量的方程和表达式。

您可以分享的任何参考文献或网站将不胜感激。

0 投票
0 回答
75 浏览

r - 如何在 R 中的 nlmer 中实现自己的非线性函数?

我正在尝试实现一个新的非线性函数以在 lme4 包中的 nlmer 函数中使用。但我不确定问题是什么。这是我第一次尝试使用 nlmer,但我正在遵循我在互联网上找到的所有说明。第一个错误是关于我的数据框。

当我运行此代码时,我收到以下错误:

我在使用 lmer 函数时没有得到它(当然没有非线性函数)。这就是为什么我认为我的问题不在于我的数据框。

我无法停止思考的其他问题,固定效果中的部分:

如您所见,我的非线性函数也参与了我的固定效应公式,我不太确定如何实现它。我希望我的方法是正确的,但由于我的第一个问题,我找不到测试的机会。

0 投票
0 回答
35 浏览

matlab - 如何求解已经给定状态向量的耦合非线性 ODES 系统

我有一个非线性动力系统。我找到了运动方程,它产生了以下状态向量: 在此处输入图像描述

我使用的变量是

在此处输入图像描述

我不断收到各种错误。我究竟做错了什么。谢谢!

0 投票
1 回答
50 浏览

matlab - 我应该如何在 MATLAB 中将 syms 与求解命令一起使用?

我无法解决这组联立非线性方程来确定未知变量a, b, c, d, e

似乎我的代码进入了无限循环并且永不停止。