实验一 采样率对信号频谱的影响

实验一 采样率对信号频谱的影响
实验一 采样率对信号频谱的影响

实验一利用DFT分析信号频谱

实验一利用DFT 分析信号频谱 一、 实验目的 1. 加深对DFT 原理的理解。 2. 应用DFT 分析信号的频谱。 3. 深刻理解利用DFT 分析信号频谱的原理,分析实现过程中出现的现象及解决方法。 二、 实验设备与环境 计算机、MATLAB^件环境。 三、 实验基础理论 1. DFT 与DTFT 的关系 方法二:实际在MATLAB 十算中,上述插值运算不见得是最好的办法。 由于DFT 是DTFT 的取 样值,其相邻两个频率样本点的间距为 —,所以如果我们增加数据的长度 N,使得到的 N DFT 谱线就更加精细,其包络就越接近 DTFT 的结果,这样就可以利用 DFT 计算DTFT 如果 没有更多的数据,可以通过补零来增加数据长度。 3、利用DFT 分析连续时间函数 利用DFT 分析连续时间函数是,主要有两个处理:①抽样,②截断 对连续时间信号x a (t) 一时间T 进行抽样,截取长度为 M 则 址 ML X a (N)「-x a (t)e4dt 二「x a (nT)e jnT n=0 再进行频域抽样可得 M 4 —j 竺 n 送,T' X a (nT)e N =TX M (k) NT n =0 因此,利用DFT 分析连续时间信号的步骤如下: (1 )、确定时间间隔,抽样得到离散时间序列 x(n). (2) 、选择合适的窗函数和合适长度 M 得到M 点离散序列x M DFT 实际上是 DTFT 在单位圆上以 的抽样,数学公式表示为: N-1 _j 空 k X(k) = X(z)| 耳八 x(n)e N z” N n=0 (2 — 1) 2、利用 DFT 求DTFT 方法一:利用下列公式: 2rk X(e j )二、X(k)( ) k=0 N k= 0,1,..N - 1 (2 — 2) Sn(N ,/2) Nsin(,/2) .N A e 2为内插函数 (2— 3) (2—4) X a (r 1)|

周期信号频谱的特点

周期信号频谱的特点 在结构施工测量中,按装修工程要求将装饰施工所需要的控制点、线及时弹在墙、板上,作为装饰工程施工的控制依据。 1.地面面层测量 在四周墙身与柱身上投测出100cm水平线,作为地面面层施工标高控制线。 根据每层结构施工轴线放出各分隔墙线及门窗洞口的位置线。 2.吊顶和屋面施工测量 以1000m线为依据,用钢尺量至吊顶设计标高,并在四周墙上弹出水平控制线。对于装饰物比较复杂的吊顶,应在顶板上弹出十字分格线,十字线应将顶板均匀分格,以此为依据向四周扩展等距方格网来控制装饰物的位置。 屋面测量首先要检查各方向流水实际坡度是否符合设计要求,并实测偏差,在屋面四周弹出水平控制线及各方向流水坡度控制线。 3.墙面装饰施工测量 内墙面装饰控制线,竖直线的精度不应低于1/3000,水平线精度每3m两端高差小于±1mm,同一条水平线的标高允许误差为±3mm。外墙面装饰用铅直线法在建筑物四周吊出铅直线以控制墙面竖直度、平整度及板块出墙面的位置。 4.电梯安装测量 在结构施工中,从电梯井底层开始,以结构施工控制线为准,及时测量电梯井净空尺寸,并测定电梯井中心控制线。 测设轨道中心位置,并确定铅垂线,并分别丈量铅垂线间距,其相互偏差(全高)不应超过1mm。 每层门套两边弹竖直线,并保证电梯门坎与门前地面水平度一致。 5. 玻璃幕墙的安装测量 结构完工后,安装玻璃幕墙时,用铅垂钢丝的测法来控制竖直龙骨的竖直度,幕墙分格轴线的测量放线应以主体结构的测量放线相配合,对其误差应在分段分块内控制、分配、消化,不使其积累。幕墙与主体连接的预埋件,应按设计要求埋设,其测量放线偏差高差不大于±3mm,埋件轴线左右与前后偏差不大于10mm。 精度要求 轴线竖向投测精度不低于1/10000。平面放线量距精度不低于1/8000,标高传递精度主楼、裙房分别不超过±15mm、±10mm。 仪器选用 该工程测量选用TOPCON电子全站仪一台,2"级经纬仪两台,DS3水准仪两台,50m钢卷尺两把。激光铅直仪一台。 每次放线前,均应仔细看图,弄清楚各个轴线之见的关系。放线时要有工长配合并检查工作。放线后,质检人员要及时对所放的轴线进行检查。重要部位要报请监理进行验线,合格后方可施工。 所有验线工作均要有检查记录。 对验线成果与放线成果之间的误差处理应符合《建筑工程施工测量规程》的规定: 1. 当验线成果与放线成果之差小于1/√2 倍的限差时,放线成果可评为优良; 2. 当验线成果与放线成果之差略小于或等于√2 限差时,对放线工作评为合格(可不必改正放线成果或取两者的平均值); 3. 当验线成果与放线成果之差超过√2 限差时,原则上不予验收,尤其是重要部位,

周期矩形信号的频谱分析

