morrissey22 发表于 2007-12-24 04:34

有1stopt 2.0版本的朋友请帮个忙,平安夜求助,好人好运

要拟合个非线性公式,被催得紧,无奈买不起正版,下的1.5版SharedModel没法用,一运行就死掉了.拜托用2.0版的朋友帮忙算下下面段程序,贴个结果给我,万分万分感谢!!!


Title "OPO";

Parameters P1=4.532,P2=0.014,P3=0.25,P4=-0.022[-0.02365,-0.021],P5=0.067,P6=0.033,P7=,P8=,P9=3E-8 ,P10=1.5E-8;

Variable T1,T2,rs,r,y1,y2;

ConstStr rp=1.064 ;
ConstStr ri=1/(1/rp-1/rs) ;
ConstStr NP=(P1+(P2+P9*(T1+273.15)^2)/(rp^2-(P3+P10*(T1+273.15)^2)^2)+P5/(rp^2-P6)+P4*rp^2+P7/(rp^2-P8))^0.5 ;
ConstStr NS=(P1+(P2+P9*(T1+273.15)^2)/(rs^2-(P3+P10*(T1+273.15)^2)^2)+P5/(rs^2-P6)+P4*rs^2+P7/(rs^2-P8))^0.5 ;
ConstStr NI=(P1+(P2+P9*(T1+273.15)^2)/(ri^2-(P3+P10*(T1+273.15)^2)^2)+P5/(ri^2-P6)+P4*ri^2+P7/(ri^2-P8))^0.5 ;
ConstStr NH=(P1+(P2+P9*(T2+273.15)^2)/((r/2)^2-(P3+P10*(T2+273.15)^2)^2)+P5/((r/2)^2-P6)+P4*(r/2)^2+P7/((r/2)^2-P8))^0.5 ;
ConstStr N=(P1+(P2+P9*((T2+273.15)^2))/(r^2-(P3+P10*((T2+273.15)^2))^2)+P5/(r^2-P6)+P4*(r^2)+P7/(r^2-P8))^0.5 ;

SharedModel;
Function y1=1/(NP/rp-NS/rs-NI/ri)/(1+1.6E-5*(T1-25)+7E-9*(T1-25)^2);
y2=r/(2*(NH-N))/(1+1.6E-5*(T2-25)+7E-9*(T2-25)^2);

Data;
T1,T2,rs,r,y1,y2=
30 137 1.50929 1.064 30 7.8
50 88 1.51336 1.064 30 7.9
70 33.5 1.51765 1.064 30 8
90 49 1.52258 0.976 30 6.1
110 57.5 1.52765 1.084 30 8.3
130 137 1.5335 1.064 30 7.8
150 88 1.54 1.064 30 7.9
170 33.5 1.5469 1.064 30 8
30 49 1.52006 0.976 30.2 6.1
50 57.5 1.5244 1.084 30.2 8.3
70 137 1.52914 1.064 30.2 7.8
90 88 1.5345 1.064 30.2 7.9
110 33.5 1.54028 1.064 30.2 8
130 49 1.5466 0.976 30.2 6.1
150 57.5 1.55346 1.084 30.2 8.3
170 137 1.56297 1.064 30.2 7.8
30 88 1.533 1.064 30.4 7.9
50 33.5 1.53778 1.064 30.4 8
70 49 1.5429 0.976 30.4 6.1
90 57.5 1.5486 1.084 30.4 8.3
110 137 1.555 1.064 30.4 7.8
130 88 1.562 1.064 30.4 7.9
150 33.5 1.5696 1.064 30.4 8
170 49 1.57835 0.976 30.4 6.1
30 57.5 1.6154 1.084 31.4 8.3
50 137 1.6232 1.064 31.4 7.8
70 88 1.6325 1.064 31.4 7.9
90 33.5 1.64336 1.064 31.4 8
110 49 1.65437 0.976 31.4 6.1
130 57.5 1.6674 1.084 31.4 8.3
150 137 1.68242 1.064 31.4 7.8
170 88 1.70029 1.064 31.4 7.9
30 33.5 1.6378 1.064 31.6 8
50 49 1.6467 0.976 31.6 6.1
70 57.5 1.6571 1.084 31.6 8.3
90 137 1.6683 1.064 31.6 7.8
110 88 1.6824 1.064 31.6 7.9
130 33.5 1.698 1.064 31.6 8
150 49 1.7158 0.976 31.6 6.1
170 57.5 1.7376 1.084 31.6 8.3
30 137 1.6651 1.064 31.8 7.8
50 88 1.67647 1.064 31.8 7.9
70 33.5 1.6889 1.064 31.8 8
90 49 1.7028 0.976 31.8 6.1
110 57.5 1.7198 1.084 31.8 8.3
130 137 1.7385 1.064 31.8 7.8
150 88 1.7626 1.064 31.8 7.9
170 33.5 1.7902 1.064 31.8 8

