数字信号处理_上机作业1.docVIP

目录 数字信号处理第一次上机实验报告 1 一、实验题目 1 二、实验目的 1 三、程序 1 四、运行结果 3 五、 结论 4 数字信号处理第一次上机实验报告 姓名：石星宇 学号一、实验题目 考虑由连续时间信号，求采样后的离散时间信号及量化后的数字信号的Matlab实现问题； 输出以下三种信号： 对采样， 对采样， 编写程序计算两序列线性卷积（不允许使用Matlab自带的卷积函数） 二、实验目的 通过实验熟悉Matlab在数字信号处理中的应用 熟悉数字信号处理中常用的函数以及语句 加深对采样定理的理解 加深对线性卷积计算过程的理解 三、程序 连续时间信号，求采样后的离散时间信号及量化后的数字信号的Matlab实现 设连续时间信号为 采样时间间隔为，量化步长为。 程序如下： dt=0.001;tf=6;t=0:dt:tf; xa=sqrt(t)+cos(t); T=0.5;n=0:tf/T; x=sqrt(n*T)+cos(n*T); deltax=0.5; xq=round(x/deltax)*deltax; subplot(121) plot(t,xa,:);hold on ,grid on plot(n*T,x,o) stem(n*T,xq,*);grid on legend(连续时间信号xa,离散时间信号x,数字信号xq) subplot(122) stairs(n*T,xq),grid on legend(数字信号采样保持) set(gcf,color,w) 输出给定的三种信号 t=linspace(0,10,5000);%连续自变量 tf=10; xa=sin(2.5*pi*t);%原始连续时间信号xa fs1=2.5;fs2=4;%两个采样率 n1=0:tf*fs1; n2=0:tf*fs2;%离散自变量向量 x1=sin(2.5*pi*n1*(1/fs1));%采样率为2.5Hz的离散时间信号x1 x2=sin(2.5*pi*n2*(1/fs2));%采样率为4Hz的离散时间信号x2 figure(1) plot(t,xa,:),hold on ,grid on stem(n1/fs1,x1,ro) stem(n2/fs2,x2,m*) legend(连续时间信号xa,采样信号x1（fs=2.5Hz）,采样信号x2（fs=4Hz）) title(连续信号及其采样信号) xlabel(t),ylabel(xa(t)、x1(n)、x2(n)) 编写程序计算两序列线性卷积 function w=my_conv(x,y) len=length(x)+length(y)-1; fn=zeros(1,len); gn=zeros(1,len); fn(length(y):end)=x; for m=1:length(y) g1(m)=y(length(y)-m+1); end gn(1:length(y))=g1; for n=1:len gn1=circshift(gn,[0,n-1]); w(n)=sum(fn.*gn1); end 四、运行结果 连续时间信号、离散时间信号及量化后的数字信号如右图所示。 三种给定的信号输出如下图所示。 计算两序列的线性卷积 设输入序列如下： 用编写的程序计算结果如下，用Matlab自带的卷积函数计算，并与之对比，结果是正确的。 结论 计算机无法严格地表示连续量和模拟量，就自变量来说，只能把取的很密。将自变量时间轴按采样周期进行分割，得到离散时间自变量。通常从零开始，取顺序增加的整数，直到终点，Matlab语言写成。 第一次上机实验报告 数字信号处理 专业： 信息对抗技术 学生姓名：石星指导教师： 吕 雁 完成时间：2015年3月24日 第 3 页 第 2 页