AR模型功率谱估计

文件为AR模型功率谱估计的MATLAB程序，有详细的注释。

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