有关FindRoot的问题,为什么会出现0.i
mathematica吧
全部回复
仅看楼主
level 2
这是结果
{k -> 0.00574815\[VeryThinSpace]+ 0. I}
这是代码
seg = g -> 9.8;
se\[Rho] = \[Rho] -> 1024;
sev = v -> 0.3;
seE = capse -> 5*10^11;
seD = D -> (capse*d^3)/(12*(1 - v^2)) /. seE /. sev;
se\[CapitalGamma] = \[CapitalGamma] -> D/(\[Rho]*g) /. seD /. seg /.
se\[Rho];
Subscript[seM, e] = Subscript[M, e] -> 917*d;
se\[Sigma] = \[Sigma] -> Subscript[M, e]/(\[Rho]*g) /. se\[Rho] /.
seg /. Subscript[seM, e];
seh = h -> 1000;
seQ = Q -> 0;
se\[CapitalLambda] = \[CapitalLambda] -> Q/(\[Rho]*g) /. seQ /.
seg /. se\[Rho];
sec = c -> 1500;
Subscript[sek, 2] =
Subscript[k, 2] -> (\[Omega]^2/c^2 - k^2)^(1/2) /. sec;
sr1 = \[Omega]^2*(1 - \[Sigma]*g*Subscript[k, 2]*
Tan[Subscript[k, 2]*h]) + (1 + \[CapitalGamma]*
k^4 - \[CapitalLambda]*k^2)*g*Subscript[k, 2]*
Tan[Subscript[k, 2]*h] /. seg /. se\[CapitalGamma] /.
se\[CapitalLambda] /. seh /. se\[Sigma] /. sec /.
Subscript[sek, 2] /. d -> 10 /. sev /. se\[Rho] /.
seE /. \[Omega] -> 2*\[Pi]*(0.01 + (i - 1)*0.00002);
FindRoot[(sr1 == 0) /. seg /. se\[CapitalGamma] /.
se\[CapitalLambda] /. seh /. se\[Sigma] /. sec /.
Subscript[sek, 2] /. d -> 10 /. sev /. se\[Rho] /. seE, {k,
0.1}]
2023年03月12日 11点03分 1
吧务
level 15
你i的定义呢?
2023年04月01日 10点04分 2
1