我试着用下面的方法来求解双曲方程,以求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
我有一个任务,我需要做一个单一的定义函数,运行牛顿方法,然后能够插入其他定义的函数到它,它将解决所有这些问题。我编写了一个对有一个变量的方程有效的方法,我只需要对系统中的一个变量进行求解,但是我不知道如何在代码中解决这些问题,而不需要对所有四个变量进行求解。
我编写的运行newton方法的函数如下:
def fnewton(function, dx, x, n):
#defined the functions that need to be evaluated so that this code can be applied to any function I call
def f(x
我目前正在尝试使用方程的"rigid_body_tree.cc“文件计算库卡手臂的雅可比矩阵: Tau = J^T*F,其中Tau是库卡手臂的7个关节力矩,F是笛卡尔作用力和末端执行器处的扭矩,J^T是换位的雅可比矩阵。
drake中有一个名为transformPointsJacobian的函数,它接受缓存、点、from_body_or_frame_ind、to_body_or_frame_ind和in_terms_of_qdot。
该函数首先计算几何雅可比矩阵,输出6x7矩阵(kuka有7个关节)
然后,它采用该矩阵并使用它来确定3x7雅可比矩阵,计算如下:
J.template b
下面是一个非线性方程的例子:
x+2y=0
x^2+y^2=1
很明显,这个方程有两个根。
但是,我使用GEKKO只能得到一个根,改变x,y的初始值可以得到另一个根。
from gekko import GEKKO
m = GEKKO()
# first solution
x,y = [m.Var(1) for i in range(2)]
m.Equations([x+2*y==0,x**2+y**2==1])
m.solve(disp=False)
print(x.value,y.value)
# second solution
x.value=-1
y.value=-1
m
我有一个类似的问题,但尽管我得到了社区的帮助,但我无法解决这个问题。
我有以下问题:球向移动线移动。我想知道球什么时候和线相撞。
我尽最大努力解决这个问题:
所有变量:f(x) ... The line of the ball movement
g(x) ... The line before it moved
g1(x) ... The line after it moved
g2(x) ... The line after it moved + the radius (when the ball center hits this line, the ball border hits
我正在尝试使用lsqcurvefit将Arrhenius形式的方程拟合到一些数据点。
D = D0 * exp( -Ea / ( R * T )); % Arrhenius equation for curve fitting
D0和Ea是我要寻找的值。T是温度,代表X,D是系数,代表Y,R是气体常数。由于matlab在没有提供雅可比矩阵的情况下找不到解决方案,我计算了雅可比矩阵并包含了一个函数,正如我在早期的帖子中@m7913d所演示的那样([Matlab曲线拟合不适用于较小的值(1e-12),what can I do?](https://stackoverflow.com/q
在牛顿法中,求解非线性方程组需要找到雅可比矩阵和雅可比矩阵逆的行列式。
这是我的组件函数,
real function f1(x,y)
parameter (pi = 3.141592653589793)
f1 = log(abs(x-y**2)) - sin(x*y) - sin(pi)
end function f1
real function f2(x,y)
f2 = exp(x*y) + cos(x-y) - 2
end function f2
对于2x2的情形,我计算雅可比矩阵和雅可比矩阵的逆的行列式,
x = [2,2]
h = 0.
这可以分解成一个简单的三组方程:
a + b = 3
b + c = 5
a + c = 4
怎样才能最好地近似这些值呢?注意,在实际应用程序中,我将有更多这样的总计和变量。特别是,我想找出是否有可能根据商品清单和从杂货收据中得出的总数来有效地估算食品成本。我假设,如果我能计算出成本,将有不同的精度范围,我可以预期,因此,一个额外的将是知道有多可能的近似是正确的,在什么范围内的价格必须是。
编辑:我只是看不出这得到了一个我觉得舒服的答案,因为我一开始就没有把这个问题正确地组织起来。
在一个编程竞赛中,一个问题是:
计算方程的所有解:x + 4y + 4z = n。您将获得n,您将确定解决方案的数量。假设x,y和z是正整数。
我曾考虑使用三重循环(蛮力),但它是无效的,造成时限超过。(由于n可以是=1000 000):
int sol = 0;
for (int i = 1; i <= n; i++)
{
for (int j = 1; j <= n / 4; j++)
{
for (int k = 1; k <= n / 4; k++)
{
if (i + 4 * j + 4 * k == n)
sol
我做了一个代码(为了乐趣和练习)来平衡化学方程式。我想试着平衡一下N + A,所以它是= Z
N = 2
A = 2
Z = 6
if N + A != Z:
print('X')
balancer = ???
平衡器应该是3,所以如果我建立一个方程,例如(均衡器x,N) +A=Z,它将是真的。如果没有直接输入,我将如何使平衡器变成三。
谢谢:)
我正在用PHP创建一个CAS (计算机代数系统),但是我现在卡住了。我正在使用。
现在我写了一个记号赋予器。它将转换如下方程式:
1+2x-3*(4-5*(3x))
要这样做:
NUMBER PLUS_OPERATOR NUMBER VAR[X] MINUS_OPERATOR NUMBER MULTIPLY_OPERATOR GROUP
(其中group是另一组令牌)。我如何简化这个方程式呢?是的,我知道你可以做什么:添加X-var,但它们在子组中。我能用来处理这些令牌的最佳方法是什么?