数字信号处理实验报告2016分析.docVIP

实验四 线性时不变离散时间系统的频域分析 实验室名称：2304 实验时间:2016.10 姓 名： 学号: 专业:通信 指导教师： 成绩 教师签名： 年 月 日 一、实验目的 1. 通过MATLAB仿真一些简单的线性时不变离散时间系统，并研究它们的时域特性。 2. 通过MATLAB模拟滤波器，研究它们的性质。 3.通过使用matlab做实验来加强对传输函数的类型和频率响应和稳定性测试来强化理解概念。 二、实验内容 实验书习题:Q4.1, Q4.2 ,Q4.8，Q4.9，Q4.14,Q4.16，Q4.19，Q4.20 三、实验器材及软件 1. 笔记本电脑1台 2. MATLAB 7.0软件 实验原理 1.了解MATLAB中产生和绘制信号的基本命令。 2.序列以向量的形式储存，。 3.运用理论课知识求解问题。 五、实验步骤 根据实验题目要求进行分析； 2. 运用所学知识用MATLAB编程实现题目要求； 3. 对结果进行分析总结。 六、实验记录（数据、图表、波形、程序等） 4.1： clf; w = 0:8*pi/511:2*pi; M=input(M=); den = [1]; num=ones(1,M)/M; h=freqz(num,den,w); subplot(2,2,1) plot(w/pi,real(h)); grid title( H(e^{j\omega})的实部) xlabel(\omega /\pi); ylabel(振幅); subplot(2,2,2) plot(w/pi,imag(h)); grid title( H(e^{j\omega})的虚部) xlabel(\omega /\pi); ylabel(振幅); pause subplot(2,2,3) plot(w/pi,abs(h)); grid title(|H(e^{j\omega})|幅度谱) xlabel(\omega /\pi); ylabel(振幅); subplot(2,2,4) plot(w/pi,angle(h)); grid title(相位谱 arg[H(e^{j\omega})]) xlabel(\omega /\pi); ylabel(以弧度为单位的相位); M=5 M=9 M=20 幅度谱关于Y轴对称，相位谱关于原点对称。表示一种低通滤波器。 4.2： clf; w = 0:8*pi/511:pi; w=0:pi/511:pi; num=[0.15 0 -0.15]; den=[1 -0.5 0.7]; h=freqz(num,den,w); subplot(2,2,1) plot(w/pi,real(h)); grid title( H(e^{j\omega})的实部) xlabel(\omega /\pi); ylabel(振幅); subplot(2,2,2) plot(w/pi,imag(h)); grid title( H(e^{j\omega})的虚部) xlabel(\omega /\pi); ylabel(振幅); pause subplot(2,2,3) plot(w/pi,abs(h)); grid title(|H(e^{j\omega})|幅度谱) xlabel(\omega /\pi); ylabel(振幅); subplot(2,2,4) plot(w/pi,angle(h)); grid title(相位谱 arg[H(e^{j\omega})]) xlabel(\omega /\pi); ylabel(以弧度为单位的相位); 它表示一种带通滤波器。 4.8： clf; fc=0.45/(2*pi); n=[-8.5:1:8.5]; y=2*fc*sinc(2*fc*n); k=n+8.5; stem(k,y);title(N=20); axis([0 20 -0.2 0.6]); xlabel(时间序号?n); ylabel(振幅); grid 4.9： clf; fc=0.65/(2*pi); n=[-8.5:1:8.5]; y=2*fc*sinc(2*fc*n); k=n+8.5; stem(k,y); title(N=15); axis([0 15 -0.2 0.6]); xlabel(时间序号n); ylabel(振幅); grid 4.14： 我们得到的一阶IIR低通和高通滤波器的传递函数，分别给出了? clf; M = 2; n