信号与系统实验报告3实验3 傅里叶变换及其性质

合集下载

傅里叶实验报告

傅里叶实验报告

一、实验目的1. 了解傅里叶变换的基本原理和方法。

2. 掌握傅里叶变换在信号处理中的应用。

3. 通过实验验证傅里叶变换在信号处理中的效果。

二、实验原理傅里叶变换是一种将信号从时域转换为频域的方法,它可以将一个复杂的信号分解为一系列不同频率的正弦波和余弦波的叠加。

傅里叶变换的基本原理是:任何周期信号都可以表示为一系列不同频率的正弦波和余弦波的叠加。

三、实验仪器与材料1. 实验箱2. 信号发生器3. 示波器4. 计算机及傅里叶变换软件四、实验步骤1. 设置信号发生器,产生一个正弦信号,频率为f1,幅度为A1。

2. 将信号发生器输出的信号输入到实验箱,通过示波器观察该信号。

3. 利用傅里叶变换软件对观察到的信号进行傅里叶变换,得到频谱图。

4. 改变信号发生器的频率,分别产生频率为f2、f3、f4的正弦信号,重复步骤2-3。

5. 分析不同频率信号的频谱图,观察傅里叶变换在信号处理中的应用。

五、实验数据与结果1. 当信号发生器频率为f1时,示波器显示的信号波形如图1所示。

图1:频率为f1的正弦信号波形2. 对频率为f1的正弦信号进行傅里叶变换,得到的频谱图如图2所示。

图2:频率为f1的正弦信号的频谱图从图2可以看出,频率为f1的正弦信号在频域中只有一个频率成分,即f1。

3. 重复步骤4,分别对频率为f2、f3、f4的正弦信号进行傅里叶变换,得到的频谱图分别如图3、图4、图5所示。

图3:频率为f2的正弦信号的频谱图图4:频率为f3的正弦信号的频谱图图5:频率为f4的正弦信号的频谱图从图3、图4、图5可以看出,不同频率的正弦信号在频域中分别只有一个频率成分,即对应的f2、f3、f4。

六、实验分析与讨论1. 傅里叶变换可以将信号从时域转换为频域,方便我们分析信号的频率成分。

2. 通过傅里叶变换,我们可以得到信号的频谱图,直观地观察信号的频率成分。

3. 实验结果表明,傅里叶变换在信号处理中具有重要作用,可以应用于信号分解、滤波、调制等领域。

《信号与系统》实验三

《信号与系统》实验三
实验记录及个人小结(包括:实验源程序、注释、结果分析与讨论等)
三:
源程序:
(1):τ/T=1/4时的周期矩形脉冲的幅度谱和相位谱:
n=-20:20;
F=zeros(size(n));
forii=-20:20
F(ii+21)= sin(ii*pi/4)/(ii*pi+eps);
end
F(21)=1/4;
实验
内容
1.求图1所示周期信号( , )的傅里叶级数,用Matlab做出其前3、9、21、45项谐波的合成波形与原信号作比较,并做出其单边幅度谱和相位谱。
图1 周期为2的三角脉冲信号
2. 求图2所示的单个三角脉冲( )的傅里叶变换,并做出其幅度谱和相位谱。
图2 单个三角脉冲
3. 求不同占空比下周期矩形脉冲的幅度谱和相位谱,例如 、 。
y=1/4;
forn=1:m
y=y+4/(n*n*pi*pi)*(1-cos(n*pi/2)).*cos(n*pi.*t);
end
源代码:
t=-6:0.01:6;
d=-6:2:6;
fxx=pulstran(t,d,'tripuls');
f1=fourierseries(3,t);
f2=fourierseries(9,t);
n=1:10;
a=zeros(size(n));
fori=1:10
a(i)=angle(4/(i*i*pi*pi)*(1-cos(i*pi/2)))
end
n=0:pi:9*pi
stem(n,a,'fill','linewidth',2);
axis([0,9*pi,-0.2,0.2])

信号与系统课程设计报告傅里叶变换的对称性和时移特性

信号与系统课程设计报告傅里叶变换的对称性和时移特性

信号与系统课程设计报告--傅里叶变换的对称性和时移特性课程设计任务书2沈阳理工大学摘要本文研究的是傅里叶变换的对称性和时移特性,傅里叶变换的性质有:对称性、线性(叠加性)、奇偶虚实性、尺度变换特性、时移特性、频移特性、微分特性、积分特性、卷积特性(时域和频域);从信号与系统的角度出发,给出了激励信号的具体模型;应用Matlab软件进行仿真,将研究的信号转化成具体的函数形式,在Matlab得到最终变换结果。

使用傅里叶变换的方法、卷积的求解方法以及函数的微分等方法研究题目。

关键词: 傅里叶变换;对称性;时移特性;Matlab3沈阳理工大学目录1、Matlab介绍........................... 错误!未定义书签。

2.利用Matlab实现信号的频域分析—傅里叶变换的对称性与时移特性设计 (5)2.1.傅里叶变换的定义及其相关性质 (5)2.2.傅里叶变换的对称性验证编程设计及实现 (7)2.3.傅里叶变换的时移特性验证编程设计及实现 (11)3.总结 (13)4.参考文献 (13)4沈阳理工大学1、Matlab介绍MATLAB作为一种功能强大的工程软件,其重要功能包括数值处理、程序设计、可视化显示、图形用户界面和与外部软件的融合应用等方面。

MATLAB软件由美国Math Works公司于1984年推出,经过不断的发展和完善,如今己成为覆盖多个学科的国际公认的最优秀的数值计算仿真软件。

MATLAB具备强大的数值计算能力,许多复杂的计算问题只需短短几行代码就可在MATLAB中实现。

作为一个跨平台的软件,MATLAB已推出Unix、Windows、Linux和Mac等十多种操作系统下的版本,大大方便了在不同操作系统平台下的研究工作。

MATLAB软件具有很强的开放性和适应性。

在保持内核不变的情况下,MATLAB 可以针对不同的应用学科推出相应的工具箱(toolbox),目前己经推出了图象处理工具箱、信号处理工具箱、小波工具箱、神经网络工具箱以及通信工具箱等多个学科的专用工具箱,极大地方便了不同学科的研究工作。

信号与系统实验报告

信号与系统实验报告

信号与系统实验报告一、实验目的(1) 理解周期信号的傅里叶分解,掌握傅里叶系数的计算方法;(2)深刻理解和掌握非周期信号的傅里叶变换及其计算方法;(3) 熟悉傅里叶变换的性质,并能应用其性质实现信号的幅度调制;(4) 理解连续时间系统的频域分析原理和方法,掌握连续系统的频率响应求解方法,并画出相应的幅频、相频响应曲线。

二、实验原理、原理图及电路图(1) 周期信号的傅里叶分解设有连续时间周期信号()f t ,它的周期为T ,角频率22fT,且满足狄里赫利条件,则该周期信号可以展开成傅里叶级数,即可表示为一系列不同频率的正弦或复指数信号之和。

傅里叶级数有三角形式和指数形式两种。

1)三角形式的傅里叶级数:01212011()cos()cos(2)sin()sin(2)2cos()sin()2n n n n a f t a t a t b t b t a a n t b n t 式中系数n a ,n b 称为傅里叶系数,可由下式求得:222222()cos(),()sin()T T T T nna f t n t dtb f t n t dtTT2)指数形式的傅里叶级数:()jn tn nf t F e式中系数n F 称为傅里叶复系数,可由下式求得:221()T jn tT nF f t edtT周期信号的傅里叶分解用Matlab进行计算时,本质上是对信号进行数值积分运算。

Matlab中进行数值积分运算的函数有quad函数和int函数。

其中int函数主要用于符号运算,而quad函数(包括quad8,quadl)可以直接对信号进行积分运算。

