快捷通道

QQ登录

只需要一步，快速开始

  注册地址  找回密码

数学建模社区-数学中国

留学>>
返回列表 发新帖
查看: 8181|回复: 8
打印 上一主题 下一主题

一个非线性曲线拟合问题的请教

[复制链接]
字体大小: 正常 放大
legan 实名认证       

9

主题

4

听众

1376

积分

升级  37.6%

  • TA的每日心情
    开心
    2013-11-7 15:01

    • 签到天数: 13 天

    [LV.3]偶尔看看II

    新人进步奖

    跳转到指定楼层
    1#
    发表于 2009-8-26 15:40 |只看该作者 |正序浏览
    |招呼Ta 关注Ta
    我们的问题是想用非线性最小二乘拟合出人口增长模型中Logistic的函数  f% p, F1 |( }$ L
    x（t）=xm/(1+(xm/x0-1)*exp(-rt))求出环境允许的最大值xm和增长率r；  z. |8 w8 m$ W2 Z; @
    我们用的matlab中的lsqcurvefit和lsqnonlin函数出现了问题，我们的函数如下：
    * ]/ G3 L. x  z+ K一、用lsqcurvefit函数求解+ b% x4 n9 W5 r: S
    %定义flogistic1函数并保存
    0 n1 H' g* E# {4 R) H3 A6 R8 C' |7 F' \. O# p
    function f=flogistic1(x,tdata)
    f=x(1)/(1+(x(1)/243-1)*exp(-x(2)*t));
    %主程序  T6 g. m$ Y. ?2 U! s8 ]3 W8 ]
    tdata=2000:2009;
    ydata=[243" m7 c7 e  Y4 m4 Z2 c7 l+ k: X+ m
    289 324 372 393 417 435 549 1086/ Z" Z/ O' {+ q; I% T- }+ i
    1263];
    x0=[2000 0.5];
    x=lsqcurvefit('flogistic1',x0,tdata,ydata)
    f=flogistic1(x,tdata)
    执行后出现的问题是：
    ??? Error using ==> optim\private\lsqncommon
    User supplied function ==> flogistic1
    failed with the following error:
    Error using ==> mrdivide
    Matrix dimensions must agree.
    Error in ==> lsqcurvefit at 165
    [x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...
    二、用lsqnonlin函数求解
    %定义flogistic2函数并保存
    5 A  X) S6 \) h) R+ i
    function f=flogistic2(x)
    tdata=2000:1:2009;
    ydata=[243
    7 [* _, B: }9 ?5 Z2 H. j4 @; L, H289 324 372 393 417 435 549 1086 1263];
    f=ydata-x(1)/(1+(x(1)/243-1)*exp(-x(2)*tdata))
    %主程序& U2 ]6 T+ d  l! Z( N
    x0=[2000 0.5];
    x=lsqnonlin('flogistic2',x0)
    f=fogistic2(x)
    执行后出现的问题是
    ??? Error using ==> optim\private\lsqncommon
    User supplied function ==> flogistic2
    failed with the following error:
    Error using ==> mrdivide
    Matrix dimensions must agree.
    Error in ==> lsqnonlin at 163
    [x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    回复

    使用道具 举报

    alair009        
    头像被屏蔽

    0

    主题

    4

    听众

    361

    积分

    升级  20.33%

  • TA的每日心情
    郁闷
    2012-2-3 19:26

    • 签到天数: 5 天

    [LV.2]偶尔看看I
    9#
    发表于 2012-1-26 11:32 |只看该作者 |招呼Ta 关注Ta
    提示: 作者被禁止或删除 内容自动屏蔽
    回复

    使用道具 举报

    xieedemzitou        

    2

    主题

    4

    听众

    31

    积分

    升级  27.37%

  • TA的每日心情
    郁闷
    2012-9-24 01:55

    • 签到天数: 1 天

    [LV.1]初来乍到
    8#
    发表于 2009-9-3 12:52 |只看该作者 |招呼Ta 关注Ta
    4L有道理。。
    回复

    使用道具 举报

    杨帆 实名认证       

    4

    主题

    4

    听众

    626

    积分

    升级  6.5%

  • TA的每日心情
    奋斗
    2012-12-4 13:38

    • 签到天数: 314 天

    [LV.8]以坛为家I

    新人进步奖 最具活力勋章 发帖功臣

    群组: 数学趣味、游戏、IQ等

    群组: 数学建模

    群组: 我们一定会赢
    7#
    发表于 2009-8-30 20:16 |只看该作者 |招呼Ta 关注Ta
    谢谢楼主分享 。。。。
    回复

    使用道具 举报

    杨帆 实名认证       

    4

    主题

    4

    听众

    626

    积分

    升级  6.5%

  • TA的每日心情
    奋斗
    2012-12-4 13:38

    • 签到天数: 314 天

    [LV.8]以坛为家I

    新人进步奖 最具活力勋章 发帖功臣

    群组: 数学趣味、游戏、IQ等

    群组: 数学建模

    群组: 我们一定会赢
    6#
    发表于 2009-8-30 20:04 |只看该作者 |招呼Ta 关注Ta
    谢谢楼主分享 。。。。
    回复

    使用道具 举报

    handsome8848        

    0

    主题

    3

    听众

    95

    积分

    升级  94.74%

    该用户从未签到
    5#
    发表于 2009-8-30 19:33 |只看该作者 |招呼Ta 关注Ta
    拟合？求最大值？LZ到底要表达什么意思。。。
    # Y  [' F3 n9 L3 {2 B% A拟合还能求最大值？
    回复

    使用道具 举报

    jd_lim        

    3

    主题

    4

    听众

    52

    积分

    升级  49.47%

    该用户从未签到
    4#
    发表于 2009-8-29 00:13 |只看该作者 |招呼Ta 关注Ta
    可以用nlinfit方法直接模拟试试- O! G& {+ q& J) E4 x9 `
    非线性回归分析
    0 r: h) Q" U! R! j9 h' k3 S( a前提是把数据整理好6 k2 w9 ^( E0 [3 h" ~+ ?
    [beta,r,j]=nlinfit(x,y,'function')
    ( U5 O9 E- U3 |# I3 Lnlintool(x,y,'function',置信区间)
    回复

    使用道具 举报

    东方明珠-WDZYQ        

    16

    主题

    21

    听众

    3372

    积分

    无

    升级  45.73%

  • TA的每日心情
    开心
    2013-10-25 15:56

    • 签到天数: 4 天

    [LV.2]偶尔看看I

    自我介绍
    活泼，爱交友

    发帖功臣 新人进步奖

    群组: 数学建模

    群组: 数学趣味、游戏、IQ等

    群组: 我行我数

    群组: Mathematica研究小组

    群组: LINGO
    3#
    发表于 2009-8-26 16:46 |只看该作者 |招呼Ta 关注Ta
    查看MATLAB的help啊
    回复

    使用道具 举报

    legan 实名认证       

    9

    主题

    4

    听众

    1376

    积分

    升级  37.6%

  • TA的每日心情
    开心
    2013-11-7 15:01

    • 签到天数: 13 天

    [LV.3]偶尔看看II

    新人进步奖

    2#
    发表于 2009-8-26 15:40 |只看该作者 |招呼Ta 关注Ta
    希望好心人能帮我们解决下，感激不尽
    回复

    使用道具 举报

    返回列表 发新帖
    高级模式
    B Color Image Link Quote Code Smilies
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩

    • 电话咨询

    • 04714969085

    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问：王兆丰

    GMT+8, 2025-6-25 03:49 , Processed in 1.488812 second(s), 98 queries .

    回顶部