1.周期信号的频谱 周期信号在满足一定条件时,可以分解为无数三角信号或指数之和。这就是周期信号的傅里叶级数展开。在三角形式傅里叶级数中,各谐波分量的形式为()1cos n n A n t ω?+;在指数形式傅里叶级数中,分量的形式必定为1j n t n F e ω 与1-j -n t n F e ω 成对出现。为了把周期信号所具有的各 次谐波分量以及各谐波分量的特征(如模、相角等)形象地表示出来,通常直接画出各次谐波的组成情况,因而它属于信号的频域描述。 以周期矩形脉冲信号为lifenxi 周期信号频谱的特点。周期矩形信号在一个周期(-T/2,T/2)内的时域表达式为 ,2 0,>2 ()A t T t f t ττ ≤?=?? (2-6) 其傅里叶复数系数为 12 n n A F Sa T ωττ?? = ??? (2-7) 由于傅里叶复系数为实数,因而各谐波分量的相位为零(n F 为正)或为π±(n F 为负),因此不需要分别画出幅度频谱n F 与相位频谱n φ。可以直接画出傅里叶系数n F 的分布图。 如图2.4.1所示。该图显示了周期性矩形脉冲信号()T f t 频谱的一些性质,实际上那个也是周期性信号频谱的普遍特性: ① 离散状频谱。即谱线只画出现在1ω的整数倍频率上,两条谱线的间隔为1ω(等于2π/t )。 ② 谱线宽度的包络线按采样函数()1/2a S n ωτ的规律变化。如图2.4.2所示。但1ω 为 2π τ 时,即( )2m π ωτ =(m=1,2,……)时,包络线经过零点。在两相邻 零点之间,包络线有极值点,极值的大小分别为-0.212()2A T τ,

信号与测试实验1时率与频率

基本信号分析 一、实验目的 1.掌握基本信号的时域和频域分析方法 2.掌握信号的自相关和互相关分析,了解其应用 二、数据处理与分析 (1)幅值为1,频率为100Hz的正弦信号,上图为时域图,下图为利用快速傅里叶变换获得的频谱图。从频谱图上看出,f=100Hz时频域的幅值最大。 (2)频域为100Hz,幅值为1的方波信号,上图为时域图,下图为借助快速傅立叶变换获得的频域图。从频谱图上看出,f=100Hz时频域的幅值最大,随着频域增大,频域的幅值逐渐衰减。

(3)频率为100Hz,幅值为1的锯齿波信号图,上图为时域图,下图为借助傅立叶变换而获得的频域图。从频域图看出,在100Hz的整数倍频率上,频域幅值都出现了峰值,随着频率的增大,峰值逐渐收敛至0. (4)平均振幅为1的噪声信号,上图为时域图,下图为通过快速傅立叶变

换得出的频谱图,从频谱图可以看出,白噪声信号的频谱杂乱无章,无明显规律。 (5)由频率为50Hz、100Hz、150Hz的正弦信号组成的复合信号,上图为时域图,下图为频域图,从图中可以看出,频谱图在50、100、150Hz处出现了峰值。 (6)频率为100Hz 的正弦信号叠加噪声信号:上图为时域信号图,下图为

通过快速傅立叶变换获得的频谱图。与没有叠加噪声信号的正弦波相比,时域波形出现了毛刺,而频谱图中除了在100Hz处有峰值外,在其他频率点处也出现了一些较低的峰值。 (7)频率为100Hz的正弦信号和频率为100Hz的方波信号进行叠加,上图为时域信号,下图为频谱图。从时域图上可以看出,正弦波形叠加方波后有了明显的畸变。从频谱图上可以看出,除了100Hz处出现峰值以外,在其他频率点也出现了一些峰值。

实验:典型信号频谱分析报告