因此利用Matlab进行周期信号的傅里叶分解可以直接对信号进行运算,也可以采用符号运算方法。

quadl函数(quad系)的调用形式为:y=quadl(‘func’,a,b)或y=quadl(@myfun,a,b)。

其中func是一个字符串,表示被积函数的.m文件名(函数名);a、b分别表示定积分的下限和上限。

傅里叶变换 实验报告

傅里叶变换 实验报告

傅里叶变换实验报告傅里叶变换实验报告引言:傅里叶变换是一种重要的数学工具,广泛应用于信号处理、图像处理、物理学、工程学等领域。

本次实验旨在通过实际操作和数据分析,深入了解傅里叶变换的原理、特性以及应用。

一、实验目的本实验的目的是通过实际操作,掌握傅里叶变换的基本原理,了解其在信号处理中的应用,并能够正确进行频域分析。

二、实验仪器和材料1. 信号发生器2. 示波器3. 计算机4. 傅里叶变换软件三、实验步骤1. 将信号发生器与示波器连接,并设置合适的频率和幅度,产生一个正弦信号。

2. 通过示波器观察并记录原始信号的时域波形。

3. 将示波器输出的信号通过音频线连接到计算机的输入端口。

4. 打开傅里叶变换软件,选择输入信号源为计算机输入端口,并进行采样。

5. 在傅里叶变换软件中,通过选择合适的窗函数、采样频率和采样点数,进行傅里叶变换。

6. 观察并记录变换后的频域波形,并进行分析。

四、实验结果与分析通过实验操作和数据分析,我们得到了信号的时域波形和频域波形。

在时域波形中,我们可以清晰地看到正弦信号的周期性特征,而在频域波形中,我们可以看到信号的频率成分。

傅里叶变换将信号从时域转换到频域,通过分析频域波形,我们可以得到信号的频率成分。

在实验中,我们可以通过改变信号发生器的频率和幅度,观察频域波形的变化,进一步理解傅里叶变换的原理和特性。

此外,傅里叶变换还可以用于信号滤波。

通过观察频域波形,我们可以选择性地去除某些频率成分,从而实现信号的滤波处理。

这在音频处理、图像处理等领域中具有广泛的应用。

五、实验总结本次实验通过实际操作和数据分析,深入了解了傅里叶变换的原理、特性以及应用。

傅里叶变换作为一种重要的数学工具,在信号处理、图像处理等领域中具有广泛的应用前景。

通过本次实验,我们不仅掌握了傅里叶变换的基本原理和操作方法,还深入了解了信号的时域和频域特性。

这对于我们进一步研究和应用傅里叶变换具有重要的意义。

总之,傅里叶变换是一项重要的数学工具,通过实际操作和数据分析,我们可以更好地理解和应用傅里叶变换,为信号处理和图像处理等领域的研究和应用提供有力支持。

实验三傅里叶变换及其性质

实验三傅里叶变换及其性质

信息工程学院实验报告课程名称:信号与系统实验项目名称:实验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 ωω∞--∞==⎰,傅里叶反变换定义为:11()[()]()2j 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 的函数,而不是默认的ω,即()()jvt F v f t e dt ∞--∞=⎰。

(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章傅里叶变换

信号与系统第3章傅里叶变换

*本章要点
1.利用傅立叶级数的定义式分析周期信号的离散谱。 2.利用傅立叶积分分析非周期信号的连续谱。 3.理解信号的时域与频域间的关系。 4.用傅立叶变换的性质进行正逆变换。 5.掌握抽样信号频谱的计算及抽样定理
将信号表示为不同频率正弦分量的线性组合意义
1.从信号分析的角度 将信号表示为不同频率正弦分量的线性组合,为不同信号之 间进行比较提供了途径。
发展历史
•1822年,法国数学家傅里叶(J.Fourier,1768-1830)在研究热传导 理论时发表了“热的分析理论”,提出并证明了将周期函数展 开为正弦级数的原理,奠定了傅里叶级数的理论基础。 •泊松(Poisson)、高斯(Guass)等人把这一成果应用到电学中去, 得到广泛应用。 •19世纪末,人们制造出用于工程实际的电容器。 •进入20世纪以后,谐振电路、滤波器、正弦振荡器等一系列具 体问题的解决为正弦函数与傅里叶分析的进一步应用开辟了广 阔的前景。 •在通信与控制系统的理论研究和工程实际应用中,傅里叶变换 法具有很多的优点。 •“FFT”快速傅里叶变换为傅里叶分析法赋予了新的生命力。
一.三角函数形式的傅里叶级数
1.正交三角函数集
三角函数系1, cos x,sin x, cos 2x,sin 2x,..., cos nx,sin nx,...
在区间[-π,π]上正交,是指在三角函数系中任何不同的两个函 数的乘积在区间的积分等于零,即
cosnxdx 0(n 1,2,3,...)
傅里叶生平
1768年生于法国 1807年提出“任何周期信号
都可用正弦函数级数表示” 1829年狄里赫利第一个给出
收敛条件 拉格朗日反对发表 1822年首次发表“热的分析
理论”中

信号与系统实验报告实验三 连续时间LTI系统的频域分析报告

信号与系统实验报告实验三   连续时间LTI系统的频域分析报告

实验三 连续时间LTI 系统的频域分析一、实验目的1、掌握系统频率响应特性的概念及其物理意义;2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法,理解具有不同频率响应特性的滤波器对信号的滤波作用;3、学习和掌握幅度特性、相位特性以及群延时的物理意义;4、掌握用MA TLAB 语言进行系统频响特性分析的方法。

基本要求:掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法,深刻理解LTI 系统的频率响应特性的物理意义,理解滤波和滤波器的概念,掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。

二、实验原理及方法1 连续时间LTI 系统的频率响应所谓频率特性,也称为频率响应特性,简称频率响应(Frequency response ),是指系统在正弦信号激励下的稳态响应随频率变化的情况,包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。

上图中x(t)、y(t)分别为系统的时域激励信号和响应信号,h(t)是系统的单位冲激响应,它们三者之间的关系为:)(*)()(t h t x t y =,由傅里叶变换的时域卷积定理可得到:)()()(ωωωj H j X j Y =3.1或者: )()()(ωωωj X j Y j H =3.2)(ωj H 为系统的频域数学模型,它实际上就是系统的单位冲激响应h(t)的傅里叶变换。

即⎰∞∞--=dt e t h j H tj ωω)()(3.3由于H(j ω)实际上是系统单位冲激响应h(t)的傅里叶变换,如果h(t)是收敛的,或者说是绝对可积(Absolutly integrabel )的话,那么H(j ω)一定存在,而且H(j ω)通常是复数,因此,也可以表示成复数的不同表达形式。

在研究系统的频率响应时,更多的是把它表示成极坐标形式:)()()(ωϕωωj ej H j H = 3.4上式中,)j (ωH 称为幅度频率相应(Magnitude response ),反映信号经过系统之后,信号各频率分量的幅度发生变化的情况,)(ωϕ称为相位特性(Phase response ),反映信号经过系统后,信号各频率分量在相位上发生变换的情况。

实验三用FFT对信号作频谱分析_实验报告

实验三用FFT对信号作频谱分析_实验报告

实验三用FFT对信号作频谱分析_实验报告一、实验目的1.学习使用FFT(快速傅里叶变换)对信号进行频谱分析;2.掌握频谱分析的基本原理和方法;3.熟悉使用MATLAB进行频谱分析的操作。

二、实验原理FFT是一种基于傅里叶变换的算法,可以将时域信号转换为频域信号,并将信号的频谱特征展示出来。

在频谱分析中,我们通过分析信号的频谱可以获得信号的频率、幅值等信息,从而对信号的性质和特征进行研究。

对于一个连续信号,我们可以通过采样的方式将其转换为离散信号,再利用FFT算法对离散信号进行频谱分析。

FFT算法可以将信号从时域转换到频域,得到离散的频谱,其中包含了信号的频率分量以及对应的幅值。

MATLAB中提供了fft函数,可以方便地对信号进行FFT分析。

通过对信号进行FFT操作,可以得到信号的频谱图,并从中提取出感兴趣的频率信息。

三、实验步骤1.准备工作:(2)建立新的MATLAB脚本文件。

2.生成信号:在脚本中,我们可以通过定义一个信号的频率、幅值和时间长度来生成一个信号的波形。

例如,我们可以生成一个频率为1000Hz,幅值为1的正弦波信号,并设置信号的时间长度为1秒。

3.对信号进行FFT分析:调用MATLAB中的fft函数,对信号进行FFT分析。

通过设置采样频率和FFT长度,可以得到信号的频谱。

其中,采样频率是指在单位时间内连续采样的次数,FFT长度是指离散信号的样本点数。

4.绘制频谱图:调用MATLAB中的plot函数,并设置x轴为频率,y轴为幅值,可以绘制出信号的频谱图。

频谱图上横坐标表示信号的频率,纵坐标表示信号的幅值,通过观察可以得到信号的频率分布情况。

四、实验结果在实验过程中,我们生成了一个频率为1000Hz,幅值为1的正弦波信号,并对其进行FFT分析。

通过绘制频谱图,我们发现信号在1000Hz处有最大幅值,说明信号主要由这一频率成分组成。

五、实验总结本实验通过使用FFT对信号进行频谱分析,我们可以方便地从信号的波形中提取出频率分量的信息,并绘制出频谱图进行观察。

信号与系统第3章 傅里叶变换

信号与系统第3章  傅里叶变换

P
f
2 (t) 1 T1
t0 T1 t0
f
2 (t)d t
a0 2
1 2
n1
(an
2
bn 2 )
2
Fn _____ 帕塞瓦尔定理
n
结论:周期信号的平均功率等于傅里叶级数展开 式中基波分量及各谐波分量有效值的平方 和,即时域和频域的能量守恒。
五. 周期信f号(t)的频c0 谱 (c三n c角os函(n数1t形 式n )) n1
(1) 偶函数 f (t) f (t)
4
an T1
T1
2 0
f (t) cos(n1t)dt
Fn
Fn
an 2
bn 0
傅里叶级数中不会含有正弦项, 只可能含有直流项和余弦项。
(2) 奇函数 f (t) f (t)
a0 0 , an 0
bn
4 T1
T1
2 0
f (t) sin(n1t)d t
e j n1t
T1 n 2
画频谱图:
c0
a0
E
T1
an
2E
T1
Sa
n1
2
, n
1,2,
cn an
1)令 m
2

