AR模型算法的matlab实现和实验分析

1、 计算真实的自相关值时，采用逆Levinson-Durbin递归方法，由a、b参数得到 ， ， ， ，其中 为滤波器的阶数，再采用公式 外推得到 的自相关值； 2、 实际功率谱 ， 可调用Matlab中的FFT算法得到； 3、 自相关序列的估计值采用公式 得到； 4、 采用各种功率谱估计方法对功率谱进行估计。

function [aerr] = acm（xp）
%ACM	Find an all-pole model using the autocorrelation method
%----
%Usage: [aerr] = acm（xp） 
%
%	The input sequence x is modeled as the unit sample response of
%	a filter having a system function of the form
%		H（z） = b（0）/A（z） 
%	where the coefficients of A（z） are contained in the vector
%		a=[1 a（1） ... a（p）]
%	The input p defines the number of poles in the model.
%	The modeling error is returned in err.
%	The numerator b（0） is typically set equal to the square
%	root of err.
%
%  see also COVM PADE PRONY SHANKS
%
%---------------------------------------------------------------

x   = x（:）;
N   = length（x）;
if p>=length（x） error（‘Model order too large‘） end
X   = convm（xp+1）;
Xq  = X（1:N+p-11:p）;
a   = [1;-Xq\X（2:N+p1）];
err = abs（X（1:N+p1）‘*X*a）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      24879  2008-12-03 22:44  AR过程的线性建模与功率谱估计\程序\v1.mat

     文件        715  2008-11-18 20:35  AR过程的线性建模与功率谱估计\程序\v.mat

     文件        802  2002-11-24 02:27  AR过程的线性建模与功率谱估计\程序\acm.m

     文件        457  2002-10-08 22:43  AR过程的线性建模与功率谱估计\程序\atog.m

     文件        582  2002-10-08 22:42  AR过程的线性建模与功率谱估计\程序\ator.m

     文件        841  2008-11-19 11:27  AR过程的线性建模与功率谱估计\程序\burg.m

     文件        968  2002-10-08 22:43  AR过程的线性建模与功率谱估计\程序\convm.m

     文件        335  2002-10-08 22:44  AR过程的线性建模与功率谱估计\程序\covar.m

     文件        680  2008-11-19 10:43  AR过程的线性建模与功率谱估计\程序\covm.m

     文件        681  2002-10-08 22:45  AR过程的线性建模与功率谱估计\程序\durbin.m

     文件        112  2008-11-01 14:24  AR过程的线性建模与功率谱估计\程序\E_r.m

     文件        823  2002-10-08 22:45  AR过程的线性建模与功率谱估计\程序\fcov.m

     文件        591  2002-10-08 22:45  AR过程的线性建模与功率谱估计\程序\glev.m

     文件        435  2002-10-08 22:45  AR过程的线性建模与功率谱估计\程序\gtoa.m

     文件        666  2002-10-08 22:45  AR过程的线性建模与功率谱估计\程序\gtor.m

     文件        823  2008-11-19 10:42  AR过程的线性建模与功率谱估计\程序\mcov.m

     文件       3053  2008-11-18 20:37  AR过程的线性建模与功率谱估计\程序\r_e.mat

     文件       3053  2008-11-18 20:37  AR过程的线性建模与功率谱估计\程序\r_real.mat

     文件        638  2002-10-08 22:48  AR过程的线性建模与功率谱估计\程序\rtoa.m

     文件        608  2002-10-08 22:49  AR过程的线性建模与功率谱估计\程序\spike.m

     文件       1816  2008-12-03 20:35  AR过程的线性建模与功率谱估计\程序\test2.asv

     文件       4269  2009-01-06 00:01  AR过程的线性建模与功率谱估计\程序\test2.m

     文件       3294  2009-01-06 00:00  AR过程的线性建模与功率谱估计\程序\test2_var.m

     文件         61  2009-01-06 00:11  AR过程的线性建模与功率谱估计\程序\Read me.txt

     文件     188928  2009-02-25 16:43  AR过程的线性建模与功率谱估计\实验二.doc

     目录          0  2009-01-06 00:09  AR过程的线性建模与功率谱估计\程序

     目录          0  2009-01-05 22:47  AR过程的线性建模与功率谱估计

     文件       3543  2009-02-01 16:19  使用帮助：新手必看.htm

     文件        183  2009-02-02 09:47  Matlab中文论坛--助努力的人完成毕业设计.url

----------- ---------  ---------- -----  ----

............此处省略2个文件信息