function difn=difn(A)
difn=[];
A=diff(A);
[r,l]=size(A);
if r==1&&l>=2
for i=1:l-1
difn=[difn A(i+1)/A(i)];
end
end
if r>1
for i=1:r-1
difn=[difn;A(i+1,:)./A(i,:)];
end
end
if (difn-min(difn))<=1
return
else
warning('该指数模型比值只差不小于1')
%报错函数 error('')
end
end
上文(difn)为使用指数模型的标准,下文为函数解析计算过程
PS:如果需要使用该函数,请将两个函数皆要复制,文件名分别为difn与time_analyze
%修正指数曲线预测
function [time_an,c]=time_analyze(A)%A数据矩阵,第一列为时间,第二列相应数据
%已根据散点图判断大致模型,fun为判断出的对应函数
%y_t=L-a*b^t
B=A;
difn(B);
[r,l]=size(A);
if r>=3
if rem(r,3)~=0
A(r,:)=[];
r=r-