24小时热门版块排行榜    

查看: 1517  |  回复: 5

xulianjing

银虫 (小有名气)

[求助] 这个编号后为什么1stopt没有反应,不能执行

这是写的代码
Title "Type your title here";
Parameters a,b,c,e,f,g,,h,i;
Variable x,y;
Function y=a*exp(-x^2/b)+c*exp(-(x-170.505)^2/e)+c*exp(-(x+170.505)^2/e)+f/(x^2+g)+h/((x-170.505)^2+i)+h/((x+170.505)^2+i)        ;
Data;
550          .9289575
552.4457          .9243243
554.8914          .9289575
557.337          .9289575
559.7826          .9277992
562.2283          .9266409
564.674          .9289575
567.1196          .9289575
572.0109          .923166
574.4565          .923166
576.9022          .9243243
579.3478          .9243243
581.7935          .9254826
584.2391          .9312741
586.6848          .9335907
589.1305          .9359073
591.5761          .9382239
594.0217          .9393822
596.4674          .9428571
598.9131          .9474903
601.3587          .9486486
603.8043          .9544401
606.25          .9544401
608.6957          .9532818
611.1414          .9521235
613.587          .9521235
616.0327          .9440154
618.4783          .9405406
620.924          .9382239
623.3696          .9382239
625.8152          .9393822
628.2609          .9335907
630.7065          .9324324
633.1522          .9312741
635.5978          .9312741
638.0435          .9312741
640.4891          .9301158
642.9348          .9277992
645.3804          .9289575
647.8261          .9301158
650.2717          .9301158
652.7174          .9277992
655.1631          .9266409
657.6087          .9289575
660.0544          .9301158
662.5          .9301158
664.9457          .9266409
667.3913          .9277992
669.837          .9289575
672.2827          .9266409
677.174          .9277992
679.6196          .9289575
682.0652          .9266409
684.5109          .9289575
689.4022          .9277992
691.8478          .9254826
694.2935          .9277992
696.7391          .9277992
699.1848          .9277992
704.0761          .923166
706.5217          .9254826
708.9674          .9289575
711.4131          .9208494
713.8587          .9208494
716.3044          .923166
718.75          .9243243
721.1957          .9277992
723.6413          .9289575
726.087          .9301158
728.5327          .9266409
730.9783          .9243243
733.424          .9277992
735.8696          .9289575
738.3152          .9301158
740.7609          .9254826
743.2065          .9254826
745.6522          .923166
748.0978          .923166
750.5435          .9324324
752.9891          .9324324
755.4348          .9359073
757.8804          .9416988
760.3261          .9428571
762.7717          .9509652
765.2174          .9486486
767.6631          .9683397
770.1087          .9972973
772.5544          1.043629
775          1.089961
777.4457          1.092278
779.8913          1.096911
782.337          1.050579
784.7827          1.012355
787.2283          .9833977
789.674          .9579151
792.1196          .9428571
794.5652          .9405406
797.0109          .9359073
799.4565          .9335907
804.3478          .9289575
806.7935          .9254826
809.2391          .9243243
811.6848          .923166
814.1304          .9243243
816.5761          .923166
819.0218          .9243243
821.4674          .923166
826.3587          .9266409
828.8044          .923166
833.6957          .9289575
836.1413          .9220077
838.587          .923166
841.0327          .923166
843.4783          .923166
845.924          .9266409
848.3696          .9266409
850.8152          .9266409
853.2609          .9277992
855.7065          .9243243
858.1522          .9243243
860.5978          .9277992
863.0435          .9277992
865.4891          .9277992
867.9348          .9243243
870.3804          .9208494
872.8261          .9220077
875.2718          .9266409
877.7174          .9266409
880.1631          .9289575
882.6087          .9289575
885.0544          .9289575
887.5          .9266409
889.9457          .9254826
892.3914          .9254826
894.837          .9266409
897.2827          .9277992
899.7283          .9301158
902.174          .9301158
904.6196          .9301158
907.0652          .9277992
909.5109          .9301158
911.9565          .9335907
914.4022          .9312741
916.8478          .9254826
919.2935          .9266409
921.7391          .9277992
924.1848          .9277992
929.0761          .9335907
931.5218          .9393822
933.9674          .9440154
936.4131          .9440154
938.8587          .9428571
941.3044          .9393822
943.75          .9451737
946.1957          .9474903
948.6414          .9544401
951.087          .9555985
953.5327          .9544401
955.9783          .9486486
958.424          .9440154
960.8696          .9405406
963.3152          .9335907
965.7609          .9335907
968.2065          .9312741
970.6522          .934749
973.0978          .9324324
975.5435          .9335907
977.9891          .9324324
980.4348          .9289575
982.8804          .9243243
985.3261          .9289575
990.2174          .9301158
995.1087          .9289575
997.5544          .9277992
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
csgt0: 金币+1, 谢谢 2013-10-18 09:36:01
“Parameters a,b,c,e,f,g,,h,i;”中g和h间多了个逗号。

