数学建模社区-数学中国

标题: 一个非线性曲线拟合问题的请教 [打印本页]
作者: legan    时间: 2009-8-26 15:40
标题: 一个非线性曲线拟合问题的请教
我们的问题是想用非线性最小二乘拟合出人口增长模型中Logistic的函数8 X: E* ?: M+ T( [8 y
x（t）=xm/(1+(xm/x0-1)*exp(-rt))求出环境允许的最大值xm和增长率r；2 W6 |* h+ ]& Z* c' f1 k
我们用的matlab中的lsqcurvefit和lsqnonlin函数出现了问题，我们的函数如下：
) W4 O0 Z5 |9 C/ ~: J& ?一、用lsqcurvefit函数求解
  t7 W& x8 G3 R& l8 J%定义flogistic1函数并保存! ~* A  B. h: ^, E7 v5 z; g! g
function f=flogistic1(x,tdata)
f=x(1)/(1+(x(1)/243-1)*exp(-x(2)*t));
%主程序
5 ?* Z( T0 t' u& e3 a# P
tdata=2000:2009;
ydata=[243
0 t0 P) B. Y% M. i: S  a3 i289 324 372 393 417 435 549 10869 Z( J. y) L; X( c5 h2 }
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函数并保存' i* B7 }/ L1 g' c; C6 f
function f=flogistic2(x)
tdata=2000:1:2009;
ydata=[243
% f9 J! x6 G1 B; R289 324 372 393 417 435 549 1086 1263];
f=ydata-x(1)/(1+(x(1)/243-1)*exp(-x(2)*tdata))
%主程序# r! r+ Q& J! `( @$ X8 ^, A
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] = ...
作者: legan    时间: 2009-8-26 15:40
希望好心人能帮我们解决下，感激不尽
作者: 东方明珠-WDZYQ    时间: 2009-8-26 16:46
查看MATLAB的help啊
作者: jd_lim    时间: 2009-8-29 00:13
可以用nlinfit方法直接模拟试试1 V/ V# h2 i6 S" ?/ m" P$ l9 T3 C
非线性回归分析
" L# f3 `2 U5 r0 X4 s- Y前提是把数据整理好
. h; t3 n! J7 P[beta,r,j]=nlinfit(x,y,'function')5 |: Y6 C9 m# X% O$ B. y% x; t  V
nlintool(x,y,'function',置信区间)
作者: handsome8848    时间: 2009-8-30 19:33
拟合？求最大值？LZ到底要表达什么意思。。。
: y' T- X) P- @3 b( Y拟合还能求最大值？
作者: 杨帆    时间: 2009-8-30 20:04
谢谢楼主分享 。。。。
作者: 杨帆    时间: 2009-8-30 20:16
谢谢楼主分享 。。。。
作者: xieedemzitou    时间: 2009-9-3 12:52
4L有道理。。
作者: alair009    时间: 2012-1-26 11:32
恩，参考一下。。88450640562953597917111575064111625743336822603965473801748718741200358398380672



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5