信号与系统的基本概念 §1.1 绪言 信号与系统是一门重要的专业基础课。是许多专业(通信、信息处理、自动化、计算机、系统工程)的必修课。重要性体现在两个方面:一是我们将来从事专业技术工作的重要理论基础;二是上述各类专业硕士研究生入学考试课程。 在教学计划中起着承前启后的作用,前期课程是高数、微分方程、差分方程、工程数学中的积分变换(傅立叶变换和拉普拉斯变换),还有电路分析基础;而其本身是后续专业课(通信原理、数字信号处理)的基础。 信号 研究的主要内容:顾名思义 系统 合成:信号 系统 响应 一个典型的电系统—通信系统 信息源 转换 电信号 电信号 还原 受信者(声音、文字、图象) 发送设备 传输信道 接收设备 输入信号/激励 输出信号/响应 通信系统  系统: 控制系统 抽象为理想化的模型,讨论激励与响应的关系 经济系统  信号:时间的函数f(t),一维函数,确定信号 * 信号与系统的关系:互相依存 信号是运载消息的工具,要很好的利用信号,需经过系统的传输、处理. 系统则是为传输信号或对信号进行处理而由元器件构成的某种组合。离开了信号,系统就失去了意义. §1.2 信号 定义:信号是带有信息的(如声音、图象等)随时间(或空间)变化的物理量。 本课程主要研究电信号(电流、电压)。 信号的分类:从不同的角度 1 从函数的定义域(时间)是否连续:  连续时间信号:在连续的时间范围内有定义。t是连续的,f(t)可是,也可不是 表达方式 时间的函数(解析式),如f(t)=Asinπt 波形图表示:  上述两种表达方式,可以互换。信号和函数两个词可互相通用 离散时间信号:在一些离散的瞬间才有定义。t=kT点上有定义,其余无定义 序列f(k)=2k,k≥0 表达方式 图形表示: 序列值f(k)={0、1、2、4、8、……}  2 从信号的重复性:  周期信号:定义在(-∞,+∞)区间,每隔一定时间T重复变化 连续f(t)=f(t+mT) 离散f(k)=f(k+mK) K为整数  非周期信号:不具有周期性的信号 例:正弦序列f(k)=sinkβ β为角频率,反映周期性重复的速率, 决定序列是否具有周期性 按定义:sinkβ=sin(β·k+m·2π) β=时, =12,为整数,是周期序列,k =12 β=时,=,为有理数,是周期序列,k =31 β=时, =4π,为无理数,是非周期序列 3 实信号:物理可实现的 复信号:实际上不能产生,但理论分析重要——复指数信号 表达式:f(t)=est ,-∞<t<+∞, δ= σ+jω f(t)=e(σ+jω)t=eσ t·ejωt = eσ tcosωt+j eσ t sinωt σ>0,增幅振荡 σ<0,衰减振荡 σ=0,等幅振荡  当ω=0,f(t)= eσt为实指数信号 当σ=ω=0,f(t)=1,为直流信号 重要特性:对时间的微分和积分仍然是复指数信号。 4.从能量有限和功率有限的角度: 能量信号:也就是能量有限信号,0<E<∞(p=0),如矩形脉冲、衰减的指数 功率信号:也就是功率有限信号,0<P<∞(E —>∞),如周期信号、阶跃信号 信号f(t)的能量E  f(t)|2dt 信号f(t)的功率P|f(t)|2dt  §1.3 信号的基本运算 一 加法和乘法 f(·)=f1(·)+f2(·) 瞬时和 f(·)=f1(·)·f2(·) 瞬时积 例1.3-1 2k+0 k<-2 f1(k)+ f2(k) = 2k +2-k k=-1、-2 k+1 +2-k k≥0 0 k<-2 f1(k)× f2(k) = 1 k=-1、-2 (k+1)×2-k k≥0 二 反转和平移 反转: f(t)—>f(- t) 以纵坐标为轴反折  倒相: f(t)—>-f(t) 以横坐标为轴反折  平移:右移 f(t)—>f(t-t0) 左移 f(t)(t)—>f(t+t0)   平移与反折结合:f(t)—>f(-t-t0) 注意:先平移后反转f[-(t+t0)] 若先反转f(-t)则f(-t-t0)为左移  t 三 尺度变换(横坐标展缩) f(t)—>f(at) 若a>1,以原点(t=0)为基准,压缩1/a 若0<a<1,以原点(t=0)为基准,展宽1/a 若a<0,反转并压缩或展宽至1/|a|  四 复合运算 f(t)—>f(-at+b) 顺序:先平移f(t)—>f(t+b);再反转f(-t+b);最后尺度变换f(-at+b). 逆符合运算f(-at+b)—>f(t) 顺序:先尺度变换 f(-t+b);再反转f(t+b);最后平移f(t) 例:已知f(5-2t)的波形如图,试画出f(t)的波形 解题思路:f(5-2t)f(5-2×2t)= f(5-t) f(5+t)f(5+t-5)= f(t)   §1.4 阶跃函数和冲激函数 重要性:完成信号的时域分解 f(t)可分解为不同时刻、不同幅度阶跃函数的连续和 f(t)可分解为不同时刻、不同幅度冲激函数的连续和 可使信号的分析、尤其是系统的分析更加简单、灵活 必要性:不是普通函数,而是奇异函数,有许多特殊的性质 重 点:引入两个函数的概念,讨论的性质 一 阶跃函数和冲激函数的定义 0, t<0 1 阶跃函数 rn(t)= , t=0 波形: 1, t>0  2 冲激函数 pn(t) 幅度—>∞ 宽度—>0 强度始终为1 波形:  表达式: 0 t< rn(t)= t -<t< —> [条件:n—>∞,斜率无限大,区间 1 t> (-,)—>0] 0 t<- rn‘(t)= pn(t)=  -<t< —> [条件:n—>∞,幅度无限大,宽度—>0] 0 t> pn(t)的强度始终为1 3 与的关系 =  =(注意积分上、下限)= 0,t>0 1,t<0 4 冲激函数的另2种定义  荻拉克给出 =0, t≠0 ,函数波形下的面积为1 该定义物理概念较明确,最易理解  的广义函数定义 (严格的数学定义) 检验函数:连续的,具有任意阶导数,且及其各阶导数在无限远处急速下降(|t|—>∞,比1/|t|m下降更快)的普通函数(如e-| t2|等) 按广义函数理论,意义为: 作用于的效果是给它赋值 5 的移位及强度表示: : t=0处的冲激 : t= t1处的冲激 : t=t2处,强度为A : t=t3处,强度为-A  二 的性质 1 取样性质(筛选性质): 2 与普通函数的乘积:f(t)= f(0) ∵[ f(t)]dt=[f(t)]dt=f(0) 又[f(0)]dt= f(0)dt= f(0) ∴ 按广义函数相等的原理可得: f(t)= f(0) f(t)dt=f(0)dt= f(0) 注意公式成立的条件:f(t)、也必须属于急降的检验函数。 例:t·=0·=0 e-αt·= e-α0·= e-3t-1dt=e-3·0-1·1=e-1(其中1为强度) 3 移位 dt= 对普通函数f(t),也有 f(t)= f(t1) f(t)dt= f(t1) 分段连续函数在区间(-∞,∞ )的导数。 跳跃度 Ji=f(ti+)-f(ti-) 广义函数概念:ti处导数为:Ji ∴ f‘(t)= f‘c(t)+Ji 例:1.4-2 求f‘(t) -∞<t<∞ 解:f(t) = 0 t<0,t>3 2+t 0<t<3  方法一:直接用上述结论 两个间断点 t1=0,J1= f(0+)-f(0-)=2 t2=3,J2= f(3+)-f(3-)=0-4=-4 ∴ f‘(t)=[-]+2-4 方法二:从函数求导 f(t)=(2+t)[-] f‘(t) =(2+t)‘[-] +(2+t)[-]‘ =[-]+(2+t)[-] =[-]+2-4 所求得的f‘(t)如下图  尺度变换:= 实际是强度变化,而不是展缩 推导:从dt研究: 若a>0,|a|=a,令x=at ∵ dt== 而dt= ∴ = 若a<0,同理可证 5 奇偶性:是偶函数 取a=-1,= 三 的导数和积分 1 导数定义:: dt= -dt = - 推导:分步积分 =dt =|-dt = 0-dt = - n阶导数::dt =(-1)n 2 导数的性质: 与f(t)的乘积:f(t)= f(0)- f`(0) 移位: f(t)= f(t1)- f`(t1) f(t)dt = - f`(t1) 尺度变换:=· =· 奇偶性:取a = - 1,=(-1)n 当n为偶数时,有= 是偶函数 当n为奇数时,有= - 是奇函数 3 的积分: 积分的区间为(-∞,t)时,区间为(-∞,+∞)时 =dx dt=1 =dx dt=0 非普通,仅是表达形式 r(t)=dx =1·dx =t· 普通积分 * 有关信号的几个概念: 1.无时限信号:在t(-∞,+∞)内均有f(t)≠0 2.有始信号:t<t1时f(t)= 0, t>t1时f(t)≠0 3.有终信号:t<t2时f(t)≠0,t>t2时f(t)= 0 4.因果信号:t<0时f(t)=0; t>0时f(t)≠0 , f(t)·v(t)表示 5.反因果信号:t≥ 0时f(t)=0; t<0时f(t)≠0 , f(t)·v(-t) 6.时限信号:在(t1,t2)内,f(t)≠0 *.抽样信号:f(t)==Sa(t) -∞<t<∞ 性质:(1)是t 的偶函数。 (2)f(t)= f(0)=1 (3)当t=kπ(k=……)时,f(t)=0 (4)f(t)dt=dt=π (5)=0  例:写出f(t)的时域表达式,并画出波形,求f(t)、f‘(t)、f‘‘(t) f(t) = sint[v(t)-v(t-π)] f`‘(t)= cost[v(t)-v(t-π)] + sint[-] = cost[v(t)-v(t-π)] f‘‘(t)= - sint[v(t)-v(t-π)] + cost[-] = - sint[v(t)-v(t-π)] + +  § 1.5 系 统 系统分析:实际物理问题→数学模型→求出解答→结果的物理解释。 主要讨论:??即时系统(无记忆系统):响应仅取决于激励,即电阻组成,用代数方程描述 ??动态系统(记忆系统):相应与激励有关,而且与过去历史状态有关(初始条件) 。含有记忆元件(电容、电感),由微分方程描述。 系统的描述: ??数学模型? ??框图表示 两种描述可互换。 1.系统的数学模型 ??连续系统?— 微分方程 RLC串联电路  由KVL: uL(t)+ uR(t)+ uC(t)= uS(t) 由各元件端口电压与电流的关系:i(t)=C*uC’(t) uR(t)=R*i(t)=R*C*uC’(t) uL(t)=L*i’(t)=L*C*uC”(t) 整理:uC”(t)+R/L* uC’(t)+1/L/C * uC(t)=1/L/C * uS(t) 二阶线性微分方程求解:需已知初始条件 uC(0), uC’(0). 结论:有以上数例可见,虽然系统的具体内容各不相同,但描述各系统的数学模型都是微分方程,因此在系统分析中,常抽去系统的物理含义,而作为一般意义下的系统来研究,以便于揭示系统的一般特性。 ??离散系统---差分方程 例1:人口问题 y(k) = y(k-1) + a*y(k-1) – b*y(k-1) + f(k) ↓ ↓ ↓ ↓ ↓ 第k年人口 第(k-1)年人口 出生 死亡 迁入 整理:y(k)-(1+a-b)*y(k-1)=f(k) ??一阶差分方程 ??结论:有以上数例可见,虽然系统的内容各不相同,但描述这些离散时间系统的数学模型都是差分方程,因而也能用相同的数学方法来分析。 2.系统的框图表示 ??连续系统:基本单元有三个:积分器、加法器、数乘器??   例1.5-2、已知框图表示,写出微分方程。 解:设右方积分器的输出为x(t) 左输出:x”(t)=f(t)-a0*x(t)-a1*x’(t) →f(t)=x”(t)+ a1*x’(t)+ a0*x(t) (1) 右  输 出:y(t)= b2*x”(t)+ b1*x’(t)+ b0*x(t) (2) 为求y(t)与f(t)的关系,消去中间变量x(t)及其导数。 由(2): a0*y= b2*( a0*x”)+ b1*( a0*x’)+ b0*( a0*x) a1*y’= b2*( a1 *x”)’+ b1*( a1*x’)’+ b0*( a1*x)’ y”= b2*(x”)”+ b1*(x’)”+ b0*(x)” 相加:y”+ a1*y’+ a0*y= b2*[x”+ a1*x’+ a0*x]” +b1*[x”+ a1*x’+ a0*x]’+ b0*[x”+ a1*x’+ a0*x] ∴ y”(t)+ a1*y’(t)+ a0*y(t)= b2*f”(t)+ b1*f’(t)+ b0*f(t) ?离散系统:延迟单元,加法器、数乘器 ??f(k)?→  → y(k)=f(k-1) 例1.5-3、 已知离散系统框图,写出差分方程。  解:2个延迟单元为二阶系统,设左边延迟单元输入为x(k) 左:x(k)=f(k)-a0*x(k-2)- a1*x(k-1)→ x(k)+ a1*x(k-1)+ a0*x(k-2)=f(k) (1) 右: y(k)= b2*x(k)- b0*x(k-2) (2) 为消去x(k),将y(k)按(1)式移位。 a1*y(k-1)= b2* a1*x(k-1)+ b0* a1*x(k-3) (3) a0*y(k-2)= b2* a0*x(k-2)-b0* a0*x(k-4) (4) (2)、(3)、(4)三式相加:y(k)+ a1*y(k-1)+ a0*y(k-2)= b2*[x(k)+ a1*x(k-1)+a0*x(k-2)]- b0*[x(k-2)+a1*x(k-3)+a0*x(k-4)] ∴ y(k)+ a1*y(k-1)+ a0*y(k-2)= b2*f(k)- b0*f(k-2)═>差分方程 ??结论:已知框图,写方程的步骤。 选中间变量x(.)。 写出个加法器输出信号的方程。 消去中间变量。 ??动态系统是否为线性系统: 例1、y(t)= f(t)*x(0)+ ∫t0f(x)dx 非线性 分解特性:??当 f(t)=0时,应得到yx (t)←x(0) , 实际 y(t)=0. 当x(t)=0时,得到yf (t)= ∫t0f(x)dx ∴不满足。 例2: y(t)= t*x(t)+sint*f(t) 分解特性: 令f(t)=0, yx(t)= t*x(t). 令x(0)=0,yf (t)= sint*f(t) 满足。 零输入线性;T[a1*x1(0)+ a2*x2 (0)]= t*[ a1* x1(0)+ a2* x2 (0)] = a1*t*x1(0)+ a2*t* x2 (0) = a1*T[x1 (0)]+ a2*T[x2 (0)] 零状态线性:T[a1* f1(t)+a2* f2 (t)] = sint*[ a1* f1 (t)+ a2* f2 (t)] = a1*sint* f1 (t)+ a2*sint* f2 (t) 均满足,为线性。 §1.6 系 统 的 性 质 系统可分为: ? 线性系统和非线性系统 ? 时变系统和时不变系统 ?? 因果系统和非因果系统 ???稳定系统和非稳定系统 线性时不变系统(Linear Time Invariant)LTI系统。 激励与响应关系的表示:  线性: 齐次性: T[a*f(t)= a*T[f(t)]= a*y(t) 输入增大a倍,响应也增大a倍。 可加性: T[f1 (t)+ f2 (t)]=T[f1 (t)]+T[f2 (t)] = y1 (t)+ y2 (t) 等于激励和响应等于响应之和。 线性:既是齐次的又是可加的。 T[a1* f1 (t)+ a2* f2(t)]= a1*T[f1 (t)]+ a2*T[f2 (t)] 例;1、积分器 f(t)→→y(t)= T[f(t)]=∫t-∞ f(x)dx a1* f1 (t)+ a2* f2 (t) →→ T[a1* f1 (t)+ a2* f2 (t)] = ∫t-∞[a1* f1 (x)+ a2* a2 (x)]dx = a1*∫t-∞f1 (x)dx+ a2*∫t-∞f2 (x)dx = a1*T[f1 (t)]+ a2*T[f2 (t)] ∴ 此系统是线性的。 2、平方运算 f(x) → 平方 → y(t)=T[f(t)]=f2(t) a*f(t)→平方→T[a*f(t)]=[a*f(t)]2=a2*f2(t) = a2*T[f(t)]≠ a*T[f(t)] f1(x)+ f2(x)→平方→T [f1 (x)+ f2 (x)]=[ f1 (x)+ f2 (x)]2≠f1 2(x)+ f2 2(x) ∴ 此系统是非线性。 3、f(k)→系统→y(k)=f(k)*f(k-1) a*f(k)→系统→T[a*f(k)]= a*f(k)*a*f(k-1) = a2*f(k)*f(k-1) ≠ a*f(k)*f(k-1) f1 (k)+ f2 (k) → 系统→T[f1 (k)+ f2 (k)] =[ f1 (k)+ f2 (k)]*[ f1 (k-1)+ f2 (k-1)]≠f1 (k)* f1 (k-1)+ f2 (k)* f2 (k-1) ∴ 此系统是非线性的。 ??动态系统中线性性质的应用:  → y(t) = T[{x(0)},{f(t)}] 线性T[{x(0)},0]+T[0,{f(t)}] = yx (t)+ yf (t) 分解特性: 零输入响应:yx (t)=T[{x(0)},0],对多个x(0),满足零输入线性。 零状态响应:yf (t)=T[0,{f(t)}],对多个f(t),满足零状态线性。 现性:满足分解特性,又具有零状态线性和零输入线性的系统。 二、时不变性 定义:系统的参数都是常数,称为时不变(非时变)系统。 性质:激励延迟一定的时间td,则它引起的零状态响应也延迟td. 设T[{0},f(t)]= yf (t) T[{0},f(t- td )]= yf (t- td ) 例 1、yf (t)= e-t * f(t) ═> 时变 2、y(k)+(k-1)* y(k-1) = f(k) ═> 时变 3、T[{0},f(t- td )] = e-t * f(t- td ) ═> 非时变 4、yf (t- td ) = e-(t-td) * f(t- td ) ═>非时变 三、本书只讨论线性时不变系统。 即LTI(Linear Time Invariant)的特性。 LTI—常系数线性微分方程(差分方程)描述。 如:y’(t) + 2 * y (t) = f’(t)- 2 * f(t) y’(t) + sint * y(t) = f(t) 线性系统 y”(t) + [y(t)]2 = f(t) 非线性系统 证:y1”(t)+[ y1”(t)]2 = f1 (t) y2”(t)+[ y2 (t)]2 = f2 (t) f1 (t)+ f2 (t)=[ y1 (t)+ y2 (t)]”+[ y1 (t)]2+[ y2 (t)]2 2、LTI连续系统具有微分特性: 若 T[{0},f(t)]= yf (t) 则 T[{0},f’(t)]= yf’(t) 证:由时不变 T[{0},f(t-△t)]= yf (t-△T) 由线性T[{0},[f(t)-f(t-△t)]/ △t] =[ yf (t)- yf (t-△T)]/ △T. △t→0, T[{0},f’(t)]= yf’(t). 3、LTI连续系统具有积分特性 若T[{0},f(t)]= yf (t),且f(-∞)=0, yf (-∞)=0 则 T[{0},∫t-∞f(x)dx]=∫t-∞yf (t)dt 例:LTI系统,在零状态条件下激励f1 (t)与响应y1 (t)如图,求f2 (t)时的y2 (t). 解:因有f2 (t)= ∫t-∞f1 (τ)dτ 故有y2 (t)= ∫t-∞y1 (τ)dτ     ??LTI连续系统具有微分特性:T[{0},f(t)]= yf (t), →T[[0],df(t)/dt]=d[yf(t)]/dt ??LTI连续系统具有积分特性: T[{0},f(t)]= yf (t), →T[{0},∫t-∞f(x)dx]=∫t-∞yf (x)dx 条件是f(-∞)=0, yf (-∞)=0. 例:1.6-1 [P31} 解:(1)首先符合分解特性: ??yx (t) = a * x(0) 零输入线性 ??yf (t) = b *∫0tf(τ)dτ t≥0 零状态线性 ∴是线性的。 时不变性针对零状态响应yf (t),判断yf (t-t0)是否等于T[f(t-t0)]. ∵yf (t-t0)= b *∫0(t-t0) f(τ)dτ t-t0≥0 T[f(t-t0)]= b *∫0tf(τ-t0)dτ t≥t0 = b *∫-t0(t-t0)f(x)dx t≥t0 (x=τ-t0) = b *∫0-t0f(x)dx+b* ∫0t-t0f(x)dx = b *∫0t-t0f(x)dx t≥t0 {条件:t=0时,接入f(t).即t<0时,f(t)= 0} ∴ T[f(t-t0)]= yf (t-t0)是时不变系统。  (2)同上。 三、因果性、 因果系统:响应(零状态响应)不出现于激励之前的系统。 定义:如果f(t)= 0,t<t0则yf(t)= T[{0},f(t)]= 0 t<t0 因果:yf(t)= 3 * f(t-1)   非因果:yf(t)= 3 * f(t+1) yf(t)=f(2*t) yf(1)=f(2) yf(k)=f(k2) yf(2)=f(4)   四、稳定性:对有界激励,其零状态响应也是有界的。 定义:如果|f(t)|<∞ ,则| yf (t)|<∞ ∵yf(k)= f(k) + f(k-1) ∴系统是稳定的。 ∵yf(t)=∫0t f(x)dx取f(t)=ε(t)有界, 但yf(t)= t*ε(t)无限。 ∴系统是不稳定的。 §1.7 LTI系统分析方法概述 系统分析的任务: 1、建立描述系统的数学方程式(微分/差分方程)。 2、给定初始状态和f(t),求响应y(t) ???描述系统的方法: ( 建立数学方程式 ) 输入输出法(外部法):建立f(t)与y(t)的直接关系, 考虑内部状态x(t),微分/差分方程,适用于单输入/出系统。 2、 状态变量法(内部法):状态方程和输出方程,适用于多输入/出系统。 ??输入输出分析法:(求解方程的方法) 时域分析 连续系统(2):求解微分方程;卷积积分. 离散系统(3):求解差分方程;卷积积分. 二、变换域分析: 连续系统傅里叶变换(4);频谱,频域分析。 卷积 连续系统拉式变换(5); 微分方程→ 代数方程。 ↓ 离散系统z变换(6); 差分方程→ 代数方程。 乘积 系统函数(7);决定系统的特性。 ??出发点:把复杂信号分解为基本信号之和。 ??基本信号:δ(t)、ε(t)、ejwt 、sinωt、est 等。 ??学习本课程的原则: 1、物理语言描述与数学语言描述并重; 2、信号分析与系统分析并重; 3、时域分析法与变换域分析法并重; 4、连续时间系统与离散时间系统并重;