matlab指数增长和阻滞增长拟合代码

代码：T=1:13X=3.04，6.63，11.13，28.43，49.22，87.85，134.11，200.00，208.69，303.03，410.73，580.03C=0.40，0.97，2.31，4.45，7.63，10.11，19.86，33.44，38.03，75.07，112.08，152.39D=2.06，2.35，3.17，8.06，15.22，28.82，43.31，54.67，55.83，75.26，102.01，174.05，E=0.54，1.18，1.98，5.06，8.77，15.65，23.89，36.63，37.17，53.98，73.16，123.31，123.31%指数增长模型y=log(x)；cc=log(c)；DD=log(d)；ee=log(e)；A=polyfit(t，y，1)；H=polyfit(t，cc，1)；I=polyfit(t，dd，1)；K=polyfit(t，ee，1)；r=a(1)；R2=h(1)；R3=I(1)；R4=k(1)；x0=exp(a(2)；x1=x0 * exp(r * t)；c0=exp(h(2)；C1=c0 * exp(R2 * t)；d0=exp(I(2)；D1=d0 * exp(R3 * t)；E0=exp(k(2)；E1=E0 * exp(R4 * t)；%块增长模型F= (a，t) a (1)。/(1(a(1)/x(1)-1)* exp(-a(2)*(t-t(1)；F= (h，t) h (1)。/(1(h(1)/c(1)-1)* exp(-h(2)*(t-t(1)；F= (I，t) I (1)。/(1(I(1)/d(1)-1)* exp(-I(2)*(t-t(1)；F= (k，t) k (1)。/(1(k(1)/e(1)-1)* exp(-k(2)*(t-t(1)；A=lsqcurvefit(f，800 1，t，x)；H=lsqcurvefit(f，800 1，t，c)；I=lsqcurvefit(f，800 1，t，d)；K=lsqcurvefit(f，800 1，t，e)；figure(1)；Subplot(2，1，1)；Plot(t、x、o、t、x1、r:)。) :X2=f(a，t)；Plot (t、x、o、t、x1、r :)。t、x2、g *-)Legend(原始数据，指数增长模型，块增长模型，2)Xlabel期间Ylabel公共编号合计(仅限于)Title(信息类)；Figure(1)Subplot(2，1，2)；Plot(t、c、o、t、c1、r:)。) :C2=f(h，t)；Plot (t、c、o、t、C1、r :)。t、C2、g *-)Legend(原始数据，指数增长模型，块增长模型，2)Xlabel期间Ylabel公共编号合计(仅限于)标题(生活流)；Figure(2)Subplot(2，1，1)；Plot(t、d、o、t、d1、r:)。) :D2=f(i，t)；Plot (t、d、o、t、D1、r :)。t、D2、g *-)Legend(原始数据，指数增长模型，块增长模型，2)Xlabel期间Ylabel公共编号合计(仅限于)Title(娱乐类)；Figure(2)Subplot(2，1，2)；Plot(t、e、o、t、e1、r:)。) :E2=f(k，t)；Plot (t、e、o、t、E1、r :)。t、E2、g *-)Legend(原始数据，指数增长模型，块增长模型，2)Xlabel期间Ylabel公共编号合计(仅限于)Title(其他)；参数信息类生活流娱乐类其他类普通弧总量程序代码：x=1336033626；Y=3.04，6.63，11.13，28.43，49.22，87.85，134.11，200.00，208.69，303.03，410.73，580.03B=2.04，2.14，3.68，10.86，17.60，33.27，47.05，75.26，77.66，88.72，123.48，130.25C=0.40，0.97，2.31，4.45，7.63，10.11，19.86，33.44，38.03，75.07，112.08，152.39D=2.06，2.35，3.17，8.06，15.22，28.82，43.31，54.67，55.83，75.26，102.01，174.05，E=0.54，1.18，1.98，5.06，8.77，15.65，23.89，36.63，37.17，53.98，73.16，123.31，123.31ST _=500 30 0.2；ft _=fittype(A0/(1b0 * exp(-k0 *(x-1)、Dependent、y、independent、x、Coefficients，a0，b0，k0 )；Cf_=fit(x、y、ft_、Startpoint、st_)Plot(cf_，fit，0.95)；hold on，plot (x，y，ro)f1t _=fittype(a1/(1b1 * exp(-k1 *(x-1)、Dependent、b、independent、x、Coefficients，a1，b1，k1 )；C1f _=fit (x、b、f1t _、startpoint、ST _)Plot(c1f_，fit，0.95)；plot，plot (x，b，g *)f2t _=fittype(a2/(1b2 * exp(-k2 *(x-1)、Dependent、c、independent、x、Coefficients，a2，b2，k2 )；C2f _=fit (x、c、f2t _、startpoint、ST _)Plot(c2f_，fit，0.95)；plot，plot (x，c，b)f3t _=fittype(a3/(1b3 * exp(-k3 *(x-1)、Dependent、d、independent、x、Coefficients，a3，b3，k3 )；C3f _=fit (x、d、f3t _、startpoint、ST _)Plot(c3f_，fit，0.95)；hold on，plot (x，d，*)f4t _=fittype(a4/(1b4 * exp(-k