怎么用matlab求解Logistic模型中的三个参数?
题目
怎么用matlab求解Logistic模型中的三个参数?
将Logistic模型的方程变为y=b/(1+a*exp(-kt)),已知t=[0,5,10,24,33,48,57,72,96,120,144,168,192,216];y=[0,0.028,0.103,0.336,0.450,0.597,0.716,0.778,0.835,0.849,0.816,0.839,0.811,0.816].怎么求出参数a,b,k的值,要具体的程序,
答案
建立m函数文件存为logistic1
function f=logistic1(b)
t=[0,5,10,24,33,48,57,72,96,120,144,168,192,216];y=[0,0.028,0.103,0.336,0.450,0.597,0.716,0.778,0.835,0.849,0.816,0.839,0.811,0.816];
f = y-b(1)./(1+b(2).*exp(-b(3).*t));
b0=[10,2,2];
>> b=leastsq('logistic1',b0)
b =
0.8221 13.9173 0.0818
或者cftool
General model:
f(x) = b/(1+a*exp(-k*x))
Coefficients (with 95% confidence bounds):
a = 13.92 (6.301,21.53)
b = 0.822 (0.7911,0.853)
k = 0.08184 (0.06479,0.0989)
Goodness of fit:
SSE:0.01404
R-square:0.9898
Adjusted R-square:0.9879
RMSE:0.03572
举一反三
已知函数f(x)=x,g(x)=alnx,a∈R.若曲线y=f(x)与曲线y=g(x)相交,且在交点处有相同的切线,求a的值和该切线方程.
我想写一篇关于奥巴马的演讲的文章,写哪一篇好呢?为什么好
最新试题
- 《村晚》描写了什么季节的风景
- FFT对信号进行频谱分析时,信号的频率的分辨率与什么有关?能否给出其数学关系?
- 如图,在平面直角坐标系xOy中,已知直线l1、l2,都经过点A(-4,0),
- 40乘55怎么用简便运算
- 能帮我把这几个英语单词连成句子吗?computer,the,help,summer camp,needs,with ..
- 已知△ABC为等腰三角形,底边BC=1,腰AB=2,求腰上的中线BM的长
- 求一篇猜测动物是什么的英语作文50词
- 克拉伯龙方程PV=nRT的字母各代表啥
- 若x2-7x+1=0,求x2+1/x2的值,我女儿的数学题,
- “兼山,艮,君子以思不出其位”什么意思?
热门考点
- 当X趋于0时,{sin[ln(1+2x)]-sin[ln(1-x)]}/x的极限等于?请求具体过程
- 连词成句1we,the,run,should,corridor,in,not,2,a,medicine,day.three,this,times,take
- He received a letter from his brother yesterday HE(两个单词的空)his brother yesterday 同义转换
- 双控开关有几种接法
- 3的2007次方乘以3的2008次方等于多少
- 坚持党的基本路线一百年不动摇与四项基本原则之间的关系 []]
- 设地球质量为m,地球半径为R,地球表面的重力加速度为g,问:当人造卫星离海平面3R时,它所受到的地球引力是在地面时的多少倍?
- 数学公式(初中)全集,急需!
- 这样的词语有吗力大如牛这样的词语有吗力大如牛
- 楚庄王好猎 翻译