|
|
这是针对实数的情况，从这个结果看，看起来象是第3个条件不好满足
第3个约束要是 =0 可能看起来更象些 .(取alpha=2情况，前面在复数区域也搜索过，好不到哪里去)
当然也可能是程序搜索的随机点还不够多(一般应该不会吧)
Clear["Global`*"];
V = \[Theta] +
2 \[Alpha]*
Log[Sqrt[y + b/a] \[Gamma]*
BesselJ[1/3, ((y + b/a)^2)^(1/3) 2/3 Sqrt[-a]]];
eq1 = (V /. y -> -1);
eq2 = (V /. y -> 1) - 1;
eq3 = (D[V, {y, 2}] /. y -> -1) - 1;
eq4 = (D[V, {y, 3}] /. y -> 1);
eq = Abs[eq1] + Abs[eq2] + 1*Abs[eq3] + 1*Abs[eq4];
r1 = NMinimize[eq /. \[Alpha] -> 2, {\[Theta], a, b, \[Gamma]},
Method -> {"RandomSearch", "SearchPoints" -> 200}]
(*r2=NMinimize[eq/.\[Alpha]->2,{\[Theta],a,b,\[Gamma]},Method->{\
"RandomSearch",Method->"InteriorPoint", "SearchPoints"->20}]*)
eq1 /. r1[[2]] /. \[Alpha] -> 2
eq2 /. r1[[2]] /. \[Alpha] -> 2
eq3 /. r1[[2]] /. \[Alpha] -> 2
eq4 /. r1[[2]] /. \[Alpha] -> 2 |
|
回复此楼
机械专硕调剂
+3