使用 fzero 求多项式的实根

Question

我只想找到方程的真正根；

4*sqrt((1-(z^2/f1^2))*(1-z^2))-(2-z^2)^2-(m*z^4*sqrt(1-z^2/f1^2)/ ...
  sqrt(1-((z^2/f1^2)/y^2)))

我知道方程包含复根，但我不想看到它们。此外，我的代码失败并说；

使用 fzero 时出错（第 242 行）区间端点处的函数值 必须是有限且实数。

scholte 中的错误（第 21 行）x=fzero(fun,x0)

这是我的代码；

rho2 = 1000; %kg/m3    
rho1 = 2700; %kg/m3    
cl2 = 1481; %m/s    
cl1 = 5919; %m/s    
m = rho2/rho1;    
y = cl2/cl1;    
poi = 0.25;

f1 = (sqrt((1-2*poi)/(2*(1-poi))))^-1;

fun = @(z) 4*sqrt((1-(z^2/f1^2))*(1-z^2))-(2-z^2)^2- ...
      (m*z^4*sqrt(1-z^2/f1^2)/sqrt(1-((z^2/f1^2)/y^2)));

x0 = [1 10];

x = fzero(fun, x0)

我多次更改 x0 间隔，但显示相同的错误。如何修复我的代码？

Answer 1

正如 Matlab 告诉你的那样，你的问题是 Function values at interval endpoints must be finite and real，在你的情况下它们不是真实的：

fun(x0(1))

ans =

  -1.0000 + 0.1454i

您的函数可能过于复杂，fzero 无法处理。但是我不是专家，让我们看看是否有比我知识更多的人可以为您指出解决该方程式的正确方向。