2
m
即在
2
m,m为整数处有零点。
2)
2
2
T1
T1
零点间谱线个数
3) c n值为正,相位为0,值为负,相位为π
4)谱线间隔为 1 带宽
2
T1
,第一个过零点带宽定义为信号的
1 3
s in31t
1 4
sin41t
E
1 n1

频域特性分析实验报告

频域特性分析实验报告

一、实验目的1. 理解频域分析在信号与系统分析中的重要性。

2. 掌握使用MATLAB进行频域分析的基本方法。

3. 通过实验,分析典型信号和系统的频域特性。

4. 熟悉并运用傅里叶变换、拉普拉斯变换等频域分析方法。

二、实验原理频域分析是信号与系统分析的重要方法之一,它将时域信号转换到频域进行分析,从而揭示信号的频率组成和系统对信号的频率响应特性。

主要分析方法包括傅里叶变换、拉普拉斯变换等。

三、实验步骤1. 实验一:傅里叶变换(1)选择一个典型信号,如正弦波、方波等。

(2)使用MATLAB的傅里叶变换函数进行变换。

(3)观察并分析信号的频谱图,包括频率、幅度等特性。

2. 实验二:拉普拉斯变换(1)选择一个典型信号,如指数函数、指数衰减函数等。

(2)使用MATLAB的拉普拉斯变换函数进行变换。

(3)观察并分析信号的复频域特性,包括极点、零点等。

3. 实验三:系统频率响应分析(1)设计一个典型系统,如滤波器、控制器等。

(2)使用MATLAB的系统函数和频率响应函数进行频率响应分析。

(3)观察并分析系统的幅频响应、相频响应等特性。

四、实验结果与分析1. 实验一:傅里叶变换以正弦波为例,进行傅里叶变换实验。

- 正弦波时域波形如图1所示。

- 正弦波的频谱图如图2所示。

图1:正弦波时域波形图2:正弦波频谱图从图2可以看出,正弦波的频谱只有一个频率成分,即正弦波本身的频率。

2. 实验二:拉普拉斯变换以指数函数为例,进行拉普拉斯变换实验。

- 指数函数时域波形如图3所示。

- 指数函数的复频域特性如图4所示。

图3:指数函数时域波形图4:指数函数复频域特性从图4可以看出,指数函数的拉普拉斯变换具有一个极点,表示信号在复频域中的位置。

3. 实验三:系统频率响应分析以一阶低通滤波器为例,进行频率响应分析实验。

- 滤波器的传递函数为:H(s) = 1 / (1 + s)- 使用MATLAB的系统函数和频率响应函数进行频率响应分析。

信号与系统实验报告

信号与系统实验报告

