习题一 绘制典型信号及其频谱图

习题一绘制典型信号及其频谱图

(1) 绘制单边指数信号及其频谱图的MATLAB程序如下：

close all;

E=1;a=1;

t=0:0.01:4;

w=-30:0.01:30;

f=E*exp(-a*t);

F=E./(a+j*w);

plot(t,f);xlabel('t');ylabel('f(t)');

figure;

plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|');

figure;

plot(w,20*log10(abs(F)));xlabel('\omega');ylabel('|F(\omega)| in dB');

figure;

plot(w,angle(F));xlabel('\omega');ylabel('\phi(\omega)');

请更改参数，调试此程序，绘制单边指数信号的波形图和频谱图。观察参数a对信号波形及其频谱的影响。

上述代码（E=1；a=1）的图形如下所示：

现改变参数再绘制图形：① E=1;a=2;

图形如下所示：

② E=2;a=1; 图形如下所示：

③ E=2;a=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

对正弦信号的采样频谱分析.doc

H a r b i n I n s t i t u t e o f T e c h n o l o g y 课程设计 课程名称：课程设计2 设计题目：对正弦信号的抽样频谱分析院系：电子与信息工程学院 班级：0805203 设计者：褚天琦 学号：1080520314 指导教师：郑薇 设计时间：2011-10-15 哈尔滨工业大学

一、题目要求： 给定采样频率fs，两个正弦信号相加，两信号幅度不同、频率不同。要求给定正弦信号频率的选择与采样频率成整数关系和非整数关系两种情况，信号持续时间选择多种情况分别进行频谱分析。 二、题目原理与分析： 本题目要对正弦信号进行抽样，并使用fft对采样信号进行频谱分析。因此首先对连续正弦信号进行离散处理。实际操作中通过对连续信号间隔相同的抽样周期取值来达到离散化的目的。根据抽样定理，如果信号带宽小于奈奎斯特频率（即采样频率的二分之一），那么此时这些离散的采样点能够完全表示原信号。高于或处于奈奎斯特频率的频率分量会导致混叠现象。设抽样周期为TS(抽样角频率为ωS)，则 可见抽样后的频谱是原信号频谱的周期性重复，当信号带宽小于奈奎斯特频率的二分之一时不会产生频谱混叠现象。 因此，我们对采样频率的选择采取fs>2fo,fs=2fo,fs<2fo三种情况进行分析。对信号采样后，使用fft函数对其进行频谱分析。为了使频谱图像更加清楚，更能准确反映实际情况并接近理想情况，我们采用512点fft。取512点fft不仅可以加快计算速度，而且可以使频谱图更加精确。若取的点数较少，则会造成频谱较大的失真。 三、实验程序： 本实验采用matlab编写程序，实验中取原信号为 ft=sin(2πfXt)+2sin(10πfXt)，取频率f=1kHz，实验程序如下： f=1000;fs=20000;Um=1; N=512;T=1/fs; t=0:1/fs:0.01; ft=Um*sin(2*pi*f*t)+2*Um*sin(10*pi*f*t); subplot(3,1,1); plot(t,ft);grid on; axis([0 0.01 1.1*min(ft) 1.1*max(ft)]); xlabel('t'),ylabel('ft'); title('抽样信号的连续形式'); subplot(3,1,2); stem(t,ft);grid on; axis([0 0.01 1.1*min(ft) 1.1*max(ft)]); xlabel('t'),ylabel('ft');

实验一利用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)|

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 对信号作频谱分析 一、 实验目的与要求 学习用FFT 对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便正确应用FFT 。 二、 实验原理 用FFT 对信号作频分析是学习数字信号处理的重要内容，经常需要进行分析的信号是模拟信号的时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D 和分析误差。频谱分辨率直接和FFT 的变换区间N 有关，因为FFT 能够实现的频率分辨率是2π/N ，因此要求2π/N 小于等于D 。可以根据此式选择FFT 的变换区间N 。误差主要来自于用FFT 作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N 较大时，离散谱的包络才能逼近连续谱，因此N 要适当选择大一些。 三、 实验步骤及内容（含结果分析） （1）对以下序列进行FFT 分析： x 1(n)=R 4(n) x 2(n)= x 3(n)= 选择FFT 的变换区间N 为8和16两种情况进行频谱分析，分别打印出幅频特性曲线，并进行讨论、分析与比较。 【实验结果如下】： n+1 0≤n ≤3 8-n 4≤n ≤7 0 其它n 4-n 0≤n ≤3 n-3 4≤n ≤7 0 其它 n

