| 查看: 1457 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
在几个条件下求解
|
|||
| 我把mathematica的文件上传了，懂的帮忙看一下，就是在条件求a3,a4，我不太会用mathematica,谢了。 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Solvea3a4.nb
2013-06-21 04:34:27, 5.83 K
» 猜你喜欢
347求调剂
已经有4人回复
-
一志愿西交288化工专硕求调剂
已经有3人回复
-
一志愿中国科学院大学265求调剂
已经有4人回复
-
291求调剂
已经有13人回复
-
070300化学求调剂
已经有3人回复
-
322求调剂
已经有7人回复
-
0703化学
已经有15人回复
-
材料学硕333求调剂
已经有15人回复
-
313求调剂
已经有3人回复
-
083000学硕274求调剂
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有没有黑暗条件下光催化剂催化甲基橙脱色这种情况？
已经有16人回复
- 数学物理方程求解！（带边界条件的偏微分方程）
已经有13人回复
- 求解担保人条件第五条
已经有3人回复
- 希夫碱什么条件会分解
已经有18人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 求解什么是Biomass!!!!
已经有16人回复
- 求解 DIFMN02_ALERT_2_A
已经有9人回复
- 太原理工大学、盐城工学院、中科院上海研究所、中科院重庆研究所？改选哪一个？
已经有32人回复
- 【求助】初学有限元，弱弱问一句，如果已知边界都是位移边界条件，怎么求解？
已经有11人回复
- 原料药制胶囊时吸湿性大，该如何解决？加什么辅料可以解决问题？
已经有4人回复
- 【求助】求问一个寻找最优解问题，是泛函问题么
已经有22人回复
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jollage(xzhdty代发): 金币+10, 代楼主发放金币 2013-06-25 07:13:49
jollage(xzhdty代发): 金币+10, 代楼主发放金币 2013-06-25 07:13:49
|
可能你给的表达式有点差别 (或者这里有个零点的问题 ) 我试了下 你给的值比上面的值 用附件的表达式要大些 下面是代码 可能效率比较低 算一下大概3分钟(全局最小值) 代码给了2个结果 我觉得第一个结果更合理 （这里面比较奇怪的是 第一条件其实应该比第二个强 但是得到的最小值却比第二个小－－ 这点我也还没想明白 细致之处需要进一步考虑） Clear["Global`*"]; Dk = 16 a3 + 6 (1/16 - 2 a3 - 4 a4) + 30 a4 + 2 (-9 a3 - 4 (1/16 - 2 a3 - 4 a4) - 16 a4) Cos[k] + 2 (1/16 - 2 a3 - 4 a4) Cos[2 k] + 2 a3 Cos[3 k] + 2 a4 Cos[4 k] // FullSimplify; Dkp = 5 (16 a3 + 6 (1/16 - 2 a3 - 4 a4) + 30 a4 + 2 (-9 a3 - 4 (1/16 - 2 a3 - 4 a4) - 16 a4) Cos[k] + 2 (1/16 - 2 a3 - 4 a4) Cos[2 k] + 2 a3 Cos[3 k] + 2 a4 Cos[4 k]) + k (-2 (-9 a3 - 4 (1/16 - 2 a3 - 4 a4) - 16 a4) Sin[k] - 4 (1/16 - 2 a3 - 4 a4) Sin[2 k] - 6 a3 Sin[3 k] - 8 a4 Sin[4 k]) // FullSimplify; Ob = 0.08904499121480983 + 3.616287830531132 a3 + 9.423861054539543 a4; f1[at3_?NumericQ, at4_?NumericQ] := NMinimize[Dk /. {a3 -> at3, a4 -> at4}, {k}][[1]]; f2[at3_?NumericQ, at4_?NumericQ] := NMinimize[Dkp /. {a3 -> at3, a4 -> at4}, {k}][[1]]; re1 = NMinimize[{Ob, f1[a3, a4] > 0, f2[a3, a4] > 0}, {a3, a4}, MaxIterations -> 1000] // AbsoluteTiming \[Epsilon] = 10^-20; f3[at3_?NumericQ, at4_?NumericQ] := NMinimize[{Dk /. {a3 -> at3, a4 -> at4}, Pi - \[Epsilon] > k > \[Epsilon]}, {k}][[1]]; f4[at3_?NumericQ, at4_?NumericQ] := NMinimize[{Dkp /. {a3 -> at3, a4 -> at4}, Pi - \[Epsilon] > k > \[Epsilon]}, {k}][[1]]; re2 = NMinimize[{Ob, f3[a3, a4] > 0, f4[a3, a4] > 0}, {a3, a4}, MaxIterations -> 1000] // AbsoluteTiming |
8楼2013-06-23 17:44:25
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
2楼2013-06-21 16:40:31
3楼2013-06-22 03:05:19
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
|
奇怪，这次试了下，发现结果报告是无下限。 同时，直接取k=1,代入进去求，也是无下限， 而且 从表达式上看，Dk,Dkp,OB 都是a3,a4的线性叠加， 而且系数比(a3系数/a4系数)分别是0.324和0.348之间 而OB的表达式中比例系数比0.38多点，这导致全局无上下限 （因为这种情况下，可以取a3/a4=0.35,a3＝-Inifinity,满足约束，但Ob负无穷大.） 取k为其他数时，相互间的比例系数变化很大，不存在一个比值落在所有的Ob 表达式 我觉得可能是你OB表示式中a3,a4前面系数的比值不对 或者是a3,a4还有些限制 在k=1的时候， 把Ob改成=0.33a3+a4 （常数和比例因子扔了） 这时候有极小值..... 问题本身是不是求全局的极值 还是求局部的？ 本身对a3,a4是不是还有约束条件? 画图的 |
4楼2013-06-22 07:55:27