实验三常见信号的MATLAB 表示及运算一、实验目的1.熟悉常见信号的意义、特性及波形2.学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法二、实验原理根据MATLAB 的数值计算功能和符号运算功能,在MATLAB 中,信号有两种表示方法,一种是用向量来表示,另一种则是用符号运算的方法;在采用适当的MATLAB 语句表示出信号后,就可以利用MATLAB 中的绘图命令绘制出直观的信号波形了;1.连续时间信号从严格意义上讲,MATLAB 并不能处理连续信号;在MATLAB 中,是用连续信号在等时间间隔点上的样值来近似表示的,当取样时间间隔足够小时,这些离散的样值就能较好地近似出连续信号;在MATLAB 中连续信号可用向量或符号运算功能来表示; ⑴ 向量表示法对于连续时间信号()f t ,可以用两个行向量f 和t 来表示,其中向量t 是用形如12::t t p t =的命令定义的时间范围向量,其中,1t 为信号起始时间,2t 为终止时间,p 为时间间隔;向量f 为连续信号()f t 在向量t 所定义的时间点上的样值; ⑵ 符号运算表示法如果一个信号或函数可以用符号表达式来表示,那么我们就可以用前面介绍的符号函数专用绘图命令ezplot 等函数来绘出信号的波形; ⑶ 常见信号的MATLAB 表示 单位阶跃信号单位阶跃信号的定义为:10()0t u t t >⎧=⎨<⎩方法一: 调用Heavisidet 函数首先定义函数Heavisidet 的m 函数文件,该文件名应与函数名同名即;%定义函数文件,函数名为Heaviside,输入变量为x,输出变量为y function y= Heavisidety=t>0; %定义函数体,即函数所执行指令%此处定义t>0时y=1,t<=0时y=0,注意与实际的阶跃信号定义的区别;方法二:数值计算法在MATLAB 中,有一个专门用于表示单位阶跃信号的函数,即stepfun 函数,它是用数值计算法表示的单位阶跃函数()u t ;其调用格式为:stepfunt,t0其中,t 是以向量形式表示的变量,t0表示信号发生突变的时刻,在t0以前,函数值小于零,t0以后函数值大于零;有趣的是它同时还可以表示单位阶跃序列()u k ,这只要将自变量以及取样间隔设定为整数即可; 符号函数符号函数的定义为:10sgn()1t t t >⎧=⎨-<⎩在MATLAB 中有专门用于表示符号函数的函数sign ,由于单位阶跃信号 t 和符号函数两者之间存在以下关系:1122()sgn()t t ε=+,因此,利用这个函数就可以很容易地生成单位阶跃信号;2.离散时间信号离散时间信号又叫离散时间序列,一般用()f k 表示,其中变量k 为整数,代表离散的采样时间点采样次数;在MATLAB 中,离散信号的表示方法与连续信号不同,它无法用符号运算法来表示,而只能采用数值计算法表示,由于MATLAB 中元素的个数是有限的,因此,MATLAB 无法表示无限序列;另外,在绘制离散信号时必须使用专门绘制离散数据的命令,即stem 函数,而不能用plot 函数; 单位序列()k δ单位序列()k δ的定义为10()0k k k δ=⎧=⎨≠⎩单位阶跃序列()u k单位阶跃序列()u k 的定义为10()0k u k k ≥⎧=⎨<⎩3.卷积积分两个信号的卷积定义为:MATLAB 中是利用conv 函数来实现卷积的;功能:实现两个函数1()f t 和2()f t 的卷积;格式:g=convf1,f2说明:f1=f 1t,f2=f 2t 表示两个函数,g=gt 表示两个函数的卷积结果;三、实验内容1.分别用MATLAB 的向量表示法和符号运算功能,表示并绘出下列连续时间信号的波形: ⑴ 2()(2)()tf t e u t -=- ⑵[]()cos()()(4)2tf t u t u t π=--1 t=-1::10;t1=-1::; t2=0::10;f1=zeros1,lengtht1,ones1,lengtht2;f=2-exp-2t.f1; plott,faxis-1,10,0, syms t;f=sym'2-exp-2theavisidet'; ezplotf,-1,10;2t=-2::8;f=0.t<0+cospit/2.t>0&t<4+0.t>4; plott,f syms t;f=sym'cospit/2heavisidet-heavisidet-4 '; ezplotf,-2,8;2.分别用MATLAB 表示并绘出下列离散时间信号的波形: ⑵ []()()(8)f t k u k u k =-- ⑶()sin()()4k f k u k π= 2 t=0:8; t1=-10:15;f=zeros1,10,t,zeros1,7; stemt1,faxis-10,15,0,10; 3 t=0:50; t1=-10:50;f=zeros1,10,sintpi/4; stemt1,faxis-10,50,-2,23.已知两信号1()(1)()f t u t u t =+-,2()()(1)f t u t u t =--,求卷积积分12()()()g t f t f t =*,并与例题比较;t1=-1::0; t2=0::1; t3=-1::1;f1=onessizet1; f2=onessizet2; g=convf1,f2;subplot3,1,1,plott1,f1; subplot3,1,2,plott2,f2; subplot3,1,3,plott3,g;与例题相比较,gt 的定义域不同,最大值对应的横坐标也不同;4.已知{}{}12()1,1,1,2,()1,2,3,4,5f k f k ==,求两序列的卷积和 ;N=4; M=5; L=N+M-1; f1=1,1,1,2;f2=1,2,3,4,5; g=convf1,f2; kf1=0:N-1; kf2=0:M-1; kg=0:L-1;subplot1,3,1,stemkf1,f1,'k';xlabel'k'; ylabel'f1k';grid onsubplot1,3,2,stemkf2,f2,'k';xlabel'k'; ylabel'f2k';grid onsubplot1,3,3;stemkg,g,'k';xlabel'k'; ylabel'gk';grid on 实验心得:第一次接触Mutlab 这个绘图软件,觉得挺新奇的,同时 ,由于之前不太学信号与系统遇到一些不懂的问题,结合这些图对信号与系统有更好的了解;实验四 连续时间信号的频域分析一、实验目的1.熟悉傅里叶变换的性质 2.熟悉常见信号的傅里叶变换3.了解傅里叶变换的MATLAB 实现方法二、实验原理从已知信号()f t 求出相应的频谱函数()F j ω的数学表示为:()F j ω()j t f t e dt ω∞--∞=⎰傅里叶反变换的定义为:1()()2j t f t F j e d ωωωπ∞-∞=⎰在MATLAB 中实现傅里叶变换的方法有两种,一种是利用MATLAB 中的Symbolic Math Toolbox 提供的专用函数直接求解函数的傅里叶变换和傅里叶反变换,另一种是傅里叶变换的数值计算实现法;1.直接调用专用函数法①在MATLAB 中实现傅里叶变换的函数为:F=fourier f 对ft 进行傅里叶变换,其结果为Fw F =fourierf,v 对ft 进行傅里叶变换,其结果为Fv F=fourier f,u,v 对fu 进行傅里叶变换,其结果为Fv ②傅里叶反变换f=ifourier F 对Fw 进行傅里叶反变换,其结果为fx f=ifourierF,U 对Fw 进行傅里叶反变换,其结果为fu f=ifourier F,v,u 对Fv 进行傅里叶反变换,其结果为fu 注意:1在调用函数fourier 及ifourier 之前,要用syms 命令对所有需要用到的变量如t,u,v,w 等进行说明,即要将这些变量说明成符号变量;对fourier 中的f 及ifourier 中的F 也要用符号定义符sym 将其说明为符号表达式;2采用fourier 及fourier 得到的返回函数,仍然为符号表达式;在对其作图时要用ezplot 函数,而不能用plot 函数;3fourier 及fourier 函数的应用有很多局限性,如果在返回函数中含有δω等函数,则ezplot 函数也无法作出图来;另外,在用fourier 函数对某些信号进行变换时,其返回函数如果包含一些不能直接表达的式子,则此时当然也就无法作图了;这是fourier 函数的一个局限;另一个局限是在很多场合,尽管原时间信号ft 是连续的,但却不能表示成符号表达式,此时只能应用下面介绍的数值计算法来进行傅氏变换了,当然,大多数情况下,用数值计算法所求的频谱函数只是一种近似值;2、傅里叶变换的数值计算实现法严格说来,如果不使用symbolic 工具箱,是不能分析连续时间信号的;采用数值计算方法实现连续时间信号的傅里叶变换,实质上只是借助于MATLAB 的强大数值计算功能,特别是其强大的矩阵运算能力而进行的一种近似计算;傅里叶变换的数值计算实现法的原理如下: 对于连续时间信号ft,其傅里叶变换为:其中τ为取样间隔,如果ft 是时限信号,或者当|t|大于某个给定值时,ft 的值已经衰减得很厉害,可以近似地看成是时限信号,则上式中的n 取值就是有限的,假定为N,有: 若对频率变量ω进行取样,得: 通常取:02k k k MM ωπωτ==,其中0ω是要取的频率范围,或信号的频带宽度;采用MATLAB 实现上式时,其要点是要生成ft 的N 个样本值()f n τ的向量,以及向量k j n eωτ-,两向量的内积即两矩阵的乘积,结果即完成上式的傅里叶变换的数值计算;注意:时间取样间隔τ的确定,其依据是τ必须小于奈奎斯特Nyquist 取样间隔;如果ft 不是严格的带限信号,则可以根据实际计算的精度要求来确定一个适当的频率0ω为信号的带宽;三、 实验内容1.编程实现求下列信号的幅度频谱1 求出1()(21)(21)f t u t u t =+--的频谱函数F 1jω,请将它与上面门宽为2的门函数()(1)(1)f t u t u t =+--的频谱进行比较,观察两者的特点,说明两者的关系;2 三角脉冲21||||1()0||1t t f t t -≤⎧=⎨>⎩3 单边指数信号3()()tf t e t ε-=4 高斯信号23()t f t e -=1 syms t w Gt=sym'Heaviside2t+1-Heaviside2t-1'; Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; FFP=absFFw; ezplotFFP,-10pi 10pi;grid; axis-10pi 10pi 0与()(1)(1)f t u t u t =+--的频谱比较,1()(21)(21)f t u t u t =+--的频谱函数F 1jω最大值是其的1/2; 2syms t w;Gt=sym'1+tHeavisidet+1-Heavisidet+1-tHeavisidet-Heavisidet-1'; Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; FFP=absFFw; ezplotFFP,-10pi 10pi;grid; axis-10pi 10pi 0 3syms t w Gt=sym'exp-tHeavisidet';Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; FFP=absFFw; ezplotFFP,-10pi 10pi;grid; axis-10pi 10pi -1 2 4syms t w Gt=sym'exp-t^2';Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; ezplotFFw,-30 30;grid; axis-30 30 -1 22.利用ifourier 函数求下列频谱函数的傅氏反变换122()16F j j ωωω=-+ 222()58()()65j j F j j j ωωωωω+-=++1syms t w Fw=sym'-i2w/16+w^2'; ft=ifourierFw,w,t; ft运行结果: ft =-exp4theaviside-t+exp-4theavisidet 2syms t wFw=sym'iw^2+5iw-8/iw^2+6iw+5'; ft=ifourierFw,w,t; ft运行结果: ft =diract+-3exp-t+2exp-5theavisidet实验心得matlab 不但具有数值计算能力,还能建模仿真,能帮助我们理解不同时间信号的频域分析;实验五 连续时间系统的频域分析一、实验目的1. 学习由系统函数确定系统频率特性的方法;2. 学习和掌握连续时间系统的频率特性及其幅度特性、相位特性的物理意义;3.通过本实验了解低通、高通、带通、全通滤波器的性能及特点;二、实验原理及方法频域分析法与时域分析法的不同之处主要在于信号分解的单元函数不同;在频域分析法中,信号分解成一系列不同幅度、不同频率的等幅正弦函数,通过求取对每一单元激励产生的响应,并将响应叠加,再转换到时域以得到系统的总响应;所以说,频域分析法是一种变域分析法;它把时域中求解响应的问题通过 Fourier 级数或 Fourier 变换转换成频域中的问题;在频域中求解后再转换回时域从而得到最终结果;在实际应用中,多使用另一种变域分析法:复频域分析法,即 Laplace 变换分析法;所谓频率特性,也称频率响应特性,是指系统在正弦信号激励下稳态响应随频率变化的情况,包括幅度随频率的响应和相位随频率的响应两个方面;利用系统函数也可以确定系统频率特性,公式如下:幅度响应用()ωj H 表示,相位响应用)(ωϕH 表示;本实验所研究的系统函数Hs 是有理函数形式,也就是说,分子、分母分别是m 、n 阶多项式; 要计算频率特性,可以写出为了计算出()ωj H 、)(ωϕH 的值,可以利用复数三角形式的一个重要特性: 而⎥⎦⎤⎢⎣⎡+=2sin 2cosππωωj j ,则()⎥⎦⎤⎢⎣⎡+=2sin 2cos ππωωn j n j n n利用这些公式可以化简高次幂,因此分子和分母的复数多项式就可以转化为分别对实部与虚部的实数运算,算出分子、分母的实部、虚部值后,最后就可以计算出幅度()ωj H 、相位)(ωϕH 的值了;三、实验内容a)sm m ms H )(1)(2-+=,m 取值区间 0,1,绘制一组曲线 m=,,,,; b) 绘制下列系统的幅频响应对数曲线和相频响应曲线,分析其频率特性; a %figurealpha=,,,,;colorn='r' 'g' 'b' 'y' 'k'; % r g b y m c k 红,绿,蓝,黄,品红,青,黑 for n=1:5b=0 alphan; % 分子系数向量a=alphan-alphan^2 1; % 分母系数向量 printsysb,a,'s' Hz,w=freqsb,a; w=w./pi; magh=absHz;zerosIndx=findmagh==0; maghzerosIndx=1; magh=20log10magh; maghzerosIndx=-inf; angh=angleHz;angh=unwrapangh180/pi; subplot1,2,1plotw,magh,colornn;hold onsubplot1,2,2plotw,angh,colornn;hold onendsubplot1,2,1hold offxlabel'特征角频率\times\pi rad/sample' title'幅频特性曲线 |Hw| dB';subplot1,2,2hold offxlabel'特征角频率 \times\pi rad/sample' title'相频特性曲线 \thetaw degrees';b1 %b=1,0; % 分子系数向量a=1,1; % 分母系数向量printsysb,a,'s'Hz,w=freqsb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';2 %b=0,1,0; % 分子系数向量a=1,3,2; % 分母系数向量printsysb,a,'s'Hz,w=freqsb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';3 %b=1,-1; % 分子系数向量a=1,1; % 分母系数向量printsysb,a,'s'Hz,w=freqsb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';实验心得:虽然之前用公式转换到频域上分析,但是有时会觉得挺抽象的,不太好理解;根据这些图像结合起来更进一步对信号的了解;同时,这个在编程序时,虽然遇到一些问题,但是总算解决了;实验六离散时间系统的Z域分析一、 实验目的1. 学习和掌握离散系统的频率特性及其幅度特性、相位特性的物理意义;2. 深入理解离散系统频率特性和对称性和周期性;3. 认识离散系统频率特性与系统参数之间的系统4.通过阅读、修改并调试本实验所给源程序,加强计算机编程能力; 二、 实验原理及方法对于离散时间系统,系统单位冲激响应序列)(n h 的 Fourier 变换)(ωj e H 完全反映了系统自身的频率特性,称)(ωj eH 为离散系统的频率特性,可由系统函数)(z H 求出,关系式如下:ωωj j e z z H e H ==)()( 6 – 1由于ωj e是频率的周期函数,所以系统的频率特性也是频率的周期函数,且周期为π2,因此研究系统频率特性只要在πωπ≤≤-范围内就可以了;∑∑∑∞-∞=∞-∞=∞-∞=--==n n n j j n n h j n n h en h e H )sin()()cos()()()(ωωωω6 – 2容易证明,其实部是ω的偶函数,虚部是ω的奇函数,其模ωj e H (的ω的偶函数,相位[])(arg ωj e H 是ω的奇函数;因此研究系统幅度特性)(ωj e H 、相位特性[])(arg ωj e H ,只要在πω≤≤0范围内讨论即可;综上所述,系统频率特性)(ωj eH 具有周期性和对称性,深入理解这一点是十分重要的;当离散系统的系统结构一定,它的频率特性)(ωj e H 将随参数选择的不同而不同,这表明了系统结构、参数、特性三者之间的关系,即同一结构,参数不同其特性也不同; 例如,下图所示离散系统,其数学模型由线性常系数差分方程描述:)()1()(n x n ay n y +-=系统函数:a z az z z H >-=,)(系统函数频率特性:ωωωωωsin )cos 1(1)(ja a a e e e H j j j +-=-=幅频特性:ωωcos 211)(2a a eH j -+=相频特性:[]ωωωcos 1sin arctan)(arg a a eH j --= 容易分析出,当10<<a 时系统呈低通特性,当01<<-a 时系统呈高通特性;当0=a 时系统呈全通特性;同时说明,在系统结构如图所示一定时,其频率特性随参数a 的变化而变化;三、 实验内容a 2281.011)(----=z z z H ;b 1.04.06.01.03.03.01.0)(2323+++-+-=z z z z z z z Hc 2181.011)(--+-=zz z H a %b=1,0,-1; % 分子系数向量a=1,0,; % 分母系数向量printsysb,a,'z'Hz,w=freqzb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';带通b %b=,,,; % 分子系数向量a=1,,,; % 分母系数向量printsysb,a,'z'Hz,w=freqzb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';高通c %b=1,-1,0; % 分子系数向量a=1,0,; % 分母系数向量printsysb,a,'z'Hz,w=freqzb,a;w=w./pi;magh=absHz;zerosIndx=findmagh==0;maghzerosIndx=1;magh=20log10magh; % 以分贝maghzerosIndx=-inf;angh=angleHz;angh=unwrapangh180/pi; % 角度换算figuresubplot1,2,1plotw,magh;grid onxlabel'特征角频率\times\pi rad/sample'title'幅频特性曲线 |Hw| dB';subplot1,2,2plotw,angh;grid onxlabel'特征角频率 \times\pi rad/sample'title'相频特性曲线 \thetaw degrees';带通实验心得:本来理论知识不是很强的,虽然已经编出程序得到相关图形,但是不会辨别相关通带,这让我深刻地反省;。