实验结果图形与理论分析相符。（2）对以下周期序列进行谱分析： x4(n)=cos[(π/4)*n]

x5(n)= cos[(π/4)*n]+ cos[(π/8)*n] 选择FFT的变换区间N为8和16两种情况进行频谱分析，分别打印出幅频特性曲线，并进行讨论、分析与比较。 【实验结果如下】： （3）对模拟周期信号进行频谱分析： x6(n)= cos(8πt)+ cos(16πt)+ cos(20πt) 选择采样频率Fs=64Hz，FFT的变换区间N为16、32、64三种情况进行频谱分析，分别打印出幅频特性曲线，并进行讨论、分析与比较。 【实验结果如下】：

信号与测试实验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处出现峰值以外，在其他频率点也出现了一些峰值。

信号的频谱分析及MATLAB实现

第23卷第3期湖南理工学院学报(自然科学版)Vol.23 No.3 2010年9月 Journal of Hunan Institute of Science and Technology (Natural Sciences) Sep. 2010信号的频谱分析及MATLAB实现 张登奇, 杨慧银 (湖南理工学院信息与通信工程学院, 湖南岳阳 414006) 摘 要: DFT是在时域和频域上都已离散的傅里叶变换, 适于数值计算且有快速算法, 是利用计算机实现信号频谱分析的常用数学工具. 文章介绍了利用DFT分析信号频谱的基本流程, 重点阐述了频谱分析过程中误差形成的原因及减小分析误差的主要措施, 实例列举了MATLAB环境下频谱分析的实现程序. 通过与理论分析的对比, 解释了利用DFT分析信号频谱时存在的频谱混叠、频谱泄漏及栅栏效应, 并提出了相应的改进方法. 关键词: MA TLAB; 频谱分析; 离散傅里叶变换; 频谱混叠; 频谱泄漏; 栅栏效应 中图分类号: TN911.6 文献标识码: A 文章编号: 1672-5298(2010)03-0029-05 Analysis of Signal Spectrum and Realization Based on MATLAB ZHANG Deng-qi, YANG Hui-yin (College of Information and Communication Engineering, Hunan Institute of Science and Technology, Yueyang 414006, China) Abstract:DFT is a Fourier Transform which is discrete both in time-domain and frequency-domain, it fits numerical calculation and has fast algorithm, so it is a common mathematical tool which can realize signal spectrum analysis with computer. This paper introduces the basic process of signal spectrum analysis with DFT, emphasizes the causes of error producing in spectrum analysis process and the main ways to decrease the analysis error, and lists the programs of spectrum analysis based on MATLAB. Through the comparison with the theory analysis, the problems of spectrum aliasing, spectrum leakage and picket fence effect are explained when using DFT to analyze signal spectrum, and the corresponding solution is presented. Key words:MATLAB; spectrum analysis; DFT; spectrum aliasing; spectrum leakage; picket fence effect 引言 信号的频谱分析就是利用傅里叶分析的方法, 求出与时域描述相对应的频域描述, 从中找出信号频谱的变化规律, 以达到特征提取的目的[1]. 不同信号的傅里叶分析理论与方法, 在有关专业书中都有介绍, 但实际的待分析信号一般没有解析式, 直接利用公式进行傅里叶分析非常困难. DFT是一种时域和频域均离散化的傅里叶变换, 适合数值计算且有快速算法, 是分析信号的有力工具. 本文以连续时间信号为例, 介绍利用DFT分析信号频谱的基本流程, 重点阐述频谱分析过程中可能存在的误差, 实例列出MATLAB 环境下频谱分析的实现程序. 1 分析流程 实际信号一般没有解析表达式, 不能直接利用傅里叶分析公式计算频谱, 虽然可以采用数值积分方法进行频谱分析, 但因数据量大、速度慢而无应用价值. DFT在时域和频域均实现了离散化, 适合数值计算且有快速算法, 是利用计算机分析信号频谱的首选工具. 由于DFT要求信号时域离散且数量有限, 如果是时域连续信号则必须先进行时域采样, 即使是离散信号, 如果序列很长或采样点数太多, 计算机存储和DFT计算都很困难, 通常采用加窗方法截取部分数据进行DFT运算. 对于有限长序列, 因其频谱是连续的, DFT只能描述其有限个频点数据, 故存在所谓栅栏效应. 总之, 用DFT分析实际信号的频谱, 其结果必然是近似的. 即使是对所有离散信号进行DFT变换, 也只能用有限个频谱数据近似表示连续频 收稿日期: 2010-06-09 作者简介: 张登奇(1968? ), 男, 湖南临湘人, 硕士, 湖南理工学院信息与通信工程学院副教授. 主要研究方向: 信号与信息处理

