重庆大学信号与系统实验 信号与系统2.docVIP

PAGE  信号与系统 实验报告 学生姓名：董师宁 学号 重庆大学电气工程学院 2015年6月9日 实验题目 信号的抽样与重构实验时间2015-6-3实验地点DS1407实验成绩 实验性质□验证性 □设计性 □综合性教师评价： □思路正确 □内容完整 □方案过程正确 □实验结果正确 □报告规范 其他： 一、实验目的 1.学习在Matlab软件下信号的抽样与重构过程的计算机仿真； 2.学习掌握抽样定理，通过不同抽样频率的取样与重构，加深对定理的理解； 3.学习Matlab软件在信号处理中的应用。二、实验主要内容和要求 设研究的信号如下： 分别采用1KHz 、2KHz和5KHz进行采样，作出原信号和采样信号的图形； 利用重构公式对采样信号进行重构并画出图形，然后与原信号进行比较并给出结论； 分析信号间的频谱关系并解释重构的差异。 三、实验主要步骤或过程 程序代码： close all; clear all; clc; Dt=0.00001;t=-0.001:Dt:0.001; xa=exp(-500*abs(t)); Ts1=0.0002;n1=-25:1:25; Ts2=0.001;n2=-5:1:5; Ts3=0.0005;n3=-10:1:10; x1=exp(-500*abs(n1*Ts1)); x2=exp(-500*abs(n2*Ts2)); x3=exp(-500*abs(n3*Ts3)); figure(Name,?￡?aD?o?μ?è??ù,NumberTitle,off); subplot(4,1,1);plot(t*1000,xa); xlabel(t),ylabel(·ù?μ) title(?-D?o?) subplot(4,1,2);stem(n1*Ts1*1000,x1); xlabel(t),ylabel(·ù?μ) title(5000Hz) subplot(4,1,3);stem(n2*Ts2*1000,x2); xlabel(t),ylabel(·ù?μ) title(1000Hz) subplot(4,1,4);stem(n3*Ts3*1000,x3); xlabel(t),ylabel(·ù?μ) title(2000Hz) Fs1=1/Ts1;Fs2=1/Ts2;n1Ts1=n1*Ts1;n2Ts2=n2*Ts2;Fs3=1/Ts3;n3Ts3=n3*Ts3; xa1=x1*sinc(Fs1*(ones(length(n1Ts1),1)*t-n1Ts1*ones(1,length(t)))); error1=max(abs(xa-xa1)); xa2=x2*sinc(Fs2*(ones(length(n2Ts2),1)*t-n2Ts2*ones(1,length(t)))); error2=max(abs(xa-xa2)); xa3=x3*sinc(Fs3*(ones(length(n3Ts3),1)*t-n3Ts3*ones(1,length(t)))); error3=max(abs(xa-xa3)); figure(Name,D?o?μ???11,NumberTitle,off); subplot(4,1,1);plot(t*1000,xa); title(?-oˉêy) subplot(4,1,2);plot(t*1000,xa1); title(5000Hz) subplot(4,1,3);plot(t*1000,xa2); title(1000Hz) subplot(4,1,4);plot(t*1000,xa3); title(2000Hz) Wmax=2*pi*2000; K=500;k=0:1:K;W=k*Wmax/K; Xa=xa*exp(-j*t*W)*Dt Xa=real(Xa); W=[-fliplr(W),W(2:501)]; Xa=[fliplr(Xa),Xa(2:501)]; w=pi*k/K; X1=x1*exp(-j*n1*w); X2=x2*exp(-j*n2*w); X3=x3*exp(-j*n3*w); X1=real(X1); X2=real(X2); X3=real(X3); w=[-fliplr(w),w(2:K+1)]; X1=[fliplr(X1),X1(2:K+1)]; X2=[fliplr(X2),X2(2:K+1)]; X3=[fliplr(X3),X3(2:K+1)]; figure(Name,D?o?μ??μ?×,NumberTitle,of