信号_频域分析实验报告(3篇)

信号_频域分析实验报告(3篇)

第1篇一、实验目的1. 理解信号的频域分析方法及其在信号处理中的应用。

2. 掌握傅里叶变换的基本原理和计算方法。

3. 学习使用MATLAB进行信号的频域分析。

4. 分析不同信号在频域中的特性,理解频域分析在实际问题中的应用。

二、实验原理频域分析是信号处理中一种重要的分析方法,它将信号从时域转换到频域,从而揭示信号的频率结构。

傅里叶变换是频域分析的核心工具,它可以将任何信号分解为不同频率的正弦波和余弦波的线性组合。

三、实验内容及步骤1. 信号生成与傅里叶变换- 使用MATLAB生成一个简单的正弦波信号,频率为50Hz,采样频率为1000Hz。

- 对生成的正弦波信号进行傅里叶变换,得到其频谱图。

2. 频谱分析- 分析正弦波信号的频谱图,观察其频率成分和幅度分布。

- 改变正弦波信号的频率和幅度,观察频谱图的变化,验证傅里叶变换的性质。

3. 信号叠加- 将两个不同频率的正弦波信号叠加,生成一个复合信号。

- 对复合信号进行傅里叶变换,分析其频谱图,验证频谱叠加原理。

4. 窗函数- 使用不同类型的窗函数(如矩形窗、汉宁窗、汉明窗等)对信号进行截取,观察窗函数对频谱的影响。

