- 1、本文档共4页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
功率谱估计的MATLAB实现
实 验 功率谱估计
实验目的:
1、掌握最大熵谱估计的基本原理。最终预测误差掌握周期图谱估计的基本原理。掌握传统谱估计中直接法与间接法之间的关系。
复习快速傅里叶变换与离散傅里叶变换之间关系。clc;
N=32;
SNR=30;
fs=1;
t=1:N;
t=t/fs;
y=sin(2*pi*0.175*t)+sin(2*pi*0.20*t);
x = awgn(y,SNR);
M=1;
P(M)=0;
Rx(M)=0;
for n=1:N
P(M)=P(M)+(abs(x(n)))^2;
ef(1,n)=x(n);
eb(1,n)=x(n);
end
P(M)=P(M)/N;
Rx(M)=P(M);
M=2;
A=0;
D=0;
for n=M:N
A=A+ef(M-1,n)*eb(M-1,n-1);
D=D+(abs(ef(M-1,n)))^2+(abs(eb(M-1,n-1)))^2;
end
xishu=-2*A/D;
a(M-1,M-1)=-2*A/D;
P(M)=P(M-1)*(1-(abs(xishu))^2);
FPE(M-1)=P(M)*(N+M)/(N-M);
TH=FPE(M-1);
for n=M:N
ef(M,n)=ef(M-1,n)+xishu*eb(M-1,n-1);
eb(M,n)=eb(M-1,n-1)+xishu*ef(M-1,n);
end
M=M+1;
A=0;
D=0;
for n=M:N
A=A+ef(M-1,n)*eb(M-1,n-1);
D=D+(abs(ef(M-1,n)))^2+(abs(eb(M-1,n-1)))^2;
end
xishu=-2*A/D;
a(M-1,M-1)=-2*A/D;
P(M)=P(M-1)*(1-(abs(xishu))^2);
FPE(M-1)=P(M)*(N+M)/(N-M);
for m=1:M-2
a(M-1,m)=a(M-2,m)+xishu*a(M-2,M-1-m);
end
while FPE(M-1)TH
TH=FPE(M-1);
for n=M:N
ef(M,n)=ef(M-1,n)+xishu*eb(M-1,n-1);
eb(M,n)=eb(M-1,n-1)+xishu*ef(M-1,n);
end
M=M+1;
A=0;
D=0;
for n=M:N
A=A+ef(M-1,n)*eb(M-1,n-1);
D=D+(abs(ef(M-1,n)))^2+(abs(eb(M-1,n-1)))^2;
end
xishu=-2*A/D;
a(M-1,M-1)=-2*A/D;
P(M)=P(M-1)*(1-(abs(xishu))^2);
FPE(M-1)=P(M)*(N+M)/(N-M);
for m=1:M-2
a(M-1,m)=a(M-2,m)+xishu*a(M-2,M-1-m);
end
end
T=1/fs;sum1=0;
f=0.01:0.01:0.5;
for m=1:M-1;
sum1=sum1+a(M-1,m)*exp(-j*2*pi*m*f*T);
end
s1=(abs(1+sum1)).^2;
s=P(M)*T./s1;
plot(f,10*log10(s),k);
xlabel(f/fs);
ylabel(功率谱/dB);
2、周期图谱估计
clc;clear;
N=32;
SNR=30;
fs=1;
t=1:N;
t=t/fs;
y=sin(2*pi*0.175*t)+sin(2*pi*0.20*t);
x = awgn(y,SNR);
sum1=0;
f=0.05:0.01:0.5;
for m=1:N
sum1=sum1+x(m)*exp(-j*2*pi*m*f);
end
s=(abs(sum1)).^2/N;
plot(f,10*log10(s),k);
xlabel(f/fs);
ylabel(功率谱/dB);
实验结果:
1、最大熵法估计结果:
2、周期图法估计结果:
文档评论(0)