» 本帖已获得的红花(最新10朵)

2楼2013-10-17 16:11:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xulianjing

银虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by dingd at 2013-10-17 16:11:17
“Parameters a,b,c,e,f,g,,h,i;”中g和h间多了个逗号。

但是你看这结果,最后y的计算值全为零,怎么办

"Type your title here"

迭代数: 21
计算用时(时:分:秒:微秒): 00:00:10:855
优化算法: 准牛顿法(BFGS) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.937506071323823
残差平方和(SSE): 152.931668275811
相关系数(R): 0
相关系数之平方(R^2): 0
决定系数(DC): -100
卡方系数(Chi-Square): 163.0586812
F统计(F-Statistic): 0

参数        最佳估算
----------        -------------
a                  0.900715064490214
b                  0.0980160213075578
c                  0.958751428639516
e                  0.259802284650505
f                  0.133029177086428
g                  0.494544279295951
h                  0.306610881118104
i                  0.777338797226548

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

No        实测值y        计算值y
1        0.9289575        0.0000000
2        0.9243243        0.0000000
3        0.9289575        0.0000000
4        0.9289575        0.0000000
5        0.9277992        0.0000000
6        0.9266409        0.0000000
7        0.9289575        0.0000000
8        0.9289575        0.0000000
9        0.923166        0.0000000
10        0.923166        0.0000000
11        0.9243243        0.0000000
12        0.9243243        0.0000000
13        0.9254826        0.0000000
14        0.9312741        0.0000000
15        0.9335907        0.0000000
16        0.9359073        0.0000000
17        0.9382239        0.0000000
18        0.9393822        0.0000000
19        0.9428571        0.0000000
20        0.9474903        0.0000000
21        0.9486486        0.0000000
22        0.9544401        0.0000000
23        0.9544401        0.0000000
24        0.9532818        0.0000000
25        0.9521235        0.0000000
26        0.9521235        0.0000000
27        0.9440154        0.0000000
28        0.9405406        0.0000000
29        0.9382239        0.0000000
30        0.9382239        0.0000000
31        0.9393822        0.0000000
32        0.9335907        0.0000000
33        0.9324324        0.0000000
34        0.9312741        0.0000000
35        0.9312741        0.0000000
36        0.9312741        0.0000000
37        0.9301158        0.0000000
38        0.9277992        0.0000000
39        0.9289575        0.0000000
40        0.9301158        0.0000000
41        0.9301158        0.0000000
42        0.9277992        0.0000000
43        0.9266409        0.0000000
44        0.9289575        0.0000000
45        0.9301158        0.0000000
46        0.9301158        0.0000000
47        0.9266409        0.0000000
48        0.9277992        0.0000000
49        0.9289575        0.0000000
50        0.9266409        0.0000000
51        0.9277992        0.0000000
52        0.9289575        0.0000000
53        0.9266409        0.0000000
54        0.9289575        0.0000000
55        0.9277992        0.0000000
56        0.9254826        0.0000000
57        0.9277992        0.0000000
58        0.9277992        0.0000000
59        0.9277992        0.0000000
60        0.923166        0.0000000
61        0.9254826        0.0000000
62        0.9289575        0.0000000
63        0.9208494        0.0000000
64        0.9208494        0.0000000
65        0.923166        0.0000000
66        0.9243243        0.0000000
67        0.9277992        0.0000000
68        0.9289575        0.0000000
69        0.9301158        0.0000000
70        0.9266409        0.0000000
71        0.9243243        0.0000000
72        0.9277992        0.0000000
73        0.9289575        0.0000000
74        0.9301158        0.0000000
75        0.9254826        0.0000000
76        0.9254826        0.0000000
77        0.923166        0.0000000
78        0.923166        0.0000000
79        0.9324324        0.0000000
80        0.9324324        0.0000000
81        0.9359073        0.0000000
82        0.9416988        0.0000000
83        0.9428571        0.0000000
84        0.9509652        0.0000000
85        0.9486486        0.0000000
86        0.9683397        0.0000000
87        0.9972973        0.0000000
88        1.043629        0.0000000
89        1.089961        0.0000000
90        1.092278        0.0000000
91        1.096911        0.0000000
92        1.050579        0.0000000
93        1.012355        0.0000000
94        0.9833977        0.0000000
95        0.9579151        0.0000000
96        0.9428571        0.0000000
97        0.9405406        0.0000000
98        0.9359073        0.0000000
99        0.9335907        0.0000000
100        0.9289575        0.0000000
101        0.9254826        0.0000000
102        0.9243243        0.0000000
103        0.923166        0.0000000
104        0.9243243        0.0000000
105        0.923166        0.0000000
106        0.9243243        0.0000000
107        0.923166        0.0000000
108        0.9266409        0.0000000
109        0.923166        0.0000000
110        0.9289575        0.0000000
111        0.9220077        0.0000000
112        0.923166        0.0000000
113        0.923166        0.0000000
114        0.923166        0.0000000
115        0.9266409        0.0000000
116        0.9266409        0.0000000
117        0.9266409        0.0000000
118        0.9277992        0.0000000
119        0.9243243        0.0000000
120        0.9243243        0.0000000
121        0.9277992        0.0000000
122        0.9277992        0.0000000
123        0.9277992        0.0000000
124        0.9243243        0.0000000
125        0.9208494        0.0000000
126        0.9220077        0.0000000
127        0.9266409        0.0000000
128        0.9266409        0.0000000
129        0.9289575        0.0000000
130        0.9289575        0.0000000
131        0.9289575        0.0000000
132        0.9266409        0.0000000
133        0.9254826        0.0000000
134        0.9254826        0.0000000
135        0.9266409        0.0000000
136        0.9277992        0.0000000
137        0.9301158        0.0000000
138        0.9301158        0.0000000
139        0.9301158        0.0000000
140        0.9277992        0.0000000
141        0.9301158        0.0000000
142        0.9335907        0.0000000
143        0.9312741        0.0000000
144        0.9254826        0.0000000
145        0.9266409        0.0000000
146        0.9277992        0.0000000
147        0.9277992        0.0000000
148        0.9335907        0.0000000
149        0.9393822        0.0000000
150        0.9440154        0.0000000
151        0.9440154        0.0000000
152        0.9428571        0.0000000
153        0.9393822        0.0000000
154        0.9451737        0.0000000
155        0.9474903        0.0000000
156        0.9544401        0.0000000
157        0.9555985        0.0000000
158        0.9544401        0.0000000
159        0.9486486        0.0000000
160        0.9440154        0.0000000
161        0.9405406        0.0000000
162        0.9335907        0.0000000
163        0.9335907        0.0000000
164        0.9312741        0.0000000
165        0.934749        0.0000000
166        0.9324324        0.0000000
167        0.9335907        0.0000000
168        0.9324324        0.0000000
169        0.9289575        0.0000000
170        0.9243243        0.0000000
171        0.9289575        0.0000000
172        0.9301158        0.0000000
173        0.9289575        0.0000000
174        0.9277992        0.0000000
这个编号后为什么1stopt没有反应,不能执行
1stopt.jpg