- 分析不同窗函数的频率分辨率和旁瓣抑制能力。

5. 信号滤波- 设计一个低通滤波器,对信号进行滤波处理,观察滤波器对信号频谱的影响。

- 分析滤波器对信号时域和频域特性的影响。

6. MATLAB工具箱- 使用MATLAB信号处理工具箱中的函数,如`fft`、`ifft`、`filter`等,进行信号的频域分析。

- 学习MATLAB工具箱中的函数调用方法和参数设置。

四、实验结果与分析1. 正弦波信号的频谱分析实验结果显示,正弦波信号的频谱图只有一个峰值,位于50Hz处,说明信号只包含一个频率成分。

2. 信号叠加的频谱分析实验结果显示,复合信号的频谱图包含两个峰值,分别对应两个正弦波信号的频率。

验证了频谱叠加原理。

3. 窗函数对频谱的影响实验结果显示,不同类型的窗函数对频谱的影响不同。

信号与系统Matlab实验作业

信号与系统Matlab实验作业

实验一典型连续时间信号和离散时间信号一、实验目的掌握利用Matlab画图函数和符号函数显示典型连续时间信号波形、典型时间离散信号、连续时间信号在时域中的自变量变换。

二、实验内容1、典型连续信号的波形表示(单边指数信号、复指数信号、抽样信号、单位阶跃信号、单位冲击信号)1)画出教材P28习题1-1(3) ()[(63)(63)]t=----的波形图。

f t e u t u t2)画出复指数信号()()j t f t e σω+=当0.4, 8σω==(0<t<10)时的实部和虚部的波形图。

t=0:0.01:10;f1='exp(0.4*t)*cos(8*t)';f2='exp(0.4*t)*sin(8*t)';figure(1)ezplot(f1,t);grid on;figure(2)ezplot(f2,t);grid on;3)画出教材P16图1-18,即抽样信号Sa(t)的波形(-20<t<20)。

t=-10:0.01:10;f='sin(t)/t';ezplot(f,t);grid on;4)用符号函数sign画出单位阶跃信号u(t-3)的波形(0<t<10)。

t=0:0.01:10;f='(sign(t-3)+1)/2';ezplot(f,t);grid on;5)单位冲击信号可看作是宽度为∆,幅度为1/∆的矩形脉冲,即t=t 1处的冲击信号为11111 ()()0 t t t x t t t otherδ∆⎧<<+∆⎪=-=∆⎨⎪⎩画出0.2∆=, t 1=1的单位冲击信号。

t=0:0.01:2;f='5*(u(t-1)-u(t-1.2))';ezplot(f,t);grid on;axis([0 2 -1 6]);2、典型离散信号的表示(单位样值序列、单位阶跃序列、实指数序列、正弦序列、复指数序列)编写函数产生下列序列:1)单位脉冲序列,起点n0,终点n f,在n s处有一单位脉冲。

傅里叶实验报告

傅里叶实验报告

傅里叶实验报告傅里叶实验报告引言傅里叶变换是一种重要的数学工具,广泛应用于信号处理、图像处理、物理学等领域。

本实验旨在通过实际操作,深入理解傅里叶变换的原理和应用。

实验设备本实验所需设备包括信号发生器、示波器、计算机等。

实验步骤1. 准备工作首先,我们需要将信号发生器连接到示波器上,以便观察信号的波形。

同时,将示波器与计算机连接,以便进行数据采集和分析。

2. 信号发生器设置将信号发生器的频率设置为50Hz,幅度设置为适当的值。

这样可以产生一个稳定的正弦信号。

3. 示波器设置将示波器的触发方式设置为外部触发,以保证观测到稳定的波形。

同时,调整示波器的水平和垂直缩放,使波形在屏幕上能够清晰显示。

4. 信号采集将示波器的输出信号通过USB接口连接到计算机上,使用相应的软件进行数据采集。

在采集过程中,需要注意保持信号的稳定性,避免干扰。

5. 数据分析将采集到的数据导入到计算机上的数据处理软件中,进行傅里叶变换。

通过傅里叶变换,我们可以将时域信号转换为频域信号,进一步分析信号的频谱特性。

实验结果通过对采集到的数据进行傅里叶变换,我们可以得到信号的频谱图。

从频谱图中,我们可以观察到信号的频率成分和强度分布情况。

通过进一步的分析,我们可以得到信号的频率、幅度、相位等信息。

实验思考傅里叶变换的应用非常广泛,例如在通信领域中,可以通过傅里叶变换将信号从时域转换为频域,从而实现信号的调制和解调。

在图像处理中,傅里叶变换可以用于图像的滤波和压缩。

在物理学中,傅里叶变换可以用于光学、声学等领域的研究。

总结通过本次实验,我们深入了解了傅里叶变换的原理和应用。

傅里叶变换是一种非常重要的数学工具,对于信号处理、图像处理、物理学等领域都具有重要意义。

通过实际操作,我们更加深入地理解了傅里叶变换的工作原理,并通过数据分析得到了实验结果。

通过实验思考,我们发现傅里叶变换在各个领域的应用都非常广泛,对于进一步研究和应用具有重要价值。

傅里叶变换实验报告

傅里叶变换实验报告

一、实验目的1. 理解傅里叶变换的基本原理及其在信号处理中的应用。

2. 掌握傅里叶变换的数学计算方法。

3. 利用MATLAB软件实现傅里叶变换,并对实验结果进行分析。

二、实验原理傅里叶变换是一种重要的信号处理方法,它可以将信号从时域转换到频域。

在频域中,信号的特征更加明显,便于分析和处理。

傅里叶变换的基本原理是将一个信号分解为不同频率的正弦波和余弦波的叠加。

傅里叶变换分为连续傅里叶变换(CFT)和离散傅里叶变换(DFT)。

CFT适用于连续信号,而DFT适用于离散信号。

在本实验中,我们将使用DFT。

三、实验步骤1. 利用MATLAB软件创建一个时域信号,如正弦波、方波或三角波。

2. 对信号进行采样,得到离散信号。

