《现代数字信号处理及其应用》仿真作业.doc

下载文档

298
1
约2.39万字
约 30页
2020-08-17 发布于浙江
举报
版权申诉
保障服务

《现代数字信号处理及其应用》仿真作业.doc

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