3楼2013-10-17 16:49:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★
xulianjing: 金币+5, 有帮助 2013-10-19 10:29:57
csgt0: 金币+1, 谢谢 2013-10-22 09:27:15
1:用5.0计算不会出现都是0的情况,但拟合效果不太好;
2:你的模型公式是固定的吗?公式中“c*exp(-(x+170.505)^2/e)”项和“h/((x-170.505)^2+i)”项分别出现了两次,确信正确吗?
3:另外一个帖子回复你用三个高斯函数拟合,效果不错。

» 本帖已获得的红花(最新10朵)

4楼2013-10-18 15:21:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xulianjing

银虫 (小有名气)

引用回帖:
4楼: Originally posted by dingd at 2013-10-18 15:21:04
1:用5.0计算不会出现都是0的情况,但拟合效果不太好;
2:你的模型公式是固定的吗?公式中“c*exp(-(x+170.505)^2/e)”项和“h/((x-170.505)^2+i)”项分别出现了两次,确信正确吗?
3:另外一个帖子回复你用三个 ...

1 能把拟合的图像和参数拟合值给我一下吗
2 “c*exp(-(x+170.505)^2/e)”项和“h/((x-170.505)^2+i)没有重复出现两次,一个是x+170.505,一个是x-170.505。
3  这是用三个函数分别拟合,进行比较的,看哪个公示拟合较好,并且参数值具有物理意义,所以要把参数值拟合出来
5楼2013-10-18 20:10:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xulianjing

银虫 (小有名气)

送红花一朵
引用回帖:
4楼: Originally posted by dingd at 2013-10-18 15:21:04
1:用5.0计算不会出现都是0的情况,但拟合效果不太好;
2:你的模型公式是固定的吗?公式中“c*exp(-(x+170.505)^2/e)”项和“h/((x-170.505)^2+i)”项分别出现了两次,确信正确吗?
3:另外一个帖子回复你用三个 ...

能把你用1stopt5得到的结果给我吗
6楼2013-10-19 15:19:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xulianjing 的主题更新
信息提示
请填处理意见