| 查看: 1260 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
boundary condition 求 free constant
|
|||
|
我不知道这个能不能用Mathematica来求，但是我用手算，太麻烦了。我过几天再试试，但是如果谁知道用mathematica来求的话，最好了。 我有个式子是 [ Last edited by wocamuchong on 2014-2-17 at 07:17 ] |
回复此楼» 猜你喜欢
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有7人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
小木虫看见有人已经知道结果了
已经有19人回复
-
2026年国自然面上资助率
已经有8人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
微信指数没变化，科研之友没阅读
已经有17人回复
-
申请2027年材料/化学类博士
已经有5人回复
-
你们的时间戳变了吗
已经有4人回复
-
今年的WR进展到哪一步了？
已经有9人回复
-
产物和副产物价值比较
已经有5人回复
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