| 查看: 1067 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
boundary condition 求 free constant
|
|||
|
我不知道这个能不能用Mathematica来求，但是我用手算，太麻烦了。我过几天再试试，但是如果谁知道用mathematica来求的话，最好了。 我有个式子是 [ Last edited by wocamuchong on 2014-2-17 at 07:17 ] |
回复此楼» 猜你喜欢
0703化学336分求调剂
已经有6人回复
-
304求调剂
已经有9人回复
-
293求调剂
已经有13人回复
-
281求调剂（0805）
已经有8人回复
-
环境领域全国重点实验室招收博士1-2名
已经有3人回复
-
材料专硕306英一数二
已经有10人回复
-
301求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有7人回复
-
302求调剂
已经有6人回复
-
26博士申请
已经有3人回复
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
|
这是针对实数的情况，从这个结果看，看起来象是第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 |
12楼2014-02-18 21:55:06
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
2楼2014-02-17 13:10:20
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
3楼2014-02-17 13:12:44
4楼2014-02-17 19:50:28
20