2

是否可以在 Matlab 中求解具有“n”个变量（2 级或 1 级）的“n”（非线性）方程组？

你能给我举个例子吗？

我尝试使用fsolve函数，但出现下一个错误：

??? Error using ==> feval Undefined function or method 'ecuaciones' for input arguments of type 'double'.

fsolve 示例

编辑：

我有一个名为“ecuaciones.m”的文件，其中包含以下内容：

function [f]=ecuaciones(x)

f(1)=x(1) + x(2) + x(3) + 2*x(4) -1,905;
f(2)=2*x(3) + 3*x(4) + x(5) + x(6) - 3,922;
f(3)=x(1) + 2*x(2) + x(4) + x(5) + x(6) - 1,961;
f(4)=1,961*x(3) + x(3)*x(5) - x(2)*x(3) - 0,02834*(x(5)^2);
f(5)=9,32845*(x(3)^2) - x(4)*x(6);
f(6)=61,4723*x(1)*x(6) - x(2)*x(5);


J = zeros (n,n);

J =

[ 1, 1, 1, 2, 0, 0]
[ 0, 0, 0, 0, 0, 0]
[ 0, 0, 2, 3, 1, 1]
[ 1, 2, 0, 1, 1, 1]
[ 0, -x3, x5 - x2 + 1961/1000, 0, x3 - (1417*x5)/25000, 0]
[ 0, 0, (2625725246496519*x3)/140737488355328, -x6, 0, -x4]
[ (8651457105425229*x6)/140737488355328, -x5, 0, 0, -x2, (8651457105425229*x1)/140737488355328];

然后，在 Matlab 控制台中，我尝试：

>> x0=[1,1,1,1,1,1];
>> [x,feval,flag]=fsolve('ecuaciones',x0);
4

1 回答 1

2

使用函数句柄

 >> x = fsolve( @ecuaciones, x0 );

您需要确保ecuaciones.m在 matlab 的路径中。

更新：
为了更好地收敛，您可能还想计算雅可比行列式。
在您的情况下，雅可比矩阵是一个nbyn矩阵，它是wrt在该点J(l,k)的偏导数。 您需要修改以返回两个输出，第一个是，第二个是f(l)x(k)x
ecuaciones.mfJ

function [f J] = ecuaciones( x )

f(1)=x(1) + x(2) + x(3) + 2*x(4) -1,905;
f(2)=2*x(3) + 3*x(4) + x(5) + x(6) - 3,922;
f(3)=x(1) + 2*x(2) + x(4) + x(5) + x(6) - 1,961;
f(4)=1,961*x(3) + x(3)*x(5) - x(2)*x(3) - 0,02834*(x(5)^2);
f(5)=9,32845*(x(3)^2) - x(4)*x(6);
f(6)=61,4723*x(1)*x(6) - x(2)*x(5);

J = [...
      1, 1, 1, 2, 0, 0;...
      0, 0, 0, 0, 0, 0;...
      0, 0, 2, 3, 1, 1;...
      1, 2, 0, 1, 1, 1;...
      0, -x3, x5 - x2 + 1961/1000, 0, x3 - (1417*x5)/25000, 0;...
      0, 0, (2625725246496519*x3)/140737488355328, -x6, 0, -x4;...
      (8651457105425229*x6)/140737488355328, -x5, 0, 0, -x2, (8651457105425229*x1)/140737488355328];

你可以打电话fsolve

>> x = fsolve( @ecuaciones, x0, struct('Jacobian','on') );
于 2013-03-17T13:46:37.373 回答