声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2364|回复: 2

[1stopt] 求必过某几个点都拟合方法!!!!!!!着急!!!!

[复制链接]
发表于 2009-12-1 11:53 | 显示全部楼层 |阅读模式

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

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

x
我用1stopt拟合的程序,是网上都1.5版本,程序如下:
Title "fdhjku";
Parameters a3,a4,a5,a6,a7;
Constant pLG=159.0058,VG=169.4644,VL=30.69681,T=620,R=83.1451;
Variable V,y;
Function   y=(1/V-1/60*1/V^2*(-60*pLG*VL^7*VG^6+60*ln(pLG*VL/R/T)*R*T*VL^6*VG^6-30*a3*R*T*VL^4*VG^6-20*a4*R*T*VL^3*VG^6-15*a5*R*T*VL^2*VG^6-12*a6*R*T*VL*VG^6-10*a7*R*T*VG^6+60*pLG*VG^7*VL^6-60*ln(pLG*VG/R/T)*R*T*VL^6*VG^6+30*a3*R*T*VL^6*VG^4+20*a4*R*T*VL^6*VG^3+15*a5*R*T*VL^6*VG^2+12*a6*R*T*VL^6*VG+10*a7*R*T*VL^6)/R/T/VL^5/VG^5/(-VG+VL)+1/V^3*a3+1/V^4*a4+1/V^5*a5+1/V^6*a6+1/V^7*a7)*R*T;
Data;
  V,   y=
515399.8    0.1
171739.1    0.3
103006.9    0.5
73550.26    0.7
57185.43    0.9
51457.73    1
25682.99    2




1069.903    44.05
1065.897    44.2
1061.917    44.35

请问,高人,如果要求必过红色字都两个点,该如何写程序呢,按照教程写了一个,运行的结果好像不行呀,哪一位高人能否帮我写写程序呢,不胜感激!!!!!!!
回复
分享到:

使用道具 举报

发表于 2009-12-1 22:57 | 显示全部楼层
运行了下,似乎每个点都在线上。

均方差(RMSE): 6.18479832992445E-6
残差平方和(RSS): 3.82517303818363E-10
相关系数(R): 0.999999999999961
相关系数之平方(R^2): 0.999999999999922
决定系数(DC): 0.999999999999903
卡方系数(Chi-Square): 1.57528795825961E-10
F统计(F-Statistic): 12942074647558.3

参数        最佳估算
----------        -------------
a3                 -94.9883443817995
a4                 106174.072223172
a5                 -966.700930901553
a6                 -10313.6208991701
a7                 3242.41108545838

====== 结果输出 =====

No        实测值y        计算值y
1        0.1        0.100001561363804
2        0.3        0.300004019572309
3        0.5        0.500005754202415
4        0.7        0.700006389496424
5        0.9        0.900006311666702
6        1        1.00000604285477
7        2        1.9999925261177
8        44.05        44.0500058980491
9        44.2        44.1999897900596
10        44.35        44.3500043294297
 楼主| 发表于 2009-12-2 22:31 | 显示全部楼层
楼上的朋友,很是感谢。
只是我只是列了一部分点,我只是想知道怎么能够保证过某个点或者几个点。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-10 15:50 , Processed in 0.057419 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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