毕业设计求达人帮忙用matlab编写个verhulst模型的程序 [复制链接]

本人想要做一个有关于verhulst模型对于位移量测预报的程序，用matlab编写了相关的程序 7<*,O&![|  
自己编了一个 如下 Qp;FVUw9  
如输入x0=[.................................] V2SHF  
[#SiwhF|  
1{<r~  

Lng@'Yr  
function verhulst(x0)                   l`A
e&nc6  
n=length(x0);                     `< 82"cAT{  
x1=zeros(1,n);                     Hva/C{Y  
x1(1)=x0(1); =p[a Cb i  
for i=2:n                         2xBYJoF(  
    x1(i)=x1(i-1)+x0(i); z[qi~&7:v  
end :CsrcT=  
for  i=2:n                         M2nWvU$  
     y(i-1)=x0(i); 6 u1|pX8  
end goOw.~dZ'  
for i=1:n-1                       3ESrd"W=  
    c(i)=0.5*(x1(i+1)+x1(i)); YCq:]  
end Nwg?(h#  
for i=1:n-1                       3=("vR`!  
    d(i)=-0.25*(x1(i+1)+x1(i))*(x1(i+1)+x1(i)); n}dLfg*  
end XFv^jSF  
B=[c' d'];                         } 4ZWAzH  
au=inv(B'*B)*B'*y';             cR=o!2O  
for i=1:n+1;                       cAn_:^
 
    ago(i)=(au(1)/au(2))/(1+((au(1)/(au(2)*x0(1))-1))*exp(-au(1)*(i-1))); LvCX(yjZ*  
end iEx4va-j  
yc(1)=ago(1); 25;(`Td5  
for i=1:n                         y7Nd3\v [\  
    yc(i+1)=ago(i+1)-ago(i); ]wUH*\(y  
end 3$9V4v@2  
但是 运行起来 感觉总是不理想 偏差很大 希望MATLAB高人指点下 程序有没有问题   XG E.*aI  
如果 有 有建设性的意见的 一定掌声 鲜花 土币 感谢  急啊

可惜啊 我不会 敢问楼主这个问题解决了没呢？