信号分析与处理答案(第二版)

第二章习题参考解答

2.1 求下列系统的阶跃响应和冲激响应。

(1)

解当激励为时，响应为，即：

由于方程简单，可利用迭代法求解：

，，

…，

由此可归纳出

的表达式：

利用阶跃响应和冲激响应的关系，可以求得阶跃响应：

(2)

解 (a)求冲激响应

，当时，。

特征方程，解得特征根为

。所以：

…(2.1.2.1)

通过原方程迭代知，，，代入式(2.1.2.1)中得：

解得，代入式(2.1.2.1)：

…(2.1.2.2)

可验证满足式(2.1.2.2)，所以：

(b)求阶跃响应通解为

特解形式为

，，代入原方程有

，即

完全解为

通过原方程迭代之，，由此可得

解得，。所以阶跃响应为：

(3)

解

(4)

解

当t>0时，原方程变为：。

…(2.1.3.1)

…(2.1.3.2) 将(2.1.3.1)、(2.1.3.2)式代入原方程，比较两边的系数得：

阶跃响应：2.2 求下列离散序列的卷积和。

(1)

解用表格法求解

(2)

解用表格法求解

(3)和如题图2.2.3所示

解用表格法求解

(4)

解

(5)

解

(6)

解参见右图。

当时：

当时：

当时：

当时：

当时：

(7) ，

解参见右图：

当时：

当时：

当时：

当时：

当时：

(8) ，

解参见右图

当时：

当时：

当时：

当时：

(9)

，

解

(10) ，

解

或写作：

2.3 求下列连续信号的卷积。

(1) ，解参见右图：

当时：

当时：

当时：

当

时：

当

时：

当时：

(2) 和

如图2.3.2所示

解当

时：

当

时：

当时：

当时：

当时：

(3) ，

解

(4) ，

解

(5) ，

解参见右图。当时：

当时：

当时：

当时：

(6) ，

解

(7) ，

解

(8) ，

解

(9) ，

解

2.4 试求题图2.4示系统的总冲激响应表达式。

解

2.5 已知系统的微分方程及初始状态如下，试求系统的零输入响应。

(1) ；

解，，

(2) ；，

解，，

，，可定出

(3) ；，

解，

，，可定出

2.6 某一阶电路如题图2.6所示，电路达到稳定状态后，开关S于时闭合，试求输出响应。

解由于电容器二端的电压在ｔ＝０时不会发生突变，所以。

根据电路可以立出ｔ＞０时的微分方程：

，整理得

齐次解：

非齐次特解：设代入原方程可定出Ｂ＝２

则：

，

2.7 积分电路如题图2.7所示，已知激励信号为，试求零状态响应。解根据电路可建立微分方程：

当时：

由可定出，

根据系统的时不变性知，当时：

当时：

2.8 求下列离散系统的零输入响应。

(1) ；，

解

由，

可定出，

，

(2) ；，

解

由，，可定出

.

(3) ；，，

解特征方程，

，

由

可定出

2.9 求下列离散系统的完全响应。

(1) ；

解齐次方程通解：

非齐次方程特解：代入原方程得：

由可定出

(2) ；，

解齐次方程通解：

非齐次方程特解：代入原方程定出

由可定出

2.10 试判断下列系统的稳定性和因果性。

(1)

解因果的；稳定的。

(2)

解因为冲激响应不满足绝对可和条件，所以是不稳定的；非因果的。

(3)

解稳定的，非因果的。

(4)

解不稳定的，因果的。

(5)

解不稳定的，因果的。

(6) (为实数)

解时：不稳定的，因果的；

时：稳定的，因果的；

时：不稳定的，因果的。

(7)

解不稳定的，非因果的。

(8)

解稳定的，非因果的。

2.11 用方框图表示下列系统。

(1)

(2)

(3)

*2.12 根据系统的差分方程求系统的单位脉冲响应。

(1)

解

当时：,

由原方程知当时：，由此可定出

(2)

解

当时：齐次方程的通解为，

由原方程迭代求解可得为：

由此可以定出

*2.13 根据系统的微分方程求系统的单位冲激响应。

(1)

解

当时：，，代入原方程可确定

(2)

解

当时：

代入原方程，比较两边系数得：

*2.14 试求下列系统的零输入响应、零状态响应、强迫响应、自由响应。(1) ；，，

解(a)求强迫响应：假设特解为：

代入原方程，可定出；则强迫响应

(a)求自由响应：

利用冲激平衡法可知：

可定出；所以

完全解形式：，由定出

即完全响应为：

所以自由响应为：

(b)求强迫响应：假设特解为：

