电子信息工程系实验报告
课程名称： 数字信号处理
实验项目名称：实验1 离散时间信号与系统的傅里叶分析 时间：
2012-3-17
班级：电信092 姓名：XXX 学号：910706201
实 验 目 的:
用傅里叶变换对离散时间信号和系统进行频域分析。
实 验 环 境:
计算机、MATLAB 软件
实 验 原 理：
对信号进行频域分析即对信号进行傅里叶变换。对系统进行频域分析即对其单位脉冲响应进行傅里叶变
换，得到系统的传输函数；也可由差分方程经过傅里叶变换直接求其传输函数，传输函数代表的就是频率响应特性。而传输函数是w 的连续函数，计算机只能计算出有限个离散频率点的传输函数值，故可在0～2∏之间取许多点，计算这些点的传输函数的值，并取它们的包络，所得包络即所需的频率特性。
实 验 内 容 和 步 骤:
1、已知系统用下面差分方程描述：y (n )=x (n )+ay (n -1)，试在a =0.95和a =0.5 两种情况下用傅立叶变换分析系统的频率特性。要求写出系统的传输函数，并打印|H (e j ω)|~ω曲线。
解：B=1;A=[1,-0.95]; [H,w]=freqz(B,A,'whole');
subplot(1,3,1);plot(w/pi,abs(H),'linewidth',2);grid on;
xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('幅频响应特性');
axis([0,2,0,2.5]);
B=1;A=[1,-0.5];[H,w]=freqz(B,A,'whole');
subplot(1,3,3);plot(w/pi,abs(H),'linewidth',2);grid on;
xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('幅频响应特性');
axis([0,2,0,2.5]);
图形如下图1、2所示：
图1 a=0.95时的幅频响应特性 图2 a=0.5时的幅频响应特性
2、已知两系统分别用下面差分方程描述: y 1(n )=x (n )+x (n -1) y 2(n )=x (n )-x (n -1)
试分别写出它们的传输函数，并分别打印|H (e j ω)| ~ω曲线。
解：B=[1,1];A=1;[H,w]=freqz(B,A,'whole');
subplot(1,2,1);plot(w/pi,abs(H),'linewidth',2);grid on;
成 绩：
指导教师（签名）：
xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('幅频响应特性');
axis([0,2,0,2.2]);
B=[1,-1];A=1;[H,w]=freqz(B,A,'whole');
subplot(1,2,2);plot(w/pi,abs(H),'linewidth',2);grid on;
xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('幅频响应特性');
axis([0,2,0,2.2]);
图形如下图3、4所示：
图3 y1(n)的幅频特性图4 y2(n)的幅频特性图5 x(n)=R3(n)的幅频响应特性
3、已知信号x(n)=R3(n)，试分析它的频域特性，要求打印|X(e jω)|~ω曲线。
解：B=[1,0,0,-1];A=[1,-1];[H,w]=freqz(B,A,'whole');
plot(w/pi,abs(H),'linewidth',2);grid on;
xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('幅频响应特性');
axis([0,2,0,1.3]);
图形如上图5所示：
4、假设x(n)=δ(n)，将x(n)以2为周期进行周期延拓，得到x'(n)，试分析它的频率特性，并画出其幅频特性。解：n=-20:20; x=[ones(1,1),zeros(1,1)];
xtide=x(mod(n,2)+1); grid on;
subplot(1,2,1),stem(n,xtide,'.')
xlabel('n'),ylabel('xtide(n)');title('延拓后周期为2的周期序列')
B=1;A=1;[H,w]=freqz(B,A,'whole');
subplot(1,2,2);plot(w/pi,abs(H),'linewidth',2);grid on;
xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('幅频响应特性');axis([0,2,0,2]);
图形如下图6、7所示：
图6 x(n)以2为周期进行周期延拓所得序列图7 x(n) 周期延拓后的幅频响应特性
学生实验报告 （理工类） 课程名称：信号与线性系统专业班级：M11通信工程 学生学号：1121413017 学生姓名：王金龙 所属院部：龙蟠学院指导教师：杨娟
20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制 实验报告书写要求 实验报告原则上要求学生手写，要求书写工整。若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 （1）细致观察，及时、准确、如实记录。 （2）准确说明，层次清晰。 （3）尽量采用专用术语来说明事物。 （4）外文、符号、公式要准确，应使用统一规定的名词和符号。 （5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细，一律用红色笔批改。实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。 实验报告装订要求
实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。
实验项目名称：常用连续信号的表示 实验学时： 2学时 同组学生姓名： 无 实验地点： A207 实验日期： 11.12.6 实验成绩： 批改教师： 杨娟 批改时间： 一、实验目的和要求 熟悉MATLAB 软件；利用MATLAB 软件，绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机，装有MATLAB7.0软件 三、实验过程 1. 绘制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?=； 2. 绘制指数信号at Ae t (f =），其中A=1，0.4a -=； 3. 绘制矩形脉冲信号，脉冲宽度为2； 4. 绘制三角波脉冲信号，脉冲宽度为4；斜度为0.5； 5. 对上题三角波脉冲信号进行尺度变换，分别得出)2t (f ，)2t 2(f -； 6. 绘制抽样函数Sa （t ）,t 取值在-3π到+3π之间； 7. 绘制周期矩形脉冲信号，参数自定； 8. 绘制周期三角脉冲信号，参数自定。 四、实验结果与分析 1．制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?= 实验代码： A=1;
快速傅里叶变换实验报告
————————————————————————————————作者：————————————————————————————————日期: ?
快速傅里叶变换实验报告 机械34班 刘攀 2０1３0105５8 一、 基本信号（函数)的FF Ｔ变换 1. 000()sin()sin 2cos36x t t t t π ωωω=+++ 1) 采样频率08s f f =,截断长度N ＝１６； 取02ωπ=ｒad/ｓ，则0f =1Hz ，s f =8Ｈz ，频率分辨率 f ?=s f f N ?= =0.5Hz 。 最高频率c f ＝30f =３Hz ，s f ＞2c f ,故满足采样定理,不会发生混叠现象。 截断长度02T T =,整周期截取,不会发生栅栏效应。理论上有一定的泄漏,但在整周期 截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图如下:
幅值误差0A ?=,相位误差0??=。 2) 采样频率08s f f =，截断长度N=3２; 取02ωπ=rad/s ，则0f =１Hz,s f ＝８Hz ，频率分辨率f ?=s f f N ?==０.2５Hz 。 最高频率c f =30f =3H ｚ,s f >２c f ,故满足采样定理,不会发生混叠现象。 截断长度04T T =,整周期截取,不会发生栅栏效应。理论上有一定的泄漏，但在整周期 截取的情况下,旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图如下：
幅值误差0A ?=，相位误差0??=。 2. 00()sin()sin116x t t t π ωω=++ 1) 采样频率08s f f =,截断长度N=16； 取02ωπ=ra ｄ/s,则0f =１Hz ，s f =８Hz,频率分辨率f ?=s f f N ?==0．5H ｚ。 最高频率c f =110f =11H ｚ,s f <２c f ，故不满足采样定理，会发生混叠现象。 截断长度02T T =,整周期截取，不会发生栅栏效应。理论上有一定的泄漏,但在整周期 截取的情况下，旁瓣上的采样都约为 0，泄漏现象没有体现出来。 频谱图：
实验一信号与系统的时域分析 一、实验目的 1、熟悉与掌握常用的用于信号与系统时域仿真分析的MA TLAB函数; 2、掌握连续时间与离散时间信号的MA TLAB产生,掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MA TLAB编程; 3、牢固掌握系统的单位冲激响应的概念,掌握LTI系统的卷积表达式及其物理意义,掌握卷积的计算方法、卷积的基本性质; 4、掌握利用MA TLAB计算卷积的编程方法,并利用所编写的MA TLAB程序验证卷积的常用基本性质; 掌握MA TLAB描述LTI系统的常用方法及有关函数,并学会利用MATLAB求解LTI系统响应,绘制相应曲线。 基本要求:掌握用MA TLAB描述连续时间信号与离散时间信号的方法,能够编写MATLAB程序,实现各种信号的时域变换与运算,并且以图形的方式再现各种信号的波形。掌握线性时不变连续系统的时域数学模型用MA TLAB描述的方法,掌握卷积运算、线性常系数微分方程的求解编程。 二、实验原理 信号(Signal)一般都就是随某一个或某几个独立变量的变化而变化的,例如,温度、压力、 声音,还有股票市场的日收盘指数等,这些信号都就是随时间的变化而变化的,还有一些信号,例如在研究地球结构时,地下某处的密度就就是随着海拔高度的变化而变化的。一幅图片中的每一个象素点的位置取决于两个坐标轴,即横轴与纵轴,因此,图像信号具有两个或两个以上的独立变量。 在《信号与系统》课程中,我们只关注这种只有一个独立变量(Independent variable)的信号,并且把这个独立变量统称为时间变量(Time variable),不管这个独立变量就是否就是时间变量。 在自然界中,大多数信号的时间变量都就是连续变化的,因此这种信号被称为连续时间信号(Continuous-Time Signals)或模拟信号(Analog Signals),例如前面提到的温度、压力与声音 信号就就是连续时间信号的例子。但就是,还有一些信号的独立时间变量就是离散变化的,这种信号称为离散时间信号。前面提到的股票市场的日收盘指数,由于相邻两个交易日的日收盘指数相隔24小时,这意味着日收盘指数的时间变量就是不连续的,因此日收盘指数就是离散时间信号。 而系统则用于对信号进行运算或处理,或者从信号中提取有用的信息,或者滤出信号中某些无用的成分,如滤波,从而产生人们所希望的新的信号。系统通常就是由若干部件或单元组成的一个整体(Entity)。系统可分为很多不同的类型,例如,根据系统所处理的信号的不同,系统可分为连续时间系统(Continuous-time system)与离散时间系统(Discrete-time system),根
1 / 7 信息工程学院实验报告 课程名称：信号与系统 实验项目名称：实验 3 傅里叶变换及其性质实验时间： 2013-11-29 班级： 姓名：学号： 一、实验目的: 1、学会运用MATLAB 求连续时间信号的傅里叶（Fourier ）变换； 2、学会运用MATLAB 求连续时间信号的频谱图； 3、学会运用MATLAB 分析连续时间信号的傅里叶变换的性质。 二、实验环境: 1、硬件：在windows 7 操作环境下； 2、软件：Matlab 版本7.1 三、实验原理： 3.1傅里叶变换的实现 信号()f t 的傅里叶变换定义为：() [()] ()j t F F f t f t e dt ， 傅里叶反变换定义为： 1 1()[()] ()2 j t f t F F f e d 。 信号的傅里叶变换主要包括MATLAB 符号运算和MATLAB 数值分析两种方法，下面分别加以探讨。同时， 学习连续时间信号的频谱图。 3.1.1 MATLAB 符号运算求解法 MATLAB 符号数学工具箱提供了直接求解傅里叶变换与傅里叶反变换的函数fourier( )和ifourier( )。 Fourier 变换的语句格式分为三种。 （1）F=fourier(f)：它是符号函数 f 的Fourier 变换，默认返回是关于的函数。 （2）F=fourier(f,v) ：它返回函数F 是关于符号对象 v 的函数，而不是默认的 ，即 ()()j v t Fv fte d t 。 （3）F=fourier(f,u,v)：是对关于u 的函数f 进行变换，返回函数F 是关于v 的函数，即 ()()jvu F v f t e du 。 傅里叶反变换的语句格式也分为三种。（1）f=ifourier(F)：它是符号函数F 的Fourier 反变换，独立变量默认为 ，默认返回是关于 x 的函数。 （2）f=ifourier(F,u)：它返回函数 f 是u 的函数，而不是默认的 x 。 （3）f=ifourier(F,u,v) ：是对关于v 的函数F 进行反变换，返回关于 u 的函数f 。 成 绩： 指导教师（签名）：
电子信息工程系实验报告 课程名称： 数字信号处理 实验项目名称：实验1 离散时间信号与系统的傅里叶分析 时间： 2012-3-17 班级：电信092 姓名：XXX 学号：910706201 实 验 目 的: 用傅里叶变换对离散时间信号和系统进行频域分析。 实 验 环 境: 计算机、MATLAB 软件 实 验 原 理： 对信号进行频域分析即对信号进行傅里叶变换。对系统进行频域分析即对其单位脉冲响应进行傅里叶变 换，得到系统的传输函数；也可由差分方程经过傅里叶变换直接求其传输函数，传输函数代表的就是频率响应特性。而传输函数是w 的连续函数，计算机只能计算出有限个离散频率点的传输函数值，故可在0～2∏之间取许多点，计算这些点的传输函数的值，并取它们的包络，所得包络即所需的频率特性。 实 验 内 容 和 步 骤: 1、已知系统用下面差分方程描述：y (n )=x (n )+ay (n -1)，试在a =0.95和a =0.5 两种情况下用傅立叶变换分析系统的频率特性。要求写出系统的传输函数，并打印|H (e j ω)|~ω曲线。 解：B=1;A=[1,-0.95]; [H,w]=freqz(B,A,'whole'); subplot(1,3,1);plot(w/pi,abs(H),'linewidth',2);grid on; xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('幅频响应特性'); axis([0,2,0,2.5]); B=1;A=[1,-0.5];[H,w]=freqz(B,A,'whole'); subplot(1,3,3);plot(w/pi,abs(H),'linewidth',2);grid on; xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|');title('幅频响应特性'); axis([0,2,0,2.5]); 图形如下图1、2所示： 图1 a=0.95时的幅频响应特性 图2 a=0.5时的幅频响应特性 2、已知两系统分别用下面差分方程描述: y 1(n )=x (n )+x (n -1) y 2(n )=x (n )-x (n -1) 试分别写出它们的传输函数，并分别打印|H (e j ω)| ~ω曲线。 解：B=[1,1];A=1;[H,w]=freqz(B,A,'whole'); subplot(1,2,1);plot(w/pi,abs(H),'linewidth',2);grid on; 成 绩： 指导教师（签名）：
1.请用MATLAB编写程序，实现任意两个有限长度序列的卷积和。要求用图 形显示两个序列及卷积结果。 解：y(n)=∑x(i)h(n-i) 假设x(n)={1,2,3,4,5}; h(n)={3,6,7,2,1,6}; y(n)=x(n)*h(n) 验证：y[n]=[1,12,28,46,65,72,58,32,29,30] 【程序】 N=5 M=6 L=N+M-1 x=[1,2,3,4,5] h=[3,6,7,2,1,6] y=conv(x,h) nx=0:N-1 nh=0:M-1 ny=0:L-1 subplot(131);stem(nx,x,'*b');xlabel('n');ylabel('x(n)');grid on subplot(132);stem(nh,h,'*b');xlabel('n');ylabel('h(h)');grid on subplot(133);stem(ny,y,'*r');xlabel('n');ylabel('y(h)');grid on 【运行结果】
2.已知两个序列x[n]=cos(n*pi/2), y[n]=e j*pi*n/4x[n],请编写程序绘制 X(e jw)和Y(e jw)和幅度和相角,说明它们的频移关系。 –提示：用abs函数求幅度，用angle求相角。 【程序】 n=0:15; x=cos(n*pi/2); y=exp(j*pi*n/4).*x; X=fft(x); Y=fft(y); magX=abs(X); angX=angle(X); magY=abs(Y); angY=angle(Y); subplot(221);stem(n,magX,'*r');xlabel('频率');ylabel('幅度');grid on; subplot(222);stem(n,angX,'*b');xlabel('频率');ylabel('相位');grid on; subplot(223);stem(n,magY,'*r');xlabel('频率');ylabel('幅度');grid on; subplot(224);stem(n,angY,'*b');xlabel('频率');ylabel('相位');grid on;
实验一：信号的时域分析 一、实验目的 1.观察常用信号的波形特点及产生方法 2.学会使用示波器对常用波形参数的测量 二、实验仪器 1.信号与系统试验箱一台（型号ZH5004） 2.40MHz双踪示波器一台 3.DDS信号源一台 三、实验原理 对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(t)信号、钟形信号、脉冲信号等。 1、信号：指数信号可表示为f(t)=Ke at。对于不同的a取值，其波形表现为不同的形式，如下图所示： 图1―1 指数信号 2、信号：其表达式为f（t）=Ksin（ωt＋θ），其信号的参数：振幅K、角频率ω、与初始相位θ。其波形如下图所示：
图1－2 正弦信号 3、指数衰减正弦信号：其表达式为其波形如下图： 图1－3 指数衰减正弦信号 4、Sa(t)信号：其表达式为：。Sa(t)是一个偶函数，t= ±π,±2π,…,±nπ时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：
图1－4 Sa（t）信号 5、钟形信号（高斯函数）：其表达式为：其信号如下图所示： 图1－5 钟形信号 6、脉冲信号：其表达式为f(t)=u(t)-u(t-T)，其中u(t)为单位阶跃函数。其信号如下图所示： 7、方波信号：信号为周期为T，前T/2期间信号为正电平信号，后T/2期间信号为负电平信号，其信号如下图所示 U(t)
实验四：离散傅里叶变换 实验原理： DFT的快速算法FFT利用了的三个固有特性：（1）对称性（2）周期性（3）可约性。FFT算法基本上可以分为两大类，即按时间抽选法(DIT,Decimation-In-Time)和按频率抽选法(DIF,Decimation-In-frequency)。 MATLAB中提供了进行快速傅里叶变换的fft函数： X=fft(x)，基2时间抽取FFT算法，x是表示离散信号的向量；X是系数向量； X=fft(x,N),补零或截断的N点DFT，当x得长度小于N时，对补零使其长度为N，当x的长度大于N时，对x截断使其长度为N。 实验内容： =60; n=[0:1:k/2]; xa1=2*sin(10*pi*n/k)+cos(18*pi*n/k); subplot(321) stem(n,xa1) xlabel('N');ylabel('x(n)'); xk1=fft(xa1);xk1=abs(xk1) subplot(322) stem(n,xk1) xlabel('k');ylabel('X(k)'); n=[0:1:k*]; xa1=2*sin(10*pi*n/k)+cos(18*pi*n/k); subplot(323) stem(n,xa1) xlabel('N');ylabel('x(n)'); xk1=fft(xa1);xk1=abs(xk1) subplot(324) stem(n,xk1) xlabel('k');ylabel('X(k)'); n=[0:1:k*2]; xa1=2*sin(10*pi*n/k)+cos(18*pi*n/k); subplot(325) stem(n,xa1) xlabel('N');ylabel('x(n)'); xk1=fft(xa1);xk1=abs(xk1) subplot(326) stem(n,xk1) xlabel('k');ylabel('X(k)');
信号与系统实验总结及心得体会 2011211204 刘梦颉2011210960 信号与系统是电子信息类专业的一门重要的专业核心基础课程，该课程核心的基本概念、基本理论和分析方法都非常重要，而且系统性、理论性很强，是将学生从电路分析领域引入信号处理与传输领域的关键性课程，为此开设必要的实验对我们加强理解深入掌握基本理论和分析方法，以及对抽象的概念具体化有极大的好处，而且为后续专业课程的学习提供了理论和大量实验知识储备，对以后的学术科研和创新工作都是十分重要的。下面我将从实验总结、心得体会、意见与建议等三方面作以总结。 一．实验总结 本学期我们一共做了四次实验，分别为：信号的分类与观察、非正弦周期信号的频谱分析、信号的抽样与恢复（PAM）和模拟滤波器实验。 1.信号的分类与观察 主要目的是：观察常用信号的波形特点以及产生方法，学会用示波器对常用波形参数进行测量。主要内容是：利用实验箱中的S8模块分别产生正弦信号、指数信号和指数衰减正弦信号，并用示波器观察输出信号的波形，测量信号的各项参数，根据测量值计算信号的表达式，并且与理论值进行比较。 2.非正弦信号的频谱分析 主要目的是：掌握频谱仪的基本工作原理和正确使用方法，掌握非正弦周期信好的测试方法，理解非正弦周期信号频谱的离散性、谐波性欲收敛性。主要内
容是：通过频谱仪观察占空比为50%的方波脉冲的频谱，和占空比为20%的矩形波的频谱，并用坐标纸画图。 3.信号的抽样与恢复 主要目的是：验证抽样定理，观察了解PAM信号的形成过程。主要内容是：通过矩形脉冲对正弦信号进行抽样，再把它恢复还原过来，最后用还原后的图形与原图形进行对比，分析实验并总结。 4.模拟滤波器实验 主要目的是：了解RC无源和有源滤波器的种类、基本结构及其特性，比较无源和有源滤波器的滤波特性，比较不同阶数的滤波器的滤波效果。主要内容：利用点频法通过测试无源低通、高通、带通和有源带阻，以及有源带通滤波器的幅频特性，通过描点画图形象地把它们的特点表现出来。 通过对信号与实验课程的学习，我掌握了一些基本仪器的使用方法，DDS 信号源、实验箱、示波器、频谱仪等四种实验仪器。初步了解了对信号的测试与分析方法对以前在书本上看到的常见信号有了更加具体的认识，使得书本上的知识不再那么抽象。 DDS信号源，也就是函数发生器，可以产生固定波形，如正弦波、方波或三角波，频率和幅度可以调节。实验箱是很多个信号实验装置的集合，可谓集多种功能于一身，其中包括函数发生器、模拟滤波器、函数信号的产生与测量、信号的抽样与恢复等模块。示波器能把抽象的电信号转换成具体的图像，便于人们研究各种电现象的变化过程。利用示波器能观察各种不同的信号幅度随时间变化的波形曲线，还可以用它测试各种不同的电量，如电压、电流、频率、相位差、
离散信号的变换（MATLAB 实验） 一、实验目的 掌握用Z 变换判断离散系统的稳定与否的方法，掌握离散傅立叶变换及其基本性质和特点，了解快速傅立叶变换。 二、实验内容 1、已经系统函数为 5147.13418.217.098.2250 5)(2342-++--+=z z z z z z Z H (1) 画出零极点分布图，判断系统是否稳定； (2)检查系统是否稳定； (3) 如果系统稳定，求出系统对于u(n)的稳态输出和稳定时间b=[0,0,1,5,-50];a=[2,-2.98,0.17,2.3418,-1.5147]; subplot(2,1,1);zplane(b,a);title('零极点分布图'); z=roots(a); magz=abs(z) magz = 0.9000 0.9220 0.9220 0.9900 n=[0:1000]; x=stepseq(0,0,1000); s=filter(b,a,x); subplot(2,1,2);stem(n,s);title('稳态输出'); （1）因为极点都在单位园内，所以系统是稳定的。 （2）因为根的幅值（magz ）都小于1，所以这个系统是稳定的。 （3）稳定时间为570。 2、综合运用上述命令，完成下列任务。 (1) 已知)(n x 是一个6点序列： ???≤≤=其它,050,1)(n n x
计算该序列的离散时间傅立叶变换，并绘出它们的幅度和相位。 要求：离散时间傅立叶变换在[-2π,2π]之间的两个周期内取401个等分频率上进行数值求值。 n=0:5;x=ones(1,6); k=-200:200;w=(pi/100)*k; X=x*(exp(-j*pi/100)).^(n'*k); magX=abs(X);angX=angle(X); subplot(2,1,1);plot(w/pi,magX);grid;title('幅度'); subplot(2,1,2);plot(w/pi,angX);grid;title('相位'); (2) 已知下列序列： a. ,1000),52.0cos()48.0cos()(≤≤+=n n n n x ππ； b ．)4sin()(πn n x =是一个N ＝32的有限序列； 试绘制)(n x 及它的离散傅立叶变换 )(k X 的图像。 a ． n=[0:1:100];x=cos(0.48*pi*n)+cos(0.52*pi*n); subplot(2,1,1);plot(n,x);title('x(n)的图像'); X=dft(x,101); magX=abs(X); subplot(2,1,2);plot(n,magX);title('丨X(k)丨的图像');
中南大学 信号与系统试验报告 姓名: 学号: 专业班级:自动化 实验一 基本信号的生成 1．实验目的 ● 学会使用MATLAB 产生各种常见的连续时间信号与离散时间信号； ● 通过MATLAB 中的绘图工具对产生的信号进行观察，加深对常用信号的 理解； ● 熟悉MATLAB 的基本操作，以及一些基本函数的使用，为以后的实验奠 定基础。 2．实验内容 ⑴ 运行以上九个例子程序，掌握一些常用基本信号的特点及其MATLAB 实现方法；改变有关参数，进一步观察信号波形的变化。 ⑵ 在 k [10:10]=- 范围内产生并画出以下信号： a) 1f [k][k]δ=； b) 2f [k][k+2]δ=； c) 3f [k][k-4]δ=； d) 4f [k]2[k+2][k-4]δδ=-。
源程序： k=-10:10; f1k=[zeros(1,10),1,zeros(1,10)]; subplot(2,2,1) stem(k,f1k) title('f1[k]') f2k=[zeros(1,8),1,zeros(1,12)]; subplot(2,2,2) stem(k,f2k) title('f2[k]') f3k=[zeros(1,14),1,zeros(1,6)]; subplot(2,2,3) stem(k,f3k) title('f3[k]') f4k=2*f2k-f3k; subplot(2,2,4) stem(k,f4k) title('f4[k]') ⑶ 在 k [0:31]=范围内产生并画出以下信号： a) ()()k k 144f [k]sin cos π π=； b) ()2k 24f [k]cos π =； c) ()()k k 348f [k]sin cos π π=。 请问这三个信号的基波周期分别是多少？ 源程序： k=0:31; f1k=sin(pi/4*k).*cos(pi/4*k); subplot(3,1,1) stem(k,f1k) title('f1[k]') f2k=(cos(pi/4*k)).^2; subplot(3,1,2) stem(k,f2k) title('f2[k]') f3k=sin(pi/4*k).*cos(pi/8*k); subplot(3,1,3) stem(k,f3k) title('f3[k]') 其中f1[k]的基波周期是4, f2[k]的基波周期是4, f3[k]的基波周期是16。
戶幵,戈丿、弟实验报告 课程名称：彳 _____________ 指导老师 _____________ 成绩： ____________________ 实验名称：离散傅里叶变换和快速傅里叶变换 实验类型: _________________ 同组学生姓名： 一、实验目的和要求(必填) 二、实验内容和原理(必填) 三、主要仪器设备(必填) 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析(必填) 七、讨论、心得 一、实验目的和要求 1. 掌握DFT 的原理和实现 2. 掌握FFT 的原理和实现，掌握用 FFT 对连续信号和离散信号进行谱分析的方法。 二、实验内容和原理 2.1 DTFT 和 DFT N 1 如果x(n)为因果有限长序列,n=0,1,...,N-1，则x(n)的DTFT 表示为：X(e j ) x(n)e n 0 序列的N 点DFT 是DTFT 在［0,2 n 上的N 点等间隔采样，采样间隔为2 d N 。通过DFT , 可以完成由一组有限个信号采样值 x(n)直接计算得到一组有限个频谱采样值 X(k)。X(k)的幅 度谱为X(k) v 'x R (k ) X |2(k ) , X R (k)和X i (k)分别为X(k)的实部和虚部。X(k)的相位谱 为(k) 列吩 序列x(n)的离散事件傅里叶变换(DTFT )表示为: X(e j ) x( n)e x(n)的离散傅里叶变换(DFT )表达式为: X(k) x(n)e n 0 j^nk N (k 0,1,…,N 1)
IDFT )定义为 x(n)丄 N \(k)e j_Nnk (n 0,1，…，N 1) N n 0 2.2 FFT 快速傅里叶变换(FFT )是DFT 的快速算法，它减少了 DFT 的运算量，使数字信号的处理 速度大大提高。 三、主要仪器设备 PC 一台，matlab 软件 四、实验内容 4.1第一题 x(n)的离散时间 傅里叶变换(DTFT ) X(e j Q )并绘图。 0 其2他n 2; (2)已知 x(n) 2n 0 n 10。 0其他 4.1.1理论分析 1) 由DTFT 计算式, X (Q)是实数，可以直接作出它的图像。 离散傅里叶反变换 求有限长离散时间信号 (1)已知 x(n) X( ) x(n)e j n e 2j 1 5j e 1 e j e 2? e 2? 0.5j e 0.5 j e sin(2.5 )
离散时间信号的傅里叶变换和离散傅里叶变换 摘要 本文主要介绍了离散时间信号的离散时间傅里叶变换及离散傅里叶变换，说明其在频域的具体表示和分析，并通过定义的方法和矩阵形式的表示来给出其具体的计算方法。同时还介绍了与离散时间傅里叶变换（DTFT ）和离散傅里叶变换（DFT ）相关的线性卷积与圆周卷积，并讲述它们之间的联系，从而给出了用圆周卷积计算线性卷积的方法，即用离散傅里叶变换实现线性卷积。 1. 离散时间傅里叶变换 1.1离散时间傅里叶变换及其逆变换 离散时间傅里叶变换为离散时间序列x[n]的傅里叶变换，是以复指数序列｛n j e ω-｝的序列来表示的（可对应于三角函数序列），相当于傅里叶级数的展开，为离散时间信号和线性时不变系统提供了一种频域表示，其中ω是实频率变量。时间序列x[n]的离散时间傅里叶变换)(ωj e X 定义如下： ∑∞ -∞ =-= n n j j e n x e X ωω ][)( （1.1） 通常)(ωj e X 是实变量ω的复数函数同时也是周期为π2的周期函数，并且)(ωj e X 的幅度函数和实部是ω的偶函数，而其相位函数和虚部是ω的奇函数。这是由于： ) ()()(tan ) ()()() (sin )()()(cos )()(2 22 ωωωωωωωωωωθωθωθj re j im j im j re j j j im j j re e X e X e X e X e X e X e X e X e X = +=== （1.2） 由于式（1.1）中的傅里叶系数x[n]可以用下面给出的傅里叶积分从)(ωj e X 中算出： ωπ ωπ πω d e e X n x n j j )(21 ][?- = （1.3）
实验四离散傅里叶变换的性质 一、实验目的 1. 熟悉matlab软件中离散傅里叶变换的实现方法及FFT函数的使用方法； 2. 通过软件仿真，加深对离散傅里叶变换性质的理解。 二、实验内容 1. 验证离散傅里叶变换的线性性质； 2. 掌握用matlab实现圆周移位的方法； 3. 验证圆周卷积与线性卷积的关系。 三、实验步骤 1. 验证线性性质 设两个有限长序列分别为xn1=[3,1,-2,2,3,4]，xn2=[1,1,1,1]，做4DFT[xn1]+2DFT[xn2],及DFT[4xn1+2xn2]的运算，比较它们的结果。 代码如下： clear,N=20;n=[0:1:N-1]; xn1=[3,1,-2,2,3,4];n1=0:length(xn1)-1; %定义序列xn1 xn2=[1,1,1,1];n2=0:length(xn2)-1; %定义序列xn2 yn1=4*xn1;yn2=2*xn2;[yn,ny]=seqadd(yn1,n1,yn2,n2); %计算4xn1+2xn2 xk1=fft(xn1,N);xk2=fft(xn2,N); %分别求DFT[xn1] 和DFT[xn2] yk0=4*xk1+2*xk2; %计算4DFT[xn1]+2DFT[xn2] yk=fft(yn,N); %计算DFT[4xn1+2xn2] subplot(2,1,1);stem(n,yk0);title('傅里叶变换之和') %显示4DFT[xn1]+2DFT[xn2] subplot(2,1,2);stem(n,yk);title('序列和之傅里叶变换') %显示DFT[4xn1+2xn2] 运行结果如图1所示，从图中可知，用两种方法计算的DFT完全相等，所以离散傅里叶变换的线性性质得到验证。
实验四 离散信号的频域分析 一、 实验目的 1. 掌握序列的傅里叶变换、离散傅里叶级数、离散傅里叶变换、快速傅里叶变换的Matlab 实现； 2. 学习用FFT 对连续信号和离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT 。 二、 实验内容及步骤 1. 计算序列的DTFT 和DFT ，观察栅栏效应 设)()(4n R n x =，要求用MATLAB 实现： （1）计算)(n x 的傅里叶变换)(ωj e X ，并绘出其幅度谱； （2）分别计算)(n x 的4点DFT 和8点DFT ，绘出其幅度谱。并说明它们和)(ωj e X 的关系。 （提示：DFT 变换可用MA TLAB 提供的函数fft 实现，也可以自己用C 语言或matlab 编写） 2.计算序列的FFT ，观察频谱泄漏 已知周期为16的信号)1612cos()1610cos()(n n n x π π +=。 （1） 截取一个周期长度M=16点，计算其16点FFT ，并绘出其幅度谱； （2） 截取序列长度M=10点，计算其16点FFT ，绘出其幅度谱，并与（1）的结果进行比 较，观察频谱泄漏现象，说明产生频谱泄漏的原因。 三、 实验报告要求 1. 结合实验中所得给定典型序列幅频特性曲线，与理论结果比较，并分析说明误差产生的原因以及用FFT 作谱分析时有关参数的选择方法。 2. 总结实验所得主要结论。 1. 计算序列的DTFT 和DFT ，观察栅栏效应 设)()(4n R n x =，要求用MATLAB 实现： （1）计算)(n x 的傅里叶变换)(ωj e X ，并绘出其幅度谱； （2）分别计算)(n x 的4点DFT 和8点DFT ，绘出其幅度谱。并说明它们和)(ωj e X 的关系。 （1）代码： n=0:3; M=10;
电 子 科 技 大 学 实 验 报 告 学生姓名：项阳 学 号： 2010231060011 指导教师：邓建 一、实验项目名称：离散时间傅里叶变换 二、实验目的： 熟悉序列的傅立叶变换、傅立叶变换的性质、连续信号经理想采样后进行重建，加深对时域采样定理的理解。 三、实验内容： 1. 求下列序列的离散时间傅里叶变换 (a) ()(0.5)()n x n u n = (b) (){1,2,3,4,5}x n = 2. 设/3()(0.9),010,j n x n e n π=≤≤画出()j X e ω并观察其周期性。 3. 设()(0.9),1010,n x n n =--≤≤画出()j X e ω并观察其共轭对称性。 4. 验证离散时间傅里叶变换的线性、时移、频移、反转（翻褶）性质。 5. 已知连续时间信号为t a e t x 1000)(-=，求： (a) )(t x a 的傅里叶变换)(Ωj X a ； (b) 采样频率为5000Hz ，绘出1()j X e ω,用理想内插函数sinc()x 重建)(t x a ，并对结果进行讨论； (c) 采样频率为1000Hz ，绘出2()j X e ω，用理想内插函数sinc()x 重建)(t x a ，并对结果进行讨论。 四、实验原理：
1. 离散时间傅里叶变换(DTFT)的定义： 2．周期性：()j X e ?是周期为2π的函数 (2)()()j j X e X e ??π+= 3．对称性：对于实值序列()x n ，()j X e ?是共轭对称函数。 *()() Re[()]Re[()] Im[()]Im[()]()() ()() j j j j j j j j j j X e X e X e X e X e X e X e X e X e X e ??????????-----===-=∠=-∠ 4．线性：对于任何12,,(),()x n x n αβ，有 1212[()()][()][()]F x n x n F x n F x n αβαβ+=+ 5．时移 [()][()]()j k j j k F x n k F x n e X e e ωωω---== 6．频移 00()[()]()j n j F x n e X e ωωω-= 7．反转（翻褶） [()]()j F x n X e ω--= 五、实验器材（设备、元器件）： PC 机、Windows XP 、MatLab 7.1 六、实验步骤： 本实验要求学生运用MATLAB 编程产生一些基本的离散时间信号，并通过MATLAB 的几种绘图指令画出这些图形，以加深对相关教学内容的理解，同时也通过这些简单的函数练习了MATLAB 的使用。 [()]()()(), ()j j jn z e n n F x n X e X z x n e x n ωωω∞-==-∞∞=-∞===<∞∑∑收敛条件为：
班级: 姓名: 学号: 实验日期: 一、实验名称脉搏、语音及图像信号的傅里叶分析 二、实验目的 1、了解常用周期信号的傅里叶级数表示。 2、了解周期脉搏信号、语音信号及图像信号的傅里叶分析过程 3、理解体会傅里叶分析的理论及现实意义 三、实验仪器 脉搏语音实验仪器，数字信号发生器，示波器 四、实验原理 1、周期信号傅里叶分析的数学基础 任意一个周期为T 的函数f(t)都可以表示为傅里叶级数： 0001 0000000001()(cos sin )21()()1()cos()() 1()sin()()n n n n n f t a a n t b n t a f t d t a f t n t d t b f t n t d t ππ πππ πωωωωπωωωπωωωπ∞=---=++== =∑??? 其中0ω为角频率，称为基频，0a 为常数，n a 和n b 称为第n 次谐波的幅 值。任何周期性非简谐交变信号均可用上述傅里叶级数进行展开，即分解为一系列不同次谐波的叠加。 对于如图1所示的方波，一个周期内的函数表达式为: (0t<)2() (-t 0)2h f t h ππ?≤??=??-≤
其傅里叶级数展开为： 0100041()()sin(21)21411(sin sin 3 sin 5)35n h f t n t n h t t t ωπωωωπ∞==--=+++∑L 同理：对于如图2所示的三角波，函数表达式为： 4t (-t<)44()232(1) (t )44h T T f t t T T h T π?≤??=??-≤
实验<编号> 学号姓名分工 11350023 韦能龙编写代码 11350024 熊栗问题分析1.问题描述 实验二信号的合成与分解
2. 问题分析 此次主要是考察傅里叶的合成与分解，运用分解公式求出系数，运用合成公式合成函数，三角波和矩形波是很典型的连个列子，这个大作业只要分解出系数还有用合成公式，基本上就解决了问题了。 3. 实验代码与实验结果 (1)周期性矩形波的系数表示 ,.....7,5,3,1),2 sin(2==n npi kpi a k 代码： t = -3:0.001:3; M = 1;%M =1,7,29,99 T = 2; W = 2*pi/T; f1 = 0*ones(1,length(t)); for n= -M:2:M a = 2/(n*pi)*sin(n*pi/2); f1 = f1+a*exp(j*n*W*t); end plot(t,f1) xlabel('t') ylabel('f(t)') title('M=1，7,29,99时的方波') ylim([-1.5 1.5]); hold on plot(t , zeros(1,length(t))) hold off 图像： M =1时：
M= 7： M = 29
M = 99 （2）三角波的系数表示：
??--==101)()(1dt e t x dt e t x T a jkwt T jkwt k )2 (sin 42 1 2 2 20npi pi n a a n == 代码： t = -3:0.001:3; M = 1;%M =1,7,29,99 T = 1; W = 2*pi/T; G1= 0*ones(1,length(t)); for n= -M:M if n==0 a =1/2; else a = 4/(n^2*pi^2)*(sin(n*pi/2)^2) ; end G1 = G1+a*exp(j*n*W*t); end G1 = G1-0.5; plot(t,G1) xlabel('t') ylabel('G(t)') title('M=1时的三角波') ylim([-1.5 1.5]); hold on plot(t , zeros(1,length(t))) hold off M=1 时
实验二离散时间傅里叶变换和离散傅里叶变换 一．实验目的 1.深刻理解离散时间信号傅里叶变换的定义，与连续傅里叶变换之间的关系； 2.深刻理解序列频谱的性质（连续的、周期的等）； 3.能用MATLAB编程实现序列的DTFT，并能显示频谱幅频、相频曲线； 4.深刻理解DFT的定义、DFT谱的物理意义、DFT与DTFT之间的关系； 5.能用MATLAB编程实现有限长序列的DFT； 6.熟悉循环卷积的过程，能用MA TLAB编程实现循环卷积运算。 二．实验原理 1.离散时间信号的频谱和图示化 2.离散傅里叶变换的定义和图示化 三．实验结果 w=[0:2:500]*pi*2/500; h=(1+0.9*exp(-j*w))./(1-0.9*exp(-j*w)); magh=abs(h); plot(w/pi,magh);grid;xlabel('f');ylabel('|H(w)|'); n=[0:127]; m=[0:127]; x=exp(j*2*pi/128*m.*n); [xk]=dft(x,128); stem(n,xk);xlabel('n');ylabel('xk');
n=[0:127]; m=[0:127]; x=cos(2*pi/128*m.*n); [xk]=dft(x,128); stem(n,xk);xlabel('n');ylabel('xk'); n=[0:127]; m=[0:127]; [xk]=dft(x,128); stem(n,xk);xlabel('n');ylabel('xk'); n=[0:127];m=[0,127]; x=sin(n); [xk]=dft(x,128); stem(n,xk);xlabel('n');ylabel('xk');