原创力文档- 1、本文档共30页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
第三章仿真作业
3.17 (1)
N=32;
noise=(randn(1,N)+j*randn(1,N))/sqrt(2);
f1=0.15;
f2=0.17;
f3=0.26;
SNR1=30;
SNR2=30;
SNR3=27;
A1=10^(SNR1/20);
A2=10^(SNR2/20);
A3=10^(SNR3/20);
signal1=A1*exp(1i*2*pi*f1*(0:N-1));
signal2=A2*exp(1i*2*pi*f2*(0:N-1));
signal3=A3*exp(1i*2*pi*f3*(0:N-1));
un=signal1+signal2+signal3+noise;
Uk=fft(un,2*N);
Sk=(1/N)*abs(Uk).^2;
r0=ifft(Sk);
r1=[r0(N+2:2*N),r0(1:N)];
r=xcorr(un,N-1,biased);
r11=real(r1);
r12=imag(r1);
r1=real(r);
r2=imag(r);
m=1-N:N-1;
subplot(2,2,1);
stem(m,r11,o);
xlabel(m);
ylabel(实部);
title(基于FFT的自相关函数快速计算);
subplot(2,2,2);
stem(m,r12,o);
xlabel(m);
ylabel(虚部);
subplot(2,2,3);
stem(m,r1);
xlabel(m);
ylabel(实部);
title(教材式（3.1.2）估计的自相关函数);
subplot(2,2,4);
stem(m,r2);
xlabel(m);
ylabel(虚部);
（2）
N=256;
noise=(randn(1,N)+j*randn(1,N))/sqrt(2);
f1=0.15;
f2=0.17;
f3=0.26;
SNR1=30;
SNR2=30;
SNR3=27;
A1=10^(SNR1/20);
A2=10^(SNR2/20);
A3=10^(SNR3/20);
signal1=A1*exp(1i*2*pi*f1*(0:N-1));
signal2=A2*exp(1i*2*pi*f2*(0:N-1));
signal3=A3*exp(1i*2*pi*f3*(0:N-1));
un=signal1+signal2+signal3+noise;
NF=1024;
Spr=fftshift((1/NF)*abs(fft(un,NF)).^2);
Sprmax=max(Spr);
Spr=Spr/Sprmax;
f=(-(NF/2)+1:(NF/2))/NF;
plot(f,20*log(Spr));
xlabel(f);
ylabel(归一化功率谱/dB);
title(周期图法);
N=256;
noise=(randn(1,N)+j*randn(1,N))/sqrt(2);
f1=0.15;
f2=0.17;
f3=0.26;
SNR1=30;
SNR2=30;
SNR3=27;
A1=10^(SNR1/20);
A2=10^(SNR2/20);
A3=10^(SNR3/20);
signal1=A1*exp(1i*2*pi*f1*(0:N-1));
signal2=A2*exp(1i*2*pi*f2*(0:N-1));
signal3=A3*exp(1i*2*pi*f3*(0:N-1));
un=signal1+signal2+signal3+noise;
M=64;
r=xcorr(un,N-1,biased);
NF=1024;
BT=fftshift(fft(r,NF));
BTmax=max(BT);
BT=BT/BTmax;%归一化幅度
f=(-(NF/2)+1:(NF/2))/NF;
plot(f,20*log(BT));
xlabel(f);
ylabel(归一化功率谱/dB);
title(BT法);
（3）
N=256;
noise=(randn(1,N)+j*randn(1,N))/sqrt(2);
f1=0.15;
f2=0.17;
f3=0.26;
SNR1=30;
SNR2=30;
SNR3=27;
A1=10^(SNR1/20);
A2=10^(SNR2/20);
A3=10^(SNR3/20);
signal1=A1*exp(1i*2*pi*f1*(0:N-1));
signal2=A2*exp(1i*2*pi*f2*(0:N-1));
signal3=A3*exp(1i*2*pi*f3*(0:N-1));
un=signal1+signal2+signal3+nois
文档评论（0）