一、填空题（每空1分, 共10分）
1．序列()sin(3/5)x n n π=的周期为 10 。
2．线性时不变系统的性质有交换律、结合律、分配律。
3．因果序列x(n)，在Z →∞时，X(Z)=x(0)。
4．抽样序列的Z 变换与离散傅里叶变换DFT 的关系为k N j e Z π2=。
5．序列x(n)=(1，-2，0，3；n=0，1，2，3), 圆周左移2位得到的序列为{0，3，1，-2; n=0,1,2,3}。
6．设LTI 系统输入为x(n) ，系统单位序列响应为h(n)，则系统零状态输出()()()y n x n h n =*。
7．对4()()x n R n =的Z 变换为：，其收敛域为： 4
11,01z z z --->-
二、单项选择题（每题2分, 共20分）
1．δ(n)的Z 变换是
A.1
B.δ(ω)
C.2πδ(ω)
D.2π
2．序列x 1（n ）的长度为4，序列x 2（n ）的长度为3，则它们线性卷积的长度是
A. 3
B. 4
C. 6
D. 7
3．LTI 系统，输入x （n ）时，输出y （n ）；输入为3x （n-2），输出为
A. y （n-2）
B.3y （n-2）
C.3y （n ）
D.y （n ）
4．下面描述中最适合离散傅立叶变换DFT 的是
A.时域为离散序列，频域为连续信号
B.时域为离散周期序列，频域也为离散周期序列
C.时域为离散无限长序列，频域为连续周期信号
D.时域为离散有限长序列，频域也为离散有限长序列
5．若一模拟信号为带限，且对其抽样满足奈奎斯特条件，理想条件下将抽样信号通过 即可完
全不失真恢复原信号
A.理想低通滤波器
B.理想高通滤波器
C.理想带通滤波器
D.理想带阻滤波器
6．下列哪一个系统是因果系统
A.y(n)=x (n+2)
B. y(n)= cos(n+1)x (n)
C. y(n)=x (2n)
D.y(n)=x (- n)
7．一个线性时不变离散系统稳定的充要条件是其系统函数的收敛域包括
A. 实轴
B.原点
C.单位圆
D.虚轴
8．已知序列Z 变换的收敛域为｜z ｜>2，则该序列为
A.有限长序列
B.无限长序列
C.反因果序列
D.因果序列
9．若序列的长度为M ，要能够由频域抽样信号X(k)恢复原序列，而不发生时域混叠现象，则频
域抽样点数N 需满足的条件是
A.N≥M
B.N≤M
C.N≤2M
D.N≥2M
10．设因果稳定的LTI系统的单位抽样响应h(n)，在n<0时，h(n)=
A.0
B.∞
C. -∞
D.1
三、判断题（本题共10个小题，每小题1分，共10分）
1．序列的傅立叶变换是频率ω的周期函数，周期是2π。（1）2．x(n)= sin（ω0n)所代表的序列不一定是周期的。（1）3．FIR离散系统的系统函数是z的多项式形式。（1）4．y(n)=cos[x(n)]所代表的系统是非线性系统。（1）5．FIR滤波器较IIR滤波器的最大优点是可以方便地实现线性相位。（1）6．用双线性变换法设计IIR滤波器，模拟角频转换为数字角频是线性转换。（0）7．对正弦信号进行采样得到的正弦序列一定是周期序列。（0）8．常系数差分方程表示的系统为线性移不变系统。（0）9．FIR离散系统都具有严格的线性相位。（0）10．在时域对连续信号进行抽样，在频域中，所得频谱是原信号频谱的周期延拓。（0）四、简答题（每题5分，共20分）
1．用DFT对连续信号进行谱分析的误差问题有哪些？
答：混叠失真；截断效应（频谱泄漏）；栅栏效应
2．画出模拟信号数字化处理框图，并简要说明框图中每一部分的功能作用。
答：
第1部分：滤除模拟信号高频部分；第2部分：模拟信号经抽样变为离散信号；第3部分：按照预制要求对数字信号处理加工；第4部分：数字信号变为模拟信号；第5部分：滤除高频部分，平滑模拟信号。
3．简述用双线性法设计IIR数字低通滤波器设计的步骤。
答：确定数字滤波器的技术指标；将数字滤波器的技术指标转变成模拟滤波器的技术指标；按模拟滤波器的技术指标设计模拟低通滤波器；将模拟低通滤波器转换成数字低通滤波器。
4．8点序列的按时间抽取的（DIT）基-2 FFT如何表示？
五、计算题 （共40分）
1．已知2
(),2(1)(2)
z X z z z z =>+-，求x(n)。（6分）
1．解：由题部分分式展开 ()(1)(2)12
F z z A B z z z z z ==++-+- 求系数得 A=1/3 ， B=2/3
所以 2
32131)(-++=z z z z z F 收敛域?z ?>2，故上式第一项为因果序列象函数，第二项为反因果序列象函数， 则 12()(1)()(2)()33
k k f k k k εε=
-+
2．写出差分方程表示系统的直接型和级联．．型结构。（8分） )1(3
1)()2(81)1(43)(-+=-+--n x n x n y n y n y
3．计算下面序列的N 点DFT 。
（1）)0()
()(N m m n n x <<-=δ（4分） （2）)
0()(2N m e n x mn N j <<=π （4分） 解：（1） kn N W k X =)(
（2）???≠==m
k m k N k X ,0,)( 4．设序列x(n)={1，3，2，1；n=0,1,2,3 }，另一序列h(n) ={1，2，1，2；n=0,1,2,3}，
（1）求两序列的线性卷积 y L (n)； （4分）
（2）求两序列的6点循环卷积y C (n)。 （4分）
（3）说明循环卷积能代替线性卷积的条件。（2分）
解：（1） y L (n)={1，5，9，10，10，5，2；n=0,1,2…6}
（2） y C (n)= {3，5，9，10，10，5；n=0,1,2,4,5}
（3）c ≥L 1+L 2-1
5．设系统由下面差分方程描述：
)1()2()1()(--+-=n x n y n y n y
（1）求系统函数H （z ）；（2分）
（2）限定系统稳定．．
，写出H （z ）的收敛域，并求出其单位脉冲响应h(n)。（6分） 解：（1） 1
)(2--=z z z z H （2）511522
z -+<< )1()251(51)()251(51)(--+---
=n u n u n h n n
实验一 离散时间系统的时域分析 一、实验目的 １. 运用MATLAB 仿真一些简单的离散时间系统，并研究它们的时域特性。 ２. 运用MATLAB 中的卷积运算计算系统的输出序列，加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。 二、实验原理 离散时间系统其输入、输出关系可用以下差分方程描述： ∑=∑=-=-M k k N k k k n x p k n y d 00] [][ 当输入信号为冲激信号时，系统的输出记为系统单位冲激响应 ][][n h n →δ，则系统响应为如下的卷积计算式： ∑∞ -∞=-= *=m m n h m x n h n x n y ][][][][][ 当h[n]是有限长度的（n ：[0，M]）时，称系统为FIR 系统；反之，称系统为IIR 系统。在MA TLAB 中，可以用函数y=Filter(p,d,x) 求解差分方程，也可以用函数 y=Conv(x,h)计算卷积。 例1 clf; n=0:40; a=1;b=2; x1= 0.1*n; x2=sin(2*pi*n); x=a*x1+b*x2; num=[1, 0.5,3]; den=[2 -3 0.1]; ic=[0 0]; %设置零初始条件 y1=filter(num,den,x1,ic); %计算输入为x1(n)时的输出y1(n) y2=filter(num,den,x2,ic); %计算输入为x2(n)时的输出y2(n) y=filter(num,den,x,ic); %计算输入为x (n)时的输出y(n) yt= a*y1+b*y2; %画出输出信号 subplot(2,1,1) stem(n,y); ylabel(‘振幅’)； title(‘加权输入a*x1+b*x2的输出’)；
测控1005班齐伟0121004931725 （18号）实验一差分方程、卷积、z变换 一、实验目的 通过该实验熟悉 matlab软件的基本操作指令，掌握matlab软件的使用方法，掌握数字信号处理中的基本原理、方法以及matlab函数的调用。 二、实验设备 1、微型计算机1台； 2、matlab软件1套 三、实验原理 Matlab 软件是由mathworks公司于1984年推出的一套科学计算软件，分为总包和若干个工具箱，其中包含用于信号分析与处理的sptool工具箱和用于滤波器设计的fdatool工具箱。它具有强大的矩阵计算和数据可视化能力，是广泛应用于信号分析与处理中的功能强大且使用简单方便的成熟软件。Matlab软件中已有大量的关于数字信号处理的运算函数可供调用，本实验主要是针对数字信号处理中的差分方程、卷积、z变换等基本运算的matlab函数的熟悉和应用。 差分方程（difference equation)可用来描述线性时不变、因果数字滤波器。用x表示滤波器的输入，用y表示滤波器的输出。 a0y[n]+a1y[n-1]+…+a N y[n-N]=b0x[n]+b1x[n-1]+…+b M x[n-M] (1) ak,bk 为权系数，称为滤波器系数。 N为所需过去输出的个数，M 为所需输入的个数卷积是滤波器另一种实现方法。 y[n]= ∑x[k] h[n-k] = x[n]*h[n] (2) 等式定义了数字卷积，*是卷积运算符。输出y[n] 取决于输入x[n] 和系统的脉冲响应h[n]。 传输函数H(z)是滤波器的第三种实现方法。 H(z)=输出/输入= Y(z)/X(z) （3）即分别对滤波器的输入和输出信号求z变换，二者的比值就是数字滤波器的传输函数。 序列x[n]的z变换定义为 X (z)=∑x[n]z-n (4) 把序列x[n] 的z 变换记为Z{x[n]} = X(z)。
测试技术与信号处理实验报告 机械转子底座的振动测量和分析 一、实验目的 1．掌握磁电式速度传感器的工作原理、特点和应用。
2．掌握振动的测量和数据分析。 二、实验内容和要求 先利用光电式转速传感器测量出电机的转速；然后利用磁电式速度传感器测量机械转子底座在该电机转速下的振动速度；对测量出的振动速度信号进行频谱分析；找出振动信号的主频与电机转速之间的关系。 三、实验步骤 1．启动实验程序“机械转子系统的振动测量.exe”; 输入个人信息，也可以启动之后通过单击“修改”按钮修改个人信息。 2．单击“采样设置”按钮，输入采集卡连接磁电速度传感器的采样通道号，批量采样频率（建议设为10KHz）、批量采样点数（建议设为10000）。 3．打开转子电机的电源，单击“单点采样”。 4．旋转调节旋钮改变转子的转速,观察图形区显示的磁电速度传感器采集到的转子底座振动信号；如果振动信号比较小，可适当提高转子的转速。 5．转子转速的测量： (1) 单击“采样设置”按钮，输入采集卡连接光电转速传感器的 采样通道号、批量采样频率（建议值为10KHz）、批量采样点 数（建议值为10000）。 (2) 单击“批量采样”按钮，开始采样；采样完成之后，采集到 的波形信号会显示在图形窗口，系统会自动计算出转子的速度
并显示出来。记录下此时的转子的转速（单位：r/s）。 (3) 再重复步骤(2)测量2次。以三次测量的平均值作为此时转子 的转速。 转速的测量结果 单点采样采集通道6,测量3组数据 6．振动信号的测量和频谱分析： (1) 单击“采样设置”按钮，输入采集卡连接磁电速度传感器的 采样通道号、批量采样频率（建议设为10KHz）、批量采样点 数（建议设为10000）。 (2) 单击“批量采样”按钮，开始采样；采样完成之后，采集到 的波形信号会显示在图形窗口。如果信号不正常，重复点击“批 量采样”按钮 (3) 单击“保存”按钮，将采集到的磁电传感器的信号数据保存 为文本文件。文件必须保存到“C:\ExperiData\”目录下。可单 击“保存设置”更改文件名。 (4) 打开刚保存的文本文件，文件前面几行保存了个人信息、采 样频率、采样通道、保存的数据个数等信息。文件中共有四列 数据，第一列为数据的序号，第二列为磁电传感器检测到的数 据。
2.1 有一个理想采样系统，其采样角频率Ωs =6π，采样后经理想低通滤波器H a (j Ω)还原，其中 ?? ???≥Ω<Ω=Ωππ 3032 1 )(，，j H a 现有两个输入，x 1(t )=cos2πt ，x 2(t )=cos5πt 。试问输出信号y 1(t )，y 2(t )有无失真？为什么？ 分析：要想时域采样后能不失真地还原出原信号，则采样角频率Ωs 必须大于等于信号谱最高角频率Ωh 的2倍，即满足Ωs ≥2Ωh 。 解：已知采样角频率Ωs =6π，则由香农采样定理，可得 因为x 1(t )=cos2πt ，而频谱中最高角频率ππ π32621=< =Ωh ，所以y 1(t )无失真； 因为x 2(t )=cos5πt ，而频谱中最高角频率ππ π32 652=>=Ωh ，所以y 2(t )失真。 3.2 设x (n )的傅里叶变换为X (e j ω)，试利用X (e j ω )表示下列序列的傅里叶变换： （1） )1()1()(1n x n x n x --+-= （2） )]()([2 1 )(2n x n x n x -+= * 分析：利用序列翻褶后的时移性质和线性性质来求解，即 )()(ωj e X n x ?，)()(ωj e X n x -?- )()(ωωj m j e X e n m x --?- 解：（1）由于)()]([ω j e X n x DTFT =，)()]([ωj e X n x DTFT -=-，则 )()]1([ωωj j e X e n x DTFT --=- )()]1([ωωj j e X e n x DTFT -=-- 故ωωωωω cos )(2])[()]([1j j j j e X e e e X n x DTFT ---=+= （2）由于)()]([ω j e X n x DTFT * * =- 故)](Re[2 ) ()()]([2ωωωj j j e X e X e X n x DTFT =+= * 3.7 试求下列有限长序列的N 点离散傅里叶变换（闭合形式表达式）：
第二章 信号的描述与分析 补充题2-1-1 求正弦信号0()sin()x t x ωt φ=+的均值x μ、均方值2 x ψ和概率密度函数 p (x )。 解答： (1)0 00 11lim ()d sin()d 0T T x T μx t t x ωt φt T T →∞== +=? ? ，式中02π T ω = —正弦信号周期 (2) 2 222 2 2 0000 1 1 1cos 2() lim ()d sin ()d d 22 T T T x T x x ωt φψx t t x ωt φt t T T T →∞-+== += = ? ? ? (3)在一个周期内 012ΔΔ2Δx T t t t =+= 000 2Δ[()Δ]lim x x T T T t P x x t x x T T T →∞<≤+=== Δ0Δ000 [()Δ]2Δ2d ()lim lim ΔΔd x x P x x t x x t t p x x T x T x →→<≤+==== 正弦信号 x
2-8 求余弦信号0()sin x t x ωt 的绝对均值x μ和均方根值rms x 。 2-1 求图示2.36所示锯齿波信号的傅里叶级数展开。
2-4周期性三角波信号如图2.37所示，求信号的直流分量、基波有效值、信号有效值及信号的平均功率。
2-1 求图示2.36所示锯齿波信号的傅里叶级数展开。 补充题2-1-2 求周期方波（见图1-4）的傅里叶级数（复指数函数形式），划出|c n|–ω和φn–ω
图，并与表1-1对比。 解答：在一个周期的表达式为 00 (0)2 () (0) 2 T A t x t T A t ? --≤
《测试信号分析与处理》 实验一差分方程、卷积、z变换 一、实验目的 通过该实验熟悉 matlab软件的基本操作指令，掌握matlab软件的使用方法，掌握数字信号处理中的基本原理、方法以及matlab函数的调用。 二、实验设备 1、微型计算机1台； 2、matlab软件1套 三、实验原理 Matlab 软件是由mathworks公司于1984年推出的一套科学计算软件，分为总包和若干个工具箱，其中包含用于信号分析与处理的sptool工具箱和用于滤波器设计的fdatool工具箱。它具有强大的矩阵计算和数据可视化能力，是广泛应用于信号分析与处理中的功能强大且使用简单方便的成熟软件。Matlab软件中已有大量的关于数字信号处理的运算函数可供调用，本实验主要是针对数字信号处理中的差分方程、卷积、z变换等基本运算的matlab函数的熟悉和应用。 差分方程（difference equation)可用来描述线性时不变、因果数字滤波器。用x表示滤波器的输入，用y表示滤波器的输出。 a0y[n]+a1y[n-1]+…+a N y[n-N]=b0x[n]+b1x[n-1]+…+b M x[n-M] (1) ak,bk 为权系数，称为滤波器系数。
N为所需过去输出的个数，M 为所需输入的个数卷积是滤波器另一种实现方法。 y[n]= ∑x[k] h[n-k] = x[n]*h[n] (2)等式定义了数字卷积，*是卷积运算符。输出y[n] 取决于输入x[n] 和系统的脉冲响应h[n]。 传输函数H(z)是滤波器的第三种实现方法。 H(z)=输出/输入= Y(z)/X(z) （3） 即分别对滤波器的输入和输出信号求z变换，二者的比值就是数字滤波器的传输函数。 序列x[n]的z变换定义为 X (z)=∑x[n]z-n (4)把序列x[n] 的z 变换记为Z{x[n]} = X(z)。 由X(z) 计算x[n] 进行z 的逆变换x[n] = Z-1{X(z)}。 Z 变换是Z-1的幂级数，只有当此级数收敛，Z 变换才有意义，而且同一个Z 变换等式，收敛域不同，可以代表不同序列的Z 变换函数。 这三种数字滤波器的表示方法之间可以进行相互转换。 四、实验步骤 1、熟悉matlab软件基本操作指令。读懂下列matlab程序指令，键入程序并 运行，观察运行结果。 Conv.m% 计算两个序列的线性卷积； %-----------------------------------------------------------------
《测试技术与信号处理》习题答案 第二章 信号分析基础 1、请判断下列信号是功率信号还是能量信号： （1）)()(10cos 2 ∞<<-∞=t e t x t π （2）)()(||10∞<<-∞=-t e t x t 【解】（1）该信号为周期信号，其能量无穷大，但一个周期内的平均功率有限，属功率信号。 （2）信号能量：? ∞ ∞ -= =10 1 )(2dt t x E ，属于能量信号。 2、请判断下列序列是否具有周期性，若是周期性的，请求其周期。)8 ()(π-=n j e n x 【解】设周期为N ，则有：8 )8 8()()(N j N n j e n x e N n x ?==+-+π 若满足)()(n x N n x =+，则有1)8/sin()8/cos(8/=-=-N j N e jN 即：k N π28/=，k N π16=，k = 0，1，2，3，… N 不是有理数，故序列不是周期性的。 3、已知矩形单脉冲信号x 0(t)的频谱为X 0(ω)=A τsinc(ωτ/2) ，试求图示三脉冲信号的频谱。 【解】三脉冲信号的时域表达式为：)()()()(000T t x t x T t x t x -+++= 根据Fourier 变换的时移特性和叠加特性，可得其频谱： )]cos(21)[2 ( sin )()()()(000T c A e X X e X X T j T j ωωτ τωωωωωω+=++=- 4、请求周期性三角波（周期为T ，幅值为0—A ）的概率分布函数F(x)与概率密度函数p(x) 。 【解】在一个周期T 内，变量x (t )小于某一特定值x 的时间间隔平均值为：T A x t i = ? 取n 个周期计算平均值，当∞→n 时，可有概率分布函数：A x nT t n x F i n =?=∞→lim )( 概率密度函数：A dx x dF x p 1 )()(== t -τ/2 0 τ/2 -T T
信号分析与处理课后习题答案 第五章快速傅里叶变换 1.如果一台通用计算机的速度为平均每次复乘需要50us ，每次复加需要10us ，用来就散N=1024点的DFT ，问： （1）直接计算需要多少时间？用FFT 计算呢？ （2）照这样计算，用FFT 计算快速卷积对信号进行处理是，估计可实现实时处理的信号最高频率？ 解： 分析：直接利用DFT 计算：复乘次数为N 2，复加次数为N(N-1)； 利用FFT 计算：复乘次数为20.5log N N ，复加次数为2log N N ； （1） 直接DFT 计算： 复乘所需时间2215010245052.4288T N us us s =?=?= 复加所需时间2(1)101024(10241)1010.47552T N N us us s =-?=-?= 所以总时间1262.90432DFT T T T s =+= FFT 计算： 复乘所需时间3220.5log 500.51024log 1024500.256T N N us us s =?=???= 复加所需时间422log 101024log 1024100.1024T N N us us s =?=??= 所以总时间为340.3584FFT T T T s =+= （2） 假设计算两个N 长序列1()x n 和2()x n 的卷积 计算过程为如下： 第一步：求1()X k ，2()X k ；所需时间为2FFT T ? 第二步：计算12()()()X k X k X k =?，共需要N 次复乘运算 所需时间为501024500.0512To N us us s =?=?= 第三步：计算(())IFFT X k ，所需时间为FFT T 所以总时间为230.35840.0512 1.1264FFT T T To s s s =?+=?+= 容许计算信号频率为N/T=911.3Hz 2.设x(n)是长度为2N 的有限长实序列，()X k 为x(n)的2N 点得DFT 。
机械工程测试技术基础(第三版)试卷集. 一、填空题 １、周期信号的频谱是离散的，而非周期信号的频谱是的。 ２、均方值Ψx2表示的是信号的强度，它与均值μx、方差σx2的关系是。 ３、测试信号调理电路主要有、、。 ４、测试系统的静态特性指标有、、。 ５、灵敏度表示系统输出与输入之间的比值，是定度曲线的。 ６、传感器按信号变换特性可分为、。 ７、当时，可变磁阻式电感传感器的输出和输入成近似线性关系，其灵敏度S趋于。 ８、和差特性的主要内容是相临、相反两臂间阻值的变化量符合、的变化，才能使输出有最大值。 ９、信号分析的过程主要包括：、。 10、系统动态特性在时域可用来描述，在复数域可用来描述，在频域可用来描述。 11、高输入阻抗测量放大电路具有高的共模抑制比，即对共模信号有抑制作用，对信号有放大作用。 12、动态应变仪上同时设有电阻和电容平衡旋钮，原因是导线间存在。 13、压控振荡器的输出电压是方波信号，其与输入的控制电压成线性关系。 14、调频波的解调又称，其解调电路称为。 15、滤波器的通频带宽和响应时间成关系。 16、滤波器的频率分辨力主要由其决定。 17、对于理想滤波器，滤波器因数λ＝。 18、带通滤波器可由低通滤波器（f c2）和高通滤波器（f c1）而成（f c2> f c1）。 19、测试系统的线性度和滞后度是由误差引起的；而重复性误差是 由误差引起的。 二、问答题（共30分） １、什么是测试？说明测试系统的构成及各组成部分的作用。（10分） ２、说明电阻丝应变片和半导体应变片的异同点，各有何优点？（10分） ３、选用传感器的原则是什么？（10分） 三、计算题（共55分） １、已知信号x(t)=e-t (t≥0)， (1) 求x(t)的频谱函数X(f)，并绘制幅频谱、相频谱。 (2) 求x(t)的自相关函数R x (τ) 。（15分） ２、二阶系统的阻尼比ξ=0.2，求ω=ωn时的幅值误差和相位误差，如果使幅值误差不大于10％，应取多大阻尼比？。（10分）３、一电容传感器，其圆形极板r = 4mm，工作初始间隙δ0 =0.3mm， （1）工作时如果传感器的工作间隙变化Δδ=±2μm，求电容的变化量。 （2）如果测量电路灵敏度S1=100mv/pF，读数仪表灵敏度S2=5格/mv，在 Δδ=±2μm时，读数仪表的指示值变化多少格？ （ε0 = 8.85×10-12 F/m）(8分) ４、已知RC低通滤波器的R=1KΩ，C=1MF，当输入信号μx= 100sin1000t时， 求输出信号μy 。（7分） ５、（1）在下图中写出动态应变仪所包含的各个电路环节。 （2）如被测量x(t) = sinωt，载波y(t)=sin6ωt，画出各环节信号的波形图。(15分 一、填空题： 1、连续 2、￠x2=H x2+óx2 3、电桥、放大、调制解调电路 4、非线性度、灵敏度、回程误差 5、斜率 6、组合型、一体化型 7、Δó〈〈ó0定位8、相邻相反相对相同9、信号分析、信号处理 10、传递函数、频率函数、脉冲响应函数11、差模12、分布电容13、频率14、鉴频、鉴频器15、反比16、带宽B 17、1 18、串联19、
第七章信号处理与分析 ６．１概述 数字信号在我们周围无所不在。因为数字信号具有高保真、低噪声和便于信号处理的优点，所以得到了广泛的应用，例如电话公司使用数字信号传输语音，广播、电视和高保真音响系统也都在逐渐数字化。太空中的卫星将测得数据以数字信号的形式发送到地面接收站。对遥远星球和外部空间拍摄的照片也是采用数字方法处理，去除干扰，获得有用的信息。经济数据、人口普查结果、股票市场价格都可以采用数字信号的形式获得。因为数字信号处理具有这么多优点，在用计算机对模拟信号进行处理之前也常把它们先转换成数字信号。本章将介绍数字信号处理的基本知识，并介绍由上百个数字信号处理和分析的VI构成的LabVIEW分析软件库。 目前，对于实时分析系统，高速浮点运算和数字信号处理已经变得越来越重要。这些系统被广泛应用到生物医学数据处理、语音识别、数字音频和图像处理等各种领域。数据分析的重要性在于，无法从刚刚采集的数据立刻得到有用的信息，如下图所示。必须消除噪音干扰、纠正设备故障而破坏的数据，或者补偿环境影响，如温度和湿度等。 通过分析和处理数字信号，可以从噪声中分离出有用的信 息，并用比原始数据更全面的表格显示这些信息。下图显示的是 经过处理的数据曲线。
用于测量的虚拟仪器(VI) 用于测量的虚拟仪器(VI)执行的典型的测量任务有： ●计算信号中存在的总的谐波失真。 ●决定系统的脉冲响应或传递函数。 ●估计系统的动态响应参数，例如上升时间、超调量等等。 ●计算信号的幅频特性和相频特性。 ●估计信号中含有的交流成分和直流成分。 在过去，这些计算工作需要通过特定的实验工作台来进行，而用于测量的虚拟仪器可以使这些测量工作通过LabVIEW程序语言在台式机上进行。这些用于测量的虚拟仪器是建立在数据采集和数字信号处理的基础之上，有如下的特性： ●输入的时域信号被假定为实数值。 ●输出数据中包含大小、相位，并且用合适的单位进行了刻度，可用来直接进行 图形的绘制。 ●计算出来的频谱是单边的（single_sided），范围从直流分量到Nyquist频率(二 分之一取样频率)。（即没有负频率出现） ●需要时可以使用窗函数，窗是经过刻度地，因此每个窗提供相同的频谱幅度峰 值，可以精确地限制信号的幅值。 一般情况下，可以将数据采集VI的输出直接连接到测量VI的输入端。测量VI的输出又可以连接到绘图VI以得到可视的显示。 有些测量VI用来进行时域到频域的转换，例如计算幅频特性和相频特性、功率谱、网路的传递函数等等。另一些测量VI可以刻度时域窗和对功率和频率进行估算。 本章我们将介绍测量VI中常用的一些数字信号处理函数。 LabVIEW的流程图编程方法和分析VI库的扩展工具箱使得分析软件的开发变得更加简单。LabVIEW 分析VI通过一些可以互相连接的VI，提供了最先进的数据分析技术。你不必像在普通编程语言中那样关心分析步骤的具体细节，而可以集中注意力解决信号处理与分析方面的问题。LabVIEW 6i版本中，有两个子模板涉及信号处理和数学，分别是Analyze 子模板和Methematics子模板。这里主要涉及前者。 进入Functions模板Analyze》Signal Processing子模板。 其中共有6个分析VI库。其中包括： ①．Signal Generation（信号发生）：用于产生数字特性曲线和波形。 ②．Time Domain（时域分析）：用于进行频域转换、频域分析等。 ③．Frequency Domain（频域分析）： ④．Measurement（测量函数）：用于执行各种测量功能，例如单边FFT、频谱、比例加窗以及泄漏频谱、能量的估算。
信号分析与处理 第一章绪论：测试信号分析与处理的主要内容、应用；信号的分类，信号分析与信号处理、测试信号的描述，信号与系统。 测试技术的目的是信息获取、处理和利用。 测试过程是针对被测对象的特点，利用相应传感器，将被测物理量转变为电信号，然后，按一定的目的对信号进行分析和处理，从而探明被测对象内在规律的过程。 信号分析与处理是测试技术的重要研究内容。 信号分析与处理技术可以分成模拟信号分析与处理和数字信号分析与处理技术。 一切物体运动和状态的变化，都是一种信号，传递不同的信息。 信号常常表示为时间的函数，函数表示和图形表示信号。 信号是信息的载体，但信号不是信息，只有对信号进行分析和处理后，才能从信号中提取信息。 信号可以分为确定信号与随机信号；周期信号与非周期信号；连续时间信号与离散时间信号；能量信号与功率信号；奇异信号； 周期信号无穷的含义，连续信号、模拟信号、量化信号，抽样信号、数字信号 在频域里进行信号的频谱分析是信号分析中一种最基本的方法：将频率作为信号的自变量，在频域里进行信号的频谱分析； 信号分析是研究信号本身的特征，信号处理是对信号进行某种运算。 信号处理包括时域处理和频域处理。时域处理中最典型的是波形分析，滤波是信号分析中的重要研究内容； 测试信号是指被测对象的运动或状态信息，表示测试信号可以用数学表达式、图形、图表等进行描述。 常用基本信号（函数）复指数信号、抽样函数、单位阶跃函数单位、冲激函数（抽样特性和偶函数）离散序列用图形、数列表示，常见序列单位抽样序列、单位阶跃序列、斜变序列、正弦序列、复指数序列。 系统是指由一些相互联系、相互制约的事物组成的具有某种功能的整体。被测系统和测试系统统称为系统。输入信号和输出信号统称为测试信号。系统分为连续时间系统和离散时间系统。
《信号分析与处理》(第二版)-徐科军、黄云志-课后答案
————————————————————————————————作者：————————————————————————————————日期：
Chap1. 1.4 ()()()()()()()()()()()() ()()()()()()()121 2 122 12112 2 121 2 2 2y 11102 y 0.5111 y 0.5 1.513y 0 13 013 y 0.5111 0.5 1.513t t t t t x t x t x x t d x x t x x t d t d t t t x x t d t d t t t t t or t t or t t t t t t t τττ ττττ τττττττττττ+∞ -∞ ----=*=-=-≤≤???=≤≤??=-= -=+-<≤=-= -=-++<<=≤-≥≤-≥??=+-<≤??-++<
()()[] ()()()[]()()()∑∞ =? ? ? ???Ω-Ω-+=- =-= =??? ??<≤<≤-=1002212 2 01cos cos cos 1cos 141cos 1cos 1 5 .0202 20 (a)n n n t n n n t n n n t x n n b n n a a T t t T t T t x πππππ πππ 代入公式得： ()() ()()() ()[] ()()[]()()∑∞ =Ω-? ? ? ???Ω-Ω-+=- =-= ==Ω=Ω-=1002222 2 012 212cos 1cos cos 11411cos 11 5.0cos 2 (b)n n n T jn t n n t n n n t x n b n n a a n n X e n X T t x t x πππππππ得到：根据时移性质： ()() ()()()[]()()[]() ∑?∑∞ =-∞ =Ω-+=-=Ω==Ω+=102232 20 2 0201 00 3cos cos 12 21cos 12cos 41 cos 2 (c)n T n n n t n n n t x n n dt t n t x T a a t n a a t x ππ ππ偶对称， 1.12 ()()dt e t x j X t j ?+∞ ∞ -Ω-=Ω频谱密度函数：
结课作业 课程名称测试信号处理与分析学生专业测控技术与仪器 学生学号912101170116 学生姓名陈昊飞 任课教师吴健 成绩
一、（20分）用标准数字电压表在标准条件下，对 被测的10 V 直流电压信号进行了10次独立测量，测量值如表1所列。由该数字电压表的检定证书给出，其示值误差按3倍标准差计算为3.5×10-6V 。同时在进行电压测量前，对数字电压表进行了24h 的校准，在10 V 点测量时，24h 的示值稳定度不超过士15μV 。试分析评定对该10V 直流电压的测量结果。 答：此次测量为静态测量，只考虑静态误差，不涉及动态误差。 在不考虑系统误差的情况下，对此10次测量进行标准不确定度的A 类评定，其平均值0001043.10_ =x ，其标准差 6 10982.8-?=σ，平均值的实验标准差6_ 1084.2)(-?=x s ，单次实验的测量结果表示为 )]([_ _ x s x ±，为61084.20001043.10-?±。 根据示值误差的判定应用σ3准则，不含粗大误差的测量值范围为（10.000077~10.000131），判断此次测量不含有粗大误差。 实际值=测量值-示值误差，所以实际值为10.0001043-3.56 10-?=10.0001008，修正后的 结果为6 1084.20001008 .10-?±。 15μV=156 6 1084.210--?>?V ，测量A 类不确定度没有超过示值稳定度，其结果是可靠的。 综上所述，最终的结果为6 1084.20001008 .10-?±。 二、（20分）测量某半导体的两参量x 和y 所得数据如表2所示。试分析x ， y 之间的关系。（要求给出详细分析过程和MATlab 源程序） 答：在未对x ，y 做任何处理时对（xi ，yi ）做多项式拟合，参考书50页程序得到： MATLAB 程序如下： clear
2013—2014学年第二学期 《测试信号处理与分析》 实训报告 学院：机械与汽车工程学院 专业：测控技术与仪器 班级：11级测控二班 姓名： 学号： 指导教师
【摘要】：现代信号分析处理技术发展的非常迅速，各种信号专业处理软件也出现在了人们的视野中，这些软件给人们带来了极大方便。本次实训，我们主要学习了INV1612型软件以及DASP信号分析处理系统，切实的感受到了方便。本次实训为期两周，包括在实验室做简支梁的振动信号测试及柔性转子的共振试验，在创新实验室测铣床的振动信号，和在圆楼三楼的AutoCAD机房对铣床的振动信号进行分析等等。在此过程中要基本掌握简支梁的震动信号的测试方法和数据分析，INV1612型多功能柔性转子测试系统、INV1601型振动与控制教学实验系统及MATLAB软件对信号的采集和处理的方法，同时在试验中遇到的问题及我们一起解决的过程。在本次实训中，充分要求了动手能力，实训中的每一项数据都要求自己动手去采集处理，从中我学会了很多知识与方法。 【关键词】：测试信号实训软件知识方法 一、简支梁 1、简支梁的概念 一种简易的支架，包括两个在一平面上可交叉扣合的条形支架，所述每个条形支架的两端为一端高一端低的结构，所述低的一端为钩状结构，钩状体与条形支架主体之间可伸缩的连接，使得每个支架针对不同大小的支撑物在长度方向上可调，且本支架结构简单，节省材料。简支梁就是承载两端竖向荷载，而不提供扭矩的支撑结构。只有两端支撑在柱子上的梁，主要承受正弯矩，一般为静定结构。体系温变、混凝土收缩徐变、张拉预应力、支座移动等都不会在梁中产生附加内力，受力简单，简支梁为力学简化模型。对于简支梁来说，梁的两端搭在两个支撑物上，两端铰接，现实看是只有两端支撑在柱子上的梁，主要承受弯距的单跨结构.一般为静定结构。 2、用“双踪示波比较法”测量简谐振动的频率 用“双踪示波比较法”测量简谐振动的频率，实验仪器有：INV1601B 型振动教学实验仪、INV1601T 型振动教学实验台、速度传感器、调速电机及调压器。
1．在系统特性测量中常用白噪声信号作为输入信号，然后测量系统的输出，并将输出信号的频谱作为系统频率特性。请用卷积分定理解释这样做的道理。 答：白噪声是指功率谱密度在整个频域内均匀分布的噪声，所有频率具有相同能量的随机噪声称为白噪声。在其频谱上是一条直线。系统频率特性：传递函数的一种特殊情况，是定义在复平面虚轴上的传递函数。时域卷积分定理：两个时间函数的卷积的频谱等于各个时间函数的乘积，即在时域中两信号的卷积等效于在频域中频谱相乘。频域卷积分定理：两个时间函数的频谱的卷积等效于时域中两个时间函数的乘积。y(t)=h(t)*x(t),对y(t)作付式变换，转到相应的频域下Y(f)=H(f)X(f)，由于x(t)是白噪声，付式变换转到频域下为一定值，假定X(f)=1，则有Y(f)=H(f)，此时就是传递函数。 2．用1000Hz的采样频率对200Hz的正弦信号和周期三角波信号进行采样，请问两个信号采样后是否产生混叠？为什么？ 采样频率ωs（2π／Ts）或fs(1／Ts)必须大于或等于信号x(t)中的最高频率ωm的两倍，即ωs＞2ωm，或fs＞2fm。 为了保证采样后的信号能真实地保留原始模拟信号的信息，采样信号的频率必须至少为原信号中最高频率成分的2倍。这是采样的基本法则，称为采样定理。 但在对信号进行采样时，满足了采样定理，只能保证不发生频率混叠，对信号的频谱作逆傅立叶变换时，可以完全变换为原时域采样信号，而不能保证此时的采样信号能真实地反映原信号。工程实际中采样频率通常大于信号中最高频率成分的3到5倍。 理论上周期三角波的频谱里包含所有奇次谐波分量，也就是说200Hz的周期三角波信号包含600Hz、1kHz、1.4kHz等等谐波，所以用1000Hz采样频率对200Hz周期三角波信号采样，会发生混叠。而对200Hz正弦信号采样不会发生混叠。 3．什么是能量泄露和栅栏效应？能量泄漏与栅栏效应之间有何关系？ 能量泄漏：将截断信号的谱XT(ω)与原始信号的谱X（ω）相比较可知，它已不是原来的两条谱线，而是两段振荡的连续谱．这表明原来的信号被截断以后，其频谱发生了畸变，原来集中在f0处的能量被分散到两个较宽的频带中去了，这种现象称之为频谱能量泄漏（Leakage）。 栅栏效应：对采样信号的频谱,为提高计算效率,通常采用FFT算法进行计算,设数据点数为N = T/dt = T.fs则计算得到的离散频率点为Xs(fi) , fi = i.fs/N , i = 0,1,2,…,N/2。这就相当于透过栅栏观赏风景,只能看到频谱的一部分,而其它频率点看不见,因此很可能使一部分有用的频率成分被漏掉,此种现象被称为栅栏效应。 频谱的离散取样造成了栅栏效应，谱峰越尖锐，产生误差的可能性就越大。例如，余弦信号的频谱为线谱。当信号频率与频谱离散取样点不等时，栅栏效应的误差为无穷大。 实际应用中，由于信号截断的原因，产生了能量泄漏，即使信号频率与频谱离散取样点不相等，也能得到该频率分量的一个近似值。从这个意义上说，能量泄漏误差不完全是有害的。如果没有信号截断产生的能量泄漏，频谱离散取样造成的栅栏效应误差将是不能接受的。 能量泄漏分主瓣泄漏和旁瓣泄漏，主瓣泄漏可以减小因栅栏效应带来的谱峰幅值估计误差，有其好的一面，而旁瓣泄漏则是完全有害的。 4．简述传递函数、频响函数和脉冲响应函数间的联系与区别。 传递函数：零初始条件下线性系统响应（即输出）量的拉普拉斯变化（或z变换）与激励（即输入）量的拉普拉斯变换之比。记作G（s）＝Y（s）／U（s），其中Y（s）、U（s)分别为输出量和输入量的拉普拉斯变换。 频响函数：(1)简谐激励时，稳态输出相量与输入相量之比。(2)瞬态激励时，输出的傅里叶变换与输入的傅里叶变换之比。(3)平稳随机激励时，输出和输入的互谱与输入的自谱之比。
数据信息采集与处理
基本内容：基于FFT的功率谱分析程序设计与应用 1.基本要求 1)对一个人为产生的信号进行采用FFT变换方法进行功率谱分析。 已知信号x(n)=80.0*COS(2*3.14*SF*n/FS) 式中： n=0,1,2 ……N-1 SF---信号频率 FS---采样频率 其FFT变换结果X(k)可用下面提供的FFT子程序求出，计算功率谱的公式为： W（k）=2(XR(k)2 +XI(k)2)/N 式中：k=0,1,2 ……N/2-1 XR(k)--- X(k)的实部 XI(k)--- X(k)的虚部 请用VB，VC或C++Builder编译器编程，或采用MATLAB计算，或采用高级语言调用MATLAB计算。处理结果为采用窗口显示时域波形和频域波形。 此信号的时域谱、频域谱、功率谱如下面图1~图3所示： 图1
图2 图3 其MATLAB代码为： FS=200; SF=10;
N=1024; n=0:N-1; t=n/FS; x=80.0*cos(2*3.14*SF*t); figure; plot(t，x); xlabel('t'); ylabel('y'); title('x=80.0*cos(2*3.14*SF*t)时域波形'); grid; y=fft(x,N); mag=abs(y); f=(0:length(y)-1)*FS/length(y);%进行对应的频率转换 figure; plot(f(1:N/2),mag(1:N/2));%做频谱图 xlabel('频率(Hz)'); ylabel('幅值'); title('x=80.0*cos(2*3.14*SF*t)幅频谱图N=1024'); grid; Py =2*(y.*conj(y))/N; %计算功率谱密度Py figure; plot(f(1:N/2),Py(1:N/2)); xlabel('频率(Hz)'); ylabel('功率谱密度'); title('x=80.0*cos(2*3.14*sf*t)功率谱密度'); grid; 2)对实验所采集的转子振动信号进行频谱分析
《测试信号分析与处理》实验指导书 实验一差分方程、卷积、z变换 一、实验目的 通过该实验熟悉 matlab软件的基本操作指令，掌握matlab软件的使用方法，掌握数字信号处理中的基本原理、方法以及matlab函数的调用。 二、实验设备 1、微型计算机1台； 2、matlab软件1套 三、实验原理 Matlab 软件是由mathworks公司于1984年推出的一套科学计算软件，分为总包和若干个工具箱，其中包含用于信号分析与处理的sptool工具箱和用于滤波器设计的fdatool工具箱。它具有强大的矩阵计算和数据可视化能力，是广泛应用于信号分析与处理中的功能强大且使用简单方便的成熟软件。Matlab软件中已有大量的关于数字信号处理的运算函数可供调用，本实验主要是针对数字信号处理中的差分方程、卷积、z变换等基本运算的matlab函数的熟悉和应用。 差分方程（difference equation)可用来描述线性时不变、因果数字滤波器。用x表示滤波器的输入，用y表示滤波器的输出。 a0y[n]+a1y[n-1]+…+a N y[n-N]=b0x[n]+b1x[n-1]+…+b M x[n-M] (1) ak,bk 为权系数，称为滤波器系数。 N为所需过去输出的个数，M 为所需输入的个数卷积是滤波器另一种实现方法。 y[n]= ∑x[k] h[n-k] = x[n]*h[n] (2) 等式定义了数字卷积，*是卷积运算符。输出y[n] 取决于输入x[n] 和系统的脉冲响应h[n]。 传输函数H(z)是滤波器的第三种实现方法。 H(z)=输出/输入= Y(z)/X(z) （3）即分别对滤波器的输入和输出信号求z变换，二者的比值就是数字滤波器的传输函数。 序列x[n]的z变换定义为 X (z)=∑x[n]z-n (4) 把序列x[n] 的z 变换记为Z{x[n]} = X(z)。
信号分析与处理课后习题答案 第五章 快速傅里叶变换 1.如果一台通用计算机的速度为平均每次复乘需要50us ，每次复加需要10us ，用来就散N=1024点的DFT ，问： （1）直接计算需要多少时间？用FFT 计算呢？ （2）照这样计算，用FFT 计算快速卷积对信号进行处理是，估计可实现实时处理的信号最高频率？ 解： 分析：直接利用DFT 计算：复乘次数为N 2，复加次数为N(N-1)； 利用FFT 计算：复乘次数为20.5log N N ，复加次数为2log N N ； （1） 直接DFT 计算： 复乘所需时间2215010245052.4288T N us us s =?=?= 复加所需时间2(1)101024(10241)1010.47552T N N us us s =-?=-?= 所以总时间1262.90432DFT T T T s =+= FFT 计算： 复乘所需时间3220.5log 500.51024log 1024500.256T N N us us s =?=???= 复加所需时间422log 101024log 1024100.1024T N N us us s =?=??= 所以总时间为340.3584FFT T T T s =+= （2） 假设计算两个N 长序列1()x n 和2()x n 的卷积 计算过程为如下： 第一步：求1()X k ，2()X k ；所需时间为2FFT T ? 第二步：计算12()()()X k X k X k =?，共需要N 次复乘运算 所需时间为501024500.0512To N us us s =?=?= 第三步：计算(())IFFT X k ，所需时间为FFT T 所以总时间为230.35840.0512 1.1264FFT T T To s s s =?+=?+= 容许计算信号频率为N/T=911.3Hz 2.设x(n)是长度为2N 的有限长实序列，()X k 为x(n)的2N 点得DFT 。