时间序列分析之修正指数模型

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-