声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2504|回复: 3

[1stopt] 求教高人1stopt的拟合问题,总是提示出错

[复制链接]
发表于 2010-1-5 17:55 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
学习1stopt不久,自己也算了一些东西,很方便,但是现在这个程序不行了,请大家帮我看看是什么问题呢?程序如下:

Title "Type your title here";
Parameters p51,p52,p53,p54,p55,p61,p62,p63,p64,p65,p71,p72,p73,p74,p75;
Constant R=83.1451,Tc=647.096;
Variable V,pLG,VG,VL,T,y;
Function   y=R*T*(1/V-(1/10)*(60*ln(pLG*VL/(R*T))*R*T*VL^6*VG^6-60*ln(pLG*VG/(R*T))*R*T*VL^6*VG^6
+10*VL^5*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)*VG^2
+10*VL^5*VG*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
+10*VL^5*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T+p55*Tc^2/T^2)*VG^3
-10*VL^2*VG^4*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
-10*VL^3*VG^4*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)
+10*VL^4*VG^3*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)
+10*VL^4*VG^2*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
-10*VL*VG^5*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
-10*VL^3*VG^5*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T+p55*Tc^2/T^2)
-10*VL^2*VG^5*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)
-10*VL^5*pLG*VG^8-80*pLG*VL^7*VG^6+80*pLG*VG^7*VL^6+10*VL^8*VG^5*pLG
+10*VL^9*VG^4*pLG-10*VL^4*pLG*VG^9+10*VL^5*R*T*VG^7-10*VL^7*R*T*VG^5
-10*VL^8*R*T*VG^4+10*VL^4*R*T*VG^8-(8*(p61+p62*T/Tc+p63*T^2/Tc^2
+p64*Tc/T+p65*Tc^2/T^2))*R*T*VL^6*VG-(5*(p51+p52*T/Tc+p53*T^2/Tc^2
+p54*Tc/T+p55*Tc^2/T^2))*R*T*VL^6*VG^2+(8*(p61+p62*T/Tc+p63*T^2/Tc^2
+p64*Tc/T+p65*Tc^2/T^2))*R*T*VL*VG^6+(5*(p51+p52*T/Tc+p53*T^2/Tc^2
+p54*Tc/T+p55*Tc^2/T^2))*R*T*VL^2*VG^6+(10*(p71+p72*T/Tc+p73*T^2/Tc^2
+p74*Tc/T+p75*Tc^2/T^2))*R*T*VG^6-(10*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T
+p75*Tc^2/T^2))*R*T*VL^6)/(VG^4*VL^4*R*T*(-VL^3+3*VL^2*VG+VG^3-3*VL*VG^2)*V^2)
+(1/10)*(60*ln(pLG*VL/(R*T))*R*T*VL^6*VG^7-60*ln(pLG*VG/(R*T))*R*T*VL^6*VG^7
+60*ln(pLG*VL/(R*T))*R*T*VL^7*VG^6-60*ln(pLG*VG/(R*T))*R*T*VL^7*VG^6
-10*VG*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)*R*T*VL^6
+40*VL^5*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)*VG^3
+40*VL^5*VG^2*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
+40*VL^5*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T+p55*Tc^2/T^2)*VG^4
+10*VL*VG^6*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
+5*VL^3*VG^6*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T+p55*Tc^2/T^2)
+8*VL^2*VG^6*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)
-40*VL^2*VG^5*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T+p75*Tc^2/T^2)
-40*VL^4*VG^5*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T+p55*Tc^2/T^2)
-40*VL^3*VG^5*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)
+80*pLG*VG^8*VL^6-80*VL^8*VG^6*pLG+40*VL^9*VG^5*pLG-40*VL^5*pLG*VG^9
+40*VL^5*R*T*VG^8-40*VL^8*R*T*VG^5+(8*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T
+p65*Tc^2/T^2))*R*T*VL*VG^7-(8*(p61+p62*T/Tc+p63*T^2/Tc^2
+p64*Tc/T+p65*Tc^2/T^2))*R*T*VL^6*VG^2-(5*(p51+p52*T/Tc+p53*T^2/Tc^2
+p54*Tc/T+p55*Tc^2/T^2))*R*T*VL^6*VG^3+(10*(p71+p72*T/Tc+p73*T^2/Tc^2
+p74*Tc/T+p75*Tc^2/T^2))*R*T*VG^7-(8*(p61+p62*T/Tc+p63*T^2/Tc^2
+p64*Tc/T+p65*Tc^2/T^2))*R*T*VL^7*VG-(5*(p51+p52*T/Tc+p53*T^2/Tc^2
+p54*Tc/T+p55*Tc^2/T^2))*R*T*VL^7*VG^2+(5*(p51+p52*T/Tc+p53*T^2/Tc^2
+p54*Tc/T+p55*Tc^2/T^2))*R*T*VL^2*VG^7-(10*(p71+p72*T/Tc+p73*T^2/Tc^2
+p74*Tc/T+p75*Tc^2/T^2))*R*T*VL^7)/(VL^4*VG^4*R*T*(-VL^3+3*VL^2*VG+VG^3
-3*VL*VG^2)*V^3)-(1/10)*((15*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T
+p55*Tc^2/T^2))*R*T*VL^2*VG^6+(18*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T
+p65*Tc^2/T^2))*R*T*VL*VG^6+(20*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T
+p75*Tc^2/T^2))*R*T*VG^6-(15*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T
+p55*Tc^2/T^2))*R*T*VL^6*VG^2-(18*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T
+p65*Tc^2/T^2))*R*T*VL^6*VG-(20*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T
+p75*Tc^2/T^2))*R*T*VL^6+60*ln(pLG*VL/(R*T))*R*T*VL^6*VG^6
-60*ln(pLG*VG/(R*T))*R*T*VL^6*VG^6-30*VL^5*pLG*VG^8-90*pLG*VL^7*VG^6
+90*pLG*VG^7*VL^6+30*VL^5*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2+p64*Tc/T
+p65*Tc^2/T^2)*VG^2+30*VL^5*VG*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2+p74*Tc/T
+p75*Tc^2/T^2)+30*VL^5*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T
+p55*Tc^2/T^2)*VG^3-30*VL*VG^5*R*T*(p71+p72*T/Tc+p73*T^2/Tc^2
+p74*Tc/T+p75*Tc^2/T^2)-30*VL^3*VG^5*R*T*(p51+p52*T/Tc+p53*T^2/Tc^2
+p54*Tc/T+p55*Tc^2/T^2)-30*VL^2*VG^5*R*T*(p61+p62*T/Tc+p63*T^2/Tc^2
+p64*Tc/T+p65*Tc^2/T^2)+30*VL^5*R*T*VG^7-30*VL^7*R*T*VG^5
+30*VL^8*VG^5*pLG)/(R*T*VL^3*VG^3*(-VL^3+3*VL^2*VG+VG^3-3*VL*VG^2)*V^4)
+(p51+p52*T/Tc+p53*T^2/Tc^2+p54*Tc/T+p55*Tc^2/T^2)/V^5+(p61+p62*T/Tc
+p63*T^2/Tc^2+p64*Tc/T+p65*Tc^2/T^2)/V^6+(p71+p72*T/Tc+p73*T^2/Tc^2
+p74*Tc/T+p75*Tc^2/T^2)/V^7);  方程项很多
Data;
V,pLG,VG,VL,T,y=
227113342.7        0.006116548        3710987.122        18.01900658        273.16        0.0001
113555658.6        0.006116548        3710987.122        18.01900658        273.16        0.0002
75703097.02        0.006116548        3710987.122        18.01900658        273.16        0.0003
56776816.07        0.006116548        3710987.122        18.01900658        273.16        0.0004
45421047.33        0.006116548        3710987.122        18.01900658        273.16        0.0005
37850534.66        0.006116548        3710987.122        18.01900658        273.16        0.0006
32443025.44        0.006116548        3710987.122        18.01900658        273.16        0.0007
28387393.35        0.006116548        3710987.122        18.01900658        273.16        0.0008
25233012.66        0.006116548        3710987.122        18.01900658        273.16        0.0009
22709507.93        0.006116548        3710987.122        18.01900658        273.16        0.001
20644822.06        0.006116548        3710987.122        18.01900658        273.16        0.0011
18924250.32        0.006116548        3710987.122        18.01900658        273.16        0.0012
(数据很多,这里就略了)

一运行,就提示“读数据出错”,都急死了,不知道到底是什么问题,初步估计是函数的问题,可是折腾了半天,还是没有找到点上。

请高人帮我找找,到底是什么问题。非常感谢!!!!!!!!!!

[ 本帖最后由 wlxydgq 于 2010-1-5 18:03 编辑 ]
回复
分享到:

使用道具 举报

发表于 2010-1-5 20:54 | 显示全部楼层
把“Data;”前的汉字“方程项很多”去掉,运行一下没问题。
 楼主| 发表于 2010-1-6 10:31 | 显示全部楼层
:@L 怎么会呢?那我运行就说:读 数据有问题。
 楼主| 发表于 2010-1-7 14:55 | 显示全部楼层
问题解决了,感谢dingd!!!!!!!!!!!!!!!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-10 16:07 , Processed in 0.061958 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表