典型系统判定方法
线性系统
- 每一项都有输入变量 x
- 每一项输入变量的次幂都是一次
e.g. :
y(t)=2x(t)+3x(t−2) 、 y(t)=2tx(t) 、 y[n]=nx[n−3] 、$y(t)=x(t^2+3t+2)$是线性系统
但是 y(t)=2x(t)2 、 y(t)=x(t)+1 不是线性系统
时不变系统
- 时间变量 t 只能是输入变量的输入,即 t 只能在输入变量的括号里面
- 时间变量 t 的系数只能是1,而且次幂只能是1次
e.g. :
y(t)=3x(t)−x2(t) 、 y(t)=x(t)+1 等都是时不变系统
但是 y(t)=tx(t) 、 y[n]=x[t−3] 、 y[n]=x[−t] 、 y(t)=x(t2−1) 都不是时不变系统
因果系统
输出不依赖于未来输入,直接看对于某一时刻 t0 输入到系统后,对于$t\ge t0$ 的 y(t) 都有值,但是 $t<t0$ 时 y(t) 没有值(为0)
e.g. :
y(t)=x(t−3) 是因果系统
但是 y(t)=x(t+3) 、$y(t)=x(\frac{1}{2}t)$ 、 y(t)=x(3−t) 都不是因果系统
无记忆系统
一个系统的输出 y(t) 仅仅依赖当前时刻输入变量$x(t)$ 的值
e.g. :
y(t)=3x(t)−ex(t) 、 y[n]=x2[n]+3x[n]+sinx[n] 都是无记忆系统
但是 y(t)=x(t)+1 、 y(t)=dtdx(t) 、 y(t)=∫−∞tx(τ)dτ 等都是记忆系统
注意:无记忆系统必然是因果系统
可逆系统
输入变量能够唯一表达成关于输出变量的系统,即有无确定的反函数
e.g. :
y(t)=tx(t)+1 、 y[n]=∑i=−∞nx[i] 、 y(t)=∫−∞tx(τ)dτ 、 y(t)=x(21t−3) 都是可逆系统
但是 y(t)=dtdx(t) 、 y(t)=x2(t) 等这些都是不可逆系统(微分器的逆是积分,但是有个常数无法确定)
稳定系统
给一个系统输入有界的信号 x(t) 得到的输出 y(t) 也是有界的,那么该系统为稳定系统
e.g. :
y(t)=exp(x(t)) 、 y(t)=tx(t) 、 y[n]=x[n]−x[n−1]都是稳定系统
但是 y[n]=∑i=−∞nx[i] 、 y(t)=dtdx(t) 都是不稳定系统
单位冲激函数、单位延迟函数
对于冲击函数如此定义(数学定义):
这个就是数学上严格的定义,或者工程上可以看作某个矩形门函数在脉宽趋近0的逼近。
以下是 狄拉克 δ 函数(连续) 与 单位脉冲 δ[n](离散) 的完整性质汇总。
连续 δ(t) 的核心性质
| 性质 | 公式 | 说明 |
|---|
| 筛选(Sampling) | ∫−∞∞δ(t−t0)φ(t)dt=φ(t0) | φ(t) 在 t0 连续 |
| 单位元(卷积) | x(t)∗δ(t)=x(t) | 卷积的"1" |
| 时移(卷积) | x(t)∗δ(t−t0)=x(t−t0) | 搬移 t0 |
| 乘积 | f(t)δ(t−t0)=f(t0)δ(t−t0) | f 在 t0 连续 |
| 尺度 | δ(at)=∥a∥1δ(t)(a=0) | 冲激变窄则变高 |
| 对称 | δ(t)=δ(−t) | 偶函数 |
| 与 u(t) 关系 | ∫−∞tδ(τ)dτ=u(t),dtdu(t)=δ(t) | 微积分互逆 |
| 倒数(乘积) | δ′(t−t0)f(t)=f(t0)δ′(t−t0)−f′(t0)δ(t−t0) | 筛选特性 |
| 导数(冲激偶) | ∫−∞∞δ′(t−t0)φ(t)dt=−φ′(t0) | 负号来自分部积分 |
| n 阶导数 | ∫−∞∞δ(n)(t−t0)φ(t)dt=(−1)nφ(n)(t0) | |
| 卷积微分 | x(t)∗δ′(t)=x′(t) | 微分算子 |
| 复合函数 | δ(f(t))=i∑∥f′(ti)∥δ(t−ti) | ti 为 f(t)=0 的单根 |
| 能量/面积 | ∫−∞∞δ(t)dt=1 | 注意如果是 δ′(t) 积分为0 |
| 平方积分 | ∫−∞∞δ2(t)dt 无定义 | 不是普通函数,平方无意义 |
离散 δ[n] 的核心性质
δ[n] 是普通序列,不是广义函数。
| 性质 | 公式 | 说明 |
|---|
| 筛选 | n=−∞∑∞x[n]δ[n−k]=x[k] | 直接提取第 k 个样值 |
| 单位元(卷积) | x[n]∗δ[n]=x[n] | |
| 时移(卷积) | x[n]∗δ[n−k]=x[n−k] | 搬移 k 点 |
| 乘积 | x[n]δ[n−k]=x[k]δ[n−k] | |
| 与 u[n] 关系 | δ[n]=u[n]−u[n−1] | 一阶差分 |
| 累加 | u[n]=k=−∞∑nδ[k] | 累加 |
| Z 变换 | Z{δ[n]}=1 | ROC:全平面 |
| DTFT | F{δ[n]}=1 | 常数 1 |
| DFT | DFT{δ[n]}=1(对所有 k) | 全频均匀 |
以下是 u(t)(连续阶跃) 与 u[n](离散阶跃) 的完整性质汇总。
连续 u(t)
| 性质 | 公式 | 说明 |
|---|
| 定义 | u(t)={1,0,t>0t<0 | t=0 处通常定义为 1/2 或 1,视教材而定 |
| 与 δ(t) 关系 | dtdu(t)=δ(t),∫−∞tδ(τ)dτ=u(t) | 微积分互逆 |
| 积分(卷积) | x(t)∗u(t)=∫−∞tx(τ)dτ | 积分器 |
| 时移积分 | x(t)∗u(t−t0)=∫−∞t−t0x(τ)dτ | 从 t0 开始积分 |
| 乘积 | x(t)u(t−t0)={x(t),0,t>t0t<t0 | 截断/加窗 |
| 尺度 | u(at)={u(t),u(−t),a>0a<0 | 注意 a<0 时翻转 |
| 对称分解 | u(t)+u(−t)=1(t=0) | 与反阶跃互补 |
| 符号函数 | u(t)−u(−t)=sgn(t) | 或 2u(t)−1=sgn(t) |
| 拉氏变换 | L{u(t)}=s1 | Re{s}>0 |
| 双边拉氏 | L{−u(−t)}=s1 | Re{s}<0 |
| 傅里叶变换 | F{u(t)}=πδ(ω)+jω1 | 含直流冲激项 |
| u(t−t0) 的 FT | πδ(ω)+jωe−jωt0 | 时移 |
| u(t) 的导数 | u′(t)=δ(t) | 广义导数 |
| u(t) 的 n 阶导数 | u(n)(t)=δ(n−1)(t) | |
离散 u[n]
| 性质 | 公式 | 说明 |
|---|
| 定义 | u[n]={1,0,n≥0n<0 | n=0 处明确定义为 1 |
| 与 δ[n] 关系 | δ[n]=u[n]−u[n−1] | 一阶后向差分 |
| 累加 | u[n]=k=−∞∑nδ[k] | 累加器 |
| 差分(卷积) | x[n]∗δ[n]=x[n] | |
| 累加(卷积) | x[n]∗u[n]=k=−∞∑nx[k] | 累加器 |
| 时移累加 | x[n]∗u[n−k]=m=−∞∑n−kx[m] | 从 k 开始累加 |
| 乘积 | x[n]u[n−k]={x[n],0,n≥kn<k | 截断左边序列 |
| u[n] 的 Z 变换 | Z{u[n]}=1−z−11=z−1z | ∥z∥>1 |
| u[n] 的 DTFT | 1−e−jΩ1+πk=−∞∑∞δ(Ω−2πk) | 含周期冲激项 |
| −u[−n−1] 的 ZT | 1−z−11=z−1z | ∥z∥<1 |
| u[n] 的 DFT | 无闭式(无限长) | 需截断后求 |
| u[n] 的差分 | ∇u[n]=u[n]−u[n−1]=δ[n] | |
| u[n] 的 k 阶差分 | ∇ku[n] 产生有限长序列 | |
u(t) 与 u[n] 的关键区别
| 对比项 | 连续 u(t) | 离散 u[n] |
|---|
| t=0 / n=0 处 | 通常取 1/2 或 1(约定不一) | 严格为 1 |
| 导数/差分 | 广义导数为 δ(t) | 后向差分为 δ[n] |
| 积分/累加 | ∫−∞tu(τ)dτ=tu(t)=r(t)(斜坡) | ∑k=−∞nu[k]=(n+1)u[n] |
| 自卷积 | u(t)∗u(t)=tu(t) | u[n]∗u[n]=(n+1)u[n] |
| 傅里叶变换 | 含 δ(ω) | 含周期 δ(Ω−2πk) |
| 双边对应 | −u(−t)(反因果) | −u[−n−1](反因果) |
常用组合公式
| 表达式 | 结果 | 用途 |
|---|
| u(t−a)⋅u(t−b)(a<b) | u(t−b) | 取更晚的窗口 |
| u(t−a)+u(b−t)(a<b) | 在 (a,b) 上为 1,其余为 1 | 需配合理解 |
| u(t−a)−u(t−b)(a<b) | rect[a,b](t) | 门函数 |
| u[n−a]−u[n−b](a<b) | 在 n∈[a,b−1] 上为 1 | 离散门/矩形窗 |
| u[n]⋅u[n−k](k>0) | u[n−k] | 截断后右移 |
| anu[n] | 右边指数序列 | Z 变换最基础对 |
四大变换总表
连续时间
双边拉氏变换:必须标注收敛域 X(s)=∫−∞∞x(t)e−stdt
傅里叶变换:$\displaystyle X(j\omega)=\int_{-\infty}^{\infty}x(t)e^{-j\omega t}dt$ (即 s=jω 收敛域包含虚轴)
| 信号 x(t) | 傅里叶变换 X(jω) | 双边拉氏变换 X(s) | 收敛域 |
|---|
| δ(t) | 1 | 1 | 全平面 |
| δ(t−t0) | e−jωt0 | e−st0 | 全平面 |
| δ(n)(t) | (jω)n | sn | 全平面 |
| 1 | 2πδ(ω) | — | 不存在 |
| u(t) | πδ(ω)+jω1 | s1 | Re{s}>0 |
| −u(−t) | — | s1 | Re{s}<0 |
| tnu(t)(n≥0) | — | sn+1n! | Re{s}>0 |
| −tnu(−t) | — | sn+1n! | Re{s}<0 |
| e−atu(t)(a>0) | a+jω1 | s+a1 | Re{s}>−a |
| −e−atu(−t) | — | s+a1 | Re{s}<−a |
| e−a∥t∥(a>0) | a2+ω22a | a2−s22a | −a<Re{s}<a |
| ejω0t | 2πδ(ω−ω0) | — | 不存在 |
| cos(ω0t) | π[δ(ω−ω0)+δ(ω+ω0)] | — | 不存在 |
| sin(ω0t) | jπ[δ(ω−ω0)−δ(ω+ω0)] | — | 不存在 |
| cos(ω0t)u(t) | — | s2+ω02s | Re{s}>0 |
| sin(ω0t)u(t) | — | s2+ω02ω0 | Re{s}>0 |
| e−atcos(ω0t)u(t) | (a+jω)2+ω02a+jω | (s+a)2+ω02s+a | Re{s}>−a |
| e−atsin(ω0t)u(t) | (a+jω)2+ω02ω0 | (s+a)2+ω02ω0 | Re{s}>−a |
| sgn(t) | jω2 | s2 | Re{s}>0 |
| rect(τt) | τSa(2ωτ) | — | — |
| Sa(ω0t) | ω0πrect(2ω0ω) | — | — |
离散时间
DTFT:周期 2π X(ejω)=n=−∞∑∞x[n]e−jωn
z 变换:必须标注收敛域 X(z)=n=−∞∑∞x[n]z−n
关系: X(ejΩ)=X(z)z=ejΩ (要求收敛域包含单位圆)
离散时间傅里叶变换(DTFT)
| 信号 x[n] | DTFT X(ejω) |
|---|
| δ[n] | 1 |
| δ[n−k] | e−jωk |
| 1 | 2πk=−∞∑∞δ(ω−2πk) |
| u[n] | 1−e−jω1+πk=−∞∑∞δ(ω−2πk) |
| anu[n](∥a∥<1) | 1−ae−jω1 |
| −anu[−n−1](∥a∥>1) | 1−ae−jω1 |
| (n+1)anu[n](∥a∥<1) | (1−ae−jω)21 |
| ejω0n | 2πk=−∞∑∞δ(ω−ω0−2πk) |
| cos(ω0n) | πk=−∞∑∞[δ(ω−ω0−2πk)+δ(ω+ω0−2πk)] |
| sin(ω0n) | jπk=−∞∑∞[δ(ω−ω0−2πk)−δ(ω+ω0−2πk)] |
| 对称矩形窗 u[n+N1]−u[n−N1−1] | sin(ω/2)sin[ω(N1+21)] |
| πnsin(ωcn) | rect(2ωcω)(周期 2π) |
| a∥n∥(∥a∥<1) | 1−2acosω+a21−a2 |
z 变换
| 信号 x[n] | z 变换 X(z) | 收敛域 |
|---|
| δ[n] | 1 | 全平面 |
| δ[n−k] | z−k | 全平面($k>0$ 时 z=0) |
| u[n] | 1−z−11=z−1z | ∥z∥>1 |
| −u[−n−1] | 1−z−11=z−1z | ∥z∥<1 |
| anu[n] | 1−az−11=z−az | ∥z∥>∥a∥ |
| −anu[−n−1] | 1−az−11=z−az | ∥z∥<∥a∥ |
| nanu[n] | (1−az−1)2az−1=(z−a)2az | ∥z∥>∥a∥ |
| −nanu[−n−1] | (1−az−1)2az−1=(z−a)2az | ∥z∥<∥a∥ |
| n2anu[n] | (1−az−1)3az−1(1+az−1) | ∥z∥>∥a∥ |
| cos(ω0n)u[n] | 1−2cosω0z−1+z−21−cosω0z−1 | ∥z∥>1 |
| sin(ω0n)u[n] | 1−2cosω0z−1+z−2sinω0z−1 | ∥z∥>1 |
| ancos(ω0n)u[n] | 1−2acosω0z−1+a2z−21−acosω0z−1 | ∥z∥>∥a∥ |
| ansin(ω0n)u[n] | 1−2acosω0z−1+a2z−2asinω0z−1 | ∥z∥>∥a∥ |
四大变换性质对照
| 性质 | 傅里叶变换 | 双边拉氏变换 | DTFT | z 变换 |
|---|
| 线性 | ax1+bx2 | aX1(s)+bX2(s) | aX1+bX2 | aX1(z)+bX2(z) |
| 时移/位移 | x(t−t0)↔e−jωt0X | x(t−t0)↔e−st0X(s) | x[n−k]↔e−jΩkX(ejΩ) | x[n−k]↔z−kX(z) |
| 频移/指数加权 | ejω0tx(t)↔X(j(ω−ω0)) | es0tx(t)↔X(s−s0) | ejΩ0nx[n]↔X(ej(Ω−Ω0)) | anx[n]↔X(a−1z) |
| 尺度/反转 | x(at)↔∥a∥1X(jaω) | x(at)↔∥a∥1X(as) | x[−n]↔X(e−jΩ) | x[−n]↔X(z−1) |
| 共轭 | x∗(t)↔X∗(−jω) | x∗(t)↔X∗(s∗) | x∗[n]↔X∗(e−jΩ) | x∗[n]↔X∗(z∗) |
| 卷积 | x1∗x2↔X1⋅X2 | x1∗x2↔X1⋅X2 | x1∗x2↔X1⋅X2 | x1∗x2↔X1⋅X2 |
| 乘积/调制 | x1⋅x2↔2π1X1∗X2 | x1⋅x2↔2πj1X1∗X2 | x1⋅x2↔2π1X1⊛X2 | x1⋅x2↔2πj1X1∗X2 |
| 微分/差分 | dtdx↔jωX | dtdx↔sX(s) | — | x[n]−x[n−1]↔(1−z−1)X(z) |
| 频域微分 | tx(t)↔jdωdX | tx(t)↔−dsdX | nx[n]↔jdΩdX | nx[n]↔−zdzdX |
| 积分/累加 | ∫−∞txdτ↔πX(0)δ(ω)+jωX | ∫−∞txdτ↔sX(s) | ∑k=−∞nx[k]↔1−e−jΩX+πX(1)∑δ | ∑k=−∞nx[k]↔1−z−11X(z) |
| Parseval | ∫∥x∥2dt=2π1∫∥X∥2dω | — | ∑∥x∥2=2π1∫−ππ∥X∥2dΩ | ∑∥x∥2=2πj1∮∥X∥2zdz |
NOTE
- 对称性:
- 实偶 对 实偶:实偶函数变换后也是实偶的(充要条件),此时傅里叶级数不含正弦项
- 虚奇 对 实奇:纯虚奇函数变换后是纯实奇函数,此时傅里叶级数不含余弦项
- 实连续信号的傅里叶变换满足 X(−jω)=X∗(jω),幅度谱为偶函数,相位谱为奇函数
- 单边信号的导数/差分
- 连续信号:$x(t)\leftrightarrow X(s)$,则 x(n)(t)↔snX(s)−sn−1x(0−)−sn−2x′(0−)−⋯−x(n−1)(0−)
- 离散信号:$x[n-m]\leftrightarrow z^{-m}X(z)+z^{-m}\sum_{k=-m}^{-1}x[k]z^{-k}$
- 信号稳定性的规律:
- 对于连续信号,其拉氏变换的极点全部位于左半平面时,信号稳定
- 对于离散信号,其 z 变换的极点全部位于单位圆内时,信号稳定
- FT 与 DTFT 的存在条件:
- FT 存在:ROC 包含虚轴(连续)或绝对可积
- DTFT 存在:ROC 包含单位圆(离散),即 ∑∣x[n]∣<∞(绝对可和)
- u[n] 的 DTFT:含冲激项是因为 u[n] 不是绝对可和的,其 DTFT 在 ω=0 处有直流分量。
收敛域规律:
- 右边序列:收敛域在最外极点圆外(z)或最右极点右侧(s)
- 左边序列:收敛域在最内极点圆内(z)或最左极点左侧(s)
- 双边序列:收敛域为圆环(z)或带状(s)
傅里叶变换的周期延拓
我们知道单个傅里叶变换有
那么我们可以通过以下变换将离散和连续联系起来
以 T 为周期无限重复(周期延拓)
设周期延拓信号:
其中冲激串(Shah 函数)$\displaystyle \text{III}T(t)=\sum{k=-\infty}^{\infty}\delta(t-kT)$ 的傅里叶变换为:
由时域卷积 ↔ 频域乘积:
物理意义:周期信号的频谱是离散冲激串,冲激位于谐波频率 ω=kω0(ω0=2π/T) 处,权重正比于原频谱 X(jω) 在该点的采样值。
这里对于已经扩展的信号(或者本身就是周期信号)$x(t)$ 的傅里叶级数的系数
注意这里求出 ak 后,得到傅里叶级数
但是有些老师认为这样没有化简完成,根据傅里叶变换的对称性质,实偶函数的傅里叶级数不含正弦项,这里还需要用欧拉公式化简。
理想采样(时域乘积)
采样信号:
由时域乘积 ↔ 频域卷积(系数 1/2π):
物理意义:采样使频谱以 ωs=2π/T 为周期无限延拓(幅值缩放 1/T),若 X(jω) 带宽超过 π/T 就会发生混叠。
离散序列 x[n]=x(nT) 的傅里叶变换(DTFT)
定义(离散数字频率 ω,单位:rad/sample):
对比理想采样的频谱(连续角频率记为 ωc,单位:rad/s):
令 ω=ωcT(即 ωc=ω/T),立得:
总结对照表
| 操作 | 时域表达式 | 频域结果 |
|---|
| 原信号 | x(t) | X(jωc) |
| 周期延拓(周期 T ) | k=−∞∑∞x(t−kT) | T2πk=−∞∑∞X(jT2πk)δ(ωc−T2πk) |
| 理想采样(周期 T ) | k=−∞∑∞x(kT)δ(t−kT) | T1k=−∞∑∞X(j(ωc−T2πk)) |
| 离散序列 DTFT | x[n]=x(nT) | X(ejω)=T1k=−∞∑∞X(jTω−2πk) |
符号说明:表中 X(ejω) 的 ω 是离散数字频率;$X(j\cdot)$ 括号内是连续角频率(rad/s),因此出现 ω/T 的换算。
罗斯阵列
罗斯阵列(Routh-Hurwitz 判据) 是不解特征方程、不直接求极点的情况下,判断系统是否稳定(即所有极点是否都在左半平面)的代数方法。
而且对于含参的 H(s) 罗斯阵列几乎是最快的求参数何时具有稳定的方法
罗斯阵列构造
罗斯阵列的构造是机械的"两行交叉相除",核心规则只有一条:
第 1 步:摆前两行(直接抄系数)
sn 行:$an,\ a{n-2},\ a{n-4},\ a{n-6},\ \ldots$(跳一个取一个)
sn−1 行:$a{n-1},\ a{n-3},\ a{n-5},\ a{n-7},\ \ldots$(跳一个取一个)
若某位置没有系数,补 0。
第 2 步:从第 3 行开始,逐行计算
通用公式(计算新行的第 j 个元素):
例子
如有特征方程:
$$
D(s)=a4s^4+a3s^3+a2s^2+a1s+a_0=0
$$
必要条件(先快速排除):所有系数同号且不缺项,稳定。若有系数为 0 或异号,一定不稳定。
罗斯表:
| s4 | a4 | a2 | a0 |
|---|
| s3 | a3 | a1 | 0 |
| s2 | b1=a3a3a2−a4a1 | b2=a3a3a0−a4⋅0=a0 | |
| s1 | c1=b1b1a1−a3b2 | | |
| s0 | d1=c1c1b2−b1⋅0=b2 | | |
判据:
- 第一列 {a4,a3,b1,c1,d1} 全部同号(通常全为正)→ 系统稳定
- 第一列符号变化次数 = 右半平面极点个数
两种特殊情况处理
1. 某行第一列元素为 0,但其余不全为 0
用无穷小正数 ε 代替 0 继续计算,最后令 ε→0+ 看符号。
2. 某行全部为零
说明存在关于原点对称的根(共轭虚根或相反实根)。用上一行构造辅助多项式 A(s),对其求导后系数替换全零行,继续计算。辅助多项式的根也是原方程的根。