原创力文档- 1、本文档共4页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
用burg法实现功率谱估计
用burg法实现功率谱估计
一.数据为：
实现功率谱估计的matlan源代码：
clear
Fs 500;
n 0:1/Fs:0.5;
w1 200*pi;A1 5;
w2 300*pi;A2 12;
xn A1* cos w1*n +j*sin w1*n +A2* j*sin w2*n +cos w2*n +randn size n ;
%xn A1*exp jw1n +A2*exp jw2n ;
subplot 211 ;
plot n,xn ;
xlabel n ;
ylabel xn ;
title xn A1*exp jw1n +A2*exp jw2n +e n ;
ymax_xn max xn +0.2;
ymin_xn min xn -0.2;
axis [0 0.5 ymin_xn ymax_xn] ;
p floor length xn /3 +2;
nfft 1024;
[xpsd,f] pburg xn,p,nfft,Fs ;
pmax max xpsd ;
xpsd xpsd/pmax;
xpsd 10*log10 xpsd ;
subplot 212 ;
plot f,xpsd ;
title Power Spectral estimate with burg ;
ylabel Power Spectral estimate dB ;
xlabel f Hz ;
grid on;
ymin_psd min xpsd -1;
ymax_psd max xpsd +1;
axis [0 Fs/2 ymin_psd ymax_psd] ;
实验结果：
下图依次是阶数为N/3，N/2
二.数据为：
Matlab源代码：
clear
clf
clc
N 100;
Fs 500;
%产生x1信号即AR信号
vn rands 1,N ;
xn zeros 1,N ;
xn 1 vn 1 ;xn 2 vn 2 ;
a1 0.78;
a2 -0.96;
for n 3:N xn n vn n +a1*xn n-1 +a2*xn n-2
end
subplot 211 ;
plot xn ;
title xn n noise n +a1*xn n-1 +a2*xn n-2 ;
p floor length xn /5 ;
nfft 1024;
[xpsd,f] pburg xn,p,nfft,Fs ;
pmax max xpsd ;
xpsd xpsd/pmax;
xpsd 10*log10 xpsd ;
subplot 212 ;
plot f,xpsd ;
title Power Spectral estimate with burg ;
ylabel Power Spectral estimate dB ;
xlabel f Hz ;
grid on;
ymin_psd min xpsd -1;
ymax_psd max xpsd +1;
axis [0 Fs/2 ymin_psd ymax_psd] ;
文档评论（0）