在MATLAB中用多个输入求解fsolve非线性系统

Question

我正在尝试使用fsolve求解由六个非线性方程组成的系统（请参见下文）。我的方程式中还有其他一些参数，我希望能够输入到系统中或更改功能，而不必每次都进入我的函数并手动对其进行编辑。

它们是F，rho和A。现在，我将它们全部设置为1，并且我能够找到解决方案，因此至少脚本可以正确运行。

function f = pressXmanifold(x)
    F = [1,1,1];
    rho = 1;
    A = 1;
    f(1) = x(1)-x(3) - ...
        1/2*rho*(x(4)/F(1))^2*...
        A*(x(4)/x(6)*F(3)/F(1))^(-2)*...
        (1+(F(3)/F(1))^2+3*(F(3)/F(1))^2*((x(4)/x(6))^2-(x(4)/x(6))));
    f(2) = x(2)-x(3) - ...
        1/2*rho*(x(5)/F(2))^2*...
        A*(x(5)/x(6)*F(3)/F(2))^(-2)*...
        (1+(F(3)/F(2))^2+3*(F(3)/F(2))^2*((x(5)/x(6))^2-(x(5)/x(6))));
    f(3) = x(1) - x(3);
    f(4) = x(2) - x(3);
    f(5) = x(4) - x(5);
    f(6) = x(4)+x(5) - x(6);
end

对于F，我希望能够将任意1 x 3数组输入到函数中。 对于rho，我也想输入一个值。 对于A，它必须根据以下条件（以if语句给出）随迭代而改变。

if F(3)/F(1) <= 0.35 && x(3)/x(2) <= 1
    A = 1;
elseif F(3)/F(1) > 0.35 && x(3)/x(2) <= 0.4
    A = 0.9*(1-x(3)/x(2));
elseif F(3)/F(1) > 0.35 && x(3)/x(2) > 0.4
    A = 0.55;
end

目前，我仅尝试处理前两个参数。 f = pressXmanifold(x,F,rho)并删除函数中的前两行。 但是当我尝试按照MathWorks网站上的说明进行操作

options = optimoptions('fsolve','Display','none','PlotFcn',@optimplotfirstorderopt);
fun = @pressXmanifold;
x0 = [1,1,1,1,1,1];
x = fsolve(fun,x0,F,rho,options)

无法执行赋值，因为这种类型的变量不支持点索引。 createOptionFeedback中的错误（第33行）             选项。（stopTestOptions {k}）= []; prepareOptionsForSolver中的错误（第57行）     optionFeedback = createOptionFeedback（options）; fsolve错误（第157行） [options，optionFeedback] = prepareOptionsForSolver（options，'fsolve'）;

Answer 1

以F和rho作为输入重写pressXmanifold

function f = pressXmanifold(x1, x2, x3, x4, x5, x6, F, rho)
   if F(3)/F(1) <= 0.35 && x3/x2 <= 1
        A = 1;
    elseif F(3)/F(1) > 0.35 && x3/x2 <= 0.4
        A = 0.9*(1-x3/x2);
    elseif F(3)/F(1) > 0.35 && x3/x2 > 0.4
        A = 0.55;
    end
    f(1) = x1-x3 - ...
        1/2*rho*(x4/F(1))^2*...
        A*(x4/x6*F(3)/F(1))^(-2)*...
        (1+(F(3)/F(1))^2+3*(F(3)/F(1))^2*((x4/x6)^2-(x4/x6)));
    f(2) = x2-x3 - ...
        1/2*rho*(x5/F(2))^2*...
        A*(x5/x6*F(3)/F(2))^(-2)*...
        (1+(F(3)/F(2))^2+3*(F(3)/F(2))^2*((x5/x6)^2-(x5/x6)));
    f(3) = x1 - x3;
    f(4) = x2 - x3;
    f(5) = x4 - x5;
    f(6) = x4+x5 - x6;
end

然后定义F和rho

rho = 1;
F =[1,1,1];
fun = @(x)pressXmanifold(x(1), x(2), x(3), x(4), x(5), x(6), F, rho);

最后

options = optimoptions('fsolve','Display','none','PlotFcn',...
    @optimplotfirstorderopt);
x0 = [1,1,1,1,1,1];
x = fsolve(fun,x0,options)