24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

maganet

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 229.4
红花: 1
帖子: 22
在线: 56.5小时
虫号: 6953635
注册: 2017-07-15
专业: 宏观经济管理与战略

[求助] 同一个方程，在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中答案不一致，请问何故？

回复此楼