功率谱计算附matlab代码

%------------------------------------------------------------------------- % exa141202.m, for example 14.12.2; % to test aryule,arburg,arcov,armcov; %------------------------------------------------------------------------- clear all; load test x; % 用自相关法求AR模型的系数和最小预测误差能量； [a,E1]=aryule(x,10); figure(1) subplot(411) stem(abs(a),'.') ylabel(' for aryule'); % 用Burg法求AR模型的系数、反射系数和最小预测误差能量； [a,E2,k]=arburg(x,10); subplot(412) stem(abs(a),'.') ylabel(' for arburg'); % 用协方差法求AR模型的系数和最小预测误差能量； [a,E3]=arcov(x,10); subplot(413) stem(abs(a),'.') ylabel(' for cov'); % 用改进的协方差法求AR模型的系数和最小预测误差能量； [a,E4]=armcov(x,10); subplot(414) stem(abs(a),'.') ylabel(' for armcov'); figure(2) stem(abs(k),'.'); ylabel(' for arburg'); E1 E2 E3 E4