信号的频谱分析

实验三信号的频谱分析 方波信号的分解与合成实验 一、任务与目的 1. 了解方波的傅立叶级数展开和频谱特性。 2. 掌握方波信号在时域上进行分解与合成的方法。 3. 掌握方波谐波分量的幅值和相位对信号合成的影响。 二、原理（条件） PC机一台，TD－SAS系列教学实验系统一套。 1. 信号的傅立叶级数展开与频谱分析 信号的时域特性和频域特性是对信号的两种不同的描述方式。对于一个时域的周期信号f(t)，只要满足狄利克莱条件，就可以将其展开成傅立叶级数： 如果将式中同频率项合并，可以写成如下形式： 从式中可以看出，信号f(t)是由直流分量和许多余弦（或正弦）分量组成。其中第一项A0/2是常数项，它是周期信号中所包含的直流分量；式中第二项A1cos(Ωt+φ1)称为基波，它的角频率与原周期信号相同，A1是基波振幅，φ1是基波初相角；式中第三项A2cos(Ωt+φ2)称为二次谐波，它的频率是基波的二倍，A2是基波振幅，φ2是基波初相角。依此类推，还有三次、四次等高次谐波分量。 2. 方波信号的频谱 将方波信号展开成傅立叶级数为： n=1，3，5… 此公式说明，方波信号中只含有一、三、五等奇次谐波分量，并且其各奇次谐波分量的幅值逐渐减小，初相角为零。图3-1-1为一个周期方波信号的组成情况，由图可见，当它包含的分量越多时，波形越接近于原来的方波信号，还可以看出频率较低的谐波分量振幅较大，它们组成方波的主体，而频率较高的谐波分量振幅较小，它们主要影响波形的细节。

（a）基波（b）基波＋三次谐波 （c）基波＋三次谐波＋五次谐波 （d）基波＋三次谐波＋五次谐波＋七次谐波 （e）基波＋三次谐波＋五次谐波＋七次谐波＋九次谐波 图3-1-1方波的合成 3. 方波信号的分解 方波信号的分解的基本工作原理是采用多个带通滤波器，把它们的中心频率分别调到被测信号的各个频率分量上，当被测信号同时加到多路滤波器上，中心频率与信号所包含的某次谐波分量频率一致的滤波器便有输出。在被测信号发生的实际时间内可以同时测得信号所包含的各频率分量。本实验便是采用此方法，实验中共有5路滤波器，分别对应方波的一、三、五、七、九次分量。 4. 信号的合成 本实验将分解出的1路基波分量和4路谐波分量通过一个加法器，合成为原输入的方波信号，信号合成电路图如图3-1-2所示。 图3-1-2 三、内容与步骤 本实验在方波信号的分解与合成单元完成。 1. 使信号发生器输出频率为100Hz、幅值为4V的方波信号，接入IN端。 2. 用示波器同时测量IN和OUT1端，调节该通路所对应的幅值调节电位器，使该通路输出方波的基波分量，基波分量的幅值为方波信号幅值的4/π倍，频率于方波相同并且没有相位差.（注意：出厂时波形调节电位器已调到最佳位置，其波形基本不失真，基本没有相位差。若实验中发现存在波形失真或有相位差的现象，请适当调节波形调节电位器，使波形恢复正常。） 3. 用同样的方法分别在OUT3、OUT5、OUT7、OUT9端得到方波的三、五、七、九此谐波分量（注意其他谐波分量各参数应当满足式3-1-1所示）。 4. 完成信号的分解后，先后将OUT1与IN1、OUT3与IN2、OUT5与IN3、OUT7与IN4、OUT9与IN5连接起来，即进行谐波叠加（信号合成），分别测量（1）基波与三次谐波；（2）基波、三次谐波与五次谐波；（3）基波、三次谐波、五次谐波与七次谐波；（4）基波、三次谐波、五次谐波、七次谐波与九次谐波合成后的波形。并分别保

