一、实验目的 1在理论学习的基础上,通过本实验加深对快速傅立叶变换的理解; 2熟悉并掌握按时间抽取FFT算法的程序; 3了解应用FFT进行信号频谱分析过程中可能出现的问题,例如混淆、泄漏、栅栏效应等,以便在实际中正确应用FFT。 二、实验内容 1仔细分析教材第六章‘时间抽取法FFT ’的算法结构,编制出相应的用FFT 进行信号分析的C语言(或MATLAB 语言)程序; 用MATLAB语言编写的FFT源程序如下: %% 输入数据f、N、T及是否补零 clc; clear; f=input('输入信号频率f:'); N=input('输入采样点数N:'); T=input('输入采样间隔T:'); C=input('信号是否补零(补零输入1,不补零输入0):'); %补零则输入1,不补则输入0 if(C==0) t=0:T:(N-1)*T; x=sin(2*pi*f*t); b=0; e lse b=input('输入补零的个数:'); while(log2(N+b)~=fix(log2(N+b))) b=input('输入错误,请重新输入补零的个数:'); end t=0:T:(N+b-1)*T; x=sin(2*pi*f*t).*(t<=(N-1)*T); end %% fft算法的实现 A=bitrevorder(x); % 将序列按二进制倒序 N=N+b; M=log2(N); % M为蝶形算法的层数 W=exp(-j*2*pi/N); for L=1:1:M %第L层蝶形算法 B=2^L/2; % B为每层蝶形算法进行加减运算的两个数的间隔 K=N/(2^L); % K为每层蝶形算法中独立模块的个数 for k=0:1:K-1 for J=0:1:B-1
傅里叶分析之掐死教程(完整版)更新于2014.06.06 Heinrich · 6 个月前 作者:韩昊知乎:Heinrich 微博:@花生油工人知乎专栏:与时间无关的故事 谨以此文献给大连海事大学的吴楠老师,柳晓鸣老师,王新年老师以及张晶泊老师。 转载的同学请保留上面这句话,谢谢。如果还能保留文章来源就更感激不尽了。 我保证这篇文章和你以前看过的所有文章都不同,这是12年还在果壳的时候写的,但是当时没有来得及写完就出国了……于是拖了两年,嗯,我是拖延症患者…… 这篇文章的核心思想就是: 要让读者在不看任何数学公式的情况下理解傅里叶分析。 傅里叶分析不仅仅是一个数学工具,更是一种可以彻底颠覆一个人以前世界观的思维模式。但不幸的是,傅里叶分析的公式看起来太复杂了,所以很多大一新生
上来就懵圈并从此对它深恶痛绝。老实说,这么有意思的东西居然成了大学里的杀手课程,不得不归咎于编教材的人实在是太严肃了。(您把教材写得好玩一点会死吗?会死吗?)所以我一直想写一个有意思的文章来解释傅里叶分析,有可能的话高中生都能看懂的那种。所以,不管读到这里的您从事何种工作,我保证您都能看懂,并且一定将体会到通过傅里叶分析看到世界另一个样子时的快感。至于对于已经有一定基础的朋友,也希望不要看到会的地方就急忙往后翻,仔细读一定会有新的发现。 ——————————————以上是定场诗—————————————— 下面进入正题: 抱歉,还是要啰嗦一句:其实学习本来就不是易事,我写这篇文章的初衷也是希望大家学习起来更加轻松,充满乐趣。但是千万!千万不要把这篇文章收藏起来,或是存下地址,心里想着:以后有时间再看。这样的例子太多了,也许几年后你都没有再打开这个页面。无论如何,耐下心,读下去。这篇文章要比读课本要轻松、开心得多…… p.s.本文无论是cos还是sin,都统一用“正弦波”(Sine Wave)一词来代表简谐波。 一、什么是频域 从我们出生,我们看到的世界都以时间贯穿,股票的走势、人的身高、汽车的轨迹都会随着时间发生改变。这种以时间作为参照来观察动态世界的方法我们称其为时域分析。而我们也想当然的认为,世间万物都在随着时间不停的改变,并且永远不会静止下来。但如果我告诉你,用另一种方法来观察世界的话,你会发现世界是永恒不变的,你会不会觉得我疯了?我没有疯,这个静止的世界就叫做频域。 先举一个公式上并非很恰当,但意义上再贴切不过的例子: 在你的理解中,一段音乐是什么呢?
实验报告 课程名称信号与线性系统分析 实验名称周期信号的傅里叶级数和频谱分析实验类型验证(验证、综合、设计、创新) 3日实验四、周期信号的傅里叶级数和频谱分析1实验目的 1)学会利用MATLAB分析傅里叶级数展开,并理解傅里叶级数的物理含义; 2)学会利用MATLAB分析周期信号的频谱特性。 2实验原理及实例分析 周期信号可以再函数的区间里展成在完备正交信号空间中的无穷级数。如果完备的正交函数集是三角函数集或指数函数集,那么,周期信号所展开的无穷级数就分别成为“三角型傅里叶级数”或“指数型傅里叶级数”,统称为傅里叶级数。
2.1周期信号的傅里叶级数 (基本原理请参阅教材第四章的4.1节和4.2节。) 例1:周期方波信号)(t f 如图1所示,试求出该信号的傅里叶级数,利用MATLAB 编程实现其各次谐波的叠加,并验证Gibbs 现象。 图1 周期方波信号)(t f 的波形图 解:从理论分析可知,周期方波信号)(t f 的傅里叶级数展开式为 )9sin 9 17sin 715sin 513sin 31(sin 4)(00000 +++++=t t t t t t f ωωωωωπ 其中,ππω220== T 。则可分别求出1、3、5、9、19、39、79、159项傅里叶级数求和的结果,其MATLAB 程序如下,产生的图形如图2所示。 close all;clear all; clc t = -2:0.0001:2; omega = 2 * pi; y = square(2 * pi * t,50); n_max = [1 3 5 9 19 39 79 159]; N = length(n_max); for k = 1:N fk = zeros(1,length(t)); for n = 1:2:n_max(k) bn = 4 / (pi * n); fk = fk + bn * sin(n * omega * t); end figure;plot(t,y,t,fk,'Linewidth',2); xlabel('t(sec)');ylabel('部分和的波形'); f(t) t(sec)
MATLAB实验傅里叶分析
实验七 傅里叶变换 一、实验目的 傅里叶变换是通信系统、图像处理、数字信号处理以及物理学等领域内的一种重要的数学分析工具。通过傅里叶变换技术可以将时域上的波形分 布变换为频域上的分布,从而获得信号的频谱特性。MATLAB 提供了专门的函数fft 、ifft 、fft2(即2维快速傅里叶变换)、ifft2以及fftshift 用于实现对信号的傅里叶变换。本次实验的目的就是练习使用fft 、ifft 以及fftshift 函数,对一些简单的信号处理问题能够获取其频谱特性(包括幅频和相频特性)。 二、实验预备知识 1. 离散傅里叶变换(DFT)以及快速傅里叶变换(FFT)简介 设x (t )是给定的时域上的一个波形,则其傅里叶变换为 2()() (1)j ft X f x t e dt π∞--∞=? 显然X ( f )代表频域上的一种分布(波形),一般来说X ( f )是复数。而傅里叶逆变换定义为: 2()() (2)j ft x t X f e df π∞-∞ =?
因此傅里叶变换将时域上的波形变换为频域上的波形,反之,傅里叶逆变换则将频域上的波形变换为时域上的波形。 由于傅里叶变换的广泛应用,人们自然希望能够使用计算机实现傅里叶变换,这就需要对傅里叶变换(即(1)式)做离散化处理,使 之符合电脑计算的特征。另外,当 把傅里叶变换应用于实验数据的分 析和处理时,由于处理的对象具有 离散性,因此也需要对傅里叶变换 进行离散化处理。而要想将傅里叶 变换离散化,首先要对时域上的波 形x (t )进行离散化处理。采用一个 时域上的采样脉冲序列: δ (t -nT ), n = 0, 1, 2, …, N -1; 可以实现上述目的,如图所示。其中N 为采样点数,T 为采样周期;f s = 1/T 是采样频率。注意采样时,采样频率f s 必须大于两倍的信号频率(实际是截止频率),才能避免混迭效应。 接下来对离散后的时域波形()()()(x t x t t n T x n T δ= -=的傅里叶变换()X f 进行离散处理。与上述做法类 似,采用频域上的δ脉冲序列: x (t δ x (t )δ t t t
实验二傅里叶分析及应用 姓名学号班级 一、实验目的 (一)掌握使用Matlab进行周期信号傅里叶级数展开和频谱分析 1、学会使用Matlab分析傅里叶级数展开,深入理解傅里叶级数的物理含义 2、学会使用Matlab分析周期信号的频谱特性 (二)掌握使用Matlab求解信号的傅里叶变换并分析傅里叶变换的性质 1、学会运用Matlab求连续时间信号的傅里叶变换 2、学会运用Matlab求连续时间信号的频谱图 3、学会运用Matlab分析连续时间信号的傅里叶变换的性质 (三)掌握使用Matlab完成信号抽样并验证抽样定理 1、学会运用MATLAB完成信号抽样以及对抽样信号的频谱进行分析 2、学会运用MATLAB改变抽样时间间隔,观察抽样后信号的频谱变化 3、学会运用MATLAB对抽样后的信号进行重建 二、实验条件 需要一台PC机和一定的matlab编程能力 三、实验内容 2、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT,并画出其频谱图(包括幅度谱和相位谱)[注:图中时间单位为:毫秒(ms)]。
符号运算法: Ft= sym('t*(Heaviside(t+2)-Heaviside(t+1))+Heaviside(t+1)-Heaviside(t-1)+(-t)*(Heavi side(t-1)-Heaviside(t-2))'); Fw = fourier(Ft); ezplot(abs(Fw)),grid on; phase = atan(imag(Fw)/real(Fw)); ezplot(phase);grid on; title('|F|'); title('phase'); 3、试用Matlab 命令求ω ωωj 54 -j 310)F(j ++= 的傅里叶反变换,并绘出其时域信号图。
实验二傅里叶分析及应用 、实验目的 (一)掌握使用Matlab 进行周期信号傅里叶级数展开和频谱分析 1、学会使用Matlab 分析傅里叶级数展开,深入理解傅里叶级数的物理含义 2、学会使用Matlab 分析周期信号的频谱特性 二)掌握使用Matlab 求解信号的傅里叶变换并分析傅里叶变换的性质 1、学会运用Matlab 求连续时间信号的傅里叶变换 2、学会运用Matlab 求连续时间信号的频谱图 3、学会运用Matlab 分析连续时间信号的傅里叶变换的性质 三)掌握使用Matlab 完成信号抽样并验证抽样定理 1、学会运用MATLAB 完成信号抽样以及对抽样信号的频谱进行分析 2、学会运用MATLAB 改变抽样时间间隔,观察抽样后信号的频谱变化 3、学会运用MATLAB 对抽样后的信号进行重建 、实验条件 Win7系统,MATLAB R2015a 三、实验内容 1、分别利用Matlab 符号运算求解法和数值计算法求下图所示信号的FT,并画出其频谱图(包括幅度谱和相位谱)[注:图中时间单位为:毫秒(ms)]。
Code: ft = sym( ' (t+2)*(heaviside(t+2)-heavisi de(t+1))+(heaviside(t+1)-heav iside(t- 1))+(2-t)*(heaviside( t-1)-heaviside(t- 2))' ); fw = simplify(fourier(ft)); subplot(2, 1, 1); ezplot(abs(fw)); grid on; title( 'amp spectrum' ); phi = atan(imag(fw) / real(fw)); subplot(2, 1, 2); ezplot(phi); grid on ; title( 'phase spectrum' ); 符号运算法 Code: dt = 0.01; t = -2: dt: 2; ft (t+2).*(uCT(t+2)- uCT(t+1))+(u CT(t+1)-uCT(t- 1))+(2-t).*(uCT (t-1)- uCT(t-2)); N = 2000; k = -N: N; w = pi * k / (N*dt); fw = dt*ft*exp(-i*t'*w); fw = abs(fw); plot(w, fw), grid on; axis([-2*pi 2*pi -1 3.5]); 数值运算法
一、傅立叶变化的原理; (1)原理 正交级数的展开是其理论基础!将一个在时域收敛的函数展开成一系列不同频率谐波的叠加,从而达到解决周期函数问题的目的。在此基础上进行推广,从而可以对一个非周期函数进行时频变换。 从分析的角度看,他是用简单的函数去逼近(或代替)复杂函数,从几何的角度看,它是以一族正交函数为基向量,将函数空间进行正交分解,相应的系数即为坐标。从变幻的角度的看,他建立了周期函数与序列之间的对应关系;而从物理意义上看,他将信号分解为一些列的简谐波的复合,从而建立了频谱理论。 当然Fourier积分建立在傅氏积分基础上,一个函数除了要满足狄氏条件外, 一般来说还要在积分域上绝对可积,才有古典意义下的傅氏变换。引入衰减因子e^(-st),从而有了Laplace变换。(好像走远了)。 (2)计算方法 连续傅里叶变换将平方可积的函数f(t)表示成复指数函数的积分或级数形式。 这是将频率域的函数F(ω)表示为时间域的函数f(t)的积分形式。 为 连续傅里叶变换的逆变换 (inverse Fourier transform) 即将时间域的函数f(t)表示为频率域的函数F(ω)的积分。 一般可称函数f(t)为原函数,而称函数F(ω)为傅里叶变换的像函数,原函数和像函数构成一个傅里叶变换对(transform pair)。 二、傅立叶变换的应用; DFT在诸多多领域中有着重要应用,下面仅是颉取的几个例子。需要指出 的是,所有DFT的实际应用都依赖于计算离散傅里叶变换及其逆变换的快速算
法,即快速傅里叶变换(快速傅里叶变换(即FFT )是计算离散傅里叶变换及其逆变换的快速算法。)。(1)、频谱分析DFT 是连续傅里叶变换的近似。因此可以对连续信号x(t)均匀采样并截断以得到有限长的离散序列,对这一序列作离散傅里叶变换,可以分析连续信号x(t)频谱的性质。前面还提到DFT 应用于频谱分析需要注意的两个问题:即采样可能导致信号混叠和截断信号引起的频谱泄漏。可以通过选择适当的采样频率(见奈奎斯特频率)消减混叠。选择适当的序列长度并加窗可以抑制频谱泄漏。(2)、数据压缩由于人类感官的分辨能力存在极限,因此很多有损压缩算法利用这一点将语音、音频、图像、视频等信号的高频部分除去。高频信号对应于信号的细节,滤除高频信号可以在人类感官可以接受的范围内获得很高的压缩比。这一去除高频分量的处理就是通过离散傅里叶变换完成的。将时域或空域的信号转换到频域,仅储存或传输较低频率上的系数,在解压缩端采用逆变换即可重建信号。(3)、OFDM OFDM (正交频分复用)在宽带无线通信中有重要的应用。这种技术将带宽为N 个等间隔的子载波,可以证明这些子载波相互正交。尤其重要的是,OFDM 调制可以由IDFT 实现,而解调可以由DFT 实现。OFDM 还利用DFT 的移位性质,在每个帧头部加上循环前缀(Cyclic Prefix ),使得只要信道延时小于循环前缀的长度,就能消除信道延时对传输的影响。三、傅里叶变换的本质; 傅里叶变换的公式为dt e t f F t j ?+∞∞--=ωω)()(可以把傅里叶变换也成另外一种形式: t j e t f F ωπ ω),(21)(=可以看出,傅里叶变换的本质是内积,三角函数是完备的正交函数集,不同频率的三 角函数的之间的内积为0,只有频率相等的三角函数做内积时,才不为0。)(2,21)(2121Ω-Ω==?Ω-ΩΩΩπδdt e e e t j t j t j
实验三周期信号的傅里叶级数分析及MATLAB实现 一、实验目的: 1.利用MATLAB实现周期信号的分解与合成,并图示仿真结果; 2.用MATLAB实现周期信号的频谱,画图观察和分析周期信号的频谱; 3.通过MATLAB对周期信号频谱的仿真,进一步加深对周期信号频谱理论知识的理解。 二、实验内容 9.1(a):程序: display('Please input the value of m(傅里叶级数展开项数)'); m=input('m='); t=-3*pi:0.01:3*pi; n=round(length(t)/4); f=cos(t).*(heaviside(t+2.5*pi)-heaviside(t+1.5*pi)+heaviside(t+0.5*pi)-heaviside(t-0.5 *pi)+heaviside(t-1.5*pi)-heaviside(t-2.5*pi)); y=zeros(m+1,max(size(t))); y(m+1,:)=f'; figure(1); plot(t/pi,y(m+1,:)); grid; axis([-3 3 -1 1.5]); title('半波余弦'); xlabel('单位:pi','Fontsize',8); x=zeros(size(t)); kk='1'; syms tx n T=2*pi; fx=sym('cos(tx)'); Nn=30; An=zeros(m+1,1); Bn=zeros(m+1,1); a0=2*int(fx,tx,-T/4,T/4)/T an=2*int(fx*cos(2*pi*(n+eps/2)*tx/T),tx,-T/4,T/4)/T bn=2*int(fx*sin(2*pi*(n+eps/2)*tx/T),tx,-T/4,T/4)/T An(1)=double(vpa(a0,Nn)); An(2)=0.5; for k=2:m An(k+1)=double(vpa(subs(an,n,k),Nn)); Bn(k+1)=double(vpa(subs(bn,n,k),Nn));
实验七 傅里叶变换 一、实验目的 傅里叶变换是通信系统、图像处理、数字信号处理以及物理学等领域内的一种重要的数学分析工具。通过傅里叶变换技术可以将时域上的波形分 布变换为频域上的分布,从而获得信号的频谱特性。MA TLAB 提供了专门的函数fft 、ifft 、fft2(即2维快速傅里叶变换)、ifft2以及fftshift 用于实现对信号的傅里叶变换。本次实验的目的就是练习使用fft 、ifft 以及fftshift 函数,对一些简单的信号处理问题能够获取其频谱特性(包括幅频和相频特性)。 二、实验预备知识 1. 离散傅里叶变换(DFT)以及快速傅里叶变换(FFT)简介 设x (t )是给定的时域上的一个波形,则其傅里叶变换为 2()() (1)j ft X f x t e dt π∞ --∞ =? 显然X ( f )代表频域上的一种分布(波形),一般来说X ( f )是复数。而傅里叶逆变换定义为: 2()() (2)j ft x t X f e df π∞ -∞ =? 因此傅里叶变换将时域上的波形变换为频域上的波形,反之,傅里叶逆变换则将频域上的波形变换为时域上的波形。 由于傅里叶变换的广泛应用,人们自然希望能够使用计算机实现傅里叶变换,这就需要对傅里叶变换(即(1)式)做离散化处理,使之符合电脑计算的特征。另外,当把傅里叶变换应用于实验数据的分析和处理时,由于处理的对象具有离散性,因此也需要对傅里叶变换进行离散化处理。而要想将傅里叶变换离散化,首先要对时域上的波形x (t )进行离散化处理。采用一个时域上的采样脉冲序列: δ (t -nT ), n = 0, 1, 2, …, N -1; 可以实现上述目的,如图所示。其中N 为采样点数,T 为采样周期;f s = 1/T 是采样频率。注意采样时,采样频率f s 必须大于两倍的信号频率(实际是截止频率),才能避免混迭效应。 接下来对离散后的时域波形()()()()x t x t t nT x nT δ=-=的傅里叶变换()X f 进行离散处理。与上述做法类似,采用频域上的δ脉冲序列: δ ( f -n/T 0), n = 0, 1, 2, …, N -1;T 0= NT 为总采样时间 可以实现傅里叶变换()X f 的离散化,如下图示。不难看出,离散后的傅里叶变换其频率间隔(频率轴上离散点的间隔,即频域分辨率) x (t ) δ 脉冲序列 x (t )δ (t -nT ) t t t
实验二傅里叶分析及应用 -、实验目的 (一)掌握使用Matlab进行周期信号傅里叶级数展开和频谱分析 1、学会使用Matlab分析傅里叶级数展开,深入理解傅里叶级数的物理含义 2、学会使用Matlab分析周期信号的频谱特性 (二)掌握使用Matlab求解信号的傅里叶变换并分析傅里叶变换的性质 1、学会运用Matlab求连续时间信号的傅里叶变换 2、学会运用Matlab求连续时间信号的频谱图 3、学会运用Matlab分析连续时间信号的傅里叶变换的性质 (三)掌握使用Matlab完成信号抽样并验证抽样定理 1、学会运用MATLAB完成信号抽样以及对抽样信号的频谱进行分析 2、学会运用MATLAB改变抽样时间间隔,观察抽样后信号的频谱变化 3、学会运用MATLAB对抽样后的信号进行重建 、实验条件 Win7 系统,MATLAB R2015a 三、实验内容 1、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT,并画出其频
谱图(包括幅度谱和相位谱)
Code: ft = sym( ' (t+2)*(heaviside(t+2)-heavisi de(t+1))+(heaviside(t+1)-heav iside(t- 1))+(2-t)*(heaviside( t-1)-heaviside(t-2))' ); fw = simplify(fourier(ft)); subplot(2, 1, 1); ezplot(abs(fw)); gridon ; title( 'amp spectrum' ); phi = atan(imag(fw) / real(fw)); subplot(2, 1,2); ezplot(phi); grid on ; title( 'phase spectrum' ); 符号运算法 Code: dt = 0.01; t = -2: dt: 2; ft (t+2).*(uCT(t+2)-uCT(t+1))+(u CT(t+1)-uCT(t-1))+(2-t).*(uCT (t-1)-uCT(t-2)); N = 2000; k = -N: N; w = pi * k / (N*dt); fw = dt*ft*exp(-i*t'*w); fw = abs(fw); plot(w, fw), grid on; axis([-2*pi 2*pi -1 3.5]); 数值运算法amp spectrum -6-4 -2 0 2 4 6 w x10 phase spectrurri -6 -4 -2 0 2 4 6 w
第二章习题 习题2.1 设随机过程X (t )可以表示成: ()2cos(2), X t t t πθ=+-∞<<∞ 式中,θ是一个离散随机变量,它具有如下概率分布:P (θ=0)=0.5,P (θ=π/2)=0.5 试求E [X (t )]和X R (0,1)。 解:E [X (t )]=P (θ=0)2cos(2)t π+P (θ=/2)2cos(2)=cos(2)sin 22 t t t π πππ+ - cos t ω 习题2.2 设一个随机过程X (t )可以表示成: ()2cos(2), X t t t πθ=+-∞<<∞ 判断它是功率信号还是能量信号?并求出其功率谱密度或能量谱密度。 解:为功率信号。 []/2 /2/2 /21()lim ()()1lim 2cos(2)*2cos 2()T X T T T T T R X t X t dt T t t dt T ττπθπτθ→∞ -→∞-=+=+++? ? 222cos(2)j t j t e e πππτ-==+ 2222()()()(1)(1) j f j t j t j f X P f R e d e e e d f f πτπππττττδδ∞-∞---∞-∞==+=-++?? 习题2.3 设有一信号可表示为: 4exp() ,t 0 (){0, t<0 t X t -≥= 试问它是功率信号还是能量信号?并求出其功率谱密度或能量谱密度。 解:它是能量信号。X (t )的傅立叶变换为: (1)004 ()()441j t t j t j t X x t e dt e e dt e dt j ωωωωω +∞-+∞--+∞-+-∞====+??? 则能量谱密度 G(f)=2 ()X f =2 22 416114j f ωπ=++ 习题2.4 X (t )=12cos 2sin 2x t x t ππ-,它是一个随机过程,其中1x 和2x 是相互统计独立的高斯随机变量,数学期望均为0,方差均为2σ。试求: (1)E [X (t )],E [2()X t ];(2)X (t ) 的概率分布密度;(3)12(,)X R t t 解:(1)()[][]()[]02sin 2cos 2sin 2cos 2121=?-?=-=x E t x E t t x t x E t X E ππππ ()X P f 因为21x x 和相互独立,所以[][][]2121x E x E x x E ?=。
实验二傅里叶分析及应用 一、实验目的 (一)掌握使用Matlab进行周期信号傅里叶级数展开和频谱分析 1、学会使用Matlab分析傅里叶级数展开,深入理解傅里叶级数的物理含义 2、学会使用Matlab分析周期信号的频谱特性 (二)掌握使用Matlab求解信号的傅里叶变换并分析傅里叶变换的性质 1、学会运用Matlab求连续时间信号的傅里叶变换 2、学会运用Matlab求连续时间信号的频谱图 3、学会运用Matlab分析连续时间信号的傅里叶变换的性质 (三)掌握使用Matlab完成信号抽样并验证抽样定理 1、学会运用MATLAB完成信号抽样以及对抽样信号的频谱进行分析 2、学会运用MATLAB改变抽样时间间隔,观察抽样后信号的频谱变化 3、学会运用MATLAB对抽样后的信号进行重建 二、实验条件 Win7系统,MATLAB R2015a 三、实验内容 1、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT,并画出其频谱图(包括幅度谱和相位谱)[注:图中时间单位为:毫秒(ms)]。
符号运算法 数值运算法
t (20 π ex p(-3 t) heaviside(t) - 8 π ex p(-5 t) heaviside(t))/(2 π) 2、试用Matlab 命令求ω ω ωj 54 -j 310)F(j ++= 的傅里叶反变换,并绘出其时域信号图。 两个单边指数脉冲的叠加 3、已知门函数自身卷积为三角波信号,试用Matlab 命令验证FT 的时域卷积定理。
4、设有两个不同频率的余弦信号,频率分别为Hz f 1001=,Hz f 38002=;现在使用抽样频率Hz f s 4000=对这三个信号进行抽样,使用MATLAB 命令画出各抽样信号的波形和频谱,并分析其频率混叠现象 > > > > > > > > > >
实验二 用快速傅里叶变换对信号进行频谱分析 一、实验目的 1.理解离散傅里叶变换的意义; 2.掌握时域采样率的确定方法; 3.掌握频域采样点数的确定方法; 4.掌握离散频率与模拟频率之间的关系; 5.掌握离散傅里叶变换进行频谱分析时,各参数的影响。 二、实验原理 序列的傅里叶变换结果为序列的频率响应,但是序列的傅里叶变换是频率的连续函数,而且在采用计算机计算时,序列的长度不能无限长,为了便于计算机处理,作如下要求:序列x (n )为有限长,n 从0~N -1,再对频率ω在0~2π范围内等间隔采样,采样点数为N ,采样间隔为2π/N 。第k 个采样点对应的频率值为2πk /N 。可得离散傅里叶变换及其逆变换的定义为 ∑-=-=1 02)()(N n n N k j e n x k X π (1) ∑-==1 02)(1)(N k k N n j e k X N n x π (2) 如果把一个有限长序列看作是周期序列的一个周期,则离散傅里叶变换就是傅里叶级数。离散傅里叶变换也是周期的,周期为N 。 数字频率与模拟频率之间的关系为 s f f /2πω=,即s s T f f πωπω22== (3) 则第k 个频率点对应的模拟频率为 N kf NT k T N k f s s s k ==?=ππ212 (4) 在用快速傅里叶变换进行频谱分析时,要确定两个重要参数:采样率和频域采样点数,采样率可按奈奎斯特采样定理来确定,采样点数可根据序列长度或频率分辨率△f 来确定 f N f s ?≤,则f f N s ?≥ (5) 用快速傅里叶变换分析连续信号的频谱其步骤可总结如下: (1)根据信号的最高频率,按照采样定理的要求确定合适的采样频率f s ; (2)根据频谱分辨率的要求确定频域采样点数N ,如没有明确要求频率分辨率,则根据实际需要确定频率分辨率; (3)进行N 点的快速傅里叶变换,最好将纵坐标根据帕塞瓦尔关系式用功率来表示,
数字信号处理第2章习题解答 2.1 今对三个正弦信号1()cos(2)a x t t π=,2()cos(6)a x t t π=-,3()cos(10)a x t t π=进行理想采样,采样频率为8s πΩ=,求这三个序列输出序列,比较其结果。画出 1()a x t 、2()a x t 、3()a x t 的波形及采样点位置并解释频谱混淆现象。 解:采样周期为2184 T ππ= = 三个正弦信号采样得到的离散信号分别表示如下: 1()cos(2)cos()42a n x n n π π=?= 2()cos(6)cos()42a n x n n π π=-?=- 3()cos(10)cos()42 a n x n n π π=?= 输出序列只有一个角频率 2 π ,其中1()a x n 和3()a x n 采样序列完全相同,2()a x n 和1()a x n 、3()a x n 采样序列正好反相。 三个正弦信号波形及采样点位置图示如下: t x a 1(t )
t x a 2(t ) t x a 3(t ) 三个正弦信号的频率分别为1Hz 、3Hz 和5Hz ,而采样频率为4Hz ,采样频率大于第一个正弦信号频率的两倍,但是小于后两个正弦信号频率的两倍,因而由第一个信号的采样能够正确恢复模拟信号,而后两个信号的采样不能准确原始的模拟信号,产生频谱混叠现象。
2.3 给定一连续带限信号()a x t 其频谱当f B >时,()a X f 。求以下信号的最低采样频率。 (1)2()a x t (2)(2)a x t (3)()cos(7)a x t Bt π 解:设()a x t 的傅里叶变换为()a X j Ω (1)2 ()a x t 的傅里叶变换为 22()[()]B a a B X j X j d ππ ωωω-?Ω-? 因为22,22B B B B πωππωπ-≤≤-≤Ω-≤ 所以44B B ππ-≤Ω≤ 即2()a x t 带限于2B ,最低采样频率为4B 。 (2)(2)a x t 的傅里叶变换为 1 (/2)2 a X j Ω 2/22B B ππ-≤Ω≤,即44B B ππ-≤Ω≤ 即(2)a x t 带限于2B ,最低采样频率为4B 。 (3)()771 ()cos(7)()2 j Bt j Bt a a x t Bt x t e e πππ-= + 根据傅里叶变换的频移性质,()cos(7)a x t Bt π的傅里叶变换为 []1 ((7)((7)2 a a X j B X j B ππΩ-+Ω+ 它为一个带宽为2B 的带通信号,其通带范围为59 22 B f B ≤≤。 根据带通模拟信+号的采样定理,最小采样频率为1/4 4(1) 4.52 B B ?+=。 补充知识:带通模拟信号的采样定理 设带通模拟信号的频带限制在L f 和H f 之间,其频谱最低频率大于L f ,最高频率小于H f ,信号带宽H L B f f =-。此带通模拟信号所需最小抽样频率s f 等于 21s k f B n ??=+ ??? 式中,B 为信号带宽;n 为商( H f B )的整数部分,1,2,n = ;为商(H f B )的小数部分,01k <<。 2.5 一带通模拟信号如图所示,现用以下采样频率对其采样。 (1)25 Hz (2)50 Hz (3)100 Hz
实验报告 2010年 4月 15日 实验二用matlab实现傅立叶变换 实验目的: 1.掌握傅立叶数值实现方法(矩阵算法) 2.生成连续周期信号,掌握程序优化技巧 3.对于自定义函数参数有效性的检查 4.复习并巩固“信号与系统”相关知识内容,学习用matlab实现问题实验环境:运行于Matlab7.6环境
实验内容:本次实验参照《信号与系统》——“Matlab综合实验”55页课后练习 习题: 1.如图4.4所示锯齿波信号,分别去一个周期的抽样数据X1(t),0<=t<=1和五个周期的 数据X(t),0<=t<5,计算其傅立叶变换X1(w)和X(w),比较有和不同并解释原因。 编程如下: 方法1: %计算单位锯齿波和五个周期波形的傅立叶变换 %解法1:基本用循环实现数值的计算;对于5个周期锯齿波用内外循环来生成实现 T1=1; %一个周期锯齿波 N1=10000; t1=linspace(0,T1-T1/N1,N1)'; f1=0*t1; f1=1-2*t1; OMG=32*pi; %频率抽样区间 K1=100; %频率抽样点数 omg=linspace(-OMG/2,OMG/2-OMG/K1,K1)'; X1=0*omg; for k=1:K1 %求解五个周期函数的傅里叶变换系数for n=1:N1 X1(k)=X1(k)+T1/N1*f1(n)*exp(-j*omg(k)*t1(n)); end end fs1=0*t1; for n=1:N1 %通过傅里叶逆变换还原原函数 for k=1:K1 fs1(n)=fs1(n)+OMG/2/pi/K1*X1(k)*exp(j*omg(k)*t1(n)); end end
脉搏、语音及图像信号的傅里叶分析 一、实验目的 1、了解常用周期信号的傅里叶级数表示。 2、了解周期脉搏信号、语音信号及图像信号的傅里叶分析过程 3、理解体会傅里叶分析的理论及现实意义 二、实验原理 任意一个周期为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 3sin 5)35n h f t n t n h t t t ωπωωωπ∞==--=+++∑ 同理:对于如图2所示的三角波,函数表达式为: 4t (-t <)44()232(1) (t )44h T T f t t T T h T π?≤??=??-≤? 其傅里叶级数展开为:
1202100022281()(1)()sin(21)21 811(sin sin 3sin 5) 35n n h f t n t n h t t t ωπωωωπ∞-==---=-++∑ 图1 方波 图2 三角波 从以上各式可知,任何周期信号都可以表示为无限多次谐波的叠加,谐波次数越高,振幅越小,它对叠加波的贡献就越小,当小至一定程度时(谐波振幅小于基波振幅的5%),则高次的谐波就可以忽略而变成有限次数谐波的叠加,这对设计仪器电路是很有意义的。 三、实验仪器 脉搏语音实验仪器,数字信号发生器,示波器 四、实验内容 1.傅里叶级数的合成 (1)利用数字信号发生器产生频率分别为100Hz 、300Hz 、500Hz 的正弦信号,并使其位相相同,振幅比为:1:1/3: 1/5,将上述三个信号,分别通过加法器输入到傅里叶分析仪,观察并记录其波形。 (2)利用数字信号发生器产生方波,输入到傅里叶分析仪,并将其与上述合成后的信号相比较。两者有何差异?试分析引起的原因,应如何消除? (3)利用数字信号发生器产生频率分别为200Hz 、600Hz 、1000Hz 的正弦信号,振幅比为:1:1/32:1/52,并且保证其相位相差180°,然后通过加法器输入到傅里叶分析仪,观察并记录其波形,并与数字信号发生器产生的三角波相比较。 (4)利用傅里叶分析仪分别产生方波与三角波,进行傅里叶分析,记录各正弦波频率以及相对的幅度之间的关系,并与上述加法器输入信号相比较。 2.滤波与选频分析 对上述(4)傅里叶分析的频谱,分别选择低频段和高频段信号通过傅里叶反变换,观察它们图像并导出保存,试分析低通滤波和高通滤波图像的区别。 3.周期信号傅里叶分析的应用: (1)“脉搏信号”的傅里叶分析 1)用傅里叶分析仪软件中提供的“脉搏信号”模块和压电晶体测试自己脉搏波的信号,观察你的脉搏信号。
数字信号处理第2章习题解答 2.1 今对三个正弦信号1()cos(2)a x t t π=,2()cos(6)a x t t π=-,3()cos(10)a x t t π=进行理想采样,采样频率为8s πΩ=,求这三个序列输出序列,比较其结果。画出 1()a x t 、2()a x t 、3()a x t 的波形及采样点位置并解释频谱混淆现象。 解:采样周期为2184 T ππ= = 三个正弦信号采样得到的离散信号分别表示如下: 1()cos(2)cos()42a n x n n π π=?= 2()cos(6)cos()42a n x n n π π=-?=- 3()cos(10)cos()42 a n x n n π π=?= 输出序列只有一个角频率 2 π ,其中1()a x n 和3()a x n 采样序列完全相同,2()a x n 和1()a x n 、3()a x n 采样序列正好反相。 三个正弦信号波形及采样点位置图示如下: t x a 1(t )
t x a 2(t ) t x a 3(t ) 三个正弦信号的频率分别为1Hz 、3Hz 和5Hz ,而采样频率为4Hz ,采样频率大于第一个正弦信号频率的两倍,但是小于后两个正弦信号频率的两倍,因而由第一个信号的采样能够正确恢复模拟信号,而后两个信号的采样不能准确原始的模拟信号,产生频谱混叠现象。 2.3 给定一连续带限信号()a x t 其频谱当f B >时,()a X f 。求以下信号的最低采样频率。 (1)2()a x t (2)(2)a x t (3)()cos(7)a x t Bt π
第二章 连续时间傅里叶变换 1 周期信号的频谱分析——傅里叶级数FS (1) 狄义赫利条件:在同一个周期1T 内,间断点的个数有限;极大值和极小值 的数目有限;信号绝对可积∞ dt t f T 1 )(。 (2) 傅里叶级数:正交函数线性组合。 正交函数集可以是三角函数集}:sin ,cos ,1{11N n t n t n ∈ωω或复指数函数集 }:{1Z n e t jn ∈ω,函数周期为 T 1,角频率为1 1122T f π=π=ω。 (3) 任何满足狄义赫利条件周期函数都可展成傅里叶级数。 (4) 三角形式的FS : (i) 展开式:∑∞ =ω+ω+=1110)sin ()(n n n t n b t con a a t f (ii) 系数计算公式: (a)直流分量:?= 1 )(11 0T dt t f T a (b)n 次谐波余弦分量:N n tdt n t f T a T n ∈ω= ?,cos )(21 11 (c)n 次谐波的正弦分量:N n tdt n t f T b T n ∈ω= ?1 ,sin )(211 (iii) 系数n a 和n b 统称为三角形式的傅里叶级数系数,简称傅里叶系数。 (iv) 称11/1T f =为信号的基波、基频;1nf 为信号的n 次谐波。 (v) 合并同频率的正余弦项得: (a) ∑∞ =ψ+ω+=110)cos()(n n n t n c c t f (b) ∑∞ =θ+ω+=1 10)sin()(n n n t n d d t f n ψ和n θ分别对应合并后 n 次谐波的余弦项和正弦项的初相位。 (vi) 傅里叶系数之间的关系: (a) 000d c a == (b) n n n n n d c a θ=ψ=sin cos (c) n n n n n n d c b θ=ψ-=cos sin (d) 000a d c == (e) 2222n n n n b a d c +== (f) n n n a b arctg -=ψ (g) n n n b a arctg =θ (5) 复指数形式的FS :