代入原方程，可定出；则强迫响应

(c)求零输入响应：

由可定出

(d)求零状态响应

零状态响应＝自由响应＋强迫响应－零输入响应

＝

综上所求，有：

(2) ；，，解法一用z变换求解。方程两边进行z变换，则有：

解法二：时域解法。

求强迫响应：

当时：即为常值序列，

设特解为，代入原方程可定出

当时：仅在激励作用下，由原方程知，即：

特解在时均满足方程。

求自由响应：

完全解：

由经迭代得：

由可定出完全解中系数为：

则自由响应分量为：

零输入响应：

由可以定出：

零状态响应：

*2.15 试证明线性时不变系统具有如下性质：

(1) 若系统对激励的响应为，则系统对激励的响应为；

(2) 若系统对激励的响应为，则系统对激励的响应为。

证(1) 已知，根据系统的线性试不变性有：

；令，则有：

证(2) 已知，根据系统的线性试不变性有：

令则，

所以

证毕。

*2.16 考察题图2.16(a)所示系统，其中开平方运算取正根。

(1) 求出和之间的关系；

(2) 该系统是线性系统吗，是时不变系统吗？

(3) 若输入信号是题图2.16(b)所示的矩形脉冲(时间单位：秒)，求响应。

解(1)

由系统

框图可

信号分析与处理答案第二版完整版

信号分析与处理答案第 二版 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】

第二章习题参考解答 求下列系统的阶跃响应和冲激响应。 (1) 解当激励为时，响应为，即： 由于方程简单，可利用迭代法求解： ，， …， 由此可归纳出的表达式： 利用阶跃响应和冲激响应的关系，可以求得阶跃响应： (2) 解 (a)求冲激响应 ，当时，。 特征方程，解得特征根为。所以： …(2.1.2.1) 通过原方程迭代知，，，代入式(2.1.2.1)中得：解得，代入式(2.1.2.1)： …(2.1.2.2) 可验证满足式(2.1.2.2)，所以： (b)求阶跃响应 通解为 特解形式为，，代入原方程有，即 完全解为 通过原方程迭代之，，由此可得 解得，。所以阶跃响应为： (3)

解 (4) 解 当t>0时，原方程变为：。 …(2.1.3.1) …(2.1.3.2) 将(2.1.3.1)、式代入原方程，比较两边的系数得： 阶跃响应： 求下列离散序列的卷积和。 (1) 解用表 格法求 解 (2) 解用表 格法求 解 (3) 和 如题图2.2.3所示 解用表 格法求 解

(4) 解 (5) 解 (6) 解参见右图。 当时： 当时： 当时： 当时： 当时： (7) ， 解参见右图： 当时： 当时： 当时： 当时： 当时： (8) ，解参见右图

当时： 当时： 当时： 当时： (9) ， 解 (10) ， 解 或写作：

求下列连续信号的卷积。 (1) ， 解参见右图： 当时： 当时： 当时： 当时： 当时： 当时： (2) 和如图2.3.2所示 解当时： 当时： 当时： 当时： 当时： (3) ， 解 (4) ， 解 (5) ， 解参见右图。当时：当时： 当时：

数字信号处理期末实验 语音信号分析与处理

山东建筑大学信电学院课程设计说明书 语音信号分析与处理 摘要 用MATLAB对语音信号进行分析与处理，采集语音信号后，在MATLAB软件平台进行频谱分析;并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器滤除噪声，恢复原信号。 数字滤波器是数字信号处理的基础，用来对信号进行过滤、检测和参数估计等处理。IIR数字滤波器最大的优点是给定一组指标时，它的阶数要比相同组的FIR 滤波器的低的多。信号处理中和频谱分析最为密切的理论基础是傅立叶变换（FT）。离散傅立叶变换（DFT）和数字滤波是数字信号处理的最基本内容。 关键词：MATLAB;语音信号；加入噪声；滤波器；滤波 1. 设计目的与要求 （1）待处理的语音信号是一个在20Hz~20kHz频段的低频信号。 （2）要求MATLAB对语音信号进行分析和处理，采集语音信号后，在MATLAB平台进行频谱分析；并对所采集的语音信号加入干扰噪声，对加入噪声的信号进行频谱分析，设计合适的滤波器进行滤除噪声，恢复原信号。 1 山东建筑大学信电学院课程设计说明书

