实验四 信号的频谱分析
一．实验目的
1.掌握利用FFT 分析连续周期，非周期信号的频谱，如周期，非周期方波，正弦信号等。理解CFS ，CTFT 与DFT （FFT ）的关系。
2.利用FFT 分析离散周期，非周期信号的频谱，如周期，非周期方波，正弦信号等。理解DFS ，DTFT 与DFT （FFT ）的关系，并讨论连续信号与离散信号频谱分析方法的异同。
二．实验要求
1.编写程序完成任意信号数字谱分析算法；
2.编写实验报告。
三．实验容
1.利用FFT ，分析并画出sin(100),cos(100)t t ππ频谱，改变采样间隔与截断长度，分析混叠与泄漏对单一频率成分信号频谱的影响。
（1）sin （100*pi*t ）产生程序：
close all ;
clc;
clear;
t=0:0.0025:0.5-0.0025;
f=400*t;
w0=100*pi;
y=sin(w0*t);
a=fft(y);
b=abs(a)/200;
d=angle(a)*180/pi;
subplot(311);
plot(t,y);
title('y=sin(wt)');
xlabel('t');
ylabel('y(t)');
subplot(312);
stem(f,b);
title('振幅');
xlabel('f');
ylabel('y(t)');
subplot(313);
stem(f,d);
title('相位');
xlabel('t');
ylabel('y(t)');
混叠
close all;
clc;
clear;
t=0:0.0115:0.46-0.0115; f=(t/0.0115)*2;
w0=100*pi;
y=sin(w0*t);
a=fft(y);
b=abs(a)/40;
d=angle(a)*180/pi; subplot(311);
plot(t,y);
title('y=sin(wt)'); xlabel('t');
ylabel('y(t)');
subplot(312);
stem(f,b);
title('振幅');
xlabel('f');
ylabel('y(t)');
subplot(313);
stem(f,d);
title('相位');
xlabel('t');
ylabel('y(t)');
泄漏
close all;
clc;
clear;
t=0:0.0025:0.5-0.0075; f=800*t;
w0=100*pi;
y=sin(w0*t);
a=fft(y);
b=abs(a)/;
d=angle(a)*180/pi; subplot(311);
plot(t,y);
title('y=sin(wt)'); xlabel('t');
ylabel('y(t)');
subplot(312);
stem(f,b);
title('振幅');
xlabel('f');
ylabel('y(t)');
subplot(313);
stem(f,d);
title('相位');
xlabel('t');
ylabel('y(t)');
（2）cos(100*pi*t);
close all;
clc;
clear;
t=0:0.0025:0.5-0.0025; f=800*t;
w0=100*pi;
y=cos(w0*t);
a=fft(y);
b=abs(a)/200;
d=angle(a)*180/pi; subplot(311);
plot(t,y);
title('y=cos(wt)'); xlabel('t');
ylabel('y(t)');
grid on;
hold on;
subplot(312);
stem(f,b);
title('振幅');
xlabel('f');
ylabel('y(t)');
grid on;
hold on;
subplot(313);
stem(f,d);
title('相位');
xlabel('f');
ylabel('y(t)');
混叠
close all;
clc;
clear;
t=0:0.0115:0.46-0.0115; f=(t/0.0115)*2;
w0=100*pi;
y=cos(w0*t);
a=fft(y);
b=abs(a)/40;
d=angle(a)*180/pi; subplot(311);
plot(t,y);
title('y=cos(wt)'); xlabel('t');
ylabel('y(t)');
subplot(312);
stem(f,b);
title('振幅');
xlabel('f');
ylabel('y(t)');
subplot(313);
stem(f,d);
title('相位');
xlabel('t');
ylabel('y(t)');
泄漏
close all;
clc;
clear;
t=0:0.0025:0.5-0.0075; f=800*t;
w0=100*pi;
y=cos(w0*t);
a=fft(y);
b=abs(a)/;
d=angle(a)*180/pi; subplot(311);
plot(t,y);
title('y=cos(wt)'); xlabel('t');
ylabel('y(t)');
subplot(312);
stem(f,b);
title('振幅');
xlabel('f');
ylabel('y(t)');
subplot(313);
stem(f,d);
title('相位');
xlabel('t');
ylabel('y(t)');
实验 应用FFT 对信号进行频谱分析 一、实验目的 1、在理论学习的基础上，通过本次实验，加深对快速傅里叶变换的理解，熟悉FFT 算法及其程序的编写。 2、熟悉应用FFT 对典型信号进行频谱分析的方法。 3、了解应用FFT 进行新红啊频谱分析过程中可呢个出现的问题，以便在实际中正确应用FFT 。 二、实验原理 一个连续信号()a x t 的频谱可以用它的傅里叶变换表示为： ()()j t a a X j x t e dt +∞ -Ω-∞Ω=? （2-1） 如果对信号进行理想采样，可以得到离散傅里叶变换： ()()j n X e x n z ω +∞ --∞=∑ （2-2） 在各种信号序列中，有限长序列在数字信号处理中占有很重要的。无限长的序列往往可以用有限长序列来逼近。对于有限长的序列我们可以使用离散傅里叶变换（DFT ），这一序列可以很好的反应序列的频域特性，并且容易利用快速算法在计算机上实现当序列的长度是N 时，我们定义离散傅里叶变换为： 1 0()[()]()N kn N n X k DFT x n x n W -===∑ （2-3） DFT 是对序列傅里叶变换的灯具采样，因此可以用于序列的频谱分析。在利用DFT 进行频谱分析的时候可能有三种误差： （1）混叠现象 序列的频谱是采样信号频谱的周期延拓，周期是2/T π，因此当采样频率不满足奈奎斯特定理，即采样频率1/s f T =小于两倍的信号频率时，经过采样就会发生频谱混叠。这导致采样后的信号序列不能真实的反映原信号的频谱。 （2）泄漏现象 泄漏是不能和混叠完全分开的，因为泄漏导致频谱的扩展，从而造成混淆。为了减小混淆的影响，可以选择适当的窗函数使频谱的扩散减到最小。 （3）栅栏效应 因为DFT 是对单位圆上Z 变换的均匀采样，所以它不可能将频谱视为一个连续的函数。这样就产生了栅栏效应。减小栅栏效应的一个方法是在源序列的末端补一些零值，从而变动DFT 的点数。 三、实验内容和结果 1、观察高斯序列的时域和频域特性 （1）固定高斯序列()a x n 中的参数p=8，当q 为2，4，8时其时域和幅频特性分别如图 2.1，图2.2所示：
实验报告 实验三：用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 引言 随着科学技术的迅猛发展，电子设备和技术向集成化、数字化和高速化方向发展，而在学校特别是大学中，要想紧跟技术的发展，就要不断更新教学和实验设备。传统仪器下的高校实验教学，已严重滞后于信息时代和工程实际的需要。仪器设备很大部分陈 旧，而先进的数字仪器(如数字存储示波器)价格昂贵不可能大量采购，同时其功能较为单一，与此相对应的是大学学科分类越来越细，每一专业都需要专用的测量仪器，因此仪器设备不能实现资源共享，造成了浪费。虚拟仪器正是解决这一矛盾的最佳方案。基于PC 平台的虚拟仪器，可以充分利用学校的微机资源，完成多种仪器功能，可以组合成功能强大的专用测试系统，还可以通过软件进行升级。在通用计算机平台上，根据测试任务的需要来定义和设计仪器的测试功能，充分利用计算机来实现和扩展传统仪器功能，开发结构简单、操作方便、费用低的虚拟实验仪器，包括数字示波器、频谱分析仪、函数发生器等，既可以减少实验设备资金的投入，又为学生做创新性实验、掌握现代仪器技术提供了条件。 信号的时域分析主要是测量测试信号经滤波处理后的特征值，这些特征值以一个数值表示信号的某些时域特征，是对测试信号最简单直观的时域描述。将测试信号采集到计算机后，在测试VI 中进行信号特征值处理，并在测试VI 前面板上直观地表示出信号的特征值，可以给测试VI 的使用者提供一个了解测试信号变化的快速途径。信号的特征值分为幅值特征值、时间特征值和相位特征值。 尽管测量时采集到的信号是一个时域波形，但是由于时域分析工具较少，所以往往把问题转换到频域来处理。信号的频域分析就是根据信号的频域描述来估计和分析信号的组成和特征量。频域分析包括频谱分析、功率谱分析、相干函数分析以及频率响应函数分析。 信号在时域被抽样后,他的频谱X(j )是连续信号频谱X(j )的形状以抽样频率为间隔周期重复而得到,在重复过程中幅度被p(t)的傅里叶级数Pn加权。因为Pn只是n的函数,所以X(j )在重复的过程中不会使其形状发生变化。假定信号x(t)的频谱限制在- m~+ m的范围内, 若以间隔Ts对xa(t)进行抽样,可知抽样信号X^(t)的频谱X^(j )是以s为周期重复。显然,若在抽样的过程中s<2 m,则X^(j )将发生频谱混叠现象,只有在抽样的过程中满足s>=2 m条件,X^(j )才不会产生频谱的混叠,接收端完全可以由x^(t)恢复原连续信号xa(t),这就是低通信号抽样定理的核心内容。
频域分析实验报告 班级： 学号： 姓名：
一、实验内容： 1利用计算机作出开环系统的波特图； 2、观察记录控制系统的开环频率特性； 3、控制系统的开环频率特性分析。 二、仿真原理： 对数频率特性图（波特图）： 对数频率特性图包括了对数幅频特性图和对数相频特性图。横坐标为频率w，采用对数分度，单位为弧度/秒；纵坐标均匀分度，分别为幅值函数20lgA(w)，以dB表示；相角，以度表示。MATLAB提供了函数bode()来绘制系统的波特图，其用法如下： （1）bode(num，den)：可绘制出以连续时间多项式传递函数表示的系统的波特图。 （2）当带输出变量[mag，pha，w]或[mag，pha]引用函数时，可得到系统波特图相应的幅值mag、相角pha及角频率点w矢量或只是返回幅值与相角。相角以度为单位，幅值可转换为分贝单位：magdb=20×log10(mag) 二、实验验证 1、用Matlab作Bode图。要求:画出对应Bode图。 （1）G(S)=25/S2+4s+25 （7）G(S)=9(s2+0.2s+1)/s(s2+1.2s+9)；
图 1 图 2 (1)G(S)=25/S2+4s+25 可以看成是一个比例环节和一个振荡环节组成，所以k=1，T1=0.04,因为v=0，所以在转折频率之前都为20lgk，因为k=1所以斜率为0，经过转折频率，分段直线斜率的变化量为-40db/dec。
（7）G(S)=9(s2+0.2s+1)/s(s2+1.2s+9)； 可以看成是一个二阶微分环节和一个积分环节和一个振荡环节组成，化常数为1后，v=1，t1=1，t2=1/3，所以我们可以看到，在起始阶段是-20*vdb/dec，所以一开始斜率为-20db/dec。当经过1/3的转折频率之后分段直线的改变量为40db/dec，当经过1的转折频率之后分段直线的改变量为-40db/dec。故图像如图所示。 第二题： 典型二阶系统Gs=Wn2/s2+2ζWns+Wn2，试绘制取不同值时的Bode图。取Wn=8，ζ=0.1，0.2，0.3，,0.5，0.6； 图 3 如图所示。
周期信号频谱的特点 在结构施工测量中，按装修工程要求将装饰施工所需要的控制点、线及时弹在墙、板上，作为装饰工程施工的控制依据。 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 限差时，原则上不予验收，尤其是重要部位，
实验名称： 控制系统的频域分析 实验类型：________________同组学生姓名：__________ 一、实验目的和要求 用计算机辅助分析的方法，掌握频率分析法的三种方法，即Bode 图、Nyquist 曲线、Nichols 图。 二、实验内容和原理 （一）实验原理 1．Bode(波特)图 设已知系统的传递函数模型： 1 1211121)(+-+-+???+++???++=n n n m m m a s a s a b s b s b s H 则系统的频率响应可直接求出： 1 1211121)()()()()(+-+-+???+++???++=n n n m m m a j a j a b j b j b j H ωωωωω MATLAB 中，可利用bode 和dbode 绘制连续和离散系统的Bode 图。 2．Nyquist(奈奎斯特)曲线 Nyquist 曲线是根据开环频率特性在复平面上绘制幅相轨迹，根据开环的Nyquist 线，可判断闭环系统的稳定性。 反馈控制系统稳定的充要条件是，Nyquist 曲线按逆时针包围临界点(-1，j0)p 圈，为开环传递函数位于右半s 一平面的极点数。在MATLAB 中，可利用函数nyquist 和dnyquist 绘出连续和离散系统的乃氏曲线。 3．Nicho1s(尼柯尔斯)图 根据闭环频率特性的幅值和相位可作出Nichols 图，从而可直接得到闭环系统的频率特性。在 MATLAB 中，可利用函数nichols 和dnichols 绘出连续和离散系统的Nichols 图。 （二）实验内容 1．一系统开环传递函数为 ) 2)(5)(1(50)(-++=s s s s H 绘制系统的bode 图，判断闭环系统的稳定性，并画出闭环系统的单位冲击响应。 2．一多环系统 ) 10625.0)(125.0)(185.0(7.16)(+++=s s s s s G 其结构如图所示 试绘制Nyquist 频率曲线和Nichols 图，并判断稳定性。 （三）实验要求
声速的测量实验报告 不会写声速的测量实验报告的朋友，下面请看我给大家整理收集的声速的测量实验报告，仅供参考。 声速的测量实验报告1 实验目的:测量声音在空气中的传播速度。 实验器材:温度计、卷尺、秒表。 实验地点:平遥县状元桥东。 实验人员:爱物学理小组 实验分工:张灏、成立敬——测量时间 张海涛——发声 贾兴藩——测温 实验过程: 1 测量一段开阔地长; 2 测量人在两端准备; 3 计时员挥手致意,发声人准备发声; 4 发生人向上举手,同时发声,计时员计时(看到举手始,听到声音止) 5 多测几次,记录数据。 实验结果: 时间17∶30 温度21℃
发声时间 0.26″ 发声距离 93m 实验结论:在21℃空气中,声音传播速度为357.69m/s. 实验反思:有一定误差,卡表不够准确。 声速的测量实验报告2 实验目的： 1)探究影响声速的因素，超声波产生和接收的原理。 2)学习、掌握空气中声速的测量方法 3)了解、实践液体、固体中的声速测量方法。 4)三种声速测量方法作初步的比较研究。 实验仪器: 1)超声波发射器 2)超声波探测器 3)平移与位置显示部件。 4)信号发生器： 5)示波器 实验原理： 1)空气中： a.在理想气体中声波的传播速度为 v88 (式中8088cp cV (1) 称为质量热容比，也称"比热[容]比"，它是气体的质 量定压热容cp与质量定容热容cV的比值;M 是气体的摩尔质量，T 是绝对温度，R=8.314472(1±1.7×10-6)Jmol-1K-1为摩尔气体常量。)
标准干燥空气的平均摩尔质量为Mst =28.966�8�710-3kg/mol b.在标准状态下 (T0�8�8273.15 K,p�8�8101.3�8�8kPa)，干燥空气中的声速 为v0=331.5m/s。在室温t℃下，干燥空气中的声速为 v88v0 (2) (T0=273.15K) c.然而实际空气总会有一些水蒸气。当空气中的相对湿度为r时，若气温为t℃时饱和蒸气压为pS，则水汽分压为rps。经过对空气平均摩尔质量 M 和质量热容比8�0 的修正，在温度为t、相对湿度为r 的空气中，声速为 (在北京大气压可近似取p�8�4 101kPa;相对湿度r 可从干湿温度计上读出。温度t℃时的饱和水汽压ps可用 lgps�8�810.286�8�2 d.式(3)的计算结果与实际的超声声速真值可能有一定偏差。 引起偏差的原因有： ~状态参量的测量误差 ~理想气体理论公式的近似性 实验方法： A. 脉冲法：利用声波传播时间与传播距离计算声速 实验中用脉冲法测量，具体测量从脉冲声源(声发射器)到声探测器
信号频谱分析和测 试 返回 一、实验室名称：虚拟仪器实验室 二、实验项目名称：信号频谱分析和测试 三、实验目的 1．了解周期函数的傅立叶变换理论及虚拟频谱分析仪的工作原理； 2．熟悉典型信号的波形和频谱特征，并能够从信号频谱中读取所需的信息。 四、实验内容 1．测量典型信号（正弦波、三角波、方波）的频谱并记录； 2．用实验平台的任意波形信号源产生一个任意信号，观察其频谱。 五、实验器材（设备、元器件）： 1、计算机一台 2、SJ-8002B 电子测量实验箱一台 3、FG1617函数发生器一台 4、虚拟频谱分析仪程序 5、Q9线一条 六、实验原理 6.1 常见周期信号傅立叶展开公式与波形 1）方波 ，其中的 2）三角波 ，其中的 )7sin 715sin 513sin 31(sin 4)( +ω+ω+ ω+ωπ=t t t t A t f T π=ω2)7cos 4915sin 2513sin 91(sin 8)(2 +ω-ω+ω-ωπ=t t t t A t f T π=ω2
3）锯齿波 ，其中 6.2 信号的离散傅立叶变换（DFT ） x(t)经采样后变为x(nT ’)，T ’为采样周期，采样频率fs=1/T ’。离散信号x(nT ’)的傅里 叶变换可以表示为： ，n=0，1，…N-1 X(k)是复数，信号的频谱是它的模，为了方便显示，做归一化处理，用 来表示频谱。 频率分辨率为： FFT 是DFT 的快速算法。 6.3 虚拟频谱分析仪 数字式虚拟频谱分析仪是通过A/D 采样器件，将模拟信号转换为数字信号，传给微处 理器系统或计算机来处理.在对交流信号的测量中，根据奈奎斯特采样定理，采样速率必须 是信号频率的两倍以上，采样频率越高，时间轴上的信号分辨力就越高，所获得的信号就越 接近原始信号，在频谱上展现的频带就越宽。 本频谱分析仪采用快速傅立叶变换的方法，分析信号中所含各个频率份量的幅值。其构 成框图如图4所示： 图4频谱分析仪框图 七、实验步骤 7．1 测量典型信号（正弦波、三角波、方波）的频谱 （1） 准备工作：用Q9线连接信号发生器与实验平台的Ain1端，并用EPP 排线连接实 验平台和计算机之间的EPP 接口，最后打开电源.。信号发生器产生一个频率为10K ，峰峰 值为3V 左右的正弦波，启动实验平台配套的频谱分析软件，观察波形显示并作图。 （2）由信号源产生一个频率为10KHz ，峰值为3V 的正弦波，用数字频谱分析仪对该信 号进行频谱测量，幅度刻度方式设为线性刻度，不加窗函数，起始频率为0Hz ，结束频率为 100KHz ，Y 线性参考电压为2V ，将测量结果填入表1，并计算出频谱的理论值填入表1。 )4sin 413sin 312sin 21(sin 2)( +ω+ω+ω+ωπ+= t t t t A A t f T π=ω2()()N nk j N n e n x k X /210π--=∑=N k X )(f ?N f f s =?N kf k f f s k =??=
实验三 连续时间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 t j ωω)()( 3.3
由于H(j ω)实际上是系统单位冲激响应h(t)的傅里叶变换，如果h(t)是收敛的，或者说是绝对可积（Absolutly integrabel ）的话，那么H(j ω)一定存在，而且H(j ω)通常是复数，因此，也可以表示成复数的不同表达形式。在研究系统的频率响应时，更多的是把它表示成极坐标形式： ) ()()(ω?ωωj e j H j H = 3.4 上式中，)j (ωH 称为幅度频率相应（Magnitude response ），反映信号经过系统之后，信号各频率分量的幅度发生变化的情况，)(ω?称为相位特性（Phase response ），反映信号经过系统后，信号各频率分量在相位上发生变换的情况。)(ωj H 和)(ω?都是频率ω的函数。 对于一个系统，其频率响应为H(j ω)，其幅度响应和相位响应分别为)(ωj H 和)(ω?，如果作用于系统的信号为t j e t x 0)(ω=，则其响应信号为 t j e j H t y 0)()(0ωω= t j j e e j H 00)(0)(ωω?ω=))((000)(ω?ωω+=t j e j H 3.5 若输入信号为正弦信号，即x(t) = sin(ω0t )，则系统响应为 ))(sin(|)(|)sin()()(00000ω?ωωωω+==t j H t j H t y 3.6 可见，系统对某一频率分量的影响表现为两个方面，一是信号的幅度要被)(ωj H 加权，二是信号的相位要被)(ω?移相。 由于)(ωj H 和)(ω?都是频率ω的函数，所以，系统对不同频率的频率分量造成的幅度和相位上的影响是不同的。 2 LTI 系统的群延时 从信号频谱的观点看，信号是由无穷多个不同频率的正弦信号的加权和（Weighted sum ）所组成。正如刚才所述，信号经过LTI 系统传输与处理时，系统将会对信号中的所有频率分量造成幅度和相位上的不同影响。从相位上来看，系统对各个频率分量造成一定的相位移（Phase shifting ），相位移实际上就是延时（Time delay ）。群延时（Group delay ）的概念能够较好地反
清镇市中小学实验报告单 实验科目：科学 真题：您是市政府的工作人员，需要调查大学生村官的工作情况，面对基层群众(农民)、基层工作人员分别作一个开场白。
【题型分析】情景模拟 【参考答案】 1.面对基层群众(农民)：大叔大婶，你们好!我是咱们市政府的小王。眼看着一年又快到头了，今年的天儿不错，家里收成还好吧?……听说咱们这儿来了些大学生，到咱这儿当村官，不知道干的咋样啊?…… 2.面对基层工作人员：XX，你好!我是市政府的工作人员王XX(同时出示工作证)，这次来到这里，主要是想调查一下大学生村官的工作情况，希望得到你们的配合和协助。你们长期在基层工作，比较辛苦，对于大学生们在村里的工作情况也比较了解。知道你工作很忙，所以我会尽量抓紧时间…… 真题：金无足赤，人无完人，您有什么缺点和不足?假如您被录用，您将怎么克服您的缺点和不足? 【题型分析】认识自己 【参考答案】每一个人都有自己的缺点和不足，我也一样。我认为我的缺点是有时候有点追求完美，这样一来有可能会拖慢办事的效率，偶尔也会给合作伙伴带来压力。有缺点和不足并不可怕，重要的是能够改正缺点，弥补不足。我将从以下几个方面加强学习，进行改正。 第一，摆正心态正确认识。在遇到每一项工作的时候，我都应该准确地把握工作实质、明确目标，在思想上做好充分准备。不能因为自己的喜好而偏离工作应有的方向和要求。 第二，制定工作计划表。在开始每一项工作的时候，我都要制定好详细的工作进度控制表，将每个阶段的工作目标、内容以及时间详细列好。在开展工作的时候严格按照计划表进行，以免拖慢工作的进程。借此来养成好的工作习惯，以提高办事效率。 第三，加强与同事合作。每一项工作都不可能是自己独立一个人能够做好的，都是需要和同事一起合作才能将其做到真正的完美。这就要求我要在以后的工作中，多多与同事交流，虚心向同事学习。 假如有幸我能够通过此次考试，在以后的工作岗位上，我将严格要求自己，加强学习，发扬优点，改正缺点，竭尽所能做好领导交办的每一项任务，做一名合格的公务员。 真题：您作为单位的新进工作人员，领导让您制定一份单位的规章制度，您对单位不熟悉，您将怎么办? 【题型分析】如何做事 【参考答案】俗话说：“无规矩不成方圆”，好的规章制度能够在保证单位良好工作秩序方面发挥积极的作用。领导将这项任务交予我完成，我一定会努力做好，我将从以下几个方面入手，开展此项工作。 首先，我要请示领导，明白领导制定此规章制度的宗旨和意图。鉴于我作为新进员工对工作单位还不熟悉，我会向同事进行询问，翻查单位档案，加深对单位的人员环境、工作环境的了解。另外，要广泛收集单位同事的意见。我将采取问卷调查的方式进行意见收集，可以使用相关软件来对收回的有效问卷进行信息处理和分析，切实了解员工的意见。此外我还要刻苦学习相关管理学知识，了解国家相关法律法规，为制定一份合理的单位规章制度做好充分的理论准备。 其次，在书写规章制度时，要遵守国家法律法规的规定，制定出合法有效的规章制度。与此同时，我还要根据之前所整理的员工意见，兼顾单位的实际情况，制定一份合乎人情的规章制度。在规章制度雏形完成之后，还要让单位同事对此发表自己的意见，在充分了解意见之后对其进行反复修改。与此同时，规章制度还要做到简洁凝练，通俗易懂。在修订好初
实验一报告、用FFT 对信号作频谱分析 一、实验目的 学习用FFT 对连续信号和时域离散信号进行频谱分析的方法，了解可能出现的分析误差及其原因，以便正确应用FFT 。 二、实验内容 1．对以下序列进行频谱分析： ()() ()()4231038470n 4033 470n x n R n n n x n n n n n x n n n =+≤≤?? =-≤≤???-≤≤?? =-≤≤??? 其它其它 选择FFT 的变换区间N 为8和16两种情况进行频谱分析。分别打印其幅频特性曲线，并进行对比，分析和讨论。 2．对以下周期序列进行频谱分析： ()()45cos 4 cos cos 4 8 x n n x n n n π π π ==+ 选择FFT 的变换区间N 为8和16两种情况分别对以上序列进行频谱分析。分别打印其幅频特性曲线，并进行对比、分析和讨论。 3.对模拟信号进行频谱分析： ()8cos8cos16cos20x t t t t πππ=++ 选择采样频率64s F Hz =，对变换区间N=16，32，64 三种情况进行频谱分析。分别 打印其幅频特性，并进行分析和讨论。
三、实验程序 1.对非周期序列进行频谱分析代码： close all;clear all; x1n=[ones(1,4)]; M=8;xa=1:(M/2);xb=(M/2):-1:1;x2n=[xa,xb]; x3n=[xb,xa]; X1k8=fft(x1n,8);X1k16=fft(x1n,16); X2k8=fft(x2n,8);X2k16=fft(x2n,16); X3k8=fft(x3n,8);X3k16=fft(x3n,16); subplot(3,2,1);mstem=(X1k8);title('(1a)8点DFT[x_1(n)]'); subplot(3,2,2);mstem=(X1k16);title('(1b)16点DFT[x_1(n)]'); subplot(3,2,3);mstem=(X2k8);title('(2a)8点DFT[x_2(n)]'); subplot(3,2,4);mstem=(X2k16);title('(2b)16点DFT[x_2(n)]'); subplot(3,2,5);mstem=(X3k8);title('(3a)8点DFT[x_3(n)]'); subplot(3,2,6);mstem=(X3k16);title('(3b)16点DFT[x_3(n)]'); 2.对周期序列进行频谱分析代码： N=8;n=0:N-1; x4n=cos(pi*n/4); x5n=cos(pi*n/4)+cos(pi*n/8); X4k8=fft(x4n); X5k8=fft(x5n); N=16;n=0:N-1; x4n=cos(pi*n/4); x5n=cos(pi*n/4)+cos(pi*n/8); X4k16=fft(x4n); X5k16=fft(x5n); figure(2) subplot(2,2,1);mstem(X4k8);title('(4a)8点 DFT[x_4(n)]'); subplot(2,2,2);mstem(X4k16);title('(4b)16点DFT[x_4(n)]'); subplot(2,2,3);mstem(X5k8);title('(5a)8点DFT[x_5(n)]'); subplot(2,2,4);mstem(X5k16);title('(5a)16点DFT[x_5(n)]') 3.模拟周期信号谱分析 figure(3) Fs=64;T=1/Fs; N=16;n=0:N-1; x6nT=cos(8*pi*n*T)+cos(16*pi*n*T)+cos(20*pi*n*T); X6k16=fft(x6nT); X6k16=fftshift(X6k16);
实验二：用FFT作谱分析 一、实验目的 (1) 进一步加深DFT算法原理和基本性质的理解(因为FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质)。 (2) 熟悉FFT算法原理和FFT子程序的应用。 (3) 学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。 二、实验原理
三、实验内容 (1) 对2 中所给出的信号逐个进行谱分析。解：（1） n=0:3; xn1=[1 1 1 1]; XK18=fft(xn1,8); XK116=fft(xn1,16); n1=0:7; n2=0:15; subplot(131); stem(n,xn1); xlabel('n'); ylabel('xn1'); subplot(132); stem(n1,abs(XK18)); xlabel('n1'); ylabel('XK18'); title('xn的8点'); subplot(133); stem(n2,abs(XK116)); xlabel('n2'); ylabel('XK116'); title('xn的16点');
（2） n1=0:7; n2=0:15; xn2=[1 2 3 4 4 3 2 1]; XK28=fft(xn2,8); XK216=fft(xn2,16); subplot(131); stem(n1,xn2); xlabel('n1'); ylabel('xn2'); subplot(132); stem(n1,abs(XK28)); xlabel('n1'); ylabel('XK28'); title('xn2的8点'); subplot(133); stem(n2,abs(XK216)); xlabel('n2'); ylabel('XK216'); title('xn2的16点'); （3） n1=0:7; n2=0:15; xn3=[4 3 2 1 1 2 3 4]; XK38=fft(xn3,8);
第一章 填空： 1、信号一般分为 确定性信号 和 随机信号 两类。 2、信号的描述方法常用的有 实域描述 和 频域描述 两种。 3、周期信号用 傅立叶级数 展开到频域描述； 瞬变非周期信号用 傅立叶变换 展开到频域描述。 4、周期信号频谱的特点是：离散性、谐波性、收敛性； 瞬变非周期信号频谱的特点是 连续性 。 5、已知时域信号 x (t )，其傅里叶变换是 。 6、sin 2πf0t 的傅里叶变换 cos 2πf0t 的傅里叶变换 判断： 1、信号的时域描述与频域描述包含相同的信息量。（ √ ） 2、非周期信号的频谱一定是连续的。（ × ） 3、非周期信号幅频谱与周期信号幅值谱的量纲一样。（ × ） 简答： 1、何为信号？如何建立其模型？ 答：蕴含信息且能传输信息的物理量 在测试技术中，撇开信 2()()j ft X f x t e dt π∞--∞=?
号具体的物理性质，而是将其抽象为某个变量的函数关系，如时间的 函数x(t)、频率的函数X(f)等。这些函数就是对信号进行分析、处 理时的数学模型。 2、信号有哪些分类？ 答：确定性信号和随机信号、连续信号和离散信号、能量信号和 功率信号 3、模拟信号与数字信号如何定义？ 答：模拟信号：独立变量和幅值均取连续值的信号。 数字信号：独立变量和幅值均取离散值的信号。 4、什么是信号的时域描述和频域描述？两者有何区别？ 答：直接观察或记录到的信号，一般是以时间为独立变量，反映 的是信号幅值随时间的变化关系，因而称其为信号的时域描述。 若把信号变换成以频率为独立变量，由此来反映信号的频率结构 和各频率成分与幅值、相位之间的关系，信号的这种描述方法称之为 频域描述。 变量：时域描述以时间为变量；频域描述以频率为变量。 内容：时域描述反映信号幅值与时间的关系；频域描述反映信 号的幅值、相位与频率的关系。 用途：时域描述用于评定振动的烈度等；频域描述用于寻找振源、 故障诊断等。 5、什么是时间尺度改变特性？其对测试工作有何意义？ 答： 对测试工作的意义： ()()x t X f ?若：
频谱分析实验报告 许开龙 热能工程系2008010717 一、实验目的 通过实验，了解频谱分析的原理，掌握数据处理中的这一重要手段。 二、实验方法 1.预习实验原理，搞清程序流程和各参数的含义。 2.自己编制一个产生两个正弦波之和的程序，即, 其中A1,A2 分别为正弦波幅值，K 1=Fs/F 1 , K2=Fs/F2, Fs为采样频率,F1,F2分别为正弦波频率。将 产生的数据放入数据文件中，数据文件的格式为 T(1) , X(1) T(2) , X(2) T(3) , X(3) ……，…… T(512) , X(512) 其中T数组是正弦波采样点的时间值，X数组是正弦波采样值。 3.利用给定的频谱分析程序对信号进行分析。程序框图如下图 程序参数说明 M－FFT 的长度，应为2的幂次(64) IWIN－窗函数类型 IWIN=1，矩形窗
IWIN=2，汉明窗 L－窗长，L<=M(64) N－数据取样数(512) Fs－采样频率（一定要和对象截止频率对应） 三、实验步骤 1.调试自己编制的产生正弦波数据之和的程序，并将产生的数据放入数据文件中 2.运行频谱分析程序，画出正弦波信号的频谱图 3.改变PSDOLD程序中的M，L参数，看其对频谱的影响 四、实验结果及数据处理 1.产生正弦波数据之和程序见附件，令A1=20,A2=4,F1=60Hz,F2=200Hz,Fs=3000Hz得到的 波形如下图： 图表 1 正弦信号之和, A1=20,A2=4,F1=60Hz,F2=200Hz,Fs=3000Hz 2.频谱分析结果 图表 2 频谱分析结果F1=60Hz, F2=200Hz, Fs=3000Hz， N=512, M=256, IWIN=2, L=256
实验三 用FFT 对信号作频谱分析 姓名： 班级： 学号： 一、实验目的 学习用FFT 对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便正确应用FFT 。 二、实验原理与方法 用FFT 对信号作频谱分析是学习数字信号处理的重要内容。经常需要进行谱分析的信号是模拟信号和时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D 和分析误差。频谱分辨率直接和FFT 的变换区间N 有关，因为FFT 能够实现的频率分辨率是N /2π，因此要求D N ≤/2π。可以根据此式选择FFT 的变换区间N 。误差主要来自于用FFT 作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N 较大时离散谱的包络才能逼近于连续谱，因此N 要适当选择大一些。 周期信号的频谱是离散谱，只有用整数倍周期的长度作FFT ，得到的离散谱才能代表周期信号的频谱。如果不知道信号周期，可以尽量选择信号的观察时间长一些。 对模拟信号进行谱分析时，首先要按照采样定理将其变成时域离散信号。如果是模拟周期信号，也应该选取整数倍周期的长度，经过采样后形成周期序列，按照周期序列的谱分析进行。 三、实验内容及步骤 （1）对以下序列进行谱分析。 1423()() 1,03()8,470,4,03()3,47 0,x n R n n n x n n n n n n x n n n n =+≤≤?? =-≤≤??? -≤≤?? =-≤≤??? 其他其他 选择FFT 的变换区间N 为8和16 两种情况进行频谱分析。分别打印其幅频特性曲线。并进行对比、分析和讨论。 （2）对以下周期序列进行谱分析。
实验一、基本信号分析 一、实验目的 1.掌握基本信号的时域和频域分析方法 2.掌握信号的自相关和互相关分析，了解其应用 二、实验原理 （1）信号的时域和频域转换 目的：研究分析信号的时域特征（如持续时间、幅值、周期等）和信号的频域特征（如是否含有周期性信号、信号的频率带宽等） 转换方法：时域有限长序列 频域有限长序列： 离散傅里叶变换 （2）信号相关性 相关是用来描述一个随机过程自身在不同时刻的状态间，或者两个随机过程在某个时刻状态间线性依从关系的数字特征。 自相关函数定义为：xx 0 1()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、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法； 2、观察截短傅里叶级数而产生的“Gibbs 现象”，了解其特点以及产生的原因； 3、掌握各种典型的连续时间非周期信号的频谱特征。 二、实验容 实验前，必须首先阅读本实验原理，读懂所给出的全部例程序。实验开始时，先在计算机上运行这些例程序，观察所得到的信号的波形图。并结合例程序应该完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。 实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。 Q3-1 编写程序Q3_1，绘制下面的信号的波形图： -+-=)5cos(51 )3cos(31)cos()(000t t t t x ωωω∑∞ ==10)cos()2sin(1n t n n n ωπ 其中，ω0 = 0.5π，要求将一个图形窗口分割成四个子图，分别绘制cos(ω0t)、cos(3ω0t)、cos(5ω0t) 和x(t) 的波形图，给图形加title ，网格线和x 坐标标签，并且程序能够接受从键盘输入的和式中的项数。 抄写程序Q3_1如下： 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
信息工程学院实验报告 课程名称： 传感器原理及应用 实验项目名称： 实验四 声音传感器实验 实验时间： 班级： 姓名： 学号： 一、实 验 目 的 1. 学习 CC2530 单片机 GPIO 的使用。 2. 学习声音传感器的使用 二、实 验 原 理 1. CC2530 节点与三轴加速度传感器的硬件接口
(1). 声音传感器模块(MIC)引脚 GND：外接GND DO：数字量输出接口(0 和1) +5V：外接5V 电源 (2). 传感器模块与CC2530 模块之间的连接 2. GPIO (1). 简介 CC2530单片机具有21个数字输入/输出引脚，可以配置为通用数字I/O或外设I/O信号，配置为连接到ADC、定时器或USART外设。这些I/O口的用途可以通过一系列寄存器配置，由用户软件加以实现。 I/O端口具备如下特性： ●21个数字I/O引脚 ●可以配置为通用I/O或外部设备I/O ●输入口具备上拉或下拉能力 ●具有外部中断能力。 这21个I/O引脚都可以用作于外部中断源输入口。因此如果需要外部设备可以产生中断。外部中断功能也可以从睡眠模式唤醒设备。 (2). 寄存器简介 本次实验中主要涉及到GPIO的寄存器如下：
3. MIC 声音传感器 (1). 概述 声音传感器的作用相当于一个话筒（麦克风）。它用来接收声波，显示声音的振动图象。但不能对噪声的强度进行测量。 该传感器内置一个对声音敏感的电容式驻极体话筒。声波使话筒内的驻极体薄膜振动，导致电容的变化，而产生与之对应变化的微小电压。这一电压随后被转化成0-5V 的电压，经过比较器转换数字信号后，被数据采集器接受，并传送给计算机。 传感器特点： ●具有信号输出指示。 ●输出有效信号为低电平。 ●当有声音时输出低电平，信号灯亮。 应用范围： ●可以用于声控灯，配合光敏传感器做声光报警，以及声音控制，声音检测的场合。 (2). 使用方法 本实验利用CC2530 的GPIO 读取声音传感器模块的检测结果输出端，当检测到一定的声音时，此输出端为低电平；未检测到一定的声音时，此输出端为高电平。因此在实际应用中可以根据这种情况判断是否有声音在传感器附近产生。 4.程序流程