论坛 › 科学分析软件 › 求教高人1stopt的拟合问题，总是提示出错

wlxydgq 发表于 2010-1-5 17:55

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

学习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 编辑 ]

dingd 发表于 2010-1-5 20:54

把“Data;”前的汉字“方程项很多”去掉，运行一下没问题。

wlxydgq 发表于 2010-1-6 10:31

:@L 怎么会呢？那我运行就说：读 数据有问题。

wlxydgq 发表于 2010-1-7 14:55

问题解决了，感谢dingd！！！！！！！！！！！！！！！

查看完整版本: 求教高人1stopt的拟合问题，总是提示出错