3. 使用MATLAB的fft函数对离散信号进行傅里叶变换。

4. 分析傅里叶变换后的频谱,观察信号在不同频率下的能量分布。

5. 对频谱进行滤波处理,提取感兴趣的特征。

6. 将滤波后的频谱进行逆傅里叶变换,还原信号。

四、实验结果与分析1. 信号创建在本实验中,我们创建了一个频率为50Hz的正弦波信号,采样频率为1000Hz。

2. 傅里叶变换使用MATLAB的fft函数对信号进行傅里叶变换,得到频谱。

观察频谱,发现50Hz 处的能量最大,与信号频率一致。

3. 滤波处理对频谱进行低通滤波,保留50Hz以下的频率成分,滤除高于50Hz的频率成分。

然后对滤波后的频谱进行逆傅里叶变换,还原信号。

观察还原后的信号,发现高频噪声被滤除,信号质量得到提高。

4. 逆傅里叶变换将滤波后的频谱进行逆傅里叶变换,还原信号。

观察还原后的信号,发现其波形与原始信号基本一致,但噪声明显减少。

五、实验结论1. 通过本实验,我们掌握了傅里叶变换的基本原理和计算方法。

2. 利用MATLAB软件可以方便地实现傅里叶变换,并对实验结果进行分析。

3. 傅里叶变换在信号处理中具有广泛的应用,如信号滤波、图像处理、通信等领域。

4. 本实验验证了傅里叶变换在噪声抑制方面的有效性,有助于提高信号质量。

实验三傅里叶变换及其性质

实验三傅里叶变换及其性质

信息工程学院实验报告课程名称:信号与系统实验项目名称:实验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 ωω∞--∞==⎰,傅里叶反变换定义为:11()[()]()2j 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 的函数,而不是默认的ω,即()()jvt F v f t e dt ∞--∞=⎰。

