讨论 ,关于多变量曲线参数拟合[有具体数据与公式]

zhotihu

其中x,y是实验可得到数据，而q是实验未知变量，但是与x 有关，x一定，不知如何拟合曲线参数.

y=2*G*R*T(ln(1-q)-H*q*q/2)-R*T/2/N/k/l*(cos(F*P/2/R/T)-1)+y[0];
x=(q/(1-q))*A/exp(-H*q)/exp(F*P/R/T);

constant R=8.315, T=298, N=6.02E+23, k=3.04E-6, l=7.14E-10, F=9.96485E+4, y[0]=72;
parameter G, q, P, H,A;
variable y, x;
29.0613    0.002999024
29.3012    0.002136781
29.6696    0.001499512
29.9532    0.00109373
33.29    0.000749756
35.8587    0.00052482
39.0339    0.000374878
41.9528    0.00026241
45.1282    0.000189226
47.0045    0.000149951
49.7463    0.000112463
53.7078    7.49756E-05
59.6153    3.74878E-05
60.7799    3.36496E-05
63.844    2.38221E-05

[ 本帖最后由 zhotihu 于 2006-9-27 16:32 编辑 ]

zhotihu

版主,请删除一个,刚才电脑突然出问题,以为没有发表成功.所以重复发表了一个

shamohu

如果各参数及变量y1的取值范围不限定，将会有不止一组的解，下面给出两组：
1：
p: 4.64961423725185
k: 1.65907115664104E-8
h: 0.00346178936117054
a: -1.11657108629155
此时的y1值：
1971.1424279429
1855.96309407041
1732.83964806135
1611.44841272735
1243.57790381632
933.615290582498
637.734761126174
328.844467588354
21.8907693463606
-191.256353433324
-486.843485773292
-866.420412995473
-1460.30015273936
-1585.80513850028
-1900.9284432348

2：
p: 5.85652934148881
k: 0.00110845808741437
h: -0.621567045449482
a: -0.614460448841969

此时的y1值：
5.93699858617846
6.55877468875397
7.22154338860103
7.86488085955475
9.62006093557225
11.110821355076
12.5317571058072
14.0164993752162
15.4825529789354
16.5023201669725
17.9048489011333
19.7177528414684
22.5735855975319
23.1643504410268
24.6722553008144

zhotihu

谢谢,你做得太好了,只是我有好几个数据要处理,可否告诉我是如何处理吗? 我的数据的确要范围,

[ 本帖最后由 zhotihu 于 2006-9-26 11:23 编辑 ]

shamohu

用Origin很麻烦，1stOpt简单、好用、强大，你的问题就是用1stOpt算的。下面的链接可以参考一下，与你的问题类似。

http://www.simwe.com/forum/viewt ... =%BD%A5%BF%AA%CF%DF

zhotihu

啊,不懂啊,可是我从来没用过,要不要很久学会呢?谢谢,

shamohu

1stOpt的特点之一就是简单、容易上手。

zhotihu

但好多算法没学过，我是学化学的。呵，不过，我试试看，如不行，能否请你帮我算一下。三组数据。

zhotihu

不知有没有参考资料？如果有能否给一份：thzhou168@yahoo.com.cn

学习中，呵，

[ 本帖最后由 zhotihu 于 2006-9-26 14:21 编辑 ]

zhotihu

请shamohu 帮我一下，不知你有没有看到我发的邮件。谢谢。

shamohu

已经发过去了

nuc2005

能够给我一份，非常感谢
nuc2005@163.com

萤火虫

我要想学习一下，能给我也发一份吧    inasense@163.com