实验3.2 典型信号频谱分析 一、 实验目的 1. 在理论学习的基础上,通过本实验熟悉典型信号的波形和频谱特征,并 能够从信号频谱中读取所需的信息。 2. 了解信号频谱分析的基本方法及仪器设备。 二、 实验原理 1. 典型信号及其频谱分析的作用 正弦波、方波、三角波和白噪声信号是实际工程测试中常见的典型信号,这些信号时域、频域之间的关系很明确,并且都具有一定的特性,通过对这些典型信号的频谱进行分析,对掌握信号的特性,熟悉信号的分析方法大有益处,并且这些典型信号也可以作为实际工程信号分析时的参照资料。本次实验利用DRVI 快速可重组虚拟仪器平台可以很方便的对上述典型信号作频谱分析。 2. 频谱分析的方法及设备 信号的频谱可分为幅值谱、相位谱、功率谱、对数谱等等。对信号作频谱分析的设备主要是频谱分析仪,它把信号按数学关系作为频率的函数显示出来,其工作方式有模拟式和数字式二种。模拟式频谱分析仪以模拟滤波器为基础,从信号中选出各个频率成分的量值;数字式频谱分析仪以数字滤波器或快速傅立叶变换为基础,实现信号的时—频关系转换分析。 傅立叶变换是信号频谱分析中常用的一个工具,它把一些复杂的信号分解为无穷多个相互之间具有一定关系的正弦信号之和,并通过对各个正弦信号的研究来了解复杂信号的频率成分和幅值。 信号频谱分析是采用傅立叶变换将时域信号x(t)变换为频域信号X(f),从而帮助人们从另一个角度来了解信号的特征。时域信号x(t)的傅氏变换为: 式中X(f)为信号的频域表示,x(t)为信号的时域表示,f 为频率。 3. 周期信号的频谱分析 周期信号是经过一定时间可以重复出现的信号,满足条件: dt e t x f X ft j ?+∞ ∞--=π2)()(

习题1 绘制典型信号及其频谱图(参考模板)

习题一绘制典型信号及其频谱图 电子工程学院 202班一、单边指数信号 单边指数信号的理论表达式为 对提供的MATLAB程序作了一些说明性的补充,MATLAB程序为

figure(3); plot(w,20*log10(abs(F)));xlabel('\omega');ylabel('|F(\omega)| in dB');title(' 幅频特性/dB'); figure(4); plot(w,angle(F)*57.29577951);xlabel('\omega');ylabel('\phi(\omega)/(°) ');title('相频特性'); 调整,将a分别等于1、5、10等值,观察时域波形和频域波形。由于波形 较多,现不失代表性地将a=1和a=5时的各个波形图列表如下进行对比,其 他a值的情况类似可推知。 a15 时 域 图 像

幅频特性 幅频特性/d B 相频特性

分析: 由上表中a=1和a=5的单边指数信号的波形图和频谱图的对比可以发现,当a值增大时,信号的时域波形减小得很快,而其幅频特性的尖峰变宽,相频特性的曲线趋向平缓。 二、矩形脉冲信号 矩形脉冲信号的理论表达式为 MATLAB程序为:

clear all; E=1;%矩形脉冲幅度 width=2;%对应了时域表达式中的tao t=-4:0.01:4; w=-5:0.01:5; f=E*rectpuls(t,width); %MATLAB中的矩形脉冲函数,width即是tao,t为时间 F=E*width*sinc(w.*width/2); figure(1); plot(t,f);xlabel('t');ylabel('f(t)');title('信号时域图像'); figure(2); plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|');title('幅频特性'); figure(3); plot(w,20*log10(abs(F)));xlabel('\omega');ylabel('|F(\omega)| in dB');title(' 幅频特性/dB'); figure(4); plot(w,angle(F));xlabel('\omega');ylabel('\phi(\omega)');title('相频特性'); 调整,将分别等于1、4等值,观察时域波形和频域波形。由于波形较多,现不失代表性地将a=1和a=4时的各个波形图列表如下进行对比,其他值的情况类似可推知。 14

信号与系统实验报告-实验3--周期信号的频谱分析

信号与系统实验报告-实验3--周期信号的频谱分析

信号与系统 实验报告 实验三周期信号的频谱分析 实验三周期信号的频谱分析 实验目的: 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法; 2、观察截短傅里叶级数而产生的“Gibbs现象”,了解其特点以及产生的原因;

3、掌握各种典型的连续时间非周期信号的频谱特征。 实验内容: (1)Q3-1 编写程序Q3_1,绘制下面的信号的波形图: 其中,0 = 0.5π,要求将一个图形窗口分割成四个子图,分别绘制cos(0t)、cos(30t)、cos(50t) 和x(t) 的波形图,给图形加title,网格线和x坐标标签,并且程序能够接受从键盘输入的和式中的项数。 程序如下: clear,%Clear all variables close all,%Close all figure windows dt = 0.00001; %Specify the step of time variable t = -2:dt:4; %Specify the interval of time w0=0.5*pi; x1=cos(w0.*t); x2=cos(3*w0.*t); x3=cos(5*w0.*t); N=input('Type in the number of the harmonic components N='); x=0; for q=1:N; x=x+(sin(q*(pi/2)).*cos(q*w0*t))/q; end subplot(221) plot(t,x1)%Plot x1 axis([-2 4 -2 2]); grid on, title('signal cos(w0.*t)') subplot(222) plot(t,x2)%Plot x2 axis([-2 4 -2 2]); grid on, title('signal cos(3*w0.*t))') subplot(223) plot(t,x3)%Plot x3

周期信号的频谱分析

信号与系统 实验报告 实验三周期信号的频谱分析 实验报告评分:_______ 实验三周期信号的频谱分析 实验目的: 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法; 2、观察截短傅里叶级数而产生的“Gibbs现象”,了解其特点以及产生的原因;

3、掌握各种典型的连续时间非周期信号的频谱特征。 实验内容: (1)Q3-1 编写程序Q3_1,绘制下面的信号的波形图: 其中,0 = 0.5π,要求将一个图形窗口分割成四个子图,分别绘制cos( 0t)、cos(3 0t)、cos(5 0t)和x(t) 的波形图,给图形加title,网格线和x坐标标签,并且程序能够接受从键盘输入的和式中的项数。 程序如下: clear,%Clear all variables close all,%Close all figure windows dt = 0.00001; %Specify the step of time variable t = -2:dt:4; %Specify the interval of time w0=0.5*pi; x1=cos(w0.*t); x2=cos(3*w0.*t); x3=cos(5*w0.*t); N=input('Type in the number of the harmonic components N='); x=0; for q=1:N; x=x+(sin(q*(pi/2)).*cos(q*w0*t))/q; end subplot(221) plot(t,x1)%Plot x1 axis([-2 4 -2 2]); grid on, title('signal cos(w0.*t)') subplot(222) plot(t,x2)%Plot x2 axis([-2 4 -2 2]); grid on, title('signal cos(3*w0.*t))') subplot(223) plot(t,x3)%Plot x3 axis([-2 4 -2 2])

信号与检验测试实验一

实验一、基本信号分析 一、实验目的 1. 掌握基本信号的时域和频域分析方法 2. 掌握信号的自相关和互相关分析,了解其应用 二、实验原理 (1)信号的时域和频域转换 目的:研究分析信号的时域特征(如持续时间、幅值、周期等)和信号的频域特征(如是否含有周期性信号、信号的频率带宽等) 转换方法:时域有限长序列 频域有限长序列: 离散傅里叶变换 (2)信号相关性 相关是用来描述一个随机过程自身在不同时刻的状态间,或者两个随机过程在某个时刻状态间线性依从关系的数字特征。 自相关函数定义为: xx 01()lim ()()T T R x t x t dt T ττ→∞ =+? 互相关函数定义为: xx 0 1()lim ()()T T R x t x t dt T ττ→∞=+?

三、实验内容与步骤 (1)产生不同的周期信号,包括正弦信号、方波信号、锯齿波信号,在时域分析这些波形特征(幅值、频率(周期))。 上图为幅值为2频率为20Hz的正弦信号时域图,下图为快速傅里叶变换之后获得的频谱图。从频谱图上看出,f=20Hz时频域的幅值最大,和时域图吻合。

上图为幅值为3频率为5Hz的方波信号时域图,下图为快速傅里叶变换之后获得的频谱图。从频谱图上看出,方波信号傅里叶分解后由一个频率为5Hz 的基波和无数个高次谐波组成。以幅值衰减十倍为带宽,由图可知此方波信号带宽约为35Hz

上图为幅值为4频率为10Hz的三角波信号时域图,下图为快速傅里叶变换之后获得的频谱图。从频域图看出,在10Hz的整数倍频率上,频域幅值出现了峰值,其后有无数个谐波和基波一起组成了三角波。以幅值衰减十倍为带宽,由图可知此三角波信号带宽约为80Hz (2)在Matlab中产生随机噪声、阶跃信号(选作)、矩形脉冲(选作)

09典型信号的频谱分析

实验九 典型信号的频谱分析 一. 实验目的 1. 在理论学习的基础上,通过本实验熟悉典型信号的频谱特征,并能够从信号频谱中读取 所需的信息。 2. 了解信号频谱分析的基本原理和方法,掌握用频谱分析提取测量信号特征的方法。 二. 实验原理 信号频谱分析是采用傅里叶变换将时域信号x(t)变换为频域信号X(f),从而帮助人们从另一个角度来了解信号的特征。 图1、时域分析与频域分析的关系 信号频谱X(f)代表了信号在不同频率分量成分的大小,能够提供比时域信号波形更直观,丰富的信息。时域信号x(t)的傅氏变换为: dt e t x f X ft j ?+∞ ∞--=π2)()( (1) 式中X(f)为信号的频域表示,x(t)为信号的时域表示,f 为频率。 工程上习惯将计算结果用图形方式表示, 以频率f 为横坐标,X(f)的实部)(f a 和虚部 )(f b 为纵坐标画图,称为时频-虚频谱图; 以频率f 为横坐标,X(f)的幅值)(f A 和相位 )(f ?为纵坐标画图,则称为幅值-相位谱; 以f 为横坐标,A(f) 2为纵坐标画图,则称为 功率谱,如图所示。 频谱是构成信号的各频率分量的集合,它 完整地表示了信号的频率结构,即信号由哪些 谐波组成,各谐波分量的幅值大小及初始相 位,揭示了信号的频率信息。 图2、信号的频谱表示方法

三. 实验内容 1. 白噪声信号幅值谱特性 2. 正弦波信号幅值谱特性 3. 方波信号幅值谱特性 4. 三角波信号幅值谱特性 5. 正弦波信号+白噪声信号幅值谱特性 四. 实验仪器和设备 1. 计算机1台 2. DRVI快速可重组虚拟仪器平台1套 3. 打印机1台 五. 实验步骤 1.运行DRVI主程序,点击DRVI快捷工具条上的"联机注册"图标,选择其中的“DRVI 采集仪主卡检测”或“网络在线注册”进行软件注册。 2.在DRVI软件平台的地址信息栏中输入WEB版实验指导书的地址,在实验目录中选择 “典型信号频谱分析”,建立实验环境。 图5 典型信号的频谱分析实验环境 下面是该实验的装配图和信号流图,图中的线上的数字为连接软件芯片的软件总线数据线号,6017、6018为两个被驱动的信号发生器的名字。 图6 典型信号的频谱分析实验装配图

用FFT对信号作频谱分析

实验三:用FFT 对信号作频谱分析 一、实验原理与方法 1、用FFT 对信号作频分析是学习数字信号处理的重要内容,经常需要进行分析的信号是模拟信号的时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D 和分析误差。频谱分辨率直接和FFT 的变换区间N 有关,因为FFT 能够实现的频率分辨率是N π2,因此要求D N ≤π2。可以根据此式选择FFT 的变换区间N 。误差主要来自于用FFT 作频谱分析时,得到的是离散谱,而信号(周期信号除外)是连续谱,只有当N 较大时,离散谱的包络才能逼近连续谱,因此N 要适当选择大一些。 2、周期信号的频谱是离散谱,只有用整数倍周期的长度作FFT ,得到的离散谱才能代表周期信号的频谱。如果不知道信号周期,可以尽量选择信号的观察时间长一些。 3、对模拟信号进行谱分析时,首先要按照采样定理将其变成时域离散信号。如果是模拟周期信号,也应该选取整数倍周期长度,经过采样后形成周期序列,按照周期序列的谱分析进行。 二、实验内容 1、对以下序列进行FFT 谱分析: )()(41n R n x = ?????≤≤-≤≤+=n n n n n n x 其他0 7483 01 )(2 ?????≤≤-≤≤-=n n n n n n x 其他0 7433 04)(3 选择FFT 的变换区间N 为8和16两种情况进行频谱分析,分别打印出幅频特性曲线,并进行讨论、分析。程序见附录3.1、实验结果见图3.1。 2、对以下周期序列进行谱分析: n n x 4cos )(4π = n n n x 8cos 4cos )(5π π+= 选择FFT 的变换区间N 为8和16两种情况进行频谱分析,分别打印出幅频特性曲线,并进行讨论、分析与比较。程序见附录3.2、实验结果见图3.2。 3、对模拟周期信号进行频谱分析: t t t t x πππ20cos 16cos 8cos )(6++= 选择采样频率Fs=64Hz ,FFT 的变换区间N 为16、32、64三种情况进行频谱分析,分别打印出幅频特性曲线,并进行讨论、分析与比较。程序见附录3.3、实验结果见图3.3。

实验一离散信号的频谱分析报告

实验一离散信号的频谱分析报告 班级 姓名 学号

实验一离散信号的频谱分析报告 1 掌握采样频率的概念 2 掌握信号频谱分析方法; 3 掌握在计算机中绘制信号频谱图的方法。 ①采样频率为1000Hz,信号频率为30Hz的正弦信号y1(n) 对其进行FFT变换 ②采样频率为1000Hz,信号频率为120Hz的正弦信号y2(n)

对其进行FFT变换 ③采样频率为1000Hz, 30Hz的正弦信号和120Hz的混合信号y3(n)。 对其进行FFT变换

语音信号波形

附录程序: fs=1000;%设定采样频率 N=1024; n=0:N-1; t=n/fs; f0=30;%设定正弦信号频率 %生成正弦信号 x=sin(2*pi*f0*t); figure(1); subplot(3,2,1); plot(t,x);%作正弦信号的时域波形xlabel('t'); ylabel('y'); title('正弦信号30HZ时域波形'); grid; %进行FFT变换并做频谱图

y=fft(x,N);%进行fft变换 mag=abs(y);%求幅值 f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换figure(1); subplot(3,2,2); plot(f,mag);%做频谱图 axis([0,100,0,500]); xlabel('频率(Hz)'); ylabel('幅值'); title('正弦信号30HZ幅频谱图N=1024'); grid; %120HZ f1=120; x=sin(2*pi*f1*t); figure(1); subplot(3,2,3); plot(t,x);%作正弦信号的时域波形 xlabel('t'); ylabel('y'); title('正弦信号120HZ时域波形'); grid; %进行FFT变换并做频谱图 y=fft(x,N);%进行fft变换 mag=abs(y);%求幅值 f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换

第四章 周期信号的频谱分析

第四章 周期信号的频域分析 1. 内容提要 本章介绍连续周期信号的傅立叶级数及其基本性质;连续周期信号频谱的概念,相位谱的作用。对离散周期信号傅立叶级数和其基本性质做简单了解。 2. 学习目标 通过本章的学习,应达到以下要求: (1)掌握周期信号频谱的概念及信号频带宽度的概念。 (2)熟悉傅里叶变换的主要性质。 (3)熟悉频域分析法。 (4)了解离散傅立叶级数的概念 3. 重点难点 (1) 信号的对称性和傅立叶系数的关系 (2) 连续信号的频谱分析,包括周期信号频谱的概念,相位谱和功率谱。 4. 应用 周期信号频域分析的MATLAB 实现 5. 教案内容 4.1 连续时间信号的傅立叶变换 周期信号的定义 周期信号是定义在001/f T =(,)-∞∞区间,每隔一定的时间间隔0T ,按相同规律重复变化的信号。即对t R ?∈,存在一个大于零的0T ,使得 0()(),f t T f t t R +=?∈ 其中0T 为基波周期,002/T ωπ=为基波角频率,001/f T =为基波频率

傅立叶级数的实质 就是将复杂信号分解成为更容易处理的信号形式。 4.1.1 指数形式的傅里叶级数 连续时间信号的傅立叶级数表示为 0()jnw t n n f t C e ∞ =-∞ = ∑ 称n C 为周期信号()f t 的傅立叶系数。傅立叶系数的计算公式为 00 00 1 ()t T jn t t Cn f t e dt T ω+-= ? 4.1.2 三角形式的傅立叶级数 若函数()f t 满足狄里赫利条件,周期信号f(t) 展开成傅里叶级数。 01111212111()cos sin cos 2sin 2cos sin n n f t a a t b t a t b t a n t b n t ωωωωωω=++++++++ 0111 (cos sin )n n n a a n t b n t ωω∞ ==++∑ 式中,n 为正整数;系数0,,n n a a b 称为傅里叶系数,考虑到三角函数集是一组完备的正交函数集,因此,可得一个周期1(0,)T 的傅里叶系数: 1 11200112 11()()T T T a f t dt f t dt T T -==?? 1 10 12()cos T n a f t n tdt T ω=? 1 10 12()sin T n b f t n tdt T ω=?

实验1 信号的频谱图

大连理工大学实验报告 学院(系): 专业: 班级: _________ 姓 名: 学号: 组: ___ __________ 实验时间: 实验室: 实验台:_____________ 指导教师签字: 成绩:__________________ 实验名称 一、 实验目的和要求 1. 掌握周期信号的傅里叶级数展开; 2. 掌握周期信号的有限项傅里叶级数逼近; 3. 掌握周期信号的频谱分析; 4. 掌握连续非周期信号的傅立叶变换; 5. 掌握傅立叶变换的性质。 二、实验程序和结果 1. 已知周期三角信号如下图1-5所示,试求出该信号的傅里叶级数,利用MATLAB 编程实现其各次谐波的叠加,并验证其收敛性。 三角型号的傅里叶级数展开式如下: ) 5cos 1 3cos 1 (cos 4 2 1 )(5 3 2 2 2 ++ + += wt wt wt t f pi

clc clear t=-5:0.001:5; y=0.5*(sawtooth(pi*(t+1),0.5)+1); plot(t,y); xlabel('t'),ylabel('周期三角波') axis([-3 3 -0.5 1.5]); grid on n_max=[1 3 5 11 47]; N=length(n_max); for k=1:N n=1:2:n_max(k); b=4./(pi*pi*n.*n); x=0.5+b*cos(pi*n'*t); figure; plot(t,y,'b'); hold on; plot(t,x,'r'); hold off; xlabel('t'),ylabel(‘部分和波形'); axis([-3 3 -0.5 1.5]);grid on; title(['最高谐波次数=',num2str(n_max(k))]); end t 部分和波形 最高谐波次数 =1 t 周期三角信号

周期方波信号的频谱具有三个特点1

1.测试系统的组成:传感器+中间变换装置+显示记录装置 传感器:反映被测对象特性的物理量(如噪声、温度)检出并转换为电量; 中间变换装置:对接收到的电信号用硬件电路进行分析处理或经A/D 变换后用软件进行计算; 显示记录装置:将测量结果显示出来,提供给观察者或其它自动控制装置 2.周期方波信号的频谱具有三个特点:○ 1离散性,频谱是非周期性离散的线状频谱,成为谱线,连接个谱线顶点的曲线为频谱的包络线,它反映了各频率分量的幅度随频率的变化情况。 ○2谐波性 普线以基波频率0 ω为间隔等距离分布,任意两谐频之比都是整数或整数比,即为有理数。各次谐波的频率都是基频0ω的整数倍,相邻频率的间隔为0ω或它的整数倍。 ○ 3收敛性 周期信号的幅值频谱是收敛的。即谐波的频率越高,其幅值越小,再整个信号中所占的比重也就越小。 傅立叶变换的性质:○ 1线性叠加性○2尺度展缩性○3对称性○4时移性质○5频移性质 采样定理:信号)(t x 的傅立叶变换为)(ωX ,其频率范围为m m ωω~-,当m s ωω2 ,频谱发生混叠。采样频率s ω的选择对正确的采样是至关重要的。如果m s ωω2≥则不会发生频混关系,因此,对采样脉冲的间隔S T 须加以限制,即采样频率()s s T /2πω或()s s T f /1必须大于或等于()t x 中的最高频率m ω的两倍,这就是采样定理,其表达式为m s ωω2≥或 m s f f 2≥ 实际采样频率一般选得大于2m ω. 测试系统的静态特性 不是真测试的条件:○1系统的幅频特性在输入信号()t x 的频谱范围内为常数;○ 2系统的相频特性()ω?是过原点且具有负斜率的直线。 传感器的分类:○1按输入量分类(用它所测量的物理量来分类):测力传感器、位移传感器、 温度传感器;○ 2按其输出量分类:电路参数型传感器、发电型传感器。 参数型传感器的工作原理:将被测物理量转换为电路参数的传感器,主要有电阻式、电容式、电感式三种。 电阻式传感器是把被测量转化为电阻变化的传感器。 电阻式传感器按其工作原理可分为变阻器式和电阻应变式两类。变阻式传感器通过改变电位器触头位置实现位移到电阻的转换。 电阻应变片的工作原理基于"力→应变→电阻变化“三个基本转换环节。 半导体电阻材料应变片的工作原理主要是利用半导体材料的电阻率随应力变化,这一现象常称为压阻效应。 电容传感器是将各种被测物理量转换为电容量变化的装置。

矩形脉冲信号频谱分析

矩形脉冲信号频谱分析

小组成员: 刘鑫 龙宇 秦元成 王帅 薛冬寒 梁琼健 一、傅里叶分析方法与过程 周期信号的分解 1、三角形式 周期为T 的周期信号,满足狄里赫利(Dirichlet )条件(实际中遇到的所有周期信号都符合该条件),便可以展开为傅里叶级数的三角形式,即: ∑∑∞ =∞ =Ω+Ω+=110sin cos 21 )(n n n n t n b t n a a t f (1) ?-=Ω=2 2 ,2,1cos )(2T T n dt t n t f T a n Λ (2)

?-=Ω=2 2 ,2,1sin )(2T T n dt t n t f T b n Λ (3) 式中: T π2= Ω 为基波频率,n a 与 n b 为傅 里叶系数。 其中 n a 为n 的偶函数, n b 为n 的奇函数。 将上式中同频率项合并可写成: ∑∞ =+Ω+=++Ω++Ω+=1022110)cos 21 ... )2cos()cos(21 )(n n n t n A A t A t A A t f ???( 式中: ) arctan(... 3,2,1,2 2 0n n n n a b n b a A a A n n -==+==? (5)

n n n n n n A b A a A a ??sin cos 0 0-=== (6) 2.指数形式 由于 2 cos jx jx e e x -+= (7) 三角函数形式可以写为 t jn j n n t jn j n n t n j n t n j n e e A e e A A e e A A t f n n n n Ω--∞=Ω∞=+Ω-∞ =+Ω∑∑∑++=++=????1 10)(1)(0212121] [2 1 21)( (8) 将上式第三项中的n 用-n 代换,并考虑到 为n 的偶函数, 为n 的奇函数 则上式可写为: t jn j n n t jn j n n t jn j n n t jn j n n e e A e e A A e e A e e A A t f n n n n Ω∞ --=Ω∞=Ω--∞-=-Ω∞=∑∑∑∑++=++=-????1 101 1021 2121212121)( (9)

周期信号的频谱测试

实验一 周期信号的频谱测试 一、实验目的: 1、掌握周期信号频谱的测试方法; 2、了解典型信号频谱的特点,建立典型信号的波形与频谱之间的关系。 二、实验原理及方法: 1、信号的频谱可分为幅度谱、相位谱和功率谱,分别是 将信号的基波和各次谐波的振幅、相位和功率按频率的高低依次排列而成的图形。 2、周期连续时间信号的频谱具有离散性、谐波性、收敛性三个特点。例如正弦波、周期矩形脉冲、三角波的幅度谱分别如图1-1,1-2,1-3所示: 1 2 3 4 5 6 7 -1-0.8-0.6-0.4-0.200.2 0.40.60.81t s i n (t ) n C 1 ωω 图1-1(a) 正弦波信号 图1-1(b) 相应的幅度谱 f(t) T A 0τ/2 n C 14ω1 5ω13ω12ω1ωω 图1-2(a) 周期矩形脉冲 图1-2(b) 相应的幅度谱 因此,信号的频谱测试方法可用频谱分析仪直接测量亦可用逐点选频测量法进行测量。本实验使用GDS-806C 型号的数字存储示波器直接测试幅度谱。

用示波器直接测试,就是将其与EE1460C 函数信号发生器连好。分别输入相应频率(重复频率)和幅度的正弦波,三角波和矩形波,此时示波器将显示按频率由低到高的各输入信号的谐波分量。GDS-806C 数字存储示波器测频谱的方法,就是将MATH 键按下,F1键选择FFT(快速傅立叶转换)功能可以将一个时域信号转换成频率构成,显示器出现一条红颜色的频谱扫描线。当示波器输入了不同信号的波形时就显示它们相应的频谱, 参数的测量由调试水平(即频率)与垂直(即增益)游标获取,从而得到输入信号的频谱图。 三、实验原理图: 图1-4 实验原理图 四、实验内容及步骤: 1、测试正弦波的幅度频谱 将信号源、示波器、按图1-4连接好;信号源CH1的输出波形调为正弦波,输出频率自选,输出信号幅度自选 ,并记录幅度与频率的参数.测出前五次谐波分量.将其数据填入表一。 表一:正弦波前五次谐波的幅度谱 f/Hz 1K 2K 3K 4K 5K Cn/dB 55.2 2、测试三角波的幅度频谱 在实验步骤1的基础上将信号源CH1的输出波形调为三角波(T) ,频率自选,幅度自选.并记录幅度和周期的参数.测出前五次谐波分量。将测量数据填入表二。 表二:三角波前五次谐波的幅度谱 f/Hz 1k 3k 5k 7k 9k Cn/dB 52.8 33.6 24.8 19.2 14.4 t f(t)T1 A T1/2 n C 14ω15ω13ω12ω1ω16ω1 7ω 图1-3(a) 三角波 1-3(b) 相应的幅度谱

周期信号的频谱的特点

周期信号的频谱的特点 一、 周期信号的频谱 一个周期信号)(t f ,只要满足狄里赫利条件,则可分解为一系列谐波分量之和。其各次谐波分量可以是正弦函数或余弦函数,也可以是指数函数。不同的周期信号,其展开式组成情况也不尽相同。在实际工作中,为了表征不同信号的谐波组成情况,时常画出周期信号各次谐波的分布图形,这种图形称为信号的频谱,它是信号频域表示的一种方式。 描述各次谐波振幅与频率关系的图形称为振幅频谱,描述各次谐波相位与频率关系的图形称为相位频谱。根据周期信号展成傅里叶级数的不同形式又分为单边频谱和双边频谱。 1单边频谱 若周期信号)(t f 的傅里叶级数展开式为式(3-15),即 ∑ ∞ =+Ω+=10)cos()(n n n t n A A t f ? (3-24) 则对应的振幅频谱n A 和相位频谱n ?称为单边频谱。 例3-3 求图3-4所示周期矩形信号)(t f 的单边频谱图。 解 由)(t f 波形可知, )(t f 为偶函数,其傅里叶系数

?==2/0021)(4T dt t f T a ?=Ω=2/0)4/sin(2cos )(4T n n n tdt n t f T a ππ 0=n b 故 ∑∑∞=∞=Ω+=Ω+=110cos )4/sin(241cos 2)(n n n t n n n t n a a t f ππ 因此 410=A , ππn n A n )4/sin(2= 即 45.01=A , 32.02≈A , 15.03≈A , 04=A , 09.05≈A , 106.06≈A ┅ 单边振幅频谱如图3-5所示。 t f(t) 图 3 - 4τ τττ4 2/ 0 2/ 4--1图 3 - 50.25 0.45 0.320.150.090.106ΩΩΩΩΩΩΩ7 6 5 4 3 2 0A n 2双边频谱

第四章傅立叶变换习题

第三章傅立叶变换 第一题选择题 1.连续周期信号f (t )的频谱F(w)的特点是 D 。 A 周期连续频谱 B 周期离散频谱 C 非周期连续频谱 D 非周期离散频谱 2.满足抽样定理条件下,抽样信号f s (t)的频谱)(ωj F s 的特点是 (1) (1)周期、连续频谱; (2)周期、离散频谱; (3)连续、非周期频谱; (4)离散、非周期频谱。 3.信号的频谱是周期的连续谱,则该信号在时域中为 D 。 A 连续的周期信号 B 离散的周期信号 C 连续的非周期信号 D 离散的非周期信号 4.信号的频谱是周期的离散谱,则原时间信号为 (2) 。 (1)连续的周期信号 (2)离散的周期信号 (3)连续的非周期信号 (4)离散的非周期信号 5.已知f (t )的频带宽度为Δω,则f (2t -4)的频带宽度为( 1 ) (1)2Δω (2)ω?2 1 (3)2(Δω-4) (4)2(Δω-2) 6.若=)(1ωj F F =)()],([21ωj F t f 则F =-)]24([1t f ( 4 ) (1)ωω41)(21j e j F - (2)ωω41)2 (21j e j F -- (3)ωωj e j F --)(1 (4)ωω21)2 (21j e j F -- 7.信号f (t )=Sa (100t ),其最低取样频率f s 为( 1 ) (1)π100 (2)π 200 (3)100π (4)200 π 8.某周期奇函数,其傅立叶级数中 B 。 A 不含正弦分量 B 不含余弦分量 C 仅有奇次谐波分量 D 仅有偶次谐波分量 9.某周期偶谐函数,其傅立叶级数中 C 。 A 无正弦分量 B 无余弦分量 C 无奇次谐波分量 D 无偶次谐波分量 10.某周期奇谐函数,其傅立叶级数中 C 。 A 无正弦分量 B 无余弦分量 C 仅有基波和奇次谐波分量 D 仅有基波和偶次谐波分量 11.某周期偶函数f(t),其傅立叶级数中 A 。

周期信号频谱分析

实验名称:周期信号的频谱分析 教材名称:电工电子实验技术(下册)页码:P142 实验目的: 1、了解和掌握周期信号频谱分析的基本概念; 2、掌握Multisim软件用于频谱分析的基本方法; 3、加深理解周期信号时域参数变化对其谐波分量的影响及变化趋势。 实验任务: 1、根据9-1给定的波形和参数测量各谐波分量的幅度值。 2、根据所测数据绘制每一波形的谱线图。 设计提示: 实验电路图: 图一、分析用电路及信号发生器调整窗口 实验结果: 表9-1数据: 周期信号的频谱分析(Multisim) 0 10 20 30 40 50 60 70 80 90 100 矩形波10%-4.023 1.923 1.833 1.689 1.499 1.273 1.024 0.763 0.506 0.263 0.047 矩形波30%-2.023 5.123 3.040 0.699 0.897 1.271 0.659 0.236 0.739 0.595 0.046 矩形波50%-0.022 6.366 0.045 2.121 0.045 1.271 0.045 0.906 0.045 0.703 0.045 正弦波0 4.999 0 0 0 0 0 0 0 0 0 三角波50%0 4.053 0 0.451 0 0.162 0 0.083 0 0.050 0 三角波70%0 3.903 1.147 0.166 0.177 0.193 0.079 0.030 0.072 0.048 0 三角波90%0 3.479 1.654 1.012 0.669 0.450 0.298 0.186 0.103 0.043 0 N 0 1 2 3 4 5 6 7 8 9 10 注:谱线数取10+直流。

周期信号的频谱分析

信号与系统 实验三周期信号的频谱分析 实验报告评分:______ 实验三周期信号的频谱分析 实验目的: 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法; 2、观察截短傅里叶级数而产生的“Gibbs 现象”,了解其特点以及产生的原因; 3、掌握各种典型的连续时间非周期信号的频谱特征。 实验内容: (1)Q3-1 编写程序 Q3_1,绘制下面的信号的波形图:

其中,0 = 0.5π,要求将一个图形窗口分割成四个子图,分别绘制cos( 0t)、cos(3 0t)、cos(5 0t) 和 x(t) 的波形图,给图形加title,网格线和 x 坐标标签,并且程序能够接受从键盘输入的和式中的项数。 程序如下: clear,%Clear all variables close all,%Close all figure windows dt = 0.00001; %Specify the step of time variable t = -2:dt:4; %Specify the interval of time w0=0.5*pi; x1=cos(w0.*t); x2=cos(3*w0.*t); x3=cos(5*w0.*t); N=input('Type in the number of the harmonic components N='); x=0; for q=1:N; x=x+(sin(q*(pi/2)).*cos(q*w0*t))/q; end subplot(221) plot(t,x1)%Plot x1 axis([-2 4 -2 2]); grid on, title('signal cos(w0.*t)') subplot(222) plot(t,x2)%Plot x2 axis([-2 4 -2 2]); grid on, title('signal cos(3*w0.*t))') subplot(223) plot(t,x3)%Plot x3 axis([-2 4 -2 2]) grid on, title('signal cos(5*w0.*t))') subplot(224) plot(t,x)%Plot xt axis([-2 4 -2 2]) grid on, title('signal xt')

相关文档
最新文档