【预测模型】基于AR模型实现数据预测含Matlab源码

1 简介

​

​

2 部分代码

%% Levinson 递推算法
% 高国兴 由陆光华 p110
%
% aa 的长度为这这+1
function [aa, ee] = fLevinson(f,p)
%初始化
sig=zeros(1,p);
sig(1+0)=f(1+0);% 从0开始比公式中的下标大1
a=zeros(p,p);%下标一致
%初值
a(1,1)=-f(1+1)/f(1+0);
sig(1+1)=(1-a(1,1)*conj(a(1,1)))*f(1+0);
for(k=2:p)
sum_af=0;
for(L=1:k-1)
sum_af=sum_af+a(k-1,L)*f(1+k-L);
end
a(k,k)=-(f(1+k)+sum_af)/sig(1+k-1);
for i=1:k-1
a(k,i)=a(k-1,i)+a(k,k)*a(k-1,k-i);
end
sig(1+k)=(1+a(k,k)*conj(a(k,k)))*sig(1+k-1);
end
aa=zeros(1,p);
aa(1)=1;
for(i=2:p+1)
aa(i)=a(p,i-1);
end
ee=sig(1+p);

3 仿真结果

​

4 参考文献

[1]周渊, 王炳和. 基于AR模型法功率谱估计的Matlab实现[J]. 武警工程学院学报, 2005, 021(006):73-76.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。