如何在fsolve中确定解向量的分量？
在使用fsolve函数确定解向量的分量时，需要遵循以下步骤：
- 定义目标函数：首先，需要定义一个目标函数，该函数描述了问题的约束条件和目标。目标函数的输入参数是解向量，输出是一个向量，表示约束条件的满足程度。
- 初始化解向量：根据问题的具体要求，初始化解向量的初始值。解向量的维度应与目标函数的输入参数一致。
- 调用fsolve函数：使用fsolve函数来求解目标函数的解向量。fsolve函数是一个数值求解器，它通过迭代的方式寻找目标函数的零点，即满足约束条件的解向量。
- 设置约束条件：如果问题存在约束条件，可以通过在目标函数中引入约束条件来限制解向量的取值范围。常见的约束条件包括等式约束和不等式约束。
- 获取解向量的分量：在求解完成后，可以通过访问解向量的各个分量来获取最终的解。根据问题的要求，可能只需要解向量中的某些分量。
需要注意的是，fsolve函数是一种数值求解方法，对于复杂的问题可能存在局部最优解的情况。因此，在使用fsolve函数时，需要根据具体问题的特点进行调参和优化，以获得更准确的解。
腾讯云提供了一系列与云计算相关的产品和服务，包括云服务器、云数据库、云存储等。具体推荐的产品和产品介绍链接地址可以根据具体问题的需求和场景来选择。
相关·内容
2回答
我遇到了这个错误：
在赋值A(I) = B中，B和I中的元素数量必须相同。
yres(1)=((u - uc).^2) + ((y - yc).^2) -(d.^2)；我不明白，为什么这不会得到一个skalar?因为元素都是标量。应该如何更改才能获得标量？
诚挚的问候
编辑:感谢sloede，所有的输入都是标量的，但是我仍然得到这个错误。
在赋值A(I) = B中，B和I中的元素数量必须相同。
Error in myfun (line 7)
yres(1)=sqrt(((u - uc).^2) + ((y - yc).^2) ) -d;
Error in fsolv
浏览 0提问于2012-07-06得票数 1
回答已采纳
为什么下面的代码返回一个ValueError？
from scipy.optimize import fsolve
import numpy as np
def f(p,a=0):
x,y = p
return (np.dot(x,y)-a,np.outer(x,y)-np.ones((3,3)),x+y-np.array([1,2,3]))
x,y = fsolve(f,(np.ones(3),np.ones(3)),9)
ValueError: setting an array element with a sequence.
浏览 1提问于2015-03-03得票数 1
回答已采纳
2回答
我正在尝试将fsolve应用于一个数组：
from __future__ import division
from math import fsum
from numpy import *
from scipy.optimize import fsolve
from scipy.constants import pi
nu = 0.05
cn = [0]
cn.extend([pi*n - pi/4 for n in range(1, 5 +1)])
b = linspace(200, 600, 400)
a = fsolve(lambda a: 1/b + fsum([1/(b+cn[i]
浏览 4提问于2012-05-31得票数 1
回答已采纳
1回答
我有一个由8个方程和8个未知数组成的系统。在这组方程中没有任何数字。方程式中只有参数；其中8个是我的未知变量(以下称为未知变量)，我正在求解这些变量，还有一堆其他参数，它们不是未知值，例如alpha，beta，gamma (以下称为参数)。
我需要解方程组，这样我就可以得到参数方面的8个未知变量。
我在Maple中使用了fsolve(方程，变量)。但是，它返回一个错误，显示tau，alpha，beta，...都在方程式中，并且没有解出来。我如何才能让Maple理解这些不是我喜欢求解的变量，而是其他未知变量应该用来求解的参数呢？
下面是我的一个方程式的例子：
e1 := (L__af+L__mf
浏览 0提问于2015-09-12得票数 0
我的代码有问题。上面写着：
Attempted to access F(2); index out of bounds because numel(F)=1.
我的代码是：
function Punto3
global p n E d Vm rugo
e=0.001:0.02:1;
E=e*10^-3;
d=0.01;
Vm= 10 ;
p= 998.3;
n=1.002*10^-3;
rugo=length(E);
xo=5;
f=fsolve(@(F)ecuacion(F),xo);
end
functi
浏览 2提问于2016-08-07得票数 0
回答已采纳
4回答
当你求解一个方程组时，如何传递参数来求解？下面是一个例子：
from scipy.optimize import fsolve
def equations(x,y,a,b):
return(x+y+a+b,x-y-a-b)
x,y = fsolve(equations,[0,0],args=(1,2))
我得到：
TypeError:方程()只带4个参数(3个给定)
在这种情况下如何传递参数？
浏览 2提问于2019-02-25得票数 0
回答已采纳
2回答
我想通过在python代码中转换数学来找到python中的斜率。在数学上，斜率公式是-
m = Σ(x-x̄)*(y -ȳ)/ Σ(x-x)^2
但在python中，它是通过-
m = (np.sum(x*y) - n*x_mean * y_mean) / (np.sum(x*x) - (n * x_mean * x_mean))
我不确定我们是如何在python中导出这个公式的。你能帮帮我吗？
为了更好地理解，我在图片中附加了作为Python公式的数学公式。请也检查一下。
谢谢
浏览 3提问于2020-02-05得票数 0
回答已采纳
比如gl_FragColor = v1 * v2，我真的搞不懂它是如何乘积的，它似乎给出了向量乘矩阵的解释。
ps：v1和v2的类型都是vec4。
浏览 2提问于2012-12-16得票数 33
回答已采纳
我在matlab中有一个简单的问题。我有一个等式：，我知道a和b，我试着用这个表达式：，A*H=b，a，b，I，a，b，我试着使用这个表达式：
H=A\b;
但我得到了错误的值:示例：
A =
231 481
233 488
241 481
243 489
b =
11 31
6 20
21 31
18 22
我得到了
H =
1.1627 0.2713
-0.5396 -0.0791
所以
A*H
ans =
浏览 5提问于2012-08-20得票数 0
回答已采纳
假设我想要解一个含有N未知变量的N方程组。假设方程是以下一般形式(我已经大大简化了它，我们可以假设每一行等于零)。
x1 - const(1)
x2 - const(2)
x3 - const(3)
...
xN - const(N)
其中x1、x2、.、xN是我的变量，const是代码前面确定的常量的长度N向量。我不知道有多少方程和变量(即不能硬代码)，但我仍然希望用一般的方式来编写和求解这个系统。
在MATLAB中，我目前的解决方案是执行以下操作，其中n_vars是变量数，我的程序在前面确定了变量数。
sym_vars = sym('x',[1 n_vars]);
for
浏览 1提问于2017-04-07得票数 0
回答已采纳
1回答
我试图使用fsolve来寻找各种参数分段定义函数的交点：
max_price = 20.0
b = 0.5
consumption_func = lambda x,b : max_price if x <= b else max_price - (x-b)
func = lambda x: consumption_func(x,b)
bid_start, _, ier, msg = fsolve(lambda x: func(x) - 5.0, 0.0, xtol=0.0001, full_output=True)
但是，如果开始点位于函数常量的范围内，则求解程序无法找到这个非常简单的
浏览 2提问于2016-12-09得票数 1
回答已采纳
2回答
我正试图解决以下几个问题：
function [eta, sol] = compressible_similarity_wo
global Gamm Ma Pr omega;
Gamm = 1.4;
Ma = 2;
Pr = 0.7;
omega=0.76;
global eta_max_ode;
eta_max_ode = 20;
opt = optimset('Display','off','TolFun',1E-20);
F = fsolve(@(F) eva
浏览 4提问于2020-02-06得票数 1
回答已采纳
1回答
假设我有一个三维模型，它代表一个很小的真实世界空间，它可以被一个10“乘10”的盒子所包围。我知道现实世界的指南针方位，在盒子的中心向下看X轴。我怎么估计这个空间中其他矢量的方位呢？
浏览 0提问于2022-12-09得票数 2
我在3D空间中有两点A和B (比如Point3D类型)。B需要向A行驶，然后在到达或通过B时停下来。两个点之间的方向是已知的(可以很容易地计算)，所以速度很容易设定。如何确定点A何时到达或“通过”(给定方向)点B？
浏览 0提问于2014-10-16得票数 0
回答已采纳
2回答
运算A\B的结果是什么，其中A(1，m)和B (1，m)？
在手册中写道：
A\B returns a least-squares solution to the system of equations A*x= B.
所以它意味着x= inv (A'*A)*A'*B？然而，矩阵A'*A是奇异的.
让我们假设：
A=[1 2 3]
B=[6 7 6]
A\B
0 0 0
0 0 0
2.0000 2.3333 2.0000
如果我们使用MLS：
C = inv (A'*A) sin
浏览 0提问于2012-11-14得票数 9
回答已采纳
我有一个函数，当我给它一个输入进行测试时，它似乎工作正常，但我需要它在许多不同的变量上工作。现在在底部，它正在评估函数，初始猜测为25，指数为50。在使用它的过程中，我注意到索引实际上是随着fsolve的运行而变化的，所以我认为我对fsolve在这里的实际运行方式有一个不完整的理解，但我无法修复它。我的最终目标是能够在for循环中将索引从1循环到54，在该循环中，初始猜测将根据我已经运行的算法进行更改。所以我只需要解决这个函数输入问题。
下面是我的代码：
function [objfun] = RK_solver( RK_solver_input )
% Redlich/Kwong
浏览 0提问于2016-03-16得票数 0
我在一个m文件中有一个fsolve，它在2011b版本中工作得很好。然而，fsolve在2012a每次都会失败。会导致这种情况的功能或选项是否发生了重大变化？
浏览 3提问于2012-03-22得票数 3
回答已采纳
1回答
我想要解这个方程：
然而，当我试图用scipy fsolve解决它时，它收敛到无穷远，而不是给我解决方案。
它趋于无穷大的原因是当x趋于无穷大时，函数趋于0：
下面是示例代码：
def f(x, r): return -e ** (-r * x)
def h(r): return 2 * f(4, r) - f(2, r) - f(10, r)
x0 = np.array([1])
print(optimize.fsolve(h, x0))
通过一些其他参数，它找到了解决方案。但是，我希望代码使用不同的参数，而不仅仅是示例中的参数。我也想避免零解。
非常感谢
浏览 0提问于2019-07-25得票数 0
1回答
例如，我定义了以下函数句柄：
F = @(x, y, z)[x+y+z; x*y*z];
funcc = @(x, y)F(x, y, 0);
打电话
res = fsolve(funcc, [10; 10]);
导致错误：
Error using @(x,y)F(x,y,0)
Not enough input arguments.
Error in fsolve (line 219)
fuser = feval(funfcn{3},x,varargin{:});
Caused by:
Failure in initial user-supplied objec
浏览 3提问于2016-05-28得票数 1
回答已采纳
我有以下功能：
Eq = @(x1, x2) [x1-6, x2+3];
fsolve(Eq, [4 1])
但会得到以下错误：
??? Input argument "x2" is undefined.
Error in ==> @(x1,x2)[x1-6,x2+3]
Error in ==> fsolve at 193
fuser = feval(funfcn{3},x,varargin{:});
Error in ==> Untitled at 6
fsolve(Eq, [4, 1])
当我将函数更改为单输入函数时，它可以完美地工作。有人
浏览 0提问于2010-08-27得票数 5
回答已采纳
我相信这个问题已经被解决了好几次，但不幸的是，提出的方法都没有对我的情况有所帮助。我正在尝试使用Scipy fsolve解决一个复杂的研究，但求解器对初始条件非常敏感。
我的问题是:有没有其他的解决方案？如果没有，有没有办法调整灵敏度？
浏览 0提问于2019-08-13得票数 1
1回答
解复矩阵方程
、、、
我有一个涉及矩阵的复杂方程式：
R = expm(X)*A + (expm(X)-I)*inv(X)*B*U;
其中R、B和U是已知矩阵。
I是一个单位矩阵。
我需要解决X的问题。有没有办法在MATLAB中解决这个问题？
浏览 0提问于2019-04-28得票数 1
1回答
我有一个of3非线性方程组，我需要用matlab中的fsolve来求解它们。
function F = root2d(y)
syms b1 b2 b3 w21 w31 theta1 theta2 theta3 a1 a2 a3;
F(1) = (1+exp(-b1*(w21*y(2)+w31*y(3)-theta1)))^(-1) - a1*y(1);
F(2) = (1/(1+exp(-b2*(y(1)-theta2))))-a2*y(2);
F(3)= (1/(1+exp(-b3*(y(1)-theta3))))-a3*y(3);
这是我的函数matlab文件。我在使用ma
浏览 2提问于2017-12-10得票数 0
回答已采纳
我正在尝试为一个由隐式函数定义的模型找到优化的参数，以使用fsolve和lsqcurvefit来拟合数据集。我已经在单独的m文件中定义了3个函数:第一个是在要定义的4个参数中定义隐式函数，第二个是使用fsolve来查找所定义的隐式函数的根，第三个是使用lsqcurvefit来查找四个参数的最佳值。当然，我需要为参数定义足够好的初始值，但是在尝试了各种合理的组合之后，lsqcurvefit总是运行大约20-30次迭代(matlab在每次迭代后打印出用fsolve找到的解计算出的矢量值)，然后打印
No solution found.
fsolve stopped because the pro
浏览 2提问于2010-11-29得票数 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
我想解决Y(x):的问题
Y(x) = A - B*x + C*ln(A/Y(x))
where A, B, and C are defined
在这两个Python or mathematically中，我都不知道该怎么做。我想尝试使用fsolve或其他工具来查找根，但我不认为我可以将x and Y(x)视为彼此独立。
请帮帮忙。
浏览 3提问于2013-12-10得票数 0
2回答
问题
我有一组用小写变量表示的变量和用大写变量表示的常量。
A = a + b
B = c + d
C = a + b + c + d + e
我得到了关于熊猫DataFrame中这些方程的结构的信息，它有两列：常数和变量
例如。
df = pd.DataFrame([['A','a'],['A','b'],['B','c'],['B','d'],['C','a'],['C','b'],
[
浏览 0提问于2018-06-13得票数 3
回答已采纳
我在这里的目标是改善用户体验，这样光标就会到达用户直觉期望的位置，无论这意味着什么。
考虑一个操纵杆对每个方向都有不同的配置速度。
也许操纵杆有一个缺点，一些方向太敏感，而另一些方向不够敏感，所以你试图纠正这一点。或者你在玩FPS游戏，你很少需要向上或向下看，所以你降低了Y的敏感度。
以下是我们每个方向的最高速度：
var map = {
x: 100,
y: 200,
}
操纵杆输入给我们一个从0到1的单位向量。
现在，操纵杆向右倾斜，25%的方向倾斜，50%的方向倾斜。
joystick = (dx: 0.25, dy: -0.50)
不好意思，我不知道从这里往哪里走。
浏览 3提问于2022-07-15得票数 0
回答已采纳
3回答
在我得到的帮助之后，我一直试图在脚本中实现它，但我没有聪明地运行它。
我需要对4072x3080图像的每一个像素使用这个算法，整个过程大约需要1h30，所以我试图强迫它，但是我得到了这个错误：
ValueError Traceback (most recent call last)
<ipython-input-12-99c1f41dbba7> in <module>()
----> 1 res = scipy.optimize.fsolve(func, x0=np.ones((K.shape[0],
浏览 3提问于2019-01-20得票数 2
回答已采纳
我试图在2.5排(在UDK)中制作一个角色，在环绕树的螺旋楼梯上“爬”一棵巨大的树干。当角色向右移动时，树会旋转穿过日场序列，因此，角色似乎在移动，而在现实中，则是树在移动。
我把日场序列播放率和角色的速度连接起来，只要角色向左或向右移动就行了。然而，当它跳跃时，速度仍然影响着播放率--它不应该随着角色的上下移动，而不是左右移动。
我如何在Kismet中设置它，以便得到一个只在X平面(水平)上有速度的浮子变量？
浏览 0提问于2012-10-04得票数 1
回答已采纳
1回答
如果你能帮我弄清楚圣人的代码。我想找出这类方程的根：e^(2*pi*b*i)。有一个常数b和一个虚构的单位i。这个等式也可以写成：cos(b*2*pi)+ i*sin(b*2*pi)。你知道做这件事的命令吗？ 我正在尝试： import matplotlib.pyplot as plt
import numpy as np
import cmath
from scipy.optimize import fsolve
z=var('z')
a = cos(4*cmath.pi*z)+ cos(14*cmath.pi*z) + cos(62*cmath.pi*z)+ (I*
(
浏览 13提问于2020-01-24得票数 0
1回答
我有一个非线性方程:我在下图中画出了这个方程的左边和右边。显然，有两种解决方案。我已经尝试了scipy.optimize.fsolve，并战略性地猜测了起点，以找到两个解决方案，但这是相当不健壮的。有没有一个Python求解器，可以找到这个方程的所有解，例如，在某个有限的x值集合中？
浏览 5提问于2021-11-10得票数 0
1回答
我有一个函数f(v,u)，我定义了一个函数
solutionf(u) := fsolve(f(v,u)=v);
我需要根据solutionf(u)绘制u，但只是
plot(solutionf(u), u = 0 .. 0.4e-1)
给我一个错误
Error, (in fsolve) number of equations, 1, does not match number of variables, 2
但是，我总是可以在任何solutionf(x)上取值x。
有简单的方法来策划这个吗？或者我必须在u上创建自己的u循环，在每个点取值，并绘制插入值？
浏览 1提问于2016-02-11得票数 1
回答已采纳
我想把墙和地板分开。现在，我有一个深度图像，并在深度图像中找到表面法线。我想要分组同一向量的法线和转换图像为二值。
下面的代码是我如何找到表面法线的。
d_im = depth_image.astype("float64")
zy, zx = np.gradient(d_im)
normal = np.dstack((-zx, -zy, np.ones_like(d_im)))
n = np.linalg.norm(normal, axis=2)
normal[:, :, 0] /= n
normal[:, :
浏览 1提问于2018-11-29得票数 2
回答已采纳
1回答
我使用下面的代码来生成一个具有6个边界条件的ODE系统的数值解。
我使用初始条件来获得一个解，但是为了找到真正的解，我必须改变另外三个条件。我的职能如下：
function diff = prob5diff(M,Fx,Fy)
u0 = [pi/2 0 0]';
sSpan = [0 13];
p = @(t,u) prob5(t,u,M,Fx,Fy);
options = odeset('reltol',1e-6,'abstol',1e-6);
[s,u] = ode45(p,sSpan,u0,options);
L = length(s);
x =
浏览 2提问于2014-05-13得票数 0
回答已采纳
在网上找不到这个。使用tessalation评估着色器，我可以为vec2或vec3插值如下所示的值：
vec3 interpolate3D(vec3 v0, vec3 v1, vec3 v2)
{
return vec3(gl_TessCoord.x) * v0 + vec3(gl_TessCoord.y) * v1 +vec3(gl_TessCoord.z) * v2;
}
我如何做同样的传递和插值矩阵？对于法线映射的切线矩阵和变换矩阵，要从高度图中计算出新的顶点法线。谢谢
浏览 0提问于2019-09-10得票数 1
回答已采纳
1回答
我用C++语言编写了求解非调和势V= x^2/2 + lambda*x^4的一维薛定谔方程的程序。我的代码正在工作，我用Mathematica动画化了结果，以检查正在发生的事情。我根据已知的谐波势解测试代码(我把lambda = 0)，但是动画显示abs(Psi)是随时间变化的，我知道这对谐波势是不正确的。实际上，我看到，在某一点上，它变成常数，但在此之前，它是振荡的。
所以我知道，在时间间隔内，我需要波函数的恒定大小，但我不知道怎么做，也不知道我在哪里做错了。这是我的代码和动画的100个时间步骤和网格上的100个点。
#include <iostream>
#include &l
浏览 3提问于2018-01-09得票数 0
1回答
我有一些关于nsolve的问题，我很难找到一些函数的解决方案，给出了一些初步的猜测。然后我想尝试一下numpy/scipy解算器。
下面是一个使用渐近的程序，可以很好地解决这个问题：[0.0, -9.05567e-72, 9.42477, 3.14159]
from sympy import *
# Symbols
theta = Symbol('theta')
phi = Symbol('phi')
phi0 = Symbol('phi0')
H0 = Symbol('H0')
# Constants
phi0 = 60*pi
浏览 0提问于2012-11-14得票数 2
回答已采纳
3回答
我正在寻找一种方法来找到与游戏实体的方向成直角的向量。一个在左边，一个在右边。
如果这会以任何方式影响答案，我将使用XNA。
编辑:这是一个二维操作。我在另一个网站上看到顺时针向量是简单的-y，x和逆时针的y，-x。这似乎是纸上谈兵。
谢谢。
浏览 0提问于2011-05-24得票数 1
回答已采纳
我想开发一个用于解线性和非线性方程组的桌面应用程序。我正在考虑在Matlab中使用C#。我以前从未使用过matlab，但我正在努力为我的项目学习。
function F = myfun(x)
F = [2*x(1) - x(2) - exp(-x(1));
-x(1) + 2*x(2) - exp(-x(2))];
**x0 = [-5; -5]; % Make a starting guess at the solution**
options=optimset('Display','iter'); % Option to
浏览 1提问于2011-10-12得票数 1
我想解决一个函数f(x)=0，这样g(x)>=0。有没有可能通过matlab来实现呢？我对f(x)=0使用fsolve/fzero函数，但有时会发生g(x)<0不合适的情况。我知道我在优化中设置了约束，但当我简单地求解零时，我不确定是否可以做同样的事情。我寻找，但没有人关注它。
浏览 0提问于2019-06-22得票数 0
我已经创建了一个函数Euler.m来使用Euler方法求解一个常微分方程系统。我希望使用此函数来解决由匿名函数func=@(t) ([x(t)+4*y(t)-exp(t);x(t)+y(t)+2*exp(t)])定义的ODE系统，该函数具有y0给出的初始条件。 func=@(t) ([x(t)+4*y(t)-exp(t);x(t)+y(t)+2*exp(t)]);
y0=[4;5/4];
y_exact=@(t) [4*exp(3*t)+2*exp(-t)-2*exp(t);2*exp(3*t)-exp(-t)+exp(t)/4]; %exact solution of ODEs
a=0;
浏览 37提问于2020-04-21得票数 2
回答已采纳
我试图使用下面的代码作为指南，以解决一个非线性方程，但我继续得到错误，如“对象太深的目标数组”和“从函数调用的结果不是一个适当的浮点数组”。
from scipy.optimize import fsolve
from math import exp
def equations(vars):
x, y = vars
eq1 = x+y**2-4
eq2 = exp(x) + x*y - 3
return [eq1, eq2]
x, y = fsolve(equations, (1, 1))
print(x, y)
我的代码将发布在下面。它指出了"
浏览 0提问于2021-10-29得票数 0
1回答
我在scipy.fsolve中遇到了以下问题，但我不知道该怎么做：
U = 0.00043
ThC =1.19
Dist = 7
IncT = 0.2
pcw = 1180000
k = 1.19
B = U * pcw / (2 * k)
fugato = fsolve((((Ql/(2*math.pi* k))*math.exp(B * x)*special.kv(0, B * x))-IncT),0.01)
print fugato
我得到了错误的TypeError: 'numpy.float64'
浏览 0提问于2014-11-26得票数 1
回答已采纳
1回答
我用fsolve函数对两个方程中的两个未知数进行求解。最后，我试图解决的方程要复杂得多，但我已经很难理解下面的基本例子。
import scipy.optimize as scopt
def fun(variables) :
(x,y) = variables
eqn_1 = x ** 2 + y - 4
eqn_2 = x + y ** 2 +3
return [eqn_1, eqn_2]
result = scopt.fsolve(fun, (0.1, 1))
print(result)
这给出了结果[-2.08470396 -0.12127194]，
浏览 8提问于2022-05-27得票数 1
回答已采纳
1回答
这是我关于stackoverflow的第一个问题。抱歉，如果它有点笨拙。 我正在尝试解决描述层状土壤中的梁的静态平衡。就作用在梁上的力而言，有一些外力，还有土对外力的反作用力。显然，行动和反应必须平衡。反作用力是来自不同土层的力的总和。土的响应是位移的非线性函数。基本上，这个特殊的平衡是一个标量，因为它就是反应=作用的总和。 另一方面，梁内的内力也必须与作用在所讨论的梁单元上的局部力相平衡。假设我们有40个这样的元素。因此，我面临的问题概括如下： 我需要找到F1(x1，x2，...x40)的根，用x1，x2...x40描述桩在土壤中的位移。我传递并初步猜测了位移，Fsolve返回了一个向量长
浏览 43提问于2020-05-03得票数 0
回答已采纳
我有16个相互独立的非线性方程组，即它们不是一个方程组。一种方法是创建16个独立的子例程，并使用which来解决我通常做的事情。但我需要把子程序的数量从16个减少到1个。让我举一个例子来说明我目前正在做的事情：
u01 = .001;....u016 = .001;
options=optimset('Display','notify','MaxFunEvals',10^7,'TolX',1e-,'TolFun',1e-6,'MaxIter',10^5);
u1 = fsolve(@polsim1,
浏览 3提问于2017-07-10得票数 0
回答已采纳
我试着用解析法和狗腿法来求解一个二阶非线性方程组。我的目标函数和它的雅可比矩阵是这样的
function [F jacF]= objective(x)
F(:,1) = ((((x(:,2)./10).*k).*(x(:,1)./100)).^2).*(rZ - Rs) +(( Cmax .* ( x(:,1)./100 ) ).^2).*( w.^2.*(rZ - Rs) ) - (((x(:,2)./10).*k).*(x(:,1)./100));
F(:,2) = (x(:,2).*k).^2.*(iZ - w.*Ls) + (x(:,2).*k).^2.*x(:,
浏览 0提问于2019-09-12得票数 0
我有一个一般方程
t=tr+(ts-tr)/(1+(a*h)^n)^(1-1/n)
对于(h=0，1，2，3)，我有t=2.000，1.6300，1.2311，1.1084。因此，有4个方程有4个未知数tr，ts，a，n。
我在matlab中使用“求解”函数。
s=solve('tr+(ts-tr)/(1+(a*0)^n)^(1-1/n)=2','tr+(ts-tr)/(1+(a*1)^n)^(1-1/n)=1.63','tr+(ts-tr)/(1+(a*2)^n)^(1-1/n)=1.2311','tr+(ts-tr)/(1+(a*3)
浏览 2提问于2013-11-14得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
即时通信 IM活动推荐
腾讯云开发者
