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

本人想要做一个有关于verhulst模型对于位移量测预报的程序，用matlab编写了相关的程序 r"^P>8  
自己编了一个 如下 1--5ok h  
如输入x0=[.................................] ZM0vB% M|  
s(-$|f+s  
'(=krM9;  
*r=6bpi  
function verhulst(x0)                   }{HlY?S  
n=length(x0);                     C&r&&Pw  
x1=zeros(1,n);                     F@[l&`7  
x1(1)=x0(1); $5(co)C  
for i=2:n                         pL
Nv\M+  
    x1(i)=x1(i-1)+x0(i); u`7\o~$  
end o[aRG7C  
for  i=2:n                         L -<!,CASW  
     y(i-1)=x0(i); Lk4gjs,V  
end pFhznH{0  
for i=1:n-1                       YDO#Q= q%  
    c(i)=0.5*(x1(i+1)+x1(i)); 0jG8G
mh!  
end kk& ([xqU  
for i=1:n-1                       #TXgV0\F  
    d(i)=-0.25*(x1(i+1)+x1(i))*(x1(i+1)+x1(i)); |}=xA%)  
end 'jqkDPn  
B=[c' d'];                         5!-'~W  
au=inv(B'*B)*B'*y';             w_PnEJa9  
for i=1:n+1;                       s@V4ny9x  
    ago(i)=(au(1)/au(2))/(1+((au(1)/(au(2)*x0(1))-1))*exp(-au(1)*(i-1))); 1)w^.8f  
end ??)IPRv?yF  
yc(1)=ago(1); sBt,y_LW  
for i=1:n                         ~}+F$&  
    yc(i+1)=ago(i+1)-ago(i); ]9F$/M#  
end K8daSvc  
但是 运行起来 感觉总是不理想 偏差很大 希望MATLAB高人指点下 程序有没有问题   AKu]c-  
如果 有 有建设性的意见的 一定掌声 鲜花 土币 感谢  急啊

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