双代号时标网络计划典型例题

双代号时标网络计划典型例题 2012-07-04 11:47 来源：打印| 收藏| 字号 分享到： 例题1：下图所示的双代号网络图中，存在绘图错误的有（）。 A．循环回路 B．多个起点节点 C．多个终点节点 D．一对节点编号代表多项工作 E．节点编号顺序错误 【正确答案】CDE 【答案解析】选项C，有8、9两个终点节点，故错误；选项D，错在一对节点（1和3）表示了两项工作（B．C工作）；选项E，节点3、4编号顺序错误。参见教材P117、119 例题2：下列关于双代号网络计划绘图规则的说法，正确的有（）。 A．网络图必须正确表达各工作间的逻辑关系 B．网络图中可以出现循环回路 C．网络图中一个节点只有一条箭线引出 D．网络图中严禁出现没有箭头节点或没有箭尾节点的箭线 E．单目标网络计划只有一个起点节点和一个终点节点 【正确答案】ADE

【答案解析】双代号网络计划的绘图规则包括：①双代号网络图必须正确表达已定的逻辑关系；②双代号网络图中，严禁出现循环回路；③双代号网络图中，在节点之间严禁出现带双向箭头或无箭头的连线；④双代号网络图中，严禁出现没有箭头节点或没有箭尾节点的箭线；⑤当双代号网络图的某些节点有多条外向箭线或多条内向箭线时，为使图形简洁，可使用母线法绘制（但应满足一项工作用一条箭线和相应的一对节点表示）；⑥绘制网络图时，箭线不宜交叉。当交叉不可避免时，可用过桥法或指向法；⑦双代号网络图中应只有一个起点节点和一个终点节点（多目标网络计划除外），而其他所有节点均应是中间节点；⑧双代号网络图应条理清楚，布局合理。参见教材P117、119 例题3：某分部工程各工作之间的逻辑关系如下表所示。根据该逻辑关系表绘出的正确网络图是（）。 A．A图 B．B图 C．C图 D．D图

信号与系统实验报告-实验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

信号的频谱分析及MATLAB实现

