数字信号处理课程设计报告(2)

《数字信号处理课程设计报告(2)》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、通信与信息工程学院数字信号处理课程设计班级：电子信息科学与技术姓名：学号：指导教师：设计时间：成绩：评语：通信与信息工程学院二〇一三年题目一采样定理的验证1．设，利用filter函数求出（其中*为线性卷积）的前50个样本，并作图。（1）n=0:49;xn=0.8.^n;B=1;A=[1,-0.8];yn=filter(B,A,xn);stem(n,yn)gridon;xlabel('n');ylabel('yn=x(n)*x(n)')2．考虑模拟信号，间隔采样得到。（1）对每一个画出。（2）采用

2、内插从样本重建模拟信号，并从图中求出在中的频率（不管末端效果）。（3）采用三次样条内插从样本重建模拟信号，并从图中求出在中的频率。（1）clc;t=0:0.001:1;Ts1=0.01;N1=1/Ts1;n1=0:N1;xn1=sin(20*pi*n1*Ts1);subplot(3,1,1);stem(n1,xn1,'.');axis([0100-11]);xlabel('n');ylabel('x_1(n)');title('Ts=0.01');t=0:0.001:1;Ts2=0.05;N2=

3、1/Ts2;n2=0:N2;xn2=sin(20*pi*n2*Ts2);subplot(3,1,2);stem(n2,xn2,'.');axis([020-11]);xlabel('n');ylabel('x_2(n)');title('Ts=0.05');t=0:0.001:1;Ts3=0.1;N3=1/Ts3;n3=0:N3;xn3=sin(20*pi*n3*Ts3);subplot(3,1,3);stem(n3,xn3,'.');axis([010-11]);xlabel('n');yla

4、bel('x_3(n)');title('Ts=0.1');（2）t=0:0.001:1;T=0.1;xt=sin(20*pi*t);subplot(4,1,1);plot(t,xt);title('xt原图');xlabel('t');ylabel('xt');gridon;Ts1=0.01;n1=0:100;t1=0:Ts1:1;xn1=sin(20*pi*n1*Ts1);tn1=ones(length(n1),1)*t1-n1'*Ts1*ones(1,length(t1));%ones（M

5、,N)：生成M*N大小的全1矩阵yt1=xn1*sinc(tn1*pi/Ts1);subplot(4,1,2);plot(t1,yt1);axis([01-11]);title('sinc内插0.01恢复的xt1图');xlabel('n');ylabel('xt1');gridon;Ts2=0.05;n2=0:20;t2=0:Ts2:1;xn2=sin(20*pi*n2*Ts2);tn2=ones(length(n2),1)*t2-n2'*Ts2*ones(1,length(t2));yt2=

6、xn2*sinc(tn2*pi/Ts2);subplot(4,1,3);plot(t2,yt2);axis([01-11]);title('sinc内插0.05恢复的xt2图');xlabel('n');ylabel('xt2');gridon;Ts3=0.1;n3=0:10;t3=0:Ts3:1;xn3=sin(20*pi*n3*Ts3);tn3=ones(length(n3),1)*t3-n3'*Ts3*ones(1,length(t3));yt3=xn3*sinc(tn3*pi/Ts3);

7、subplot(4,1,4);plot(t3,yt3);axis([01-11]);title('sinc内插0.1恢复的xt3图');xlabel('n');ylabel('xt3');gridon;（3）Ts1=0.01;Fs=1/Ts1;xa1=spline(Ts1*n1,xn1,t);subplot(3,1,1);plot(t,xa1);axis([0,1,-1,1]);ylabel('x_a(t)');title('Ts=0.01');Ts2=0.05;Fs=1/Ts2;xa2=spl

8、ine(Ts2*n2,xn2,t);subplot(3,1,2);plot(t,xa2);axis([0,1,-1,1]);ylabel('x_a(t)');title('Ts=0.05');Ts3=0.1;Fs=1/Ts3;xa3=spline(Ts3*n3,xn3,t);subplot(3,1,3);plot(t,xa3);axis([0,1,-1,1]);ylabel('x_a(t)');title('Ts=0.1');结果分析与体会：采用内插从样本重建模拟信号时，0.01s和0.05s依