function?[a?E]?=?computeARpara（x?p）
%?函數computeARpara的功能是計算AR模型功率譜估計的各參數
%?輸入參數：輸入信號x；AR模型階次p
%?輸出參數：AR模型系數a；功率預測誤差E
N?=?length（x）;???????????????????????%?輸入信號長度
%%?初始值
ef?=?x;??????????????????????????????%?前向預測誤差
eb?=?x;??????????????????????????????%?后向預測誤差
a=?1;????????????????????????????????%?初始AR模型系數
E=?x‘*x/N;???????????????????????????%?初始功率預測誤差
k=?zeros（1?p）;??????????????????????%?為反射系數預分配空間，提高循環速度
E=?[E?k];????????????????????????????%?為誤差預分配空間，提高速度
%%?burg算法
for?m?=?1:p??????????????????????????%?計算1~p+1階的模型系數和功率預測誤差
?efm?=?ef（2:end）;????????????????????%?前一階次的前向預測誤差
?ebm?=?eb（1:end?-?1）;????????????????%?前一階次的后向預測誤差
?num?=?-2.*ebm‘*efm;?????????????????%?反射系數的分子項
?den?=?efm‘*efm?+?ebm‘*ebm;??????????%?反射系數的分母項
?k（m）?=?num./den;????????????????????%?當前階次的

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-05-06?22:43??AR模型功率譜估計（Burg算法）\
?????文件????????1238??2018-04-30?13:02??AR模型功率譜估計（Burg算法）\computeARpara.m
?????文件????????2528??2018-04-30?14:57??AR模型功率譜估計（Burg算法）\myBurg.m
?????文件?????????598??2018-04-30?17:28??AR模型功率譜估計（Burg算法）\myBurg_main.m