用Julia求非线性方程的零点
非线性方程的零点是指方程在某个变量取值下等于零的解。Julia是一种高性能的动态编程语言，它具有类似于Python的简洁语法和类似于C的性能。在Julia中，可以使用多种方法来求解非线性方程的零点。
一种常用的方法是使用牛顿迭代法。牛顿迭代法是一种迭代的数值方法，通过不断逼近函数的零点来求解方程。具体步骤如下：
- 定义非线性方程：首先，需要定义一个表示非线性方程的函数。例如，假设要求解方程 f(x) = 0，可以定义一个函数 f(x)。
- 初始化迭代：选择一个初始值 x0，作为迭代的起点。
- 迭代计算：根据牛顿迭代法的公式，计算下一个迭代点 xn+1 = xn - f(xn)/f'(xn)，其中 f'(xn) 表示函数 f(x) 在点 xn 处的导数。
- 判断终止条件：重复步骤3，直到满足终止条件。常见的终止条件可以是迭代次数达到一定的上限，或者两次迭代点之间的差值小于某个阈值。
- 输出结果：迭代结束后，得到一个近似的零点 xn。
非线性方程的求解是一个复杂的问题，牛顿迭代法只是其中一种方法。在实际应用中，还可以使用其他数值方法，如二分法、割线法、弦截法等。
关于Julia的使用，可以参考官方文档和教程，以及相关的开源库和包。以下是一些相关资源：
- Julia官方网站：https://julialang.org/
- Julia官方文档：https://docs.julialang.org/
- Julia中文社区：https://discourse.juliacn.com/
- Julia数值计算库：https://julialang.org/packages/
请注意，以上资源仅供参考，具体的使用方法和相关包的选择应根据实际需求和情况进行评估。
相关·内容
1回答
我想在Maple中找到以下微分方程的平衡点：
y' = h(t-1) with h(t) = 1 if t>= 0 and h(t) = 0 otherwise
我试着用这样的分段方程：
piecewise(t >= 0, h(t) = 1);
但我不知道如何继续，因为我需要求解h(t-1) = 0，而不是“h(T)=0”。
浏览 1提问于2013-10-02得票数 0
回答已采纳
1回答
我试着用下面的方法来求解双曲方程，以求x和y。我想知道在fsolve中使用相同的方程两次来找到第一个x然后是y是否有意义。我的代码如下：
from scipy.optimize import *
from numpy import *
import math
a = 1/(a_6**2)
b = 1/(b_6**2)
def function_hyper(loc):
x = loc[0]
y = loc[1]
F = empty((2))
F[0] = a*pow(x, 2) - b*pow(y, 2) - 1
浏览 6提问于2021-06-18得票数 0
我在试着解非线性方程组。 为此，我成功地使用渐近在线解决了一个由三个方程组成的系统。但对于7个等式，它不再起作用，返回给我一条"Operation timed out“消息。所以，我尝试使用标准的fsolve来解决这个问题。但是，虽然我只用一组简单的3个方程就可以得到很好的渐近结果，但fsolve给了我很差的结果。 这是我写的代码(感谢我在这个论坛上找到的)： def equations(p):
w1, w2, m = p
f0 = w1+w2-1
f1 = w1*m+w2*(1-2*m)-1/3
f2 = w1*m**2+w2*(1-2*m)**2-
浏览 27提问于2019-04-29得票数 1
回答已采纳
在MATLAB中求解约束线性和非线性方程组
我在MATLAB中用直接刚度法解决了一个有限元问题。这个问题现在被表述为一个非线性方程组：
KU=F or CF=U with C = K^{-1}.
该问题用3-Dd表示，有200多个节点。模型被用作变形估计器，力始终是变形的输入和输出。所以有很多线性方程。现在这个问题需要扩展，需要实现两种类型的约束。某些节点必须受到约束，以便它们只能围绕给定点以圆周运动方式移动。并且一些节点必须被限制在最大变形范围内。
总共有成百上千个形式为方程的方程：
a_1 X(1) + a_2 X(2) + ... - F(1) = 0
大约10个非线性形式的等式约束(也
浏览 3提问于2016-05-31得票数 0
1回答
具有多重方程的枕方程组
、、、
你好，当我有一个方程组时，我在使用fsolve函数时遇到了问题。
我的代码只让我有和我的系统中的方程一样多的变量，但是我们都知道，在实践中，你可以拥有比变量更多的方程，特别是当解不是一个具体的数字而是一系列的数字时，一个附加的方程可以帮助缩小“干草堆”的范围。
如何插入比变量更多的方程是我的问题？
假设我有以下非线性系统：
A/B=0.4583(3)
A/C=0.25
A/D=0.72(2)
B/C=0.54(54)
因此，我有以下代码：
from scipy.optimize import *
from numpy import *
def FUNC(arg):
A,B,C,D
浏览 0提问于2018-11-19得票数 1
回答已采纳
朱莉娅:当我有一个曲线图时，如何找到最佳拟合曲线/方程？我有一个用map绘制的图，但是我需要找到一个符合这个的二次方程？
浏览 62提问于2021-01-26得票数 0
我想知道并试图搜索它，但我不知道如何谷歌它，如何解决线性方程组，其中方程是平等的，彼此。我能够为两个方程编写一个具体的算法，但不能对N个方程进行动态计算。我将给出三个例子(我的方程大致是怎样的)：
C1，C2，C3，X -都不知道，但我不需要知道X的结果。
可以这样解释(最后一个方程C1 + C2 + C3 =1不在这里集成)：
请不要试图解决这个问题，我不确定这些方程是否有结果。我只是随机输入系数。但这就是我的方程式的样子。只有不同的系数。最好是在蟒蛇，但如果我知道如何用矩阵做数学的话，我会很高兴。我试着用两个未知数计算它，最后我得到了二次方程，所以在有三个未知数的情况下
浏览 10提问于2022-03-24得票数 1
回答已采纳
1回答
我会尝试描述我的问题，但最后有一些具体的问题。
我有一些困难，解决一个涡轮机问题，在Matlab。目的是用已知的物理测量和运行条件来计算涡轮的效率。为了做到这一点，我使用了科学论文中给出的方程式：
在我的例子中的问题是，要知道效率，我需要出口气体参数(温度，压力，速度…)这又取决于效率。我主要在一个名为EES (工程方程求解器)的软件中解决了这个问题，在这个软件中，方程可以按任何顺序写入，程序对它们进行分组并自动求解。由于程序中的限制，我切换到了Matlab，这不允许我扩展这个问题。
在matlab中，我首先以一个可解的顺序将所有的方程写成一个函数，除了求解问题所需的两个退出参数之外。将两个退
浏览 1提问于2012-05-31得票数 0
1回答
我有三个变量，如下所示：
float a = N * N + Q * Q - 1;
float b = 2 * M * N - 2 * N * x1 + 2 * P * Q - 2 * Q * y1 + 2 * s1 * r1;
float c = x1 * x1 + M * M - 2 * M * x1 + P * P + y1 * y1 - 2 * P * y1 - r1 * r1;
我想找出二次方程的根，因为我想用它来求圆的半径。
在Java中是否有内置函数来解决这个问题？一条好班轮就好了。
浏览 1提问于2014-09-25得票数 0
回答已采纳
1回答
我试图用numpy来求解如下形式的非线性方程组：
a(y-2.7)(1-exp(-a*z)) = (x-2.7)(1-exp(-z))
b(w-2.7)(1-exp(-b*z)) = (x-2.7)(1-exp(-z))
c(w-2.7)(1-exp(-b*z)) = (y-2.7)(1-exp(-a*z)
d([y+w]/2-2.7)(1-exp(-d*z)) = (x-2.7)(1-exp(-z))
很明显，系统中的方程和未知数一样多。a，b，c，d是上述系统的常量。这是最简单的系统，在其他情况下会有更多的方程。这些方程的解具有相似的数量级，因此我知道Levenberg-Marquar
浏览 4提问于2014-06-09得票数 0
回答已采纳
嗨，我正在尝试实现一个简单的塑性有限元代码。我不清楚如何设置通过拉格朗日乘数应用的位移控制。在线性问题的情况下，我执行了以下操作来应用显示控制，并且它起作用了。如果能量平衡是线性的，那么将得到以下方程：
然而，在非线性问题中，K是K(u)并且是非线性的。因此能量平衡方程及其导数如下所示
感觉不对劲。我该如何评估它？如何在非线性情况下修改方程的矩阵形式？那么矩阵方程是如何修改的呢？我使用添加了用latex输入的图像。如果不清楚，请让我知道。
浏览 9提问于2021-06-24得票数 0
2回答
我在R和Excel中建立了长度与干质量的幂回归方程，但系数不匹配。
我使用洪Ooi的从这个链接：。在该代码中，他们能够使用R代码从Excel复制幂方程。然而，当我尝试的时候，我得到了一些非常奇怪的系数。当用随机长度测试时，来自功率趋势线的Excel方程更准确。
守则如下：
#sample dataset of Lengths and Dry Masses
test <- structure(list(
Length = c(23, 17, 16, 25, 15, 25, 11, 22, 13, 21, 31),
DryMass = c(3.009, 1.6, 1, 4.177,
浏览 1提问于2019-08-19得票数 1
回答已采纳
1回答
我尝试用Stata来估计上述非线性模型。不幸的是，我对Stata并不满意。有没有人能帮我用Stata写上面的函数？在Stata中，如何在nl命令中写入地区虚设、时间固定效果和国家固定效果？有没有一种方法可以在Stata中写出上面的方程中的求和？或者，是否更容易估计每个单独区域的方程？
浏览 13提问于2017-10-25得票数 0
2回答
我有一道数学题，是用Maple解决的，但它必须用C++解决。问题是我找不到一个方程的数值方法(它是用内置的maple "solve“过程求解的)。我可以以某种方式将这个过程包含到我的代码中，或者找到显示Maple数值方法的完整描述的地方。这个方程是一个复杂的代数方程组。
浏览 4提问于2012-09-26得票数 1
1回答
我试图用机器学习来求解二次方程。我试图编写使用mxnet的代码，但是失败了，因为它使用线性回归，而我的问题是非线性的(据我所知)。
如果至少有一个实际根存在，则可以预测任意方程的解。
浏览 0提问于2018-01-24得票数 2
1回答
我有两个方程，一个是线性的，另一个是非线性的。我必须尽量减少第一个，同时最大化第二个。
用Z3可以实现这一点吗？它似乎无法对非线性方程组进行优化。它用非线性方程组返回“未知”。
浏览 0提问于2018-03-14得票数 0
回答已采纳
2回答
在Dymola 2021 X中显示方程关联图有一个新的特点，但关联图中的非线性方程组比平移日志中的非线性方程组多，因此我不知道如何比较这两种统计数据。如何解读这些信息？
--------------------------------------------------------------------
更新时间为12.23.2020:以下是使用Dymola2021x-HotFIX-1的结果
浏览 3提问于2020-12-02得票数 2
回答已采纳
1回答
朱莉娅说：
不符合段落内行内的大型LaTeX方程可以使用带有“语言”数学的带栅栏的代码块编写为显示方程，如下例所示。
```math
f(a) ={1}{2π}\int_{0}^{2π}(\alpha+R\cos(θ))d\theta
但是，在我的一个项目中，内联数学表达式似乎触发了“无效转义序列”错误：
"""
...
```math
Gu \sim PY(d(x，u)}，\theta{u u}，G{pi(U)})
...
"""
我是否不正确地理解了文档中显示公式的用法？
朱莉娅版本是0.7
浏览 0提问于2018-09-17得票数 1
回答已采纳
我有关于流体细丝直径(D/D0)如何随时间变薄的实验数据(t)。D/D0 (ydata)和t (xdata)是np阵列。
我想将这些数据与以下公式相吻合：
在这个方程中，我的已知变量是D/D1 = D/D0，D1=D0和t。
D/D0 = np.array([0.129 , 0.1613, 0.1935, 0.2258, 0.2581, 0.2903, 0.3226,0.3548, 0.3871, 0.4194, 0.4516, 0.4839, 0.5161, 0.5484, 0.5806, 0.6129, 0.6452, 0.6774, 0.7097, 0.7419, 0.7742
浏览 4提问于2022-03-08得票数 0
如何使用SymPy (by PyCall)求解方程，案例: solve(x^2 - 4，x)。
浏览 2提问于2016-11-15得票数 0
1回答
我有实验数据的x和y坐标，我正试着对其拟合一个椭圆。我的x数据在列表xpix中，y数据在列表ypix中。我已经用椭圆np.meshgrid和np.contour拟合了旋转的方程式。 xpix = [.....]
ypix = [.....]
x = np.linspace(-80, 80, 100)
y = np.linspace(-80, 80, 100)
X, Y = np.meshgrid(x,y)
#equation of rotated ellipse below, uses X and Y
F = (X*np.cos(63*np.pi/180) -Y*np.sin(63*np.p
浏览 23提问于2021-04-08得票数 0
我正在尝试建立一个丝杠传动的理论模型，将其作为力元素参与mbs模拟。这意味着
F=blackbox(u)
在上图中，我将系统简化为一个有代表性的模型。phi_1(t)和u_1(t)是螺母的给定运动。应计算引线的运动和返回的反作用力。它们是phi_2(t)和u_2(t)。
我写下了静态情况下的所有方程。但是我不确定我应该如何选择方程的运动来数值求解它们。
在运动学中，您可以看到有一个内部自由度。你可以选择它作为u_2，所以根据给定的phi_1和u_1，phi_2有一种可能的位置。
现在你可以根据u_2写出轴向力F和扭矩M，这样你就可以计算出影响螺纹的切向力和法向力。
在这一点上我遇到了
浏览 0提问于2013-08-23得票数 0
我有一个矩阵，它是某个参数A=A(x)的函数。我想找出这个矩阵变得奇异的点x。示例(我有一个很大的矩阵)：
syms x
A=[x sin(x); cos(x^2) 2.5];
到目前为止，我一直在象征性地计算矩阵的行列式，然后使用fzero或newtzero来找到这个特征方程的根。也就是说。
detA = det(A);
fzero(matlabFunction(detA),startingGuess)
然后我发现了这一点：，在这种情况下，人们提倡在任何情况下都不要使用行列式。
事实上，符号行列式的计算非常缓慢。然而，我尝试使用秩(A)代替在链接中的建议，它似乎不适用于符号矩阵。
有没有办法
浏览 3提问于2017-03-24得票数 2
回答已采纳
2回答
求多项式函数的所有零点
、、、
我正在使用Dymola 2013版本。我试着解决一个简单的数学问题，比如：
f= x^2 -4 ;
y=1;
f=y;
F和x被定义为实数。解决方案是2.36。但我需要用两种解决方案进行计算。所以2.36和-2.36！在我的问题中，f是一个多项式，像ax^3 + bx^2 +cx +d，y是线性的。Y= ax +b
我怎样才能得到这个问题的所有解决方案？X没有显式的值。X通常至少有两个解决方案。X应该是一个向量吗？在这种情况下，我遇到了方程的维数问题...有人能帮我吗？
浏览 2提问于2013-04-09得票数 1
在这篇文章中，，G5W给出了很好的解决方案，但是如果我们有两个未知数，那么它是一个未知成员(X)的方程呢？示例
3x - 5y = 2
如何搜索决策“？有必要将x和y值结合起来求解以下方程，从而得到2。(在这个简单方程x=4和y=2中的手工计算中)
或者如果我们和三个未知数有交易，就像这样
一般来说，这样的方程可以有无限多的解，如何推导出不超过三个解呢？
uniroot函数是如何适应这样的方程的。
浏览 0提问于2018-09-23得票数 2
回答已采纳
我有一个讨厌的表情，我在数学上玩。
(-X +(2X^2)/ (2 (-1 + X)^2 ) ((1 +2 (-1 + p) X- (-1 + p) X^2)/(-1 +X)^2)^(3/2))/X
我用平面z=0绘制了它，其中X和p都被限制在0到1之间：
Plot3D{讨厌方程在这里，0}，{p，0，1}，{X，0，1}
我决定得到这个平面交集的方程，这个方程是由讨厌的方程和z= 0生成的。所以我用了解题：
求解{那个讨厌的方程== 0}，{p，X}，Reals
输出结果甚至更糟糕，其中包含#符号(我不知道它是什么，我是数学新手)。有没有一种方法可以得到nasty方程和z=0之间的交线方程，
浏览 4提问于2013-11-18得票数 0
回答已采纳
2回答
我是Freefem++的新手，我想要解决的问题是费舍尔方程：
du/dt = d^2u/dx^2 + d^2u/dy^2 + k * u * (1-u)
du/dn = 0 - border condition
我试图用弱形式重新描述这个问题，但是Freefem在公式上出现了一个错误：
problem Fisher(uh, vh) = int2d(Th)(uh * vh/dt + Grad(uh)' * Grad(vh) ) - int2d(Th)(k * uh * vh) + int2d(Th)(uh0 * vh/dt) - int2d(Th)(k * uh * vh * uh);
浏览 4提问于2014-10-10得票数 0
回答已采纳
我想在函数$f$中计算一组向量(或元组)，但是Julia说这是不可能的。
例如:如果我有一个元组数组p=[(1,1), (1,-1), (-1,1), (-1,-1)]和一个函数f(x,y)=x+y。我想计算一下f(p[1]) = f(1,1)= 2。但是朱莉娅说这些类型是不相容的。
你能帮帮我吗？
浏览 5提问于2021-02-22得票数 5
回答已采纳
1回答
我想知道如何在Julia中使用plots (或替代方案)根据数据(而不是方程式)绘制3d表面。
我有几个包含x，y，z数据的表。下面是3个表的例子：
Table 1:
| x | y | z |
| 1.1 | 10 | 3.2 |
| 21.2 | 10 | 2.1 |
| 32.8 | 10 | 0.7 |
Table 2:
| x | y | z |
| 2.1 | 8 | 4.3 |
| 24.3 | 8 | 5.2 |
| 43.1 | 8 | 1.7 |
Table 3:
| x | y | z |
|
浏览 7提问于2021-06-22得票数 2
嗨，我正在试着用主定理解下面的方程：
T(n) =a；对于n<= 2
T(n) = T(√n) + a；否则
当我找到一个类似的方程()时，我想知道我的解是否正确。我得到T(n) = O(log log n)。对于我上面的等式，这是正确的吗？
浏览 73提问于2017-12-22得票数 1
我有一套4套政策制订和评价报告：
du/dt + A(u) * du/dx = Q(u)
其中，u是一个矩阵，包含：
u=[u1;u2;u3;u4]
A是一个4×4矩阵。A和q是u=u1;u2;u3;u4.的函数Q为4*1。
但我的问题是：
如何在MATLAB中求解上述方程？
如果用Matlab的PDE函数进行求解，能不能将其转换成一个简单的函数，而不是从Matlab的现成函数中得到的？
是否有任何方法可以显式地计算A、和Q。我的意思是，在每一个时间步骤中，我从前一个时间步骤的数据中计算出A和Q，并在使程序运行速度更快的方程中加入新的值。
浏览 8提问于2014-01-25得票数 1
1回答
Z3Py:从方程组生成抽象公式
、、、、
我的例子:方程组
伪码约束基
a = b+c
∧ e = a*c
∧ a = +2 ; some replaceable concrete values
∧ c = +18
溶液
b = -16
∧ e = -32
我想要的信息
在方程组中，我想得到以下知识：
抽象公式，我可以用它从给定的值(在约束基中)计算变量值(解决方案)。
(就像在高中，老师不仅希望看到结果，还希望看到这样一个经过改造的抽象公式。)
公式如.
b = a-c ; is an equivalent transformation from `a = b+c`
∧ e = (a-c)*c ; is
浏览 3提问于2013-02-28得票数 1
我知道阶跃传递函数是如何工作的，但是线性传递函数是如何工作的呢？你用的是什么方程？与具有两个输入和一个偏置的AND门的应答相关
浏览 3提问于2015-09-16得票数 2
在Dymola中，我经常会遇到一个非线性系统初始化失败，或者是一个在大型热流体系统中很难解决的刚性系统，但是对于一个简单的系统来说，不会出现这样的问题。我的问题是：
，，所以我想知道，求解非线性系统模型的最大能力是多少？例如，我最多可以在我的模型中包含多少个非线性方程？，在Dymola中是否有允许提高求解非线性系统能力的设置？，如何在不损害模型精度的情况下减少模型中的非线性方程的数量？
浏览 3提问于2020-04-24得票数 2
回答已采纳
1回答
求指数函数给定两点
、、、
我知道我可以用下面的命令在maxima中求解方程，但是如何解决两个不同的方程。
kill(all);
r:.5; a:1; b:-5.7; theta:theta; solve(a*e^(b*theta)=r,theta);
tex(''%);
我试着得到一个膨胀函数的方程，给出两点。我该怎么做呢。示例1在(2,12)，点2在(8,768)
浏览 3提问于2017-07-12得票数 1
回答已采纳
1回答
我想知道在下面的场景中我最好的方法是什么。我有8个未知数，然而实际上是无限数量的非线性方程，这使得系统超定。
未知数：
U M V N J S W N
方程式：
U*M + V * Catime1 – V*M – Mgtime1 = 0
J*M + W * Catime1 – W*M – Srtime1 = 0
U*N + V * Catime2 – V*N – Mgtime2 = 0
J*N + W * Catime2 – W*N – Srtime2 = 0
U*S + V * Catime3 – V*S - Mgtime3 = 0
J*S + W * Catime3 – W*S - Srt
浏览 0提问于2013-05-28得票数 0
1回答
方程的确定性和随机性部分
、、、、
我正在寻找一种既能解确定性方程又能解随机方程的数值方法。在确定性的情况下，我知道四阶RK方法是一个有价值的方法，非常有效。不幸的是，还没有成功地应用于随机方程(至少据我所知)。
现在我想知道的是，是否存在一种数值方法，可以求解这两个方程(大致上，我的意思是，与解析解相比)，在这种情况下，会是什么。例如，一个解析可解的随机方程就是Black-Scholes方程。
浏览 3提问于2018-10-31得票数 0
1回答
我想用Julia语言实现一个通用的Runge-Kutta步长函数。在Python中，我可以只传递一个函数作为这个RK4函数作为输入获取的参数之一。如果我在Julia中这样做，会有性能损失吗？
我的函数如下所示：
function uv_rk4_step(rhs,Vs,Ps)
Vs_k1 = rhs(Vs,Ps)
Vs_k1 = Ps.dt*Vs_k1
Vs_k2 = rhs((Vs+(1/2)*Vs_k1),Ps)
Vs_k2 = Ps.dt*Vs_k2
Vs_k3 = rhs((Vs+(1/2)*Vs_k2),Ps)
Vs_k3 = Ps.dt*Vs_k3
Vs_
浏览 2提问于2016-06-27得票数 2
我需要编写R代码，可靠地找到方程cos(x + b) = a的最小正解，变量a和b可以取任意值。
在0，2π范围内，方程有两种解，即abs(a) < 1。但如果我用acos求解这个方程，我只得到一个答案，这可能不是最小的正解。如图所示。
我需要做些什么才能得到最小的正解？
浏览 0提问于2019-05-11得票数 1
回答已采纳
我正在尝试实现apache commons数学优化包。我正在寻找一个"hello world“的例子来解决一个有约束的两个方程的非线性系统。我的系统是可微的，由两个独立变量组成。apache commons文档和彻底的google搜索没有给我提供任何可供我构建的代码示例--有什么建议吗？
浏览 0提问于2011-11-17得票数 3
大家好，Gekko开发人员，
我想为我的热能存储系统建立以下方程模型，
它们既是线性的，又是微分的。一些变量也是时间相关的，作为外部参数出现。从Gekko的例子中，看起来我只能添加这种格式的方程，比如mdl.Equation(T_a.dt() == mt*c_p_w*(T_b-T_a)) / d_w*c_p_w，但是对于不涉及微分的方程，我想添加为普通方程，也可以添加mdl.Equation(m[t] == mc - m[t])。
附注:我如何添加这两种样式，以及哪种模式和求解器将允许这样做。？
谢谢,
诚挚的问候,
拉胡尔
浏览 4提问于2020-11-18得票数 3
1回答
我在用朱莉娅的方法数值解一个非线性方程组。我在用Newthon方法。我唯一不知道怎么做的，就是计算雅可比矩阵。到目前为止，我还没有找到计算偏导数的函数。
我的系统：
f(x1, x2) = 2*x2^2+x1^2
g(x1, x2) = (x1-1)^2 + (x2-1/2)^2
谢谢你的支持，向你问好，斯齐蒙。
浏览 1提问于2019-01-20得票数 0
回答已采纳
我意识到S盒可以用AES实现非线性变换.然而，我不确定这如何使AES安全。例如，如果我们没有S盒，那么就可以从一组线性方程组中计算密钥：
C^1=Ax+k
C^2=AC^1+k
..。
y=AC^n+k
其中A为线性变换，k为密钥，C为中间密文，n为加密轮数，x为输入，y为最终输出。但是，如果我们加上一个S框，那么它所执行的替换就不可能表示为x，f(x)的函数，所以现在我们有：
C^1=Af(x)+k
C^2=Af(C^1)+k
..。
y=Af(C^n)+k
在我看来，这似乎也是高斯消除的牺牲品(通过将每个方程替换成下一个方程)，尽管这样一个函数用于S盒中的替换可能非
浏览 0提问于2021-06-19得票数 1
1回答
我是Julia编程的新手，我设法解决了一些一阶常微分方程，但当我想转到二阶时，我不知道如何使用求解器来实现所需的方程。 我想要解这个方程 y" + y = 0 在初始条件下 y(0) = 3
y'(0) = -0.5 我该怎么做呢？
浏览 52提问于2020-07-12得票数 3
Wolfram系统建模器有一个名为的函数，youtube视频( 20:58)介绍了“方程浏览器”，我想知道Dymola是否具有相同的功能，这将有助于调试非线性系统。
另外，它还欢迎其他一些解决方案，比如使用python来可视化方程结构。
这是“等式浏览器”的截图。
浏览 1提问于2020-05-08得票数 6
回答已采纳
1回答
我正在学习调试dymola中的模型，我可以看到，dymola将所有方程分成不同的组。例如，在dslog.txt文件中，有关于用于求解第7方程组的迭代变量的信息，但是在哪里可以找到第7组中使用的方程呢？
以下是我目前使用的设置
浏览 0提问于2019-12-11得票数 2
回答已采纳
在我的方程式中，我们有许多带有^2的表达式，以此类推。我想将"²“映射到^2，以获得类似以下内容： julia> a² == a^2
true 然而，在Julia中，上述内容并不是法律规范。有没有关于如何实现它的想法？
浏览 32提问于2019-08-28得票数 2
在Dymola建立了一个具有朗肯循环的太阳能电站模型，虽然初始化工作正常，但是在2840s的时候，由于非线性方程组的求解失败，导致仿真停止。如下面的截图所示，Dymola建议提供更好的起始值，但这没有任何意义，因为失败发生在模拟期间，而不是初始化过程中。有了更好的起始值，情况根本不会改善。我的问题是:我应该如何处理在模拟过程中无法求解非线性方程？ ?
浏览 71提问于2020-07-06得票数 1
1回答
我在Sympy (0.7.6.1)中定义了一个微分方程，例如：
>>> y = Function('y')
>>> x = Symbol('x')
>>> y0 = y(x)
>>> y1 = y0.diff(x)
>>> y2 = y1.diff(x)
>>> eq = y2 - 2*y1 + y0
>>> eq
2
d d
y(
浏览 6提问于2015-09-15得票数 1
回答已采纳
1回答
我想用Sympy解决简单的2D运动学逆问题。我知道x和y位置的正运动学方程。
x = l1*cos(theta1) + l2*cos(theta1+theta2)
y = l1*sin(theta1) + l2*sin(theta1+theta2)
如果我知道theta1和theta2这两个方程，如何用Sympy求解这两个方程？
浏览 2提问于2019-03-29得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
即时通信 IM
实时音视频活动推荐
腾讯云开发者