信号的频谱分析及MATLAB 实现（实例） 摘自：张登奇,杨慧银.信号的频谱分析及MATLAB 实现[J].湖南理工学院学报(自然科学版),2010,(03) 摘 要：DFT 是在时域和频域上都已离散的傅里叶变换，适于数值计算且有快速算法，是利用计算机实现信号频谱分析的常用数学工具。文章介绍了利用DFT 分析信号频谱的基本流程，重点阐述了频谱分析过程中误差形成的原因及减小分析误差的主要措施，实例列举了MATLAB 环境下频谱分析的实现程序。通过与理论分析的对比，解释了利用DFT 分析信号频谱时存在的频谱混叠、频谱泄漏及栅栏效应，并提出了相应的改进方法。 关键词：MATLAB ；频谱分析；离散傅里叶变换；频谱混叠；频谱泄漏；栅栏效应 3 分析实例 对信号进行频谱分析时，由于信号不同，傅里叶分析的频率单位也可能不同，频率轴有不同的定标方式。为了便于对不同信号的傅里叶分析进行对比，这里统一采用无纲量的归一化频率单位，即模拟频率对采样频率归一化；模拟角频率对采样角频率归一化；数字频率对2π归一化；DFT 的k 值对总点数归一化。同时，为了便于与理论值进行对比，理解误差的形成和大小，这里以确定信号的幅度谱分析为例进行分析说明。假设信号为：)()(t u e t x t -=，分析过程：首先利用CTFT 公式计算其模拟频谱的理论值；然后对其进行等间隔理想采样，得到)(n x 序列，利用DTFT 公式计算采样序列的数字连续频谱理论值，通过与模拟频谱的理论值对比，理解混叠误差形成的原因及减小误差的措施；接下来是对)(n x 序列进行加窗处理，得到有限长加窗序列)(n xw ，再次利用DTFT 公式计算加窗后序列)(n xw 的数字连续频谱，并与加窗前)(n x 的数字连续频谱进行对比，理解截断误差形成的原因及减小误差的措施；最后是对加窗序列进行DFT 运算，得到加窗后序列)(n xw 的DFT 值，它是对)(n xw 数字连续频谱进行等间隔采样的采样值，通过对比，理解栅栏效应及DFT 点数对栅栏效应的影响。利用MATLAB 实现上述分析过程的程序如下： clc;close all;clear; %CTFT 程序，以x(t)=exp(-t) t>=0 为例 %利用数值运算计算并绘制连续信号波形 L=4, %定义信号波形显示时间长度 fs=4,T=1/fs; %定义采样频率和采样周期 t_num=linspace(0,L,100);%取若干时点,点数决定作图精度 xt_num=exp(-1*t_num);%计算信号在各时点的数值 subplot(3,2,1);plot(t_num,xt_num),%绘信号波形 xlabel('时间(秒)'),ylabel('x(t)'),%加标签 grid,title('(a) 信号时域波形'),%加网格和标题 %利用符号运算和数值运算计算连续信号幅度谱的理论值 syms t W %定义时间和角频率符号对象 xt=exp(-1*t)*heaviside(t),%连续信号解析式 XW=fourier(xt,t,W),%用完整调用格式计算其傅氏变换 %在0两边取若干归一化频点，点数决定作图精度 w1=[linspace(-0.5,0,50),linspace(0,1.5,150)];

应用MATLAB对信号进行频谱分析

数字信号处理课程设计报告书 2011年7 月 1日 课题名称 应用MATLAB 对信号进行频谱分析 姓 名 张炜玮 学 号 20086377 院、系、部 电气系 专 业 电子信息工程 指导教师 刘鑫淼 ※※※※※※※※※ ※※ ※※ ※※ ※※ ※※※※※ ※※ 2008级数字信号处理课程设计

应用MATLAB对信号进行频谱分析 20086377 张炜玮 一、设计目的 用MATLAB语言进行编程，绘出所求波形，并且运用FFT求对连续信号进行分析。 二、设计要求 1、用Matlab产生正弦波，矩形波，并显示各自的时域波形图； 2、进行FFT变换，显示各自频谱图，其中采样率、频率、数据长度自选，要求注明； 3、绘制三种信号的均方根图谱； 4、用IFFT回复信号，并显示恢复的正弦信号时域波形图。 三、系统原理 用FFT对信号作频谱分析是学习数字信号处理的重要内容。经常需要进行频谱分析的信号是模拟信号和时域离散信号。频谱分辨率直接和FFT的变换区间N 有关，因为FFT能够实现频率分辨率是2π/N。 x(n)是一个长度为M的有限长序列，则x(n)的N点离散傅立叶变换为： X(k)=DFT[x(n)]= kn N W N n n x ∑ - = 1 ) ( ,k=0,1,...,N-1 N j e N Wπ2- = 逆变换：x(n) =IDFT[X(k)]= kn N W k X N n N - ∑ - = 1 ) ( 1 ,k=0,1,...,N-1 但FFT是一种比DFT更加快速的一种算法，提高了DFT的运算速率,为数字信号处理技术应用于各种信号处理创造了条件，大大提高了数字信号处理技术的发展。本实验就是采用FFT，IFFT对信号进行谱分析。 四、程序设计 fs=input('please input the fs:');%设定采样频率 N=input('please input the N:');%设定数据长度 t=0:0.001:1; f=100;%设定正弦信号频率 %生成正弦信号 x=sin(2*pi*f*t);

信号与检验测试实验一

实验一、基本信号分析 一、实验目的 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中产生随机噪声、阶跃信号（选作）、矩形脉冲（选作）

完整版双代号网络计划图习题

双代号网络计划图 一、选择题 1.双代号网络图中分为实工作和虚工作，据此，以下说法正确的有（A ）。 A.混凝土养护工作是实工作，用实箭线表示 B.虚工作不需要消耗时间，但需要消耗资源 C.虚工作的不需要消耗资源，但需要消耗时间 D.虚工作的作用是保证一张网络图只有一个起始节点和一个终点节点 2.已知下列双代号网络图，工作H的紧后工作有（C ）。 A.工作G、C、D B.工作G、C C.工作G、B、C D.工作G、B、C、A、D 3.已知某工程项目分部工程双代号网络计划如下图，其关键线路为（A ）。 A.①→②→③→⑤→⑦ B.①→②→③→⑤→⑦ C.①→③→⑤→⑦ D.①→②→③→⑥→⑦ 4.某工作M有A、B、C三项紧前工作，ES=5，D=3，ES=4，D=5，ES=6，CBBAA D=1，则ES为（ B ）。MC A.8 B.9 C.7 D.6 5.某工作N有A、B、C三项紧后工作，LF=10，D=3，LS=12，LF=15，D=1，CCABA则LF为（A ）。N D.10 C.14 B.0 A.7 6.某工作K的LF=20，D=3，ES=10，则TF为（D ）。KKKK A.20 B.10 C.17 D.7 7.某工作M的最早开始时间为第16天，持续时间为5天。该工作有三项紧后工作，它们的最早开始时间分别为第24天、第27天、第28天，则M工作的自由时差为（A ）。 A.3 B.6 C.7 D.5 8.有M、N两项连续工作，ES=4，ES=14，D=4，LF=30，则（A ）。MNMM A.TF=22，FF=6 B.TF=26，FF=6 MMMM C.TF=22，FF=10 D.TF=26，FF=10 MMMM9.某工作P，已知LF=15、D=2，有M、N、K三项紧前工作，ES=5、

单代号网络图和双代号网络图

单代号网络图与双代号网络图的区别 单代号网络图和双代号网络图所表达的计划内容是一致的，两者的区别仅在于绘图的符号不同。单代号网络图的箭线的含义是表示顺序关系，节点表示一项工作;而双代号网络图的箭线表示的是一项工作，节点表示联系。在双代号网络图中出现较多的虚工作，而单代号网络图没有虚工作。 单代号网络计划的表示方法 1、箭线：箭线既不占用时间，也不消耗资源。箭线仅用来表示工作之间的顺序关系。 2、节点：节点代表一项工作(节点代号、工作名称、作业时间都标注在节点圆圈或方框内)，需占用一定的时间和资源。 3、线路：从网络图的开始节点到结束节点，沿着箭线的指向所构成的若干条 '通道'即为线路。 单代号网络图的绘制 绘图的基本规则： 1)网络图必须按照已定的逻辑关系绘制－－与双代号网络图一致。 2)严禁在网络图中出现没有箭尾节点的箭线和没有箭头节点的箭线。 3)绘制网络图时，宜避免箭线交叉。当交叉不可避免时，可采用过桥法、断线法表示－－与双代号网络图一致。 4)网络图中有多项开始工作或多项结束工作时，就大网络图的两端分别设置一项虚拟的工作，作为该网络图的起点节点及终点节点。 双代号网络计划的表示方法：箭线、节点、线路

