清空筛选条件
清空筛选条件

fsolve

7 次查看（过去 30 天）
gianluca
gianluca 2011-12-12
Hi,
I'm solving a system of two non linear equations involving the same number of variables. I'm using the command fsolve, so before I've defined my functions:
function F = eqns(x)
F = [x(1).*A + (1-x(1)).*x(2).*B + (1-x(1)).*(1-x(2)).*C - b1;
x(1).*D + (1-x(1)).*x(2).*E + (1-x(1)).*(1-x(2)).*F - b2]
and save it as nle.m. A,B,C,D,E,F are constant values and b1 and b2 are my own data stored in two vectors in the workspace.
The command
x0 = [0 0];
[x,toll] = fsolve(@nle,x0)
solve my system only if I write the values of b1 and b2 in the equations defined into the m-file. Since I've thousand of b1 and b2 pairs, I would that fsolve find the zero capturing the current values of b1 and b2 from the workspace. How can I do this?
Thanks for the help

请先登录，再回答此问题。

采纳的回答

Titus Edelhofer
Titus Edelhofer 2011-12-12
Hi,
add the b1, b2 to the function as variables:
function F = eqns(x, b1, b2)
F = ...
and call pass your b1, b2 to the call for fsolve:
b1 = 42;
b2 = 1;
fun = @(x) nle(x, b1, b2);
[x,toll] = fsolve(fun, x0);
Titus
  5 个评论
显示 3更早的评论
gianluca
gianluca 2011-12-12
Hi,
this is the nle.m file:
function F = eqns(x,b1,b2)
F = [x(1).*0+(1-x(1)).*x(2).*30+(1-x(1)).*(1-x(2)).*150-b1;
x(1).*189+(1-x(1)).*x(2).*55.5+(1-x(1)).*(1-x(2)).*89.9-b2];
In the command windows:
b1 = CGR1;
b2 = S;
fun = @(x) nle(x, b1, b2);
[x,toll] = fsolve(fun, x0);
where CGR1 and S are the vector of size (28721 x 1).
I attach yuo a little part of the CGR1 and S vectors (e.g. 10 rows):
CGR1 = [37.6 38.6 39.7 40.7 41.8 42.8 43.9 44.9 46.0 47.1]
S = [67.4 67.7 67.9 68.1 68.3 68.5 68.7 69.0 69.2 77.3]
Titus Edelhofer
Titus Edelhofer 2011-12-14
Hi,
you need to put the call to fsolve into the loop as I did above.
Titus

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by