Logistic模型求解怎么用matlab求解啊？

懸賞分：100 - 解決時間：2008-11-17 23:09

已知

x=0:1:12

y=[43.65 109.86 187.21 312.67 496.58 707.65 960.25 1238.75 1560.00 1824.29 2199.00 2438.89 2737.71]

y=L/(1+a*exp(-k*x))

利用線性回歸模型所得到的a和k的估計值和L=3000作為Logistic模型的擬合初值，對Logistic模型做非線性回歸。

據(jù)說用lsqcurvefit求解，本人非專業(yè)人士不懂啊？只需要列出式子！！拜托了。

問題補充：是不是得分成每三組數(shù)一組啊？茫然！！還要擬合Gompertz模型？有會的給加分啊！！

提問者：liziqiangde - 江湖新秀五級最佳答案

x=0:1:12

y=[43.65 109.86 187.21 312.67 496.58 707.65 960.25 1238.75 1560.00 1824.29 2199.00 2438.89 2737.71]

y=L/(1+a*exp(-k*x))

利用線性回歸模型所得到的a和k的估計值和L=3000作為Logistic模型的擬合初值，對Logistic模型做非線性回歸。

%第一步,線性回歸模型得到a,k

%這里假定y=a*exp(k*x),對兩邊取ln(Matlab中，ln用log函數(shù)表示),有

%lny=lna+k*x

%即logy是x的線性函數(shù)，斜率為k*loge,截距為loga

x=0:1:12 ;

y=[43.65 109.86 187.21 312.67 496.58 707.65 960.25 1238.75 1560.00 1824.29 2199.00 2438.89 2737.71] ;

line_A=polyfit(x,log(y),1);

k=line_A(1);

a=exp(line_A(2));

plot(x,y,'*',x,a*exp(k*x))

title('線性回歸的參數(shù)曲線與已經(jīng)點的關(guān)系')

%第二步，Logistic模型

%在Matlab下輸入：edit，然后將下面兩行百分號之間的內(nèi)容，復(fù)制進去，保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function y=zhidao_liziqiangde(A,x)

%其中k=A(1),a=A(2)

k=A(1);

a=A(2);

L=3000;

y=L./(1+a*exp(-k*x));

總結(jié)

以上是生活随笔為你收集整理的Gompertz模型绘图 matlab,Logistic模型matlab求解的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。