关于fsolve求解范围的问题

qian211111 · 发表于 2023-1-13 10:33

您需要登录才可以下载或查看，没有帐号？注册

x

目标：
求解一个如下非线性方程组的解，并想将该解限定在纯虚数的范围内；
示例代码如下：

funHd=@(Zvec)func(Zvec);
x0=[13.7472i 13.7472i 0.0020i];
% options = optimoptions('fsolve','Display','none','PlotFcn',@optimplotfirstorderopt);
[zV,fv] = fsolve(funHd,x0)
function f=func(Zvec)
f=[
Zvec(3)^2 - ((- Zvec(3)^2 + Zvec(1)*Zvec(2) + 1)^2 - 4*Zvec(1)*Zvec(2))^(1/2)/(2*Zvec(2)) - Zvec(1)*Zvec(2) - (1 + 6632178539682017i/17592186044416)
Zvec(3)^2 - ((- Zvec(3)^2 + Zvec(1)*Zvec(2) + 1)^2 - (Zvec(1)/2 + Zvec(2)/2 + Zvec(3))*(2*Zvec(1) + 2*Zvec(2) - 4*Zvec(3)))^(1/2)/(Zvec(1) + Zvec(2) - 2*Zvec(3)) - Zvec(1)*Zvec(2) - (1 + 6632160806443741i/17592186044416)
Zvec(3)^2 - ((- Zvec(3)^2 + Zvec(1)*Zvec(2) + 1)^2 - 4*Zvec(1)*Zvec(2))^(1/2)/(2*Zvec(1)) - Zvec(1)*Zvec(2) - (1 + 6632178539682017i/17592186044416)
];
end2 E( X+ E' ?' o