离散时间系统的时域特性分析.docVIP

杭州电子科技大学通信工程学院 实验报告 课程名称：数字信号处理实验 实验名称：离散时间系统的时域特性分析 指导教师：何美霖 学生姓名：杨曦杰 学生学号学生班级学生专业：通信工程学院 实验日期：2022.10.9 一：实验目的 通过使用Matlab函数仿真简单得离散时间系统，研究其时域特性，以加深对离散时间系统差分方程，冲激响应和系统线性和时不变特性得理解 二：实验原理 三：预习与参考 1.所使用的主要函数 Filter函数 Conv函数 Impz函数 2．相关函数的应用实例 四：实验内容以及步骤 实验程序 clear; close all;clc B1=[0.5 0.27 0.77];A1=[1]; B2=[0.45 0.5 0.45];A2=[1 -0.53 0.46]; n=0:1:299; x=cos(20*pi*n/256)+cos(200*pi*n/256); y1=filter(B1,A1,x); y2=filter(B2,A2,x); h1=impz(B1,A1,n); h2=impz(B2,A2,n) D=5; xd=[zeros(1,D) x] yd=filter(B2,A2,xd); diff1=[zeros(1,D) y2]-yd; a=3.0;b=-2.0; x1=cos(2*pi*0.1*n); y21=filter(B2,A2,x1); x2=cos(2*pi*0.4*n); y22=filter(B2,A2,x2); x3=a*x1+b*x2; y23=filter(B2,A2,x3); y24=a*y21+b*y22; diff2=y24-y23; figure(NumberTitle, off, Name, ?μí3μ?ê?è?ê?3?2¨D?); subplot(311); stem(x,g,linewidth,1.0); title(?μí3ê?è?) subplot(312); stem(y1,b,linewidth,1.0); title(?μí31ê?3?) subplot(313); stem(y2,r,linewidth,1.0); title(?μí32ê?3?) figure(NumberTitle, off, Name, ?μí3μ?3??÷?ìó|DòáD2¨D?); subplot(211) stem(n,h1,fill,r,linewidth,1.0);hold on; title(?μí313??÷?ìó|); subplot(212) stem(n,h2,fill,b,linewidth,1.0);hold on; title(?μí323??÷?ìó|); figure(NumberTitle, off, Name, ?D???μí3ê±2?±?); subplot(311) stem(n,y2,fill,r,linewidth,1.0);hold on; title(?μí32ê?3?) subplot(312) stem(n,yd(1:300),fill,b,linewidth,1.0);hold on; title(?μí32ê±ò?ê?3?) subplot(313) stem(n,diff1(1:300),fill,g,linewidth,1.0);hold on; title(?μí32ê±ò?ê?3?2?) figure(NumberTitle, off, Name, ?D???μí3??D?); subplot(311); stem(n,y23,fill,r,linewidth,1.0);hold on; title(ê?3?1) subplot(312); stem(n,y24,fill,b,linewidth,1.0);hold on; title(ê?3?2) subplot(313); stem(n,diff2,fill,g,linewidth,1.0);hold on; title(ê?3?2?) 五：实验结果与数据处理、分析 系统2是一个线性时不变系统 根据图时不变中输入差是一条直线说明延迟后的变换等于变换后的延迟，所以系统2是时不变的。 根据图判断线性 代码中设定两个不同的输入信号，加权后分别输入系统2后相加，与加权相加后输入系统的输出作比较，可以看出输出差极小，位数为10的负十五次方，可以认为相等，说明系统2是线性 六：解答实验思考题 实验感想，这次实验虽然非常简单，但是却非常有意义，这次实验加深了我对离散时间系统差分方程，冲激响应和系统线性和时不变特性得理解，在编写的过程中要注意思路清晰，熟悉matlab相关函数的使用等 七：实验总结 这次实验让我学会了用matlab函数来表示系统差分