caizi2008 发表于 2007-12-26 18:19

不一定要用这个软件,用别的也可以解决
MATLAB
关键的问题是解决问题 不是使用软件
可以自己尝试嘛

dingd 发表于 2007-12-28 21:55

今年的问题不能拖到明年啊!不过总让人这样帮忙,毕竟不是长久之计,自己也不方便啊!

用2.5算的:

方差(RMSE): 0.043131313231802
残差平方和(SSE): 0.178589777385582
相关系数(R): 0.999998165832538
相关系数之平方(R^2): 0.999996331668439
决定系数(DC): 0.999996331584251
F统计(F-Statistic): 3165.81860689619
参数         最佳估算
-------------------- -------------
p1          4.53102163379877
p2          0.0128000212992281
p3          0.236500286993358
p4          -0.0223014256617568
p5          0.0677072161003616
p6          0.0361432354334591
p7          1.94999698881693
p8          50.0000275970557
p9          3.43047237461055E-8
p10         1.10000387160548E-8

====== 输出结果 =====
文件:Data file - 1
No 实测y1 计算y1 实测y2 计算y2
1 30 30.0019530052142 7.8 7.74793843348618
2 30 30.0103776515397 7.9 7.87330398153311
3 30 30.0166857194024 8 7.99825133255548
4 30 30.0270810812806 6.1 6.07271132508011
5 30 30.0330956974303 8.3 8.40778104725706
6 30 30.0437337972634 7.8 7.74793843348618
7 30 30.0558997754511 7.9 7.87330398153311
8 30 30.0652428929333 8 7.99825133255548
9 30.2 30.1801188136343 6.1 6.07271132508011
10 30.2 30.1874878390266 8.3 8.40778104725706
11 30.2 30.1950464563017 7.8 7.74793843348618
12 30.2 30.2053982296275 7.9 7.87330398153311
13 30.2 30.2148016092127 8 7.99825133255548
14 30.2 30.2242171520303 6.1 6.07271132508011
15 30.2 30.2327989088072 8.3 8.40778104725706
16 30.2 30.2648903719648 7.8 7.74793843348618
17 30.4 30.380095082201 7.9 7.87330398153311
18 30.4 30.387621742578 8 7.99825133255548
19 30.4 30.3937441622439 6.1 6.07271132508011
20 30.4 30.4012019980275 8.3 8.40778104725706
21 30.4 30.4106676227042 7.8 7.74793843348618
22 30.4 30.4198370625275 7.9 7.87330398153311
23 30.4 30.4277954011866 8 7.99825133255548
24 30.4 30.4396334651426 6.1 6.07271132508011
25 31.4 31.3632338970684 8.3 8.40778104725706
26 31.4 31.3655099703292 7.8 7.74793843348618
27 31.4 31.3727648946969 7.9 7.87330398153311
28 31.4 31.3834235507083 8 7.99825133255548
29 31.4 31.3849069054243 6.1 6.07271132508011
30 31.4 31.3901477245314 8.3 8.40778104725706
31 31.4 31.3960837706734 7.8 7.74793843348618
32 31.4 31.4047624186069 7.9 7.87330398153311
33 31.6 31.562204031708 8 7.99825133255548
34 31.6 31.5635160628499 6.1 6.07271132508011
35 31.6 31.5676899592855 8.3 8.40778104725706
36 31.6 31.5679766695456 7.8 7.74793843348618
37 31.6 31.5769450722216 7.9 7.87330398153311
38 31.6 31.5822830744979 8 7.99825133255548
39 31.6 31.5858695844704 6.1 6.07271132508011
40 31.6 31.5931656993047 8.3 8.40778104725706
41 31.8 31.7744407288368 7.8 7.74793843348618
42 31.8 31.7810022716229 7.9 7.87330398153311
43 31.8 31.784371998495 8 7.99825133255548
44 31.8 31.7857749313769 6.1 6.07271132508011
45 31.8 31.7920468685805 8.3 8.40778104725706
46 31.8 31.7924875459538 7.8 7.74793843348618
47 31.8 31.8004391223472 7.9 7.87330398153311
48 31.8 31.8015200206924 8 7.99825133255548

dingd 发表于 2007-12-28 22:07

这是拟合图

morrissey22 发表于 2008-1-3 07:38

楼上的朋友,非常感谢!大家good luck in 2008!:loveliness:
页: [1]
查看完整版本: 有1stopt 2.0版本的朋友请帮个忙,平安夜求助,好人好运