level 2
贴吧用户_GbJUQ5M
楼主
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
这是我的代码
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"]