本章重点
- 连续信号、离散信号的概念;周期信号与非周期信号的概念;
- 几种基本信号的理解以及信号的基本运算;
- 系统的性质;
- 连续单位阶跃信号与单位冲激信号的理解;
- 信号的基本运算;
- 系统的性质;
语音信号的特点是一维、其变化依赖于时间。它的特性有音量和频率。
系统的输入和输出都是信号,信号的种类很多。
1.1 连续时间信号和离散时间信号
信号:信号是消息的载体,一般表现为随时间变化的某种物理量。可以描述范围极其广泛的物理现象,可以分为确知信号和随机信号,或者是连续时间信号与离散时间信号。
确知信号:一个或几个自变量的函数。作为信号分析的基础。
电信号最易于传输、控制和处理。
本课程只研究确知信号。在应用和科研层面,一般是研究随机信号。现实生活中,一般是一个确知信号+一个随机信号的组合,这是因为一个确知信号总是收到随机噪声的干扰。
您可以使用Python 生成一个随机信号:
import numpy as np
import matplotlib.pyplot as plt
sample_rate = 1000
duration = 1.0
frequency = 5
t = np.linspace(0, duration, int(sample_rate * duration), endpoint=False)
random_signal = np.random.normal(0, 1, t.shape)
plt.figure(figsize=(10, 4))
plt.plot(t, random_signal)
plt.title('Continuous Random Signal')
plt.xlabel('Time (s)')
plt.ylabel('Amplitude')
plt.grid(True)
plt.show()
可以这么看:连续时间信号是函数,离散时间信号是数列。(某些时间点上的值)
- 连续信号
x1(t)→T1,x2(t)→T2若T2T1是有理数, 则Ax1(t)±Bx2(t)是周期的,其基波周期:
T0=LCM(T1,T2)
- 离散信号
x1[n]→T1,x2[n]→T2Ax1[n]±Bx2[n]一定是周期的,其基波周期:
T0=LCM(N1,N2)信号的描述:连续时间信号 x(t),x(t1,t2),离散时间信号 x(n),x(n1,n2)
在做信号处理的时候,时间可以往后退(数学意义上)。使用历史信号也是预测的基础。
连续时间信号在离散时刻点上的样本可以构成一个离散时间信号。
信号的能量与功率:
连续时间信号在[t1,t2]区间的能量定义为:
E=∫t1t2x(t)2dt连续时间信号在[t1,t2]区间的平均功率为:
P=t2−t11∫t1t2x(t)2dt需要注意的是,这个积分有可能是过零的,积分是对数,也就是出现无穷大。为了避免这种情况,我们对信号进行了分类。
离散时间信号在[n1,n2]区间的能量定义为:
E=n=n1∑n2x(n)2离散时间信号在[n1,n2]区间的平均功率为:
P=n2−n11n=n1∑n2x(n)2在无限区间上也可以定义信号的总能量:
E=T→∞lim∫−TT∣x(t)∣2 =∫−∞∞∣x(t)∣2dtE=N→∞limn=n1∑n2∣x(n)∣2我们刚刚说到最好对他们实行分类,事实上,这个问题可以按照能量的多少(有限还是无限)来分类。
- 信号的总能量有限(能量信号)
E∞<∞,P∞=0
- 信号的总能量无限,平均功率有限(功率信号)
E∞=∞,0<P∞<∞
- 信号的总能量和平均功率都无限
E∞=∞,P∞=∞判别一个信号是否是能量信号,您应该先计算它的总能量。
我们曾在高中学过函数的周期性,在周期信号上您能更直观地看出。如果一个信号是周期信号,则x(t+T)=x(t)或x(n+T)=x(n)。离散时间信号会表现出一根根火柴棒起起伏伏的样子。
值得指出的是,周期信号一定是功率信号。在无限个周期内能量越加越大,直到无穷大,但每个周期内能量是有限的。因此,从定义上可见它一定是个功率信号。
x(t)=x(t+kT)对于连续周期信号,满足上式的最小周期称为基波周期,用T0表示。
离散周期信号:一个离散信号x[n], 若对所有n均有:
x[n]=x[n+mN], m=0, 则称 x[n] 为离散周期信号。(在时间上是离散的 , 只在某些不连续的 规定瞬时 ( 本课程规定为整数值 ) 给出函数值 , 在其它时间没有定义 。)
满足上式的最小周期称为 基波周期,用N0 表示。
P∞=T1∫0T∣x(t)∣2dt或
P∞=2T1∫−TT∣x(t)∣2dtP∞=N1n=0∑N−1∣x(n)∣2或
P∞=2N+11n=−N∑N∣x(n)∣21.2 自变量变换
由于信号可视为自变量的函数,当自变量改变时,必然会使信号的特性相应地改变。
1.2.1 信号的运算与变换
x(t)=x1(t)+x2(t)指两个 ( 或多个 ) 信号之和构成另一个信号 , 它在任意瞬间的值等于两个 ( 或 ) 多个信号在同一瞬间的值的代数和。
对于离散信号有y[n]=x1[n]+x2[n],同序列号n的序列值逐项对应相加。
x(t)=x1(t)×x2(t)指两个 ( 或多个 ) 信号相乘构成另一个信号时 , 把所有相同瞬间的值一一相乘。
对于离散信号有y[n]=x1[n]×x2[n],同序列号n的序列值逐项对应相乘。
移位(time shift)
x(t±t0)x(t-t_0 ), 则是原信号向右平移t0时间;
x(t+t_0 ), 则是原信号向左平移t0时间;
对于离散信号:
x[n±n0]x[n−n0]: 向右平移n0个采样点;
x[n+n0]: 向左平移n0个采样点;
反折(time reversal)
x(−t)将x(t)沿纵坐标反折得到。
对于离散信号:
x[−n]将序列x[n]沿n=0的纵轴反折得到。
Fs = 44100;
duration = 5;
recObj = audiorecorder(Fs, 16, 1);
disp('开始录音...');
recordblocking(recObj, duration);
disp('录音完成.');
audioData = getaudiodata(recObj);
reversedAudioData = flipud(audioData);
sound(reversedAudioData, Fs);
audiowrite('reversed_audio.wav', reversedAudioData, Fs);
尺度变换(Compression and expansion)
x(at)如果将信号 x(t) 的自变量乘以正系数a为x(at) , 则:
a>1, 将波形进行压缩;
a<1, 将波形进行扩展;
若 x(t) 是已录制声音的磁带:
- x(−t) 表示磁带倒转播放产生的信号
- x(2t) 表示磁带以二倍加快播放的结果
- x(t/2) 表示原磁带放音速度降至一半产生的信号
对离散信号进行压缩时会出现什么问题?
在移位、反折、尺度变换时,一定要注意一点,那就是在对时间变量t作变换!!!
并且,要形成这样一种思路:线性变换是时移、反折、尺度打组合技。
微分(Differential)
x′(t)=dtdx(t)差分(difference)
信号的差分分为前向差分和后向差分:
一阶前向差分定义为: Δx(n)=x(n+1)−x(n)
一阶后向差分定义为:∇x(n)=x(n)−x(n−1)
积分(Integration)
x(t)=∫0tx(τ)dτ累加(Accumulation)
y[n]=k=0∑nx[k]1.2.2 信号的奇偶性
信号的奇偶分解(The odd and even decomposition of the signal)对于连续信号:
x(t)=xe(t)+xo(t)则
x(−t)=xe(−t)+xo(−t)=xe(t)−xo(t)对于连续信号,
xe(t)=2x(t)+x(−t),xo(t)=2x(t)−x(−t)在祖师爷的书里面不这么描述,使用的是Ev和Od表示奇偶。
Ev{x(t)}=2x(t)+x(−t),Od{x(t)}=2x(t)−x(−t)对于离散信号:
xe[n]=2x[n]+x[−n],xo[n]=2x[n]−x[−n]信号与其偶分量和奇分量之间还满足以下能量关系式:
连续:
∫−∞∞x2(t)dt=∫xe2(t)dt+∫xo2(t)dt离散:
n=−∞∑∞x2[n]=k=−∞∑∞xe2[n]+k=−∞∑∞xo2[n]信号的实虚分解
x(t)=α(t)+jβ(t),x∗(t)=α(t)−β(t)则:
α(t)=2x(t)+x∗(t),jβ(t)=2x(t)−x∗(t)1.3 指数信号与正弦信号
代表很广泛的一类信号。
分析很方便:其微分与积分都为指数信号。
表征很大范围信号:很多信号可以表征为复指数信号的加权和
1.3.1 连续时间复指数信号与正弦信号
连续时间复指数信号与正弦信号(continus- - time complex exponential signal and sinusoidal signals )中最常用的是复指数信号
其表示式为:
x(t)=Ceat其中:C和a一般为复数
C=α+jβ,a=r+jω0x(t)=(α+jβ)er+jω0t实指数信号 (real exponential signal)
x(t)=(α+jβ)er+jω0t若C和a为实数,则为实指数信号。
周期复指数信号 (periodic complex signal)
x(t)=Ceat令a=jω0,得x(t)=ejω0t
所有x(t)都满足x(t)=x(t+T),并且T=ω02π。因此x(t)是周期信号。
正弦信号(sinusoidal signals)
欧拉公式:
ejω0t=cos(ω0t)+jsin(ω0t)cos(ω0t)=21ejω0t+21e−jω0tsin(ω0t)=21jejω0t−21je−jω0t得
ej(ω0t+ϕ)=cos(ω0t+ϕ)+jsin(ω0t+ϕ)取实部则为正弦信号
x(t)=Acos(ω0t+ϕ)一般复指数信号
x(t)=Ceat其中,
C=α+jβ=∣C∣ejθ,a=r+jω0x(t)=∣C∣ejθertejω0t=∣C∣ertej(ω0+θ)=∣C∣ertcos(ω0t+θ)+j∣C∣ertsin(ω0t+θ)一般复指数信号
x(t)=∣C∣ertejω0t一般的复指数信号,可以将它看成是 实指数信号ert和周期复指数信号ejω0t相乘的结果;
ω0 反映了振荡信号的变化频率; r 反映了振荡信号峰值的变化趋势;
用包络线来描述信号峰值的变化趋势。
解题
x1(t)=cos6tT1=62π=3πx2(t)=cos2(4t)=21+cos(8t)直流信号不影响信号的周期。
T2=82π=4πx3(t)=cos3tT3=32π对于x1(t)+x2(t),T2T1=34为有理数,故其周期为:
T0=LCM(T1,T2)=π对于x1(t)+x3(t),T3T1=33为无理数,该周期不是周期信号。
1.3.2 离散时间复指数信号与正弦信号
复指数信号(序列):
x[n]=CeanC与a均为复数。
实指数信号
x[n]=Cean正弦信号
x[n]=cos(ω0n+ϕ)一般复指数信号
x[n]=Cean其中,
C=∣C∣ejθ,a=∣a∣ejω01.3.3 离散时间复指数序列的周期性质
连续周期复指数ejω0t具有两个性质:
ω0愈大,ejω0t振荡频率愈高;
对任何ω0, ejω0t都是周期的。
提出问题:对于离散时间 复指数信号是否也具有这两条特性呢?
问题1: ω0愈大, ejω0n振荡频率愈高???
问题2 :任何ω0, ejω0n一定是周期的?
对于有理数倍数π的ω0,ejω0n一般无法是周期的。这是因为在如此情况下的三角函数,求周期会求出无理数。
20240829作业:
基本题,P39 1.26 All,1.25 def
Details
答案
1.25(d)
x(t)=Ev{cos(4πt)u(t)}=21cos(4πt)u(t)+21cos(−4πt)u(−t)=21cos(4πt)u(t)+21cos(4πt)u(−t)=21cos(4πt)周期的,且周期为T=4π2π=21
1.25(e)
x(t)=Ev{sin(4πt)u(t)}=21sin(4πt)u(t)+21sin(−4πt)u(−t)=21sin(4πt)u(t)−21sin(4πt)u(−t)所以x(t)是非周期的。
1.25(f)
x(t)=n=−∞∑∞e−(2t−n)u(2t−n)x(t+T)=n=−∞∑∞e−(2t+(2T−n))u(2t+(2T−n))x(t)是周期的。
2T=k=1,T=211.26(a)
x[n]=sin(76πn+1)周期的,因为2πω0=76π/2π=73=Nm,N=7,满足有理数的条件。
1.26(b)
x[n]=cos(8n−π)非周期的,因为2πω0=16π1,不满足有理数条件。
1.26(c)
x[n]=cos(8πn2)不能直接像刚刚那样确定,我们不妨取x[n+kT]=cos(8π(n+kT)2)=cos(8π(n2+2kTn+k2T2)),这样看也还是看不出来。继续考虑,最后一项为常数,未知数为n,我们不妨再深入假设8πk2T2+4πkTn=2Nπ,对于任意整数n,总能找到kT使得2nkT+k2T2=16N成立,而这个成立的条件我们可以先看哈,kT必须是16的倍数。这样就符合周期性。那么它的基波周期是T=8。
1.26(d)
x[n]=cos(2πn)cos(4πn)也不能用c题的思路了。我们最好是用分解来做,
x[n]=cos(2πn)cos(4πn)=21[cos(43πn)+cos(4πn)]2πω0=83,N1=8;2πω0=81,N2=8。
所以x[n]是周期的,周期为T=8。
1.26(e)
x[n]=2cos(4πn)+sin(8πn)−2cos(2πn+6π)ω01/2π=4π/2π=81,N1=8ω02/2π=8π/2π=161,N2=16ω03/2π=2π/2π=41,N3=4故x[n]是周期的,取最小公倍数16
1.4 单位冲激函数与单位阶跃函数
1.4.1 单位阶跃函数 (Unit Step Functions)
闭合表达式(Closed Expression)
一般地,我们可以采用与阶跃函数相乘的形式来表达分段函数 , 而写出其闭合表达式 。
1.4.2 单位冲激函数
某些物理现象需要用一个 时间极短 , 但 取值很大 的函数模型来描述 , 如力学中的瞬间作用冲击力 , 电学中的雷击电闪等等 。 “ 冲激函数 ” 就是以这类实际问题为背景而引出的 。冲激函数可由不同的方式来定义 。
(a)某种函数的极限来定义: