《信号与系统——MATLAB综合实验》讲义_第二讲

中x, y表示位置，vx , vy表示速度，可知运动方程为
⎡
vx
⎤
⎢ ⎡ x ⎤ ⎢
vy
⎥ ⎥
˙
s
=
⎢ ⎢
y
⎢⎢⎢⎣vvxy
⎥ ⎥ ⎥ ⎥ ⎥⎦
=
⎢ ⎢−GM E
⎢
⎢
⎢ ⎢
−GM
E
(x2
x
3
+ y2)2
y
3
⎥ ⎥ ⎥ ⎥ ⎥ ⎥
⎢⎣
(x2 + y2 )2 ⎥⎦
其中引力常数G = 6.672×10−11，地球 质量M e = 5.97 ×1024，假设初值x(0) = −4.2×107 , y(0) = 0, vx (0) = 0, vy (0) = 4×103。
4.2 周期信号的傅里叶级数分析
z 和傅里叶变换的数值计算方法相似
例4.3 绘制周期T1 = 1，幅度E = 1的对称方波的前10项傅里叶 级数的系数（三角函数形式），并用前5项恢复原信号。
⎢ ⎣
F
(ω2
−
Δω
)
⎥ ⎦
⎣⎢e− j(ω2 −Δω )t1
e− jω1 (t1 +Δt ) e− j(ω1 +Δω )(t1 +Δt )
# e− j(ω2 −Δω )(t1 +Δt )
"
e ⎤ ⎡ − jω1(t2 −Δt) f (t1) ⎤
"
e ⎥⎥ ⎢⎢ − j(ω1+Δω )(t2 −Δt )
2008-8-13
谷源涛@清华大学电子工程系
7
3.3 零输入响应与零状态响应
z lsim 函数还可以对带有非零起始状态的LTI 系统 进行仿真，y = lsim(sys, u, t, x0)
z 注意：若用lsim 函数仿真非零起始状态响应，则 该系统必须用状态方程描述（对传递函数描述的 LTI 系统，lsim函数无法仿真非零起始状态响 应）。
+ kΔω)
=
T N
N −1 n=0
f
t1 + nΔt e− j(ω1+kΔω)(t1+nΔt)
z 同理写出逆变换
2008-8-13
∫ ∑ f (t) = 1
2π
∞ F (ω)ejωtdω
−∞
=
Ω
2π K
K −1
F (ω1
k =0
+ kΔω)e j(ω1+kΔω)t
∑ f
(t1
+
nΔt)
=
Ω
2π K
0]
⎢⎢λ2
⎥ ⎥
+
[0
⎢⎣λ3 ⎥⎦
1]
⎡⎢⎣ee12
(t (t
) )
⎤ ⎥ ⎦
2008-8-13
谷源涛@清华大学电子工程系
5
3.2 微分方程式的建立与求解
z 微分方程的解包括齐次解和特解两部分。齐次解 即系统特征方程的根，用roots函数计算。
例3.3
求微分方程
d3 dt 3
r
(t
)
+
7
d2 dt 2
z 第二种：构造一个激励信号，先保持2V足够长时间再 跳变为4V，然后即可以零初始状态一次仿真得到系统 的完全响应r1。
2008-8-13
谷源涛@清华大学电子工程系
10
知识点（2）解微分方程组
例3.6 绘制地球卫星的运行轨道。以卫星轨道为坐标平面，地球
位置为坐标原点，定义卫星的二维运动状态s = ⎡⎣x, y, vx , vy ⎤⎦T ，其
w(n) = ∑u(m)v(n +1− m) m
2008-8-13
谷源涛@清华大学电子工程系
14
3.5 卷积
例3.7 已知某系统冲激响应为h (t ) = t 2, 0 < t < 2。以- 0.5开始， 宽度为1.5秒，幅度为1的矩形脉冲e(t ) 激励该系统，求输出信 号r (t )。
2008-8-13
谷源涛@清华大学电子工程系
15
知识点（4）连续时间信号的数值 计算
z 自然界中绝大部分物理量是连续的，数字计算机 处理连续信号必须做近似，计算精度取决于算法 和数字表示的位数。
z 例如积分运算的简单数值计算方法是分段求和
∫ ∑ F (a,b) =
b a
f
(t )dt
≈
b−a N
N −1 n=0
f
⎛ ⎜⎝
2008-8-13
谷源涛@清华大学电子工程系
8
3.3 零输入响应与零状态响应
例3.5 给定下图所示电路，t < 0开关S处于1的位置，而且已经
达到稳态，将其看做起始状态，当t = 0时，S由1转向2。分别
求t > 0时i (t )的零输入响应和零状态响应。
R1 = 1Ω
i(t)
iL (t )
e(t) = 4V
的频谱恢复时域信号fs (t)，比较和原信号f (t)的差别。
z 三种方法
z 二重循环 z 循环＋矢量相乘 z 矩阵相乘
2008-8-13
谷源涛@清华大学电子工程系
20
矢量计算法
⎡ e ⎤ ⎡ − j(ω1+kΔω )t1
T
f (t1) ⎤
F (ω1
+
kΔω)
=
T N
⎢⎢e− j(ω1 +kΔω )(t1 +Δt )
例3.2 描述如下系统
状态方程
⎡⎢⎢λλ12 ⎢⎣λ3
⎤ ⎥ ⎥ ⎥⎦
=
⎡−2
⎢ ⎢
0
⎢⎣ 2
0 −3 −2
−1⎤ ⎡λ1 ⎤ ⎡1
3
⎥ ⎥
⎢⎢λ2
⎥ ⎥
+
⎢⎢0
0 ⎥⎦ ⎢⎣λ3 ⎥⎦ ⎢⎣0
0⎤ −3⎥⎥ 0 ⎥⎦
⎡⎢⎣ee12
(t (t
)⎤ )⎥⎦
输出方程
r (t) = [0
1
⎡λ1 ⎤
r
(t
)
+ 16
d dt
r
(t
)
+
12r
(t
)
=
0的齐次解。
z 特解即系统（采用微分方程表示）在给定信号激 励下的输出。用lsim函数进行仿真。
例3.4
给定微分方程
d2r (t
dt 2
)
+
2
dr (t
dt
)
+
3r
(t
)
=
de ( t
dt
)
+
e
(t
)，
如果已知(1) e(t ) = t2；(2) e(t ) = et，分别求两种情况下此方程的特解。
K −1
F (ω1
k =0
ω + k Δ )e j(ω1+kΔω )(t1+nΔt )
谷源涛@清华大学电子工程系
19
4.1 傅里叶变换
例4.2请绘制矩形脉冲
f
(t)
=
⎧⎪1 ⎨
t <1 2
⎪⎩0 otherwise
的波形(t ∈[−1,1])和频谱F (ω)(ω ∈[−8π ,8π ])，并利用计算得到
学编程 学应用
谷源涛 清华大学电子工程系
课程内容
2008-8-13
谷源涛@清华大学电子工程系
2
3. 连续时间系统的时域分析
z 引言 z 微分方程式的建立与求解 z 零输入响应与零状态响应 z 冲激响应与阶跃响应 z 卷积
2008-8-13
谷源涛@清华大学电子工程系
3
3.1 引言
z 连续时间系统的研究方法包括输入－输出法（端 口描述法）和系统状态变量分析法。
e j(ω1 +Δω )t1 e j(ω1 +Δω )(t1 +Δt )
# e j(ω1 +Δω )(t2 −Δt )
" e ⎤ ⎡ j(ω2 −Δω)t1 F (ω1) ⎤
"
e
j(ω2
−Δω
)(t1 +Δt )
⎥ ⎥
⎢ ⎢
F
(ω1
+
Δω)
⎥ ⎥
%
# ⎥⎢ # ⎥
( ) "
e
j(ω2
−Δω
)(t2
−Δt
13
3.5 卷积
z 卷积运算的数值近似
∫ f (t) =
∞ −∞
f1(τ )
f2 (t
−τ )dτ
∑ ∫ ∞
f (nT ) =
m=−∞
mT +T mT
f1(τ ) f2 (nT −τ )dτ
∑ f (nT ) ≈ T f1(mT ) f2 (nT − mT ) m
z MATLAB 提供了w = conv(u,v) 函数实现卷积和
2008-8-13
C = 1F
vC (t )
L= 1H 4
e(t) = 2V
R2
=
3Ω 2
谷源涛@清华大学电子工程系
9
3.3 零输入响应与零状态响应
z 解：采用两种方法
z 第一种：首先仿真2V电压en作用足够长时间（10秒） 后系统进入稳态，从而得到稳态的状态变量值x0，再以 其作为初始值仿真4V电压e作用下的输出rf，即是系统 的完全响应。这种方法可以得到零状态响应rzs和零输 入响应rzi。
(ω1
+
Δω
)
⎥ ⎥
⎢#⎥
( ) ⎢⎣e ⎥⎦ j(ω2 −Δω)(t1+nΔt)