- 1、本文档共10页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
现代数字信号处理仿真作业研讨
3.17
N=32;
noise=(randn(1,N)+j*rand(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(j*2*pi*f1*(0:N-1));
signal2=A2*exp(j*2*pi*f2*(0:N-1));
signal3=A3*exp(j*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);
subplot(2,2,1);
m=-31:31;
stem(m,real(r1));
title(实部);
axis([-31,31,-1500,2500]);
subplot(2,2,2);
m=-31:31;
stem(m,imag(r1));
title(虚部);
axis([-31,31,-2000,2000]);
subplot(2,2,3);
m=-31:31;
stem(m,real(r));
title(实部);
axis([-31,31,-1500,2500]);
subplot(2,2,4);
m=-31:31;
stem(m,imag(r));
title(虚部);
axis([-31,31,-2000,2000]);
N=256;
noise=(randn(1,N)+j*rand(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(j*2*pi*f1*(0:N-1));
signal2=A2*exp(j*2*pi*f2*(0:N-1));
signal3=A3*exp(j*2*pi*f3*(0:N-1));
un=signal1+signal2+signal3+noise;
NF=1024;
Pw_p=(abs(fft(un,NF)).^2)/NF;
Pw_p=10*log10(Pw_p);
Pw_p=Pw_p-max(Pw_p);
figure();
dw=2*pi/NF;
df=dw/(2*pi);
plot(-pi/(2*pi):df:pi/(2*pi)-df,fftshift(Pw_p));
title(周期图法);
xlabel(w/2π);
ylabel(归一化功率谱/dB);
M=64;
r=xcorr(un,M,biased);
NF=1024;
Pw_bt=abs(fft(r,NF));
Pw_bt=10*log10(Pw_bt);
Pw_bt=Pw_bt-max(Pw_bt);
figure();
dw=2*pi/length(Pw_bt);
df=dw/(2*pi);
plot(-pi/(2*pi):df:pi/(2*pi)-df,fftshift(Pw_bt));
title(BT法);
xlabel(w/2π);
ylabel(归一化功率谱/dB);
p=16;
r0=xcorr(un,p,biased);
r=r0(p+1:2*p+1);
a(1,1)=-r(2)/r(1);
sigma(1)=r(1)-(abs(r(2)^2))/r(1);
for m=2:p
k(m)=-(r(m+1)+sum(a(m-1,1:m-1).*r(m:-1:2)))/sigma(m-1);
a(m,m)=k(m);
for i=1:m-1
a(m,i)=a(m-1,i)+k(m)*conj(a(m-1,m-i));
end
sigma(m)=sigma(m-1)*(1-abs(k(m)^2));
end
NF=1024;
Pw_LD=sigma(p)./fftshift(abs(fft([1,a(p,:)],NF)).^2);
Pw_LD=10*log10(Pw_LD);
Pw_LD=Pw_LD-max(Pw_LD);
dw=2*pi/length(Pw_LD);
df=dw/(2*pi);
figure();
plot(-pi/(2*pi):df:pi/(2*pi)-df,Pw_LD);
titl
您可能关注的文档
- 物理沪科版八年级第二单元第二节研讨.ppt
- 物理波动光学测试题答案研讨.ppt
- 物理电路图专题训练研讨.ppt
- 物理白皮书答案汇总研讨.doc
- 物理科普阅读研讨.ppt
- 物理竞赛[热学]研讨.ppt
- 物理竞赛所有公式研讨.doc
- 物理第3章研讨.ppt
- 物理第三章教案研讨.doc
- 物理系201108031高得伟研讨.doc
- 苏少版二年级美术上册《第7课 相同图样排排排》教学设计.docx
- 人美 四年级美术上册《第13课 降落伞》教学设计.docx
- 第8课 青铜工艺之美 教学设计 冀美版(2024)美术六年级上册.docx
- 第5单元《观察物体》教学设计-2023-2024学年二年级上册数学人教版.docx
- 浙教版2023年 八年级下册 第1课走进人工智能 教学设计.docx
- 六年级下学期信息技术第7课过程编辑仔细做(教案).docx
- 部编版语文五年级下册-13.人物描写一组 第1课时.docx
- 初中校园欺凌班会教学有助于.docx
- 人音版一年级下册第7课 巧巧手 教学设计.docx
- 2024学年高中安全无小事警钟长鸣班会教学设计.docx
文档评论(0)