| 查看: 728 | 回复: 3 | ||||
afl523木虫 (小有名气)
|
[求助]
求大神帮忙编一下程序 已有1人参与
|
|
这里有8个方程，8个未知数，该怎么样解呀？ 10^(6.72-9500/1273)=v*c/w 10^(2.96-7510/1273)=b*c/x 10^(3.63-8700/1273)=v*n/y 10^(2.8-8500/1273)=b*n/z (0.02-v)/(0.26-c)=(w+y)*4.24146545/(w+x) (0.02-v)/(0.0103-n)=(w+y)*3.6359743/(y+z) (0.02-v)/(0.02-b)=(w+y)*0.548272522/(x+z) w+x+y+z=1 求解'b','c','n','v','w','x','y','z' 这8个未知数都〉0 谢谢各位了 |
回复此楼» 猜你喜欢
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有3人回复
-
关于代码变化问题，想知道的进来
已经有11人回复
-
基金中了
已经有13人回复
-
关于filecode
已经有6人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stopt 语句修改
已经有23人回复
- 求大神帮检查下一个简单的UDF
已经有13人回复
- 我的一个程序运行总是不成功，求大神帮忙
已经有3人回复
- 求助 EC50 怎么算！！！！ 求大神帮忙
已经有3人回复
- 求大神帮忙解题啊！！！
已经有4人回复
- 修改意见中的一句话，求助大神帮忙翻译一下
已经有5人回复
- 积分求助 求大神帮忙！
已经有7人回复
- 求大神们帮我解决一个关于修改系统进化树字体的问题
已经有18人回复
- 很简单的Fortran赋值出错啦，求大神们帮忙找bug...
已经有11人回复
- 哪位大神帮忙看一下MATLAB求解微分方程的问题？
已经有3人回复
- 编程菜鸟求大神引导入门。
已经有14人回复
- 求大神帮忙解答！
已经有8人回复
- fluent安装，不能用，求助大神帮忙~~~
已经有17人回复
- 求大神帮忙求KIE
已经有8人回复
- abaqus安装出错，求大神帮忙！！！
已经有3人回复
- 求大神帮忙看一下断口
已经有6人回复
- 全部身家求大神帮忙解决fluent压力脉冲怎么编写UDF的问题
已经有8人回复
- Perl 编程 请大神们帮忙！！！生物信息的牛人请进
已经有11人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
感谢参与，应助指数 +1
|
大概计算一下，解似乎不稳定唯一： 迭代数: 118 计算用时(时:分:秒:毫秒): 00:00:13:666 计算中止原因: 达到收敛判定标准 优化算法: 准牛顿法(BFGS) + 通用全局优化法 函数表达式 1: 10^(6.72-9500/1273)-(v*c/w) 2: 10^(2.96-7510/1273)-(b*c/x) 3: 10^(3.63-8700/1273)-(v*n/y) 4: 10^(2.8-8500/1273)-(b*n/z) 5: (0.02-v)/(0.26-c)-((w+y)*4.24146545/(w+x)) 6: (0.02-v)/(0.0103-n)-((w+y)*3.6359743/(y+z)) 7: (0.02-v)/(0.02-b)-((w+y)*0.548272522/(x+z)) 8: w+x+y+z-(1) 目标函数值: 1.47665524119339E-8 x: 0.96512144751794 y: 9.88405936807881E-9 z: 0.0348785344311554 w: 8.18627576923919E-9 v: 0.0199999806006773 b: 14988.6195000327 c: 7.40236045278039E-8 n: 3.08784838784678E-10 目标函数值: 2.06089097154221E-8 x: 0.914518283604818 y: 4.37575401472771E-9 z: 0.0854817107020073 w: 1.40119794298616E-9 v: 0.019999997484647 b: 82977.3109992174 c: 1.26701857774441E-8 n: 1.36701458367838E-10 目标函数值: 1.25093438027696E-5 x: 0.95744742961639 y: 4.62756870069125E-6 z: 0.042544825733514 w: 3.11708270594457E-6 v: 0.0199931835253893 b: 39.0377711429402 c: 2.81955003512332E-5 n: 1.44617589696085E-7 目标函数值: 0.000191180735154163 x: 0.876588315538741 y: 1.6916414317921E-11 z: 0.123397944417429 w: 1.37400269123709E-5 v: 0.0199958299791895 b: 8.1093114185841 c: 0.000124268650130806 n: 5.28590160184477E-13 目标函数值: 0.00932006988104486 x: 0.718323251981402 y: 0.109939780354001 z: 0.142779884301518 w: 0.0289570833630833 v: 0.0201289774015204 b: 0.0214584108070313 c: 0.260163602073569 n: 0.0103645415536109 |
2楼2015-03-09 20:36:25
afl523
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 2261.1
- 散金: 127
- 红花: 2
- 帖子: 65
- 在线: 23.2小时
- 虫号: 2062164
- 注册: 2012-10-15
- 性别: MM
- 专业: 金属结构材料
3楼2015-03-10 21:57:54
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
afl523: 金币+15, ★★★很有帮助, 谢谢啦 2015-03-12 09:25:17
afl523: 金币+15, ★★★很有帮助, 谢谢啦 2015-03-12 09:25:17
|
这是 1stopt 软件的计算结果，代码如下： Parameters x[0,],y[0,],z[0,],w[0,],v[0,],b[0,],c[0,],n[0,]; Function 10^(6.72-9500/1273)-(v*c/w) ; 10^(2.96-7510/1273)-(b*c/x) ; 10^(3.63-8700/1273)-(v*n/y) ; 10^(2.8-8500/1273)-(b*n/z) ; (0.02-v)/(0.26-c)-((w+y)*4.24146545/(w+x)); (0.02-v)/(0.0103-n)-((w+y)*3.6359743/(y+z)) ; (0.02-v)/(0.02-b)-((w+y)*0.548272522/(x+z)) ; w+x+y+z-(1); 软件获得与教程：http://muchong.com/bbs/viewthread.php?tid=4586018 |
4楼2015-03-10 22:06:27
100