2. 设计步骤 （1）选择一个语音信号或者自己录制一段语音文件作为分析对象； （2）对语音信号进行采样，并对语音信号进行FFT频谱分析，画出信号的时域波形图和频谱图； （3）利用MATLAB自带的随机函数产生噪声加入到语音信号中，对语音信号进行回放，对其进行FFT频谱分析； （4）设计合适滤波器，对带有噪声的语音信号进行滤波，画出滤波前后的时域波形图和频谱图，比较加噪前后的语音信号，分析发生的变化； （5）对语音信号进行回放，感觉声音变化。 3. 设计原理及内容 3.1 理论依据 （1）采样频率：采样频率（也称采样速度或者采样率）定义了每秒从连续信号中提取并组成离散信号的采样个数，它用赫兹（Hz）来表示。采样频率只能用 于周期性采样的采样器，对于非周期采样的采样器没有规则限制。通俗的讲，采样频率是指计算机每秒钟采集多少个声音样本，是描述声音文件的音质、音调，衡量声卡、声音文件的质量标准。采样频率越高，即采样的间隔时间越短，则在单位之间内计算机得到的声音样本数据就越多，对声音波形的表示也越精确。（2）采样位数：即采样值或取样值，用来衡量声音波动变化的参数。 （3）采样定理：在进行模拟/数字信号的的转换过程中，当采样频率f大于信s.max 号中，最高频率f的2倍时，即：f>=2f，则采样之后的数字信号完整的maxmaxs.max 保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的 5~10倍；采样频率又称乃奎斯特定理。 （4）时域信号的FFT分析：信号的频谱分析就是计算信号的傅立叶变换。连续信号与系统的傅立叶分析显然不便于直接用计算机进行计算，使其应用受到限制。而FFT是一种时域和频域均离散化的变换，适合数值计算，成为用计算机分析 离2 山东建筑大学信电学院课程设计说明书 散信号和系统的的有力工具。对连续信号和系统，可以通过时域采样，应用DFT 进行近似谱分析。

信号分析与处理习题

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 点离散傅里叶变换（闭合形式表达式）：

信号分析与处理试题

河南科技学院2006-2007学年第二学期期终考试 信号分析与处理试题 适用班级： 注意事项：1 在试卷的标封处填写院（系）、专业、班级、姓名和准考证号。 2 考试时间共100分。 一、单项选择题（本大题共10小题，每题2分，共20分） 1．下列单元属于动态系统的是（ ） A. 电容器 B.电阻器 C.数乘器 D.加法器 2．单位阶跃函数()u t 和单位冲激函数()t δ的关系是（ ） A.()/()d t dt u t δ= B.()/()du t dt t δ= C.()()u t t δ= D.()2()u t t δ= 3．()()f t t dt δ∞-∞=?（ ） A.()f t B.()t δ C.(0)f D.(0)δ 4．单位冲激函数()t δ的()F j ω=（ ） A .0 B.-1 C.1 D.2 5．设()f t 的频谱为()F j ω，则利用傅里叶变换的频移性质，0()j t f t e ω的频谱为（ ） A.0()F j ω B.()F j ω C.0[()]F j ωω+ D.0[()]F j ωω- 6．设1()f t 的频谱为1()F j ω，2()f t 的频谱为2()F j ω，利用傅里叶变换卷积定理，12()()f t f t *的频谱为（ ） A.1()F j ω B.2()F j ω C.11()()F j F j ωω* D.11()()F j F j ωω 7．序列()n m δ-的Z 变换为（ ） A.m z B.m z - C.m D.m - 8．单边指数序列()n a u n ，当（ ）时序列收敛 A.1a < B.1a ≤ C.1a > D.1a ≥ 9．取样函数()/Sa t sint t =，则(0)Sa =（ ） A.0 B.1 C.2 D.3 10．设实函数()f t 的频谱()()()F j R jX ωωω=+，下列叙述正确的是（ ）

《信号分析与处理》复习总结

