这是什么意思,怎么会出现这种情况,求助大佬
mathematica吧
全部回复
仅看楼主
level 2
RootSearch::err: The Private RootSearch function zBent had an error in branch 4 when searching for a root near Ersek`RootSearch`Private`newX. Please email a notebook containing this use of the RootSearch function to RootSearch author at ([email protected]).
这是我的代码
Needs["Ersek`RootSearch`"]
seg = g -> 9.8;
seh = h -> 4000;
Subscript[se\[Rho], 1] = Subscript[\[Rho], 1] -> 1024;(*水的密度*)
Subscript[se\[Rho], e] = Subscript[\[Rho], e] -> 917;
sec = c -> 1500;
seE = capse -> 5*10^9;
sev = v -> 0.3;
seT = T -> 0;
se\[CurlyEpsilon] = \[CurlyEpsilon] -> 1;
Subscript[se\[Lambda], n] =
Subscript[\[Lambda],
n] -> ((\[Omega]^2/c^2) - k^2)^(1/2) /. \[Omega] ->
2*\[Pi]*f /. sec;
seD = D -> (E*d^3)/(12*(1 - v^2)) /. seE /. sev /. d -> 20;
se\[Omega] = \[Omega] -> 2*\[Pi]*f;
se\[CapitalOmega] = \[CapitalOmega] -> \[Omega]*\[CurlyEpsilon] /.
se\[Omega] /. se\[CurlyEpsilon];
y1 =
((Subscript[\[Lambda], n]*h - n*\[Pi] + \[Pi] +
1/3*(Subscript[\[Lambda], n]*h -
n*\[Pi] + \[Pi])^3 +
2/15*(Subscript[\[Lambda], n]*h -
n*\[Pi] + \[Pi])^5)*Subscript[\[Lambda], n]*
Subscript[\[Rho], e]/Subscript[\[Rho], 1]*
d*(1 - (g*Subscript[\[Rho], 1]*\[CurlyEpsilon]^2)/(
d*Subscript[\[Rho], e]*\[CapitalOmega]^2))*(1 + (
D*k^4)/(Subscript[\[Rho], 1]*g) - (T*k^2)/(
Subscript[\[Rho], 1]*g)) - 1) /. seg /.
Subscript[se\[Rho], 1] /. Subscript[se\[Rho], e] /.
seT /. se\[CurlyEpsilon] /. se\[CapitalOmega] /.
Subscript[se\[Lambda], n] /. seD /. se\[CapitalOmega] /.
seh /. n -> 2;
(*\[CurlyEpsilon]->1,n->2,d->20*)
a1220 = Table[y12201 = y1 /. f -> 0.3 + (i - 1)*0.00019 /. d -> 20;
o1220 = RootSearch[y12201 == 0, {k, 0, 0.01}];
{f, k /. o1220}, {i, 10000}];
a1220
Export["root1220.csv", a1220, "CSV"]
2023年07月14日 09点07分 1
吧务
level 15
别盯着警告信息的字面不放。重要的不是警告什么意思,而是什么东西引发了警告。你好好想想以下问题:
1. 这个警告到底是在i等于几的时候冒出来的?
2. 你要怎么确定引发问题的i的取值?
3. 找到这个i之后,你应该用什么办法来分析相应的y12201?
2023年08月05日 03点08分 2
1