离散时间信号与系统的频域分析.docx

计算机与信息工程学院实验报告

专业：通信工程 年级/班级：2012级通信工程 2013—2014学年第二学期
课程名称

指导教师

本组成员
学号姓名

实验地点

实验时间

项目名称
离散时间信号与系统的频域分析
实验类型

一、 实验目的
1、掌握离散时间信号与系统的频域分析方法，从频域的角度对信号与系统的特性进行分析。
2、掌握离散时间信号傅里叶变换与傅里叶逆变换的实现方法。
3、掌握离散时间傅里叶变换的特点及应用
4、掌握离散时间傅里叶变换的数值计算方法及绘制信号频谱的方法
二、 实验仪器或设备
一台装有MATLAB的计算机
三、 实验原理
1. 离散时间系统的频率特性
在离散LTI 系统时域分析中得到系统的单位冲激响应可以完全表征系统，进而通过h[n]特性来分析系统的特性。系统单位冲激响应h[n]的傅里叶变换H () 成为LTI 系统的频率响应。与连续时间LTI 系统类似，通过系统频率响应可以分析出系统频率特性。与系统单位冲激响应h[n]一样，系统的频率响应H ( ) 反映了系统内在的固有特性，它取决于系统自身的结构及组成系统元件的参数，与外部激励无关，是描述系统特性的一个重要参数，H () 是频率的复函数可以表示为

其中，|1随频率变化的规律称为幅频特性；ϕ(ω )随频率变化的规律称为相频特性。
2. 离散时间信号傅里叶变换的数值计算方法
算法原理，
由傅里叶变换原理可知：
序列 f [n]的离散时间傅里叶变换F是ω 的连续函数。由于数据在 matlab 中以向量的形式存在，F ()只能在一个给定的离散频率的集合中计算。然而，

只有类似
形式的e− jω 的有理函数，才能计算其离散时间傅里叶变换。
四、 实验内容
1 离散时间傅里叶变换
(1) 下面参考程序是如下序列在范围−4π ≤ω ≤ 4π 的离散时间傅里叶变换

实验代码
%计算离散时间傅里叶变换的频率样本
clear all;
w=-4*pi:8*pi/511:4*pi;
num=[2 1]; den=[1 -0.6];
h=freqz(num,den,w);
subplot(2,1,1)
plot(w/pi,real(h)); grid;
title(‘实部’)
xlabel(‘omega/\pi’);
ylabel(‘振幅’);
subplot(2,1,2)
plot(w/pi, imag(h)); grid;
title(‘虚部’)
xlabel(‘omega/\pi’);
ylabel(‘振幅’);
figure;
subplot(2,1,1)
plot(w/pi, abs(h)); grid;
title(‘幅度谱’)
xlabel(‘omega/\pi’);
ylabel(‘振幅’);
subplot(2,1,2)
plot(w/pi, angle (h)); grid;
title(‘相位谱’)
xlabel(‘omega/\pi’);
ylabel(‘以弧度为单位的相位’);

修改程序，在范围 0≤ω≤π内计算如下有限长序列的离散时间傅里叶变换h[n]=[1 2 3 4 5 6 7 8 9]
clear all;
w=0*pi:8*pi/511:pi;
h=freqz(h[n]);
subplot(2,1,1)
plot(w/pi,real(h)); grid;
title(‘实部’)
xlabel(‘omega/\pi’);
ylabel(‘振幅’);
subplot(2,1,2)
plot(w/pi, imag(h)); grid;
title(‘虚部’)
xlabel(‘omega/\pi’);
ylabel(‘振幅’);
figure;
subplot(2,1,1)
plot(w/pi, abs(h)); grid;
title(‘幅度谱’)
xlabel(‘omega/\pi’);
ylabel(‘振幅’);
subplot(2,1,2)
plot(w/pi, angle (h)); grid;
title(‘相位谱’)
xlabel(‘omega/\pi’);
ylabel(‘以弧度为单位的相位’);

教师签名：

年 月 日
河南师范大学计算机与信息工程学院