信号是带有信息（如语音、音乐、图象、数据等）的随时间（和空间）变化的物理或物理现象，其图象称为信号的波形。信号是消息的表现形式，消息则是信号的具体内容。 分类：根据不同分类原则，信号可分为：连续时间信号与离散时间信号；确定信号与随机信号；周期信号和非周期信号；功率信号与能量信号等等 反因果信号：若当t ≥0时，f (t )=0;当t <0时，f (t )≠0. 系统：由若干相互作用和相互依赖的事物组合而成的具有特定功能的整体。 ???????=???≠=∞=?∞ ∞ -1)()0( 0)0( )(dt t t t t δδ()()t t δδ-= ()t δ为偶对称函数 1()d 2j t t e ωδωπ ∞-∞= ?——()t δ的逆傅立叶变 换 ()()d ()() t x t t t t x t t t δε-∞ -=-?) ()()()(000t t t x t t t x -=-δδ)(| |1 )(t a at δδ= )(t δ'是奇对称函数 ) ()(, 0)(t d d t δττδττδ='='? ? ∞ -∞ ∞ -离散时间单位: 0()(), ()()(1) m n n m n n n εδδεε+∞ ==-=--∑稳定 性 ∑? +∞-∞ =∞ +∞ -∞ <∞-=-? -z z z z n Z ε

数字信号处理期末试卷及答案

A 一、选择题（每题3分，共5题） 1、 )6 3()(π-=n j e n x ，该序列是 。 A.非周期序列 B.周期6 π = N C.周期π6=N D. 周期π2=N 2、 序列)1()(---=n u a n x n ，则)(Z X 的收敛域为 。 A.a Z < B.a Z ≤ C.a Z > D.a Z ≥ 3、 对)70() (≤≤n n x 和)190()(≤≤n n y 分别作20 点 DFT ，得 )(k X 和)(k Y ， 19,1,0),()()( =?=k k Y k X k F ，19,1,0)],([)( ==n k F IDFT n f ， n 在 范围内时，)(n f 是)(n x 和)(n y 的线性卷积。 A.70≤≤n B.197≤≤n C.1912≤≤n D.190≤≤n 4、 )()(101n R n x =，)()(72n R n x =，用DFT 计算二者的线性卷积，为使计算量尽可能的少，应使DFT 的长度N 满足 。 A.16>N B.16=N C.16

测试信号分析与处理作业实验一二

王锋 实验一：利用FFT 作快速相关估计 一、实验目的 a.掌握信号处理的一般方法，了解相关估计在信号分析与处理中的作用。 b.熟悉FFT算法程序；熟练掌握用FFT作快速相关估计的算法。 c.了解快速相关估计的谱分布的情况。 二、实验内容 a.读入实验数据[1]。 b.编写一利用FFT作相关估计的程序[2]。 c.将计算结果表示成图形的形式，给出相关谱的分布情况图。 注[1]：实验数据文件名为“Qjt.dat”。 实验数据来源：三峡前期工程 “覃家沱大桥” 实测桥梁振动数据。 实验数据采样频率：50Hz。 可从数据文件中任意截取几段数据进行分析，数据长度N 自定。 注[2]：采用Matlab 编程。 三、算法讨论及分析 算法为有偏估计，利用FFT计算相关函数 Step 1: 对原序列补N个零，得新序列x2N(n) Step2: 作FFT[x2N(n)]得到X2N(k) Step 3: 取X2N(k)的共轭，得 Step 4: 作 Step 5: 调整与的错位。 四、实验结果分析 1. 该信号可以近似为平稳信号么？ 可以近似为平稳信号，随机过程的统计特性不随样本的采样时刻而发生变化。取N=8192，分别取间隔m=500,m=700,m=1000,所得到的均值均为0.5366，方差为47369，与时间无关。

图1-1 自相关函数图 (上图表示的R0，下图为调整后的R0) 2. 该信号是否具有周期性，信噪比如何？ >> load Qjt.dat; %加载数据 N=32768; %数据长度 i=1:1:N; %提取数据 plot(i,Qjt(i)); 抛去几个极值点，从图1-2可以看出，数据具有一定的周期性，杂音比较少，说明信噪比较高。 图1-2 数据图

信号分析与处理课后习题答案

信号分析与处理课后习题答案 第五章快速傅里叶变换 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 。

信号分析与处理模拟试卷

1．具有跳变的信号在其跳变处的导数是一个 a 。 a ）强度等于跳变幅度的冲激函数 b) 幅度为无限大的冲激函数 c) 强度为无限大的冲号 d) 理想阶跃信号 2．设 x (n ) 是一个绝对可求和的信号，其有理 z 变换为 X ( z ) 。若已知 X ( z ) 在 z =0.5有一个极点，则 x (n ) 是 c 。 a ）有限长信号 b ）左边信号 c ）右边信号 d ）区间信号 3. z (t ) = 4t 2δ (2t ? 4) = b 。 a ）8δ (t ? 2) b ）16δ (t ? 2) c ）8 d ）16 4. 设两个有限长序列 x (n ) 和 h (n ) 的卷积为 y (n ) = x (n ) ? h (n ) ， y (n ) 的长度 L y 与 x (n ) 的长度L x 和 h (n ) 的长度 L h 的关系是 b 。 a ） L y = L x + L h + 1 b ） L y = L x + L h ? 1 c ） L y = L x ? L h + 1 d ） L y = L x ? L h ? 1 5. 已知 x (n ) 的 Z 变换 X ( z ) =?2.5z /(z 2 ? 1.5z ? 1), 则 X ( z ) 可能存在的收敛域是 a a )|Z|<0.5, 0.5<|Z|< 2, |Z|> 2 b) |Z|<0.5, 0.5<|Z|< 2 c) 0.5<|Z|< 2, |Z|> 2 d) |Z|> 2 二．填空题（20分，每空1分） （1）按照信号幅度和时间取值方式的不同，信号可以分为以下几种类型：连续时间信号、离散时间信号、数字信号。 （2）若一个离散时间系统满足__线性__和__时不变性则称为线性时不变系统，线性移不变系统具有因果性的充分必 要条件是系统的单位抽样响应满足下式：__h(n)=0 (当n<0时)___。 （3）快速傅里叶变换（FFT ）并不是一种新的变换形式，但它应用了系数kn N W 的_对称性__周期性__可约性__，不断地将长序列的DFT 分解成几个短序列的DFT,并减少DFT 的运算次数。其运算量是DFT 的__N 2 /[(N/2)log 2N]__倍。 （4）求积分 dt )t ()t (212-+? ∞ ∞ -δ的值为 5 。 (5)线性系统是同时具有 齐次性 和 叠加性 的系统。 (6)系统的完全响应也可以分为暂态响应和稳态响应。随着时间t 的增大而衰减为零的部分 称为系统的暂态响应 ，其余部分为系统的 稳态响应 。 （7）周期信号频谱3个典型特点：离散性、谐波性、收敛性. (8)模拟滤波器设计IIR 数字滤波器的方法有 冲激响应不变法 和 双线性变换法 。 一、判断下列说法的正误，正确请在括号里打“√”，错误请打“×”。（10分，每小题2分） 1.单位冲激函数总是满足)()(t t -=δδ （ √ ） 2.满足绝对可积条件 ∞

测试信号分析与处理作业实验五

王锋 实验五：多种功率谱估计的比较 一、实验目的 a.了解功率谱估计在信号分析中的作用； b.掌握随机信号分析的基础理论，掌握参数模型描述形式下的随机信 号的功率谱的计算方法； c.掌握在计算机上产生随机信号的方法； d.了解不同的功率谱估计方法的优缺点。 二、实验准备 有三个信号源，分别代表三种随机信号（序列）。 信号源1： 123()2cos(2)2cos(2)2cos(2)()x n f n f n f n z n πππ=+++ 其中，1230.08,=0.38,0.40f f f == z(n)是一个一阶 AR 过程，满足方程： ()(1)(1)()z n a z n e n =--+ (1)0.823321a =- e(n)是一高斯分布的实白噪声序列，方差20.1σ= 信号源2和信号源3： 都是4阶的AR 过程，它们分别是一个宽带和一个窄带过程，满足方程： ()(1)(1)(2)(2)(3)(3)(4)(4)()x n a x n a x n a x n a x n e n =--------+ e(n)是一高斯分布的实白噪声序列，方差2σ，参数如下： 三、实验内容 a. 描绘出这三个实验信号的真实功率谱波形。 b. 在计算机上分别产生这个三个信号，令所得到的数据长度 N= 256 。 注意：产生信号的时候注意避开起始瞬态点。例如，可以产生长度为512 的信号序列，然后取后面256 个点作为实验数据。 c. 分别用如下的谱估计方法，对三个信号序列进行谱估计。 1、经典谱估计 周期图法 自相关法 平均周期图法（Bartlett 法）

Welch法（可选每段64 点，重叠32 点，用Hamming 窗）2、现代谱估计 Yule - Walker方程（自相关法） 最小二乘法 注：阶次p可在3－20之间，由自己给定。 四、实验结果分析 生成的信号源

《信号分析与处理》(第二版)-徐科军、黄云志-课后标准答案

《信号分析与处理》(第二版)-徐科军、黄云志-课后答案

————————————————————————————————作者：————————————————————————————————日期：

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 ?+∞ ∞ -Ω-=Ω频谱密度函数：

吕卫阳—信号分析与处理第二次作业—北京科技大学

周期序列的频谱分析： 已知周期序列在一个周期N=4内的取值为x(n)=[0 1 2 3]采用MATLAB计算该周期序列的频谱（DTFS）。 程序： %周期序列的时域波形 x=[0 1 2 3];n=0:3; N=length(x);figure(1); stem(n,x,'*'); axis([0 4 -4 4]);grid; xlabel('n'); ylabel('x(n)'); title('周期序列时域波形'); for k=0:1:3 dk(k+1)=(x(1)*exp(-j*k*2*pi/N*0)+x(2)*exp(-j*k*2*pi/N*1)+x(3)*exp (-j*k*2*pi/N*2)... +x(4)*exp(-j*k*2*pi/N*3))/N; realdk(k+1)=real(dk(k+1)); imagdk(k+1)=imag(dk(k+1)); magnitude(k+1)=abs(dk(k+1)); phase(k+1)=angle(dk(k+1)); end %周期序列的频谱：实部和虚部 k=0:1:3; figure(2); subplot(2,1,1); stem(k,realdk(k+1),'*'); axis([0 4 -4 4]); xlabel('k'); ylabel('Real Part of d(k)');grid; subplot(2,1,2); stem(k,imagdk(k+1),'*'); axis([0 4 -4 4]); xlabel('k'); ylabel('Imaginary Part of d(k)');grid; %周期序列的频谱：幅值和相位 figure(3); subplot(2,1,1); stem(k,magnitude(k+1),'*'); axis([0 4 -4 4]);

2020年整理信号分析与处理期末试卷A.doc

学年第二学期期末考试 信号分析与处理试卷(A) 使用班级答题时间120分钟 一、判断题（本大题共10小题，每题2分，共20分） 1、单位冲激函数总是满足)t ( )t(- =δ δ。（） 2、满足绝对可积条件∞ < ?∞∞-dt)t(f的信号一定存在傅立叶变换，不满足这一条件的信号一定不存在傅立叶变换。（） 3、非周期信号的脉冲宽度越小，其频带宽度越宽。（） 4、所有周期信号的频谱都是离散谱，并且随频率的增高，幅度谱总是渐小的。（） 5、离散时间信号的频谱都是周期的。（） 6、信号()()2 7/ 8 cos+ =n n xπ是周期信号。（） 7、信号0 )4 (2= - ?∞∞-dt t δ。（） 8、因果系统时指系统在 t时刻的响应只与 t t=时刻的输入有关（） 9、线性系统是指系统同时满足叠加性和齐次性（） 10、过渡带即为通带与阻带之间的频率范围。（） 二、填空题（本大题共9小题10个空，每空2分，共20分） 1、我们把声、光、电等运载消息的物理量称为。 2、幅度有限的周期信号是信号。 3、已知}1 ,3,2{ ) ( 1 - = k f，}2,0,0,1,3{ ) ( 2 = k f，则卷积和f1(k)*f2(k)= 。 4、若信号f(t)的最高频率是2kHz，则t) f(2的乃奎斯特抽样频率为。 5、若一个离散时间系统满足_____________和____________,则称为线性时不变系统。 6、实现滤波功能的系统称为_____________。 7、 () 12 1 4 t dt δ - -= ? 8、 sin 22 t t ππ δ ???? -*+= ? ? ???? 9、周期信号频谱3个典型特点：离散性、谐波性、。 三、选择题（本大题共10小题，每题2分，共20分）

信号分析与处理课后习题答案

信号分析与处理课后习题答案 第五章 快速傅里叶变换 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 。

信号分析与处理实验

信号分析与处理实验实验一时间信号的产生 班级：自动化1101班 姓名：陈宝平 学号： 成绩：

1. 实验目的 数字信号处理系统中的信号都是以离散时间形态存在的。研究离散时间信号，首先 需要产生出各种离散时间信号。使用MATLAB 软件可以方便的产生各种常见的离散时间信号，还具有强大的绘图功能，便于用户直观地处理输出结果。 通过本实验，学习用MATLAB 产生一些常见的离散时间信号，并通过MATLAB 绘图工具对产生的信号进行观察加深对常见离散信号和信号卷积和运算的理解。 2. 实验原理 离散时间信号用x(n)来表示，自变量n 必须是整数；连续时间信号用x(t)来表示。常见的时间信号如下： (1) 单位冲激序列δ(n)=???≠=0,00,1n n ; 如图(a); 单位冲激信号δ(t)=? ??≠=0,00 ,1t t ;如图(b): (a) (b) 如果δ(n)在时间轴上延迟了k 个单位，得到δ(n-k)=? ??≠=k n k n ,0,1。 (2) 单位阶跃序列u(n)=???<≥0,00,1n n ; 如图(c)：单位阶跃信号u(t)=???<≥0 ,00 ,1t t ; 如图(d): (c) (d) 如果u(n)在时间轴上延迟了k 个单位,得到u(n-k)=???<≥k n k n ,0,1。

(3) 矩形序列R N (n)=???≥<-≤≤),0(,0) 10(,1N n n N n ,矩形序列有一个重要的参数，就是序列宽度 N 。R N (n)与u(n)之间的关系为R N (n)= u(n)- u(n-N)。如图(e): 单位矩形信号R (t)=? ??≥<≤≤),0(,0) 0(,1T t t T t , R(t)=u(t)-u(t-T),如图(f): (e) (f) （4）正弦序列x(n)=Acos(ω0n+?)。只有当 2ωπ 为有理数时，正弦序列具有周期性, 如图(g): 正弦信号x(t)=Acos(?ω+t 0),如图(h); (g) (h) (5)单边实指数序列x(n)=a n u(n)，当a>0时，该序列均取正值，当a<0时，序列在正负摆动。如图分别为：x(n)=1.2n , x(n)=(-1.2)n , x(n)=0.8n , x(n)=(-0.8)n 的图。

信号分析及处理期末考试Word版

2014-2015学年第一学期期末考试 《信号分析与处理中的数学方法》 学号： 姓名： 注意事项： 1.严禁相互抄袭，如有雷同，直接按照不及格处理； 2.试卷开卷； 3.本考试提交时间为2014年12月31日24时，逾期邮件无效； 4.考试答案以PDF 和word 形式发送到sp_exam@https://www.sodocs.net/doc/6411250406.html, 。 1、叙述卡享南—洛厄维变换，为什么该变换被称为最佳变换，何为其实用时的困难所在，举例说明其应用。 解：形为λφ(s ) = ∫C(t,s)φ(t)dt T (1-1) 的方程称为齐次佛莱德霍姆积分方程，其中φ（t ）为未知函数，λ是参数，C （t,s ）为已知的“核函数”，它定义在[0，T]×[0，T]上，我们假定它是连续的，且是对称的： (t,s)=s (s,t) (1-2) 使积分方程(1-1)有解的参数λ称为该方程的特征值，相应的解φ(t)称为该方程的特征函数。 又核函数可表示为： C(t,s)=∑λn φn (t)φn (s)∞ n=1 （1-3） 固定一个变量（例如t ），则式（1-3）表示以s 为变量的函数C(t,s)关于正交系{φn (s)} 的傅里叶级数展开，而傅里叶级数正好是λn φn (t)。 设x （t ）为一随机信号，则其协方差函数 s（t,s ）=s {[x(t)-E{x(t)}][x(s)-E{x(s)}]}是一个非随机

对称函数，而且是非负定的。为了能方便地应用式(1-3)，假定C(t,s)是正定的，在多数情况下，这是符合实际的。当然，还假定C(t,s)在[0，T]×[0，T]上连续。 现在用特征函数系{φn (t)}作为基来表示x （t ）： x(t)=∑αn φn (t)∞ n=1 (1-4) 其中 αn =∫x （t ）φn （t ）dt T 因为{φn （t ）}是归一化正交系，所以展开式（1-4）类似于傅里叶级数展开。但是因 为x （t ）是随机的，从而系数x n 也是随机的，因此这个展开式实际上并不是通常的傅里叶展开。 式(1-4)称为随机信号的卡享南-洛厄维展开。因为这种变换能使变换后的分量互不相关，而且这种展开的截断既能使均方差误差最小，又能使统计影响最小，故具有最优性。 卡享南-洛厄维变换没有固定的变换矩阵，它依赖于给定的随机向量的协方差阵。正是这种变换的特点，也是它在实际使用时的困难所在，因为它需要依照不固定的矩阵求特征值和特征向量。 卡享南-洛厄维变换应用在数据压缩技术中。按照最优化原则的数据压缩技术可以解决通讯和数据传输系统的信道容量不足和计算机存储容量不足的问题。通过对信号作正交变换，根据失真最小的原则在变换域进行压缩。卡享南-洛厄维变换被选用并不是偶然的，因为这种变换消除了原始信号x 的诸分量间的相关性，从而使数据压缩能遵循均方误差最小的准则实施。 2、最小二乘法的三种表现形式是什么？以傅里叶级数展开为例说明其各自的优缺点。 解：希尔伯特空间中线性逼近问题的求解方法称为最小二乘法。通常它有三种不同的表现形式：投影法、求导法和配方法。我们以傅里叶级数展开为例来说明。 投影法： 设X 为希尔伯特空间，{e 1,e 2,e 3……}为X 中的一组归一化正交元素，x 为X 中的某一元素。在子空间M=span{e1,e2,e3……}中求一元素m ，使得 ‖x-m0‖=min‖x?m‖m∈M (2-1) 由于M 中的元素可表示为e 1,e 2,e 3……的线性组合，那么问题就转化为求系数 α 1 ,α2……使得 ‖x-∑a k e k ‖∞ k=1=min 2-2 投影定理指出了最优系数α1,α2……应满足

信号分析与处理 杨西侠 第2章习题答案

2-1 画出下列各时间函数的波形图，注意它们的区别 1）x 1(t) = sin Ω t ·u(t ) 2）x 2(t) = sin[ Ω ( t – t 0 ) ]·u(t ) 3）x 3(t) = sin Ω t ·u ( t – t 0 ) -1

4）x2(t) = sin[ ( t – t0) ]·u( t – t0) 2-2 已知波形图如图2-76所示，试画出经下列各种运算后的波形图 （1）x ( t-2 ) （2）x ( t+2 )

（3）x (2t) （4）x ( t/2 ) （5）x (-t) （6）x (-t-2)

（7）x ( -t/2-2 ) （8）dx/dt 2-3 应用脉冲函数的抽样特性，求下列表达式的函数值 (1)?+∞ ∞--)(0t t x δ(t) dt = x(-t 0) (2)?+∞ ∞--)(0t t x δ(t) dt = x(t 0) (3)?+∞∞ --)(0t t δ u(t - 20t ) dt = u(2 t ) (4)?+∞ ∞--)(0t t δ u(t – 2t 0) dt = u(-t 0) (5)() ?+∞∞ --+t e t δ(t+2) dt = e 2-2 (6)()?+∞ ∞-+t t sin δ(t-6π ) dt = 6 π + 2 1

(7) ()()[]?+∞ ∞-Ω---dt t t t e t j 0δδ =()?+∞ ∞ -Ω-dt t e t j δ–?+∞∞ -Ω--dt t t e t j )(0δ = 1-0 t j e Ω- = 1 – cos Ωt 0 + jsin Ωt 0 2-4 求下列各函数x 1(t)与x 2(t) 之卷积，x 1(t)* x 2(t) (1) x 1(t) = u(t), x 2(t) = e -at · u(t) ( a>0 ) x 1(t)* x 2(t) =?+∞ ∞---ττττ d t u e u a )()( = ?-t a d e 0 ττ = )1(1at e a -- x 1(t)* x 2(t) =ττδτδτπ d t t u t )]1()1([)]()4 [cos(---+-+Ω?+∞ ∞- = cos[Ω(t+1)+ 4 π ]u(t+1) – cos[Ω(t-1)+ 4 π ]u(t-1) (3) x 1(t) = u(t) – u(t-1) , x 2(t) = u(t) – u(t-2) x 1(t)* x 2(t) = ? +∞ ∞ -+-----τττττd t u t u u u )]1()()][2()([ 当 t <0时，x 1(t)* x 2(t) = 0 当 0

信号分析与处理第四次作业10月31号

信号分析与处理第四次作业 （数字信号处理 王世一 P119） 2. ()x n 表示一个具有周期为N 的周期性序列，而() X k 表示它的离散傅里叶级数的系数，也是一个具有周期为N 的周期性序列。试根据 ()x n 确定() X k 的离散傅里叶级数的系数。 3．在题图3中表示了两个周期都为6的周期性序列，确定这两个序列的周期卷积的结果 3()x n ，并画出草图。 题图3 ) (~2n x 5.如果 ()x n 是一个具有周期为N 的周期性序列，它也是具有周期为2N 的周期性序列。令1() X k 表示当 ()x n 看作是具有周期为N 的周期性序列的DFS 系数。而2 () X k 表示当 ()x n 看作是具有周期为2N 的周期性序列的DFS 系数。当然1() X k 是具有周期为N 的周期性序列，而2() X k 是具有周期为2N 的周期性序列，试根据1() X k 确定2 () X k 。 7.求下列序列的DFS ： （a ）}{ 1,1,1,1-- （b ）}{ 1,,1,--j j 8．计算下列各有限长序列的离散傅里叶变换（假设长度均为N ）： （a ）()()δ=x n n （b ）0()()δ=-x n n n 00<

()x n 圆周移位两个点） ： 144()((2))()=-x n x n R n 244()(())()=-x n x n R n 题图9 n 题图10 10.在题图10中表示了两个有限长序列，试画出它们的六点圆周卷积： 11.有限长序列的DFT 对应于序列在单位圆上的z 变换的取样。例如一个10点序列()x n 的 DFT 对应于图10-1表示的10个等间隔点上()X z 的取样。我们希望找出在图10-2所示的围 线上()X z 的等间隔取样，即 2[( )()]1010 0.5() ππ +=k j z e X z 。证明图如何修改()x n 以获得一个序 列1()x n 致使1()x n 的DFT 对应于所希望的()X z 的取样。 的圆 题图11-1 12.频谱分析的模拟数据以10kHz 取样，且计算了1024个取样的离散傅里叶变换。试确定