- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
假设存在以下问题:我必须函数,其中之一由 t 参数化,并且两者都在自变量 z 上定义:f(t, z, *ags), g(z, *args)。
我想找到一对值 (t0, z0) 以便
函数及其关于 z 的导数是相同的f(t0, z0, *args) = g(z0, *args)和df/dz(t0, z0, *args) = dg/dz(z0, *args)。
我知道存在一个解决方案,并且我得到了一个合理的起点(tS,zS)。但是,至少其中一个函数仅在指定的区间 [zL .. zH] 上定义(我知道)。
我现在的问题是,这是在 python 中数值求解方程组的最佳方法。
我尝试了 scipys fsolve,但它似乎失败了,我认为是因为它无法处理有限的定义间隔。我尝试了 Differential_evolution 包,只是为了最小化复合函数,但这似乎完全是多余的。
我有所有函数及其导数的表达式(尽管它们很复杂)。
肯定有一个简单的Python寻根器,它能够求解两个非线性方程组,这两个方程仅在有限的区间内定义?!
出于某种原因,我发现的所有内容要么能够求解方程组,而不考虑限制,要么能够接受限制,但一次只能求解一个方程......
如果有人能指出我正确的方向,我将非常感激!
最佳答案
此问题需要假设 1) 存在多个解,2) 至少有一个解位于指定区间的范围内。如果不满足这些假设,您的问题就变成了优化——最小化函数之间的差异。
我生成了以下满足假设的示例:
import numpy as np
from scipy.optimize import fsolve
def f(t,z):
return t**2 + z + z**2 + np.sin(z)
def g(z):
return z**2
def dfdz(t,z):
return 1 + 2*z + np.cos(z)
def dgdz(z):
return 2*z
def solve(x):
t,z = x
#residual array
r = np.zeros(2)
#match equations
r[0] = (f(t,z) - g(z))**2
#match derivatives
r[1] = dfdz(t,z) - dgdz(z)
return r
x0 = [10,-10]
x = fsolve(solve,x0)
t,z = x
print(t,z)
在这些初始条件下,答案是[3.06998012761 -9.42477800292]。通过在超出范围时添加残差,可以使用 fsolve 来要求有界变量。这是一个相当老套的解决方案,不是很健壮。我们可以通过以下修改将 z 绑定(bind)到 [-5,0] 的范围(其中有另一个正确答案):
def solve(x,zL,zH):
t,z = x
#penalize deviation from range
if z < zL:
e = z - zL
elif z > zH:
e = z - zH
else:
e = 0
#residual array
r = np.zeros(2)
#match equations
r[0] = (f(t,z) - g(z))**2 + (e)**2
#match derivatives
r[1] = dfdz(t,z) - dgdz(z)
return r
x0 = [10,-10]
x = fsolve(solve,x0,args=(-5,0))
t,z = x
print(t,z)
新答案是[1.77245385 -3.14159263]。
如前所述,有界变量在优化问题中更常见,因此优化包可以更直观地处理它。在 python 中,有标准的 scipy.optimize.minimize,或更强大的包,例如 Pyomo或GEKKO 。这是 GEKKO 中的等效问题:
from gekko import GEKKO
#initialize a GEKKO model
m = GEKKO()
#add GEKKO variables
t = m.Var(value = -10)
z = m.Var(value = -10, lb=-5, ub=0)
#define the constraints
m.Equation(t**2 + z + z**2 + m.sin(z) == z**2)
m.Equation(1 + 2*z + m.cos(z) == 2*z)
#solve a system of non-dynamic equations
m.options.IMODE = 1
m.solve(disp=False)
print(t.value,z.value)
[-1.772454] [-3.142608]
关于python - 使用 python 在有限定义区间上数值求解非线性方程组,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/50133906/
我在 Python 中使用 matplotlib,并制作了一个带条形的直方图。现在,当直方图出现时,仅 5 的倍数出现在 x 轴上,1000 的倍数出现在 y 轴上。对于 y 轴,这完全没有问题,但对
我正在使用 JavaScript 和 jQuery。我有以下脚本每 30 秒提醒一次 hi。 $(document).ready( function() { alert("hi"); setI
已结束。此问题正在寻求书籍、工具、软件库等的推荐。它不满足Stack Overflow guidelines 。目前不接受答案。 我们不允许提出寻求书籍、工具、软件库等推荐的问题。您可以编辑问题,以便
在 Numpy(python 包)中,可以使用语法 numpy.linspace(minValue, MaxValue, numberOfSamples) 构造 float 的离散区间。 . 我看到
所以我想在 -3 到 3 的区间内制作一些数字,以便在下面绘制这些函数,所以我想要尽可能多的数字。 我这样做: double k[601]; double y[601]; for (int i = 0
我有一个 Postgresql 表,用于存储有关计划进程的信息,包括上次执行进程的时间。不同的进程对其运行频率有不同的要求。 我列出了需要重新运行的进程列表: SELECT * FROM proces
如何正确使用此类带日期间隔的查询 @SqlUpdate("delete fromlogin where created < now() - ':days days' :: interval") v
我正在尝试计算图中的间隔,我在维基百科上找到了算法的数学描述: http://en.wikipedia.org/wiki/Interval_(graph_theory) H = { n0 }
我有一个基于 Informix-SQL 的 Pawnshop 应用程序,该应用程序根据黄金的重量和纯度计算应向客户贷出多少钱。当铺的最低贷款额为 5.00 美元。当铺员工通常会借出以 5 或 0 结尾
我将 NHibernate 与代码映射一起使用,并且我有一个由此公式创建的属性。 Property(x => x.IsInOverdue, mapper => mapper .Fo
我正在尝试从头开始为 Beta 分布编写卡方拟合优度检验,而不使用任何外部函数。下面的代码报告“1”适合,即使来自 scipy.stats 的 kstest 返回零。数据是正常分布的,所以我的函数也应
如何在 C# 中将任何值四舍五入到 10 区间?例如,如果我有 11,我希望它返回 10,如果我有 136,那么我希望它返回 140。 我可以很容易地用手做 return ((int)(number
如何在 Go 中表示 PostgreSQL 区间? 我的结构看起来像这样: type Product struct { Id int Name
我想编写一个函数,将数值限制在封闭的 0,1 区间内: func clamp01(_ value:T) -> T { return value 1 ? 1 : value } 在 Swift 3
我有一个简单的表格,用于存储来自在线仪表的降水读数。这是表定义: CREATE TABLE public.precip ( gauge_id smallint,
a = y def __gt__(self, y): return not self.x > y def __eq__(self, y): return
我正在处理 pandas 数据框 D=pd.DataFrame(data=[1.0,2.0,2.0,2.0,5.0,3.0,2.0,2.0,5.0,5.0,8.0,1.0]) 我识别低于特定阈值的值
我编写了一些C++代码来解决此问题: #include #include using namespace std; unsigned int countSetBits(unsigned int n
好的,我知道之前有人用一个有限的缩放示例问过这个问题 [-1, 1]间隔 [a, b] Different intervals for Gauss-Legendre quadrature in num
我是一名优秀的程序员,十分优秀!