| 查看: 1570 | 回复: 13 | |||
[求助]
同一个方程，在matab中的答案和在1spopt中结果不一致， 已有2人参与
|
|||
|
下面是一个求解价格P的一元一次方程： (1） matlab中的代码如下： p=solve(' (-0.5)*2/0.5331*( ( p-200+( (80+0.1*120)/0.1) )*(0.3671+log((1+0.1*(0.5331-0.3671)))/0.1) - ((0.08*(200+40*0.04)+40)/0.08)/0.08 *((2.71828182845905^((0.08*(0.3671-0.04)))-1)-((0.08*(0.3671-0.04)))) -(40*0.04*0.3671)+(1/2*40*0.04^2)- ((80+0.1*120)/0.1)*0.5331 )+2*((500-0.5*p)+2)/0.5331*(0.3671+log((1+0.1*(0.5331 -0.3671)))/0.1)+1 ','p') 解得 p=607.05064790464619193396327973728 （2） 1spopt中的代码如下： Parameters p[0,] ; Function (-0.5)*2/0.5331*( ( p-200+( (80+0.1*120)/0.1) )*(0.3671+log((1+0.1*(0.5331-0.3671)))/0.1) - ((0.08* (200+40*0.04)+40)/0.08)/0.08 *((2.71828182845905^((0.08*(0.3671-0.04)))-1)-((0.08*(0.3671-0.04)))) - (40*0.04*0.3671)+(1/2*40*0.04^2)- ((80+0.1*120)/0.1)*0.5331 )+2*((500-0.5*p)+2)/0.5331*(0.3671+log((1+0.1*(0.5331 -0.3671)))/0.1)+1 =0; 解得 p= 705.803382033099 (3) matlab 先对表达式处理 w=(-0.5)*2/0.5331*( ( p-200+( (80+0.1*120)/0.1) )*(0.3671+log((1+0.1*(0.5331-0.3671)))/0.1) - ((0.08*(200+40*0.04)+40)/0.08)/0.08 *((2.71828182845905^((0.08*(0.3671-0.04)))-1)-((0.08*(0.3671-0.04)))) -(40*0.04*0.3671)+(1/2*40*0.04^2)- ((80+0.1*120)/0.1)*0.5331 )+2*((500-0.5*p)+2)/0.5331*(0.3671+log((1+0.1*(0.5331 -0.3671)))/0.1)+1 得到： w=94643453558687378453/78153286502318080 - (2993414659220048125*p)/1500543100844507136 vpa(w) 得到：1210.9977429532690068640285035303 - 1.9948874894265625590242248762822*p p=solve(' 1210.9977429532690068640285035303 - 1.9948874894265625590242248762822*p ','p') 解得 p=607.05064790464619193396327973728 （4） 再一次在 1spopt中输入的代码如下 Parameters p[0,] ; Function 1210.9977429532690068640285035303 - 1.9948874894265625590242248762822*p =0; 解得 p=607.05064790464619193396327973728 问题： 同一代码，两次在matlab 答案一致，两次在1spopt中 答案不一致，请问何故？ |
回复此楼» 本帖已获得的红花（最新10朵）
maganet» 猜你喜欢
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有10人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
需要合成515-64-0，50g，能接单的留言
已经有3人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
-
临港实验室与上科大联培博士招生1名
已经有8人回复
-
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有4人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2017-07-27 11:53:25
askuyue
铁杆木虫 (著名写手)
- 应助: 12 (小学生)
- 金币: 6810.2
- 散金: 12
- 红花: 13
- 帖子: 2024
- 在线: 188.3小时
- 虫号: 2260045
- 注册: 2013-01-25
- 性别: GG
- 专业: 数学
2楼2017-07-27 03:22:45
3楼2017-07-27 07:16:45
5楼2017-07-27 14:55:30
6楼2017-07-27 14:55:55
Mr__Right
专家顾问 (著名写手)
-
专家经验: +31 - 数学EPI: 1
- 应助: 317 (大学生)
- 金币: 14456.3
- 散金: 500
- 红花: 54
- 帖子: 2716
- 在线: 950.6小时
- 虫号: 1972612
- 注册: 2012-09-04
- 性别: GG
- 专业: 应用数学方法
- 管辖: 程序语言
7楼2017-07-27 20:26:32
Mr__Right
专家顾问 (著名写手)
-
专家经验: +31 - 数学EPI: 1
- 应助: 317 (大学生)
- 金币: 14456.3
- 散金: 500
- 红花: 54
- 帖子: 2716
- 在线: 950.6小时
- 虫号: 1972612
- 注册: 2012-09-04
- 性别: GG
- 专业: 应用数学方法
- 管辖: 程序语言
8楼2017-07-27 20:27:51
9楼2017-07-28 00:04:04
10楼2017-07-28 00:04:20