| 查看: 2452 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
tuyuting新虫 (小有名气)
|
[求助]
求助~~~二元三次方程用什么软件求极值 已有3人参与
|
||
|
方程为： z=-92.743+21.940*x+1.535*y-7.368*x^2-0.006*y^2+9.359*x*y/10000+0.803*x^3+(9.853*y^3)/1000000-(2.889*x*y^2)/10000+0.0255*y*x^2 其中x的范围为0到2.5 ； y的范围为100到200 尝试用matlab求极大值，参考网上的教程，编的代码如下 syms x1 x2; y=-92.743+21.940*x1+1.535*x2-7.368*x1^2-0.006*x2^2+9.359*x1*x2/10000+0.803*x1^3+(9.853*x2^3)/1000000-2.889*x1*x2^2/10000+0.0255*x1^2*x2 y=simple(y) dydx1=diff(y,x1) dydx2=diff(y,x2) S=solve(dydx1,dydx2) S.x1 S.x2 计算结果是一大堆公式 y = -92743/1000+1097/50*x1+307/200*x2-921/125*x1^2-3/500*x2^2+9359/10000000*x1*x2+803/1000*x1^3+9853/1000000000*x2^3-2889/10000000*x1*x2^2+51/2000*x1^2*x2 y = 803/1000*x1^3+(-921/125+51/2000*x2)*x1^2+(1097/50+9359/10000000*x2-2889/10000000*x2^2)*x1-92743/1000+307/200*x2+9853/1000000000*x2^3-3/500*x2^2 dydx1 = 2409/1000*x1^2+2*(-921/125+51/2000*x2)*x1+1097/50+9359/10000000*x2-2889/10000000*x2^2 dydx2 = 51/2000*x1^2+(9359/10000000-2889/5000000*x2)*x1+307/200+29559/1000000000*x2^2-3/250*x2 S = x1: [4x1 sym] x2: [4x1 sym] ans = 2159960412801613996538037/559496974666398295037916100*(72886429125179524/362347273980951+1/1347082716481861419470583858351152538612*6^(1/2)*((-370175839116662828260169439633164657624844729276059377616214790883620376135999116*(-5476821316928811165845458719373673857342723142053626365094440283912473672234915730406620344251339195447860696279655747108837039168333940756549124378354994324333114393516897294565622823208223590992509759772727657498758565678213186720956+10609852924969292005429057493877030663111181724710719087057168516226374251717872739658280037283020433159528825176519873227495027565825932363212*i*412301133939570303333491378185360079681893557830658491717107061914374469161463266738792864832300686314186337380891301039123524853581262522719077718180352606647106995680690230930274470027^(1/2))^(1/3)+217*。。。。省略 ans = 72886429125179524/362347273980951+1/1347082716481861419470583858351152538612*6^(1/2)*((-370175839116662828260169439633164657624844729276059377616214790883620376135999116*(-5476821316928811165845458719373673857342723142053626365094440283912473672234915730406620344251339195447860696279655747108837039168333940756549124378354994324333114393516897294565622823208223590992509759772727657498758565678213186720956+10609852924969292005429057493877030663111181724710719087057168516226374251717872739658280037283020433159528825176519873227495027565825932363212*i*412301133939570303333491378185360079681893557830658491717107061914374469161463266738792864832300686314186337380891301039123524853581262522719077718180352606647106995680690230930274470027^(1/2))^(1/3)+217*。。。省略 得不到具体的数值 麻烦各位虫友帮忙看看，代码出错在哪儿？多谢 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有294人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
8楼2015-07-15 10:02:11
WanderingHeart
铁杆木虫 (著名写手)
- 应助: 90 (初中生)
- 金币: 13890.8
- 散金: 5
- 红花: 18
- 帖子: 1885
- 在线: 467.1小时
- 虫号: 257464
- 注册: 2006-06-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2015-07-13 21:43:16
tuyuting
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 273.4
- 散金: 337
- 帖子: 245
- 在线: 433.9小时
- 虫号: 829015
- 注册: 2009-08-16
- 性别: MM
- 专业: 污染控制化学
3楼2015-07-13 21:51:01
WanderingHeart
铁杆木虫 (著名写手)
- 应助: 90 (初中生)
- 金币: 13890.8
- 散金: 5
- 红花: 18
- 帖子: 1885
- 在线: 467.1小时
- 虫号: 257464
- 注册: 2006-06-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2015-07-14 07:05:04
投票:
10