(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傅里叶变换和其性质

实验3傅里叶变换和其性质

信息工程学院实验报告学会运用MATLAB 求连续时间信号的傅里叶(Fourier )变换;学会运用MATLAB 求连续时间信号的频谱图;学会运用MA TLAB 分析连续时间信号的傅里叶变换的性质。

实 验 环 境:Windows 7 MATLAB7.1实 验 内 容 及 过 程:3.1试用MATLAB 命令求下列信号的傅里叶变换,并绘出其幅度谱和相位谱。

(1)1sin 2(1)()(1)t f t t ππ-=- (2)22sin()()t f t t ππ⎡⎤=⎢⎥⎣⎦3.2试用MATLAB 命令求下列信号的傅里叶反变换,并绘出其时域信号图。

(1)1104()35F j j ωωω=-++ (2)224()F e ωω-=3.3试用MATLAB 数值计算方法求门信号的傅里叶变换,并画出其频谱图。

门信号即1,/2()0,/2t g t t τττ⎧≤⎪=⎨>⎪⎩,其中1τ=。

3.4已知两个门信号的卷积为三角波信号,试用MA TLAB 命令验证傅里叶变换 的时域卷积定理。

3.1(1)MATLAB 源程序为:ft = sym('(sin(2*pi*(t-1)))/(pi*(t-1))'); Fw = fourier(ft); subplot(211)ezplot(abs(Fw));grid ontitle('幅度谱')phase = atan(imag(Fw)/real(Fw)); subplot(212)ezplot(phase);grid on title('相位谱')波形图如图T3-1所示:图T3-1 3.1(2)MATLAB源程序为:t = sym('((sin(pi*t))/(pi*t))^2'); Fw = fourier(ft);subplot(211)ezplot(abs(Fw));grid ontitle('幅度谱')phase = atan(imag(Fw)/real(Fw)); subplot(212)ezplot(phase);grid ontitle('相位谱')波形图如图T3-2所示:图T3-2 3.2(1)MATLAB源程序为:t=sym('t');Fw = sym('10/(3+i*w)-4/(5+i*w)'); ft = ifourier(Fw);ezplot(ft);grid on;axis([-1 3 -1 7]);波形图如图T3-3所示:图T3-3 3.2(2)MATLAB源程序为:t=sym('t');Fw = sym('exp(-4*w^2)');ft = ifourier(Fw);ezplot(ft)grid on波形图如图T3-4所示:图T3-4 3.3、MATLAB源程序为:dt = 0.001;t = -0.5:dt:0.5;ft = uCT(t+0.5)-uCT(t-0.5);N = 2000;k = -N:N;W = 2*pi*k/((2*N+1)*dt);F = dt * ft*exp(-j*t'*W);plot(W,F);grid on;axis([-pi pi -1 3]);xlabel('W'), ylabel('F(W)')title('amplitude spectrum');波形图如图T3-5所示:图T3-53.4将门函数先进行时域卷积运算,再将卷积后的结果做傅里叶变换,程序和结果如下:clear;clc;dt = 0.01; t = -2:dt:2.5;f1 = uCT(t+0.5)- uCT(t-0.5);f = conv(f1,f1)*dt;ft=sym('f');Fw = fourier(ft)Fw =2*i*pi*dirac(1,w)将一个门函数先进行傅里叶变换,再将结果与自身相乘,程序和结果如下:clear;clc;dt = 0.01; t = -2:dt:2.5;f1 = uCT(t+0.5)- uCT(t-0.5);ft=sym('f1');Fw = fourier(ft);Fw=Fw*FwFw =-4*pi^2*dirac(1,w)^2由此来验证傅里叶变换的时域卷积定理。

信号与系统实验报告-傅立叶变换

信号与系统实验报告-傅立叶变换

实验三连续时间傅里叶级数目的:本练习要检验连续时间傅里叶级数(CTFS)的性质中等题:1.满足x1(t)=x1(t+T)的最小周期T是多少?利用这个T值,用解析法求x1(t)的CTFS系数。

程序:linspace(-1,1,1000);x=sym('cos(2*pi*t)')y=sym('sin(4*pi*t)')x1=x+yezplot(x1,[-1,1]);grid波形:x1(t)程序:t=linspace(-1,1,1000);x=sym('cos(2*pi*t)')y=sym('sin(4*pi*t)')x1=x+yezplot(x1,[-1,1]);grid %可以从图形求出T=1 k=[-5:5]syms te=exp(-i*2*pi*t*k)f=x1*eFn=int(f,t,0,1);F=abs(Fn) %求Fn系数的绝对值subplot(1,2,1);Fn=double(Fn)subplot(1,2,1)stem(k,Fn)set(gca,'YLim',[-1 1.2])set(gca,'ytick',[0.1 0.2 0.3 0.4 0.5 0.6 1])set(gca,'xtick',[-5 -4 -3 -2 -1 0 1 2 3 4 5])gridF=double(F)subplot(1,2,2)stem(k,F)set(gca,'YLim',[-1 1.2])set(gca,'ytick',[0.1 0.2 0.3 0.4 0.5 0.6 1])set(gca,'xtick',[-5 -4 -3 -2 -1 0 1 2 3 4 5])grid波形:x1(t)的CTFS系数为Fn(左)以及其绝对值,比较可以得出虚部的系数。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

信息工程学院实验报告课程名称:实验项目名称:实验3 傅里叶变换及其性质 实验时间:2015/11/17 班级:通信141 姓名: 学号:一、实 验 目 的:学会运用MATLAB 求连续时间信号的傅里叶(Fourier )变换;学会运用MATLAB 求连续时间信号的频谱图;学会运用MATLAB 分析连续时间信号的傅里叶变换的性质。

二、实 验 设 备 与 器 件 软件:Matlab 2008 三、实 验 原 理 3.1傅里叶变换的实现信号()f t 的傅里叶变换定义为: ()[()]()j t F F f t f t e dt ωω∞--∞==⎰,傅里叶反变换定义为:11()[()]()2j 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 的函数,而不是默认的ω,即()()jvt F v f t e dt ∞--∞=⎰。

(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 。

值得注意的是,函数fourier( )和ifourier( )都是接受由sym 函数所定义的符号变量或者符号表达式。

3.1.2连续时间信号的频谱图信号()f t 的傅里叶变换()F ω表达了信号在ω处的频谱密度分布情况,这就是信号的傅里叶变换的物理含义。

()F ω一般是复函数,可以表示成()()()j F F eϕωωω=。

()~F ωω与()~ϕωω曲线分别称为非周期信号的幅度频谱与相位频谱,它们都是频率ω的连续函数,在形状上与相应的周期信号频谱包络线相同。

非周期信号的频谱有两个特点,密度谱和连续谱。

要注意到,采用fourier()和ifourier() 得到的返回函数,仍然是符号表达式。

若需对返回函数作图,则需应用ezplot()绘图命令。

3.1.3 MATLAB 数值计算求解法fourier( )和ifourier( )函数的一个局限性是,如果返回函数中有诸如单位冲激函数()t δ等项,则用ezplot()函数无法作图。

对某些信号求变换时,其返回函数可能包含一些不能直接用符号表达的式子,因此不能对返回函数作图。

此外,在很多实际情况中,尽管信号()f t 是连续的,但经过抽样所获得的信号则是多组离散的数值量()f n ,因此无法表示成符号表达式,此时不能应用fourier()函数对f(n)进行处理,而只能用数值计算方法来近似求解。

从傅里叶变换定义出发有0()()lim ()j tj n F f t edt f n e ωωω∞∞-∞∆→-∞--∆==∆∆∑⎰,当∆足够小时,上式的近似情况可以满足实际需要。

对于时限信号()f t ,或者在所研究的时间范围内让()f t 衰减到足够小,从而近似地看成时限信号,则对于上式可以考虑有限n 的取值。

假设是因果信号,则有1()(),01M n j n F f n e n M ωω-=-∆=∆∆≤≤-∑傅里叶变换后在ω域用MATLAB 进行求解,对上式的角频率ω进行离散化。

假设离散化后得到N 个样值,即 2,0k k k N N πω=≤≤∆-1,因此有 1()(),01M n k j n F k f n ek N ω-=-∆=∆∆≤≤-∑。

采用行向量,用矩阵表示为1*1**[()][()][]k j n T TT N M M N F k f n eω-∆=∆∆。

其要点是要正确生成()f t 的M 个样本向量[()]f n ∆与向量[]j n k eω-∆。

当∆足够小时,上式的内积运算(即相乘求和运算)结果即为所求的连续时间信号傅里叶变换的数值解。

3.2傅里叶变换的性质傅里叶变换的性质包含了丰富的物理意义,并且揭示了信号的时域和频域的关系。

熟悉这些性质成为信号分析研究工作中最重要的内容之一。

3.2.1 尺度变换特性傅里叶变换的尺度变换特性为:若()()f t F ω↔,则有1()()f at F a aω↔,其中,a 为非零实常数。

3.2.2频移特性傅里叶变换的频移特性为:若()()f t F ω↔,则有00()()j tf t eF ωωω↔-。

频移技术在通信系统中得到广泛应用,诸如调幅变频等过程都是在频谱搬移的基础上完成的。

频移的实现原理是将信号()f t 乘以载波信号0cos t ω或0sin t ω,从而完成频谱的搬移,即0000001()cos [()()]2()sin [()()]2f t t F F jf t t F F ωωωωωωωωωω↔++-↔+--四、实 验 内 容 与 步 骤4.1试用MATLAB 命令求下列信号的傅里叶变换,并绘出其幅度谱和相位谱。

(1)1sin 2(1)()(1)t f t t ππ-=- (2)22sin()()t f t t ππ⎡⎤=⎢⎥⎣⎦4.2试用MATLAB 命令求下列信号的傅里叶反变换,并绘出其时域信号图。

(1)1104()35F j j ωωω=-++ (2)224()F e ωω-=4.3试用MATLAB 数值计算方法求门信号的傅里叶变换,并画出其频谱图。

门信号即1,/2()0,/2t g t t τττ⎧≤⎪=⎨>⎪⎩,其中1τ=。

4.4已知两个门信号的卷积为三角波信号,试用MATLAB 命令验证傅里叶变换的时域卷积定理。

5.问题与思考傅里叶变换的其他性质可以用类似的方法加以验证,试举一例,说明你验证过程的思路。

解:4.1(1)MATLAB源程序为:clear;clc;ft=sym('sin(2*pi*(t-1))/(pi*(t-1))');Fw = fourier(ft);subplot(211)ezplot(abs(Fw),[-5*pi 5*pi]);grid ontitle('幅度谱');phase = atan(imag(Fw)/real(Fw));subplot(212)ezplot(phase);grid ontitle('相位谱');4.1(2)MA TLAB源程序为:clear;clc;ft = sym('(sin(pi*t)/(pi*t))^2');Fw = fourier(ft);subplot(211)ezplot(abs(Fw));grid ontitle('幅度谱');phase = atan(imag(Fw)/real(Fw));subplot(212)ezplot(phase);grid ontitle('相位谱');4.2(1)MA TLAB源程序为:clear;clc;t=sym('t');Fw= sym('10/(3+i*w)-4/(5+i*w)');ft = ifourier(Fw);ezplot(ft),grid on4.2(2)MA TLAB源程序为:clear;clc;t=sym('t');Fw = sym('exp(-4*(w^2))');ft = ifourier(Fw);ezplot(ft),grid on4.3 MA TLAB源程序为:clear;clc;ft1=sym('Heaviside(t+1/2)-Heaviside(t-1/2)');subplot(121);ezplot(ft1,[-pi pi]),grid onFw1 = simplify(fourier(ft1));subplot(122);ezplot(abs(Fw1),[-10*pi 10*pi]), grid onaxis([-10*pi 10*pi -0.2 1.2]);4.4两个门信号卷积成为三角波信号的实验程序代码:clear;clc;dt = 0.01; t = -1:dt:2.5;f1 = uCT(t+1/2)- uCT(t-1/2);f2 = uCT(t+1/2)- uCT(t-1/2);f = conv(f1,f2)*dt;n =length(f);tt = (0:n-1)*dt-2;subplot(211), plot(t,f1),grid on;axis([-1, 1, -0.2,1.2]);title('f1(t)'); xlabel('t');subplot(212), plot(tt,f),grid on;axis([-2, 2, -0.2,1.2]);title('f(t)=f1(t)*f2(t)'); xlabel('t');两个门信号卷积成为三角波信号的实验结果如图6所示:图6三角波信号傅里叶变换的实验程序代码:clear;clc;dt = 0.01;t = -4:dt:4;ft=(t+1).*uCT(t+1)-2*t.*uCT(t)+(t-1).*uCT(t-1);N = 2000;k = -N:N;W = 2*pi*k/((2*N+1)*dt);F = dt * ft*exp(-j*t'*W);plot(W,F), grid onaxis([-10*pi 10*pi -0.2 1.2]);xlabel('W'), ylabel('F(W)')title('f1(t)*f2(t)的频谱图');ft1和ft2分别傅里叶变换然后再相乘的代码:clear;clc;ft1=sym('Heaviside(t+1/2)-Heaviside(t-1/2)');Fw1=fourier(ft1);ft2=sym('Heaviside(t+1/2)-Heaviside(t-1/2)');Fw2 = fourier(ft2);Fw=Fw1.*Fw2;ezplot(Fw,[-10*pi 10*pi]);grid onaxis([-10*pi 10*pi -0.2 1.2]);三角波信号傅里叶变换的实验结果如图7所示,ft1和ft2分别傅里叶变换然后再相乘的实验结果如图8所图7 图8图7和图8几乎是一样的,所以傅里叶变换的时域卷积定理是正确的。

相关文档
最新文档