| 查看: 1293 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
boundary condition 求 free constant
|
|||
|
我不知道这个能不能用Mathematica来求，但是我用手算，太麻烦了。我过几天再试试，但是如果谁知道用mathematica来求的话，最好了。 我有个式子是 [ Last edited by wocamuchong on 2014-2-17 at 07:17 ] |
回复此楼» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有8人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有18人回复
-
感觉是下周放榜了
已经有12人回复
-
我面上完蛋了
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有25人回复
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4677.7
- 红花: 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 (幼儿园)
- 金币: 4677.7
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
2楼2014-02-17 13:10:20
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4677.7
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
3楼2014-02-17 13:12:44
4楼2014-02-17 19:50:28
40