关于fsolve求解范围的问题

qian211111

目标：
. y2 L% o: R, B; a  j: \/ U9 d2 Y: G求解一个如下非线性方程组的解，并想将该解限定在纯虚数的范围内；
* `! t: B  |2 b8 h示例代码如下：
5 C3 z3 d8 z+ F0 Z' ]

• 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)
• ];
• end
  

复制代码
3 J1 n+ h5 X: g2 `1 ?8 i  R+ C  q+ f) h* a请问有无限定解范围为纯虚数的的限定设置？
, h' \  n1 _1 m7 N# G

like_22you

带约束的非线性方程组可能对你会有帮助。

qian211111

暂时想到的方法是取变量实部的平方附加到求解函数上去，即原目标函数变为f=abs(f)+addCondtion^2;
0 B' G* _0 Y. [& m, _! G暂未找到其他方法；
0 i! q+ _2 G5 S但是这样带来的问题是，目标函数变复杂，fsolve求解结果不收敛，残差很大；
! I$ Q6 L9 R, {' V暂未解决该问题。

qian211111

like_22you 发表于 2023-1-13 11:20
0 u' n, L$ ]8 M带约束的非线性方程组可能对你会有帮助。

, J7 j3 f# K7 i  ]( D' K很有帮助，多谢；但是还是存在残差很大的问题。