24小时热门版块排行榜    

查看: 2683  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhangxiaodao

银虫 (小有名气)

[求助] 根据数据和方程的形式,求方程中的参数

RT, 数据如下列, 然后方程的形式为y=a*exp[-p/(x+q)],求方程中的p和q,谢谢各位大侠,急求啊~~
x=1.1
1.2
1.3
1.4
1.5
1.6
1.7
1.8
1.9
2
2.1
2.2
2.3
2.4
y=10363972.36
5697599.795
3269017.372
1948866.775
1202604.284
765578.5464
501320.0508
336809.4864
231640.3371
162754.7914
116618.904
85080.3473
63109.81473
47535.43773
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

CODE:
x=[1.1
1.2
1.3
1.4
1.5
1.6
1.7
1.8
1.9
2
2.1
2.2
2.3
2.4];
y=[10363972.36
5697599.795
3269017.372
1948866.775
1202604.284
765578.5464
501320.0508
336809.4864
231640.3371
162754.7914
116618.904
85080.3473
63109.81473
47535.43773];
st_ = [1000 -2 -1];
ft_ = fittype('a*exp(-p/(x+q))','dependent',{'y'},'independent',{'x'},'coefficients',{'a', 'p','q'});
[cf_,good]= fit(x,y,ft_ ,'Startpoint',st_)
h_ = plot(cf_,'fit',0.95);
legend off;  % turn off legend from plot method call
set(h_(1),'Color',[1 0 0],...
     'LineStyle','-', 'LineWidth',2,...
     'Marker','none', 'MarkerSize',6);
hold on,plot(x,y,'*')

cf_ =

     General model:
     cf_(x) = a*exp(-p/(x+q))
     Coefficients (with 95% confidence bounds):
       a =           1  (1, 1)
       p =         -42  (-42, -42)
       q =         1.5  (1.5, 1.5)

good =

           sse: 1.7672e-007
       rsquare: 1
           dfe: 11
    adjrsquare: 1
          rmse: 1.2675e-004

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

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
9楼2011-12-20 14:38:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

715211229

木虫 (正式写手)

★ ★
臭水沟(金币+2): 谢谢交流~~ 2011-12-20 11:03:15
matlab :ctool工具箱
%%%拟合效果很差

Fit computation did not converge:
Maximum number of function evaluations exceeded. Increasing
MaxFunEvals (in fit options) may allow for a better fit, or
the current equation may not be a good model for the data.

Fit found when optimization terminated:

General model:
       f(x) = a*exp(-p/(x+q))
Coefficients (with 95% confidence bounds):
       a =       11.28  (-2.505e+004, 2.507e+004)
       p =      -7.497  (-3191, 3176)
       q =     -0.4499  (-151.6, 150.7)

Goodness of fit:
  SSE: 1.309e+014
  R-square: -0.159
  Adjusted R-square: -0.3697
  RMSE: 3.45e+006


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

我是蜗牛
2楼2011-12-20 05:28:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangxiaodao

银虫 (小有名气)

送鲜花一朵
引用回帖:
楼: Originally posted by 715211229 at 2011-12-20 05:28:41:
matlab :ctool工具箱
%%%拟合效果很差

Fit computation did not converge:
Maximum number of function evaluations exceeded. Increasing
MaxFunEvals (in fit options) may allow  ...

谢谢
4楼2011-12-20 09:00:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

715211229

木虫 (正式写手)

★
臭水沟(金币+1): 谢谢交流~~ 2011-12-20 11:03:26
你看下这个帖子,我这个可能是错误的
http://muchong.com/bbs/viewthread.php?tid=3866180
我是蜗牛
5楼2011-12-20 10:48:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 14/700 2026-11-01 16:32 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 hKe7U1kad5Sg 2026-10-31 7/350 2026-11-01 15:08 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 15:01 by BCdvBAEpt3zm
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-01 14:50 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 14:50 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 12:40 by BCdvBAEpt3zm
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 12:27 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 12:11 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
请填处理意见