24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 1016  |  回复: 5

xulianjing

银虫 (小有名气)

[求助] 有1stopt的虫友给跑下程序,求!

数据x   y
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
拟合公式(1)y=a*exp(-x^2/b)+c*exp(-(x-7.578*22.5)^2/e)+c*exp(-(x+7.578*22.5)^2/e)+f/(x^2+g)+h/((x-7.578*22.5)^2+i)+h/((x+7.578*22.5)^2+i)       
(2)a*exp(-x^2/b)+c*exp(-(x-7.578*22.5)^2/e)+c*exp(-(x+7.578*22.5)^2/e)
(3)f/(x^2+g)+h/((x-7.578*22.5)^2+i)+h/((x+7.578*22.5)^2+i)       
用三个函数都分别拟合一下,得出参数


[ Last edited by fegg7502 on 2013-10-16 at 12:07 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xulianjing

银虫 (小有名气)

为什么没有回复
2楼2013-10-15 21:48:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
fegg7502: 金币+1, 3ks 2013-10-16 12:08:08
你自己先动手算一下,有问题再发上来。
3楼2013-10-16 10:44:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xulianjing

银虫 (小有名气)

引用回帖:
3楼: Originally posted by dingd at 2013-10-16 10:44:42
你自己先动手算一下,有问题再发上来。

主要是我没那软件,而matlab和origin拟合的非常不理想
4楼2013-10-16 12:08:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

★ ★
csgt0: 金币+2, 谢谢 2013-10-17 14:59:04
你的数据是3个峰,给出的三个公式似乎都很难得到好的结果。用下面3个高斯函数叠加拟合,效果不错,供参考:

y=y0+(A/(w*sqrt(PI/2)))*exp(-2*((x-xc)/w)^2)+(A1/(w*sqrt(PI/2)))*exp(-2*((x-xc1)/w1)^2)+(A2/(w*sqrt(PI/2)))*exp(-2*((x-xc2)/w2)^2);

均方差(RMSE):0.00370450035150142
残差平方和(SSE):0.0023878581766437
相关系数(R): 0.990502017066602
相关系数之平方(R^2): 0.981094245813007
确定系数(DC): 0.981094245813008
卡方系数(Chi-Square): 0.00125663151560826
F统计(F-Statistic): 963.845295035091

参数        最佳估算
----------        -------------
y0        0.927513993738507
a        2.64040820356053
w        12.4179148614587
xc        777.8268692624
a1        0.369134562035136
xc1        948.883750968836
w1        -23.9865268985633
a2        0.403125539975586
xc2        607.63176576071
w2        -22.9522975526765
有1stopt的虫友给跑下程序,求!
f2.jpg

5楼2013-10-16 13:06:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xulianjing

银虫 (小有名气)

引用回帖:
5楼: Originally posted by dingd at 2013-10-16 13:06:43
你的数据是3个峰,给出的三个公式似乎都很难得到好的结果。用下面3个高斯函数叠加拟合,效果不错,供参考:

y=y0+(A/(w*sqrt(PI/2)))*exp(-2*((x-xc)/w)^2)+(A1/(w*sqrt(PI/2)))*exp(-2*((x-xc1)/w1)^2)+(A2/(w* ...

第二个公式和第三个公式我用origin拟合了一下,现在最急需的是第一个公式的拟合,也是用matlab很不理想,老师那等着结果呢,不知咋办了
6楼2013-10-16 17:12:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xulianjing 的主题更新
信息提示
请填处理意见