声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1548|回复: 3

[编程技巧] 请教 lsqcurvefit 不能得到最优结果的问题

[复制链接]
发表于 2007-11-27 16:36 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
  我用lsqcurvefit来拟合一个公式,可是就是不能得到最优结果。
恳请各位帮帮忙。
x=[0.9         0.85        0.8          0.77        0.7];
y=[7.8240   10.0858   11.8706   12.1859   13.0324];

f=inline('log(a(1)*(1-x/a(2)))/a(3)','a','x');
[a,jm]=lsqcurvefit(f,[70 1 1],x,y)

结果如下:
Maximum number of function evaluations exceeded;
   increase options.MaxFunEvals
a =
   70.7955    0.9834    0.2253

jm =
    0.3213

选取不同的初值都是一样得不到结果,请各位帮忙看看问题在哪?

应该加上什么条件,或者是由于别的原因?

[ 本帖最后由 eight 于 2007-11-27 16:48 编辑 ]
回复
分享到:

使用道具 举报

发表于 2007-11-27 16:40 | 显示全部楼层
数据点太少?---这个问题用1stOpt比较容易解决。
发表于 2007-11-27 18:31 | 显示全部楼层
显然是初值选择的问题,如下是lingo的结果:
   Objective value:                             0.1353339
   Model Title: :Fit Problem

                       Variable           Value        Reduced Cost
                              A        434.7426            0.000000
                              B       0.9347376           0.8659544E-08
                              C       0.3568502           0.2563307E-07
                          X( 1)       0.9000000            0.000000
                          X( 2)       0.8500000            0.000000
                          X( 3)       0.8000000            0.000000
                          X( 4)       0.7700000            0.000000
                          X( 5)       0.7000000            0.000000
                          Y( 1)        7.824000            0.000000
                          Y( 2)        10.08580            0.000000
                          Y( 3)        11.87060            0.000000
                          Y( 4)        12.18500            0.000000
                          Y( 5)        13.03240            0.000000

                            Row    Slack or Surplus      Dual Price
                              1       0.1353339           -1.000000

评分

1

查看全部评分

 楼主| 发表于 2007-11-27 22:59 | 显示全部楼层
谢谢各位的帮忙。结果非常不错。
再问一下,什么是  lingo ?  能说详细一点吗?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 05:00 , Processed in 0.056539 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表