1、箭线： (1) 一根箭线表示一项工作，长短不按比例绘制。 (2)每一项工作都要消耗一定的时间和资源。 (3)虚箭线只表达施工过程的逻辑关系，不占用时间，不消耗资源。 (4) 紧靠其前面的工作称紧前工作，紧靠后面的工作叫紧后工作。 2、节点： (1)网络图中表示工作开始、结束或连接关系的圆圈称为节点。 (2)节点只是一个“瞬间”，它既不消耗时间，也不消耗资源。 (3)每个节点都要编号（编号原则是:每一个箭尾节点的号码j必须大于箭头节点的号码i）。 虚工作： 时间为零的假设工作。用虚箭线表示； 特点：不消耗时间和资源。 作用：确切表达网络图中工作之间相互制约、相互联系的逻辑关系(即紧前和紧后工作的关系）。

习题1绘制典型信号及其频谱图

习题一 绘制典型信号及其频谱图 电子工程学院202班 单边指数信号的理论表达式为 figure(4); 调整，将a 分别等于1、5、10等值，观察时域波形和频域波形。由于波形 较多，现不失代表性地将a=1和a=5时的各个波形图列表如下进行对比， 其 他a 值的情况类似可推知。 单边指数信号 信号 名称 单边 时间函数f t 频谱函数F ■ 指数 脉冲 Ee% t a 对提供的MATLAB 程序作了一些说明性的补充， MATLAB 程序为 %单边指数信号 clc; close all ; clear all ; E=1; a=1; %调整a 的值，观察不同a 的值对信号波形和频谱的影响 t=0:0.01:4; w=-30:0.01:30; f=E*exp(-a*t); F=1./(a+j*w); figure(1); plot(t,f);xlabel( 't' );ylabel( 'f(t)' );title( '信号时域图像’); figure(2); plot(w,abs(F));xlabel( '\omega' 特性'); figure (3); plot(w,20*log10(abs(F)));xlabel( );ylabel( '|F(\omega)|' ); ti tle( '幅频 '\omega' );ylabel( '|F(\omega)| in dB' );title( 幅频特性/dB'); plot(w,a ngle(F)*57.29577951);xlabel( )/ (°)' );title( '相频特性’); '\omega' );ylabel( '\phi(\omega

如何绘制简单网络图

如何绘制简单网络图 导语： 网络图(Network planning)是一种图解模型，形状如同网络，故称为网络图。网络图是由作业（箭线）、事件（又称节点）和路线三个因素组成的。绘制一张简单的网络图，其实并不难，使用专业的软件操作起来就非常轻松。 免费获取网络拓扑图软件：https://www.docsj.com/doc/8f2112106.html,/network/ 如何绘制简单网络图？ 亿图网络图绘制作软件是由亿图软件公司推出的一款专门用来绘制电脑网络计划图的软件。软件功能强大，容易上手，几乎包含所有网络图的绘制，例如基本网络图、双代号网络图、单代号网络图、Cisco网络图、机架图、网络通信图、3D网络图、AWS图等等，可以完美替代Visio。软件采用拖拽的绘图方式，界面简单明了，操作方便，用户即看机即会，无需花费多少时间学习。

软件不仅提供各种专业图库，还提供海量模板，这点是其他软件无法比拟的。强大的定制功能使得用户不仅可以自定义图形的填充和线条颜色，也可以自行绘制图库里的形状。一键导出到PDF，Word, Visio, Png 等17种文件格式，无障碍与他人分享。新版本不仅实现了跨平台，而且还支持云存储，使得团队协作更加容易。亿图网络图绘制软件是您绘制单代号网络图的不二选择。

双代号网络图的三要素： 1、工作：即箭线。是用来表示一项工作需要消耗的人力、物力以及时间的具体活动过程。也称之为工序或者作业。 2、节点：网络图中，箭线出发点和交汇处会使用圆圈，来表示该圆圈前面一项或者多项工作的结束以及允许后面一项或者多项工作的开始的时间，即为节点。 3、线路：在网络图中，从起点节点开始，沿箭头方向连续通过一系列箭线与节点，最后到达终点节点的通路称为线路。 双代号网络图的绘制规则： 1、不能出现从一个节点出发，又顺着箭头方向又回到原出发点的循环回路。 2、不能有两个或者两个以上的箭线从同一节点出发又同时指向同一节点。 3、网络图中的箭线(包括虚箭线，以下同)应保持自左向右的方向，不应出现箭头指向左方的水平箭线和箭头偏向左方的斜向箭线。 4、网络图中的节点编号应该自左向右，由小到大，确保工作的起点节点的编号小于工作终点节点的编号，而且不能出现重复编号的节点，但是可以是非连续的编号。 5、网络图中不能出现双向箭头和无箭头的连线。 6、网络图中不允许出现没有箭尾节点的箭线和没有箭头节点的箭线。 7、禁止在箭线上引入或者引出其他箭线，除非网络图的起点节点有多条箭线引出(外向箭头)或者重点节点有多条箭线引入(内向箭线)，此时为了让图形更简洁，则可以使用母线法绘图。 8、绘制网络图时，宜避免箭线交叉。但是当交叉不可避免时，也可以采用过桥法或指向法进行处理。 9、一个网络图中只能有一个起点节点和终点节点。 10、同一项工作不能在一个网络图中重复表达。

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

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

实验一离散信号的频谱分析报告 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);%进行对应的频率转换