西工大软件工程信号与系统实验五.docx

西北工业大学 《信号与系统》实验报告 学 院： 软件与微电子学院 学　　号： 姓　　名： 专 业： 实验时间： 实验地点： 软件学院实验室310 指导教师： 陈勇 西北工业大学 2017年 9月 实验目的 了解离散时间全通系统； 了解离散时间傅里叶变换性质； 了解MATLAB的模拟电话播音键； 学会使用MATLAB函数residue； 了解离散时间系统的部分分式展开。 实验要求 运用理论课程所学的离散信号的傅里叶变换的相关知识解决实验中的问题； 按照题目要求完成相应问题。 实验设备（环境） 计算机； MATLAB R2017a； 《信号与系统计算机练习:利用MATLAB》(美)巴克 著，刘树棠 译，西安交通大学出版，2000.6； 实验内容与步骤 5.1 计算DTFT的样本 基本题（a）（b）（c）（d）（e） 5.2 电话拨号音 基本题（a）（b）（c） 5.3 离散时间全通系统 基本题（a）（b） 5.4 频率采样:基于DTFT的滤波器设计 基本题（a）（b）（c） 5.5 系统辨识 基本题（a）（b）（c）（d） 5.6 离散时间系统的部分分式展开 基本题（a）（b）（c）（d） 实验结果 5.1计算DTFT的样本 (a) 由题意可得： X 包含x[n]非零样本的向量x为x=[ones(1,10)]。 (b) x=[ones(1,10)]; k=[0:99]; w=2*pi*k/100; X=(1-exp(-j*10*2*pi*k/100))./(1-exp(-j*2*pi*k/100)); subplot(2,1,1); plot(w,abs(X)); xlabel(频率); ylabel(幅值); subplot(2,1,2); plot(w,angle(X)); xlabel(幅值); ylabel(相位); (c) 代码如下： 图像如下： 图与(b)中的结果相比较，是将(b)中的图k≥N/2的样本从π≤ω2π移动到?π≤ω0的区间上构成的。 (d) 因为x[n]=u[n]-u[n-10]， 所以x[n+5]=u[n+5]-u[n-5]； 则Xrjω 显然Xrjω 画Xr对于w的图代码如下： 图像如下： 验证代码如下： 验证图像如下： 综上x[n+5]的DTFT的求法正确 (e) z=[1 2 3 4 5 4 3 2 1]; k=[0:99]; w=2*pi*k/100; w=w-pi; Z=fftshift(fft(z,100)); plot(w,real(Z)); xlabel(w); ylabel(Z(jw)); 5.2电话拨号音 (a) 创建d0到d9的代码如下： n=[0:999]; d0=sin(0.7217*n)+sin(1.0247*n); d1=sin(0.5346*n)+sin(0.9273*n); d2=sin(0.5346*n)+sin(1.0247*n); d3=sin(0.5346*n)+sin(1.1328*n); d4=sin(0.5906*n)+sin(0.9273*n); d5=sin(0.5906*n)+sin(1.0247*n); d6=sin(0.5906*n)+sin(1.1328*n); d7=sin(0.6535*n)+sin(0.9273*n); d8=sin(0.6535*n)+sin(1.0247*n); d9=sin(0.6535*n)+sin(1.1328*n); sound(d0,8192); sound(d1,8192); sound(d2,8192); sound(d3,8192); sound(d4,8192); sound(d5,8192); sound(d6,8192); sound(d7,8192); sound(d8,8192); sound(d9,8192); (b) 代码如下： 图像如下： (c) space=[zeros(1,1000)]; phone=[d8 space d5 space d7 space d6 space d4 space d0 space d2 space d8]; sound(phone,8192); 听到的声音与电话机上的一致。 5.3离散时间全通系统 (a) a1=[1]; b1=[0 0 0 1]; freqz(b1,a1,1000); (b) a2=[1 -3/4]; b2=[-3/4 1]; freqz(b2,a2,1000); 显然H1(e 5.4频率采样：基于DTFT的滤波器设计 (a) 绘图代码如下： 图像如下： (b) k=[0:8]; w=2