MATLAB信号与系统里拉普拉斯变换与极点和零点的绘制
-
什么是信号的拉普拉斯变换
我们假设输入线性时不变系统的复信号是 x ( t ) = e s 0 t , s 0 = σ + j Ω 0 x(t)=e^{s_{0}t},s_{0}=\sigma+j\varOmega_{0} x(t)=es0t,s0=σ+jΩ0,设系统的输出为 y ( t ) y(t) y(t),系统的冲激响应为 h ( t ) h(t) h(t),那么由卷积积分可以得出系统的输出为
y ( t ) = ∫ − ∞ ∞ h ( t ) x ( t − τ ) d τ = ∫ − ∞ ∞ h ( t ) e s 0 ( t − τ ) d τ y(t)=\int_{-\infty}^{\infty}h(t)x(t-\tau)d\tau=\int_{-\infty}^{\infty}h(t)e^{s_{0}(t-\tau)}d\tau y(t)=∫−∞∞h(t)x(t−τ)dτ=∫−∞∞h(t)es0(t−τ)dτ
= e s 0 ∫ − ∞ ∞ h ( τ ) e − τ s 0 d τ = x ( t ) H ( s 0 ) =e^{s_{0}}\int_{-\infty}^{\infty}h(\tau)e^{-\tau{s_{0}}}d\tau=x(t)H(s_{0}) =es0∫−∞∞h(τ)e−τs0dτ=x(t)H(s0)
我们把 x ( t ) x(t) x(t)称作线性时不变系统的特征函数。 H ( s 0 ) H(s_{0}) H(s0)通过冲激响应与系统相关联对于任意 h ( t ) h(t) h(t),我们把下式叫做 h ( t ) h(t) h(t)的拉普拉斯变换.
H ( s ) = ∫ − ∞ ∞ h ( τ ) e − τ s d τ H(s)=\int_{-\infty}^{\infty}h(\tau)e^{-\tau{s}}d\tau H(s)=∫−∞∞h(τ)e−τsdτ -
为什么要进行拉普拉斯变换
我们知道在时域中,信号经过系统的输出被定义为输入信号与冲激响应的卷积,而在时域中进行卷积非常复杂,而我们可以用拉普拉斯变换将卷积积分化为输入信号的拉普拉斯变换与冲激响应的拉普拉斯变换的乘积,进而简化了计算。即
y ( t ) = [ x ∗ h ] ( t ) = X ( s ) H ( s ) y(t)=[x*h](t)=X(s)H(s) y(t)=[x∗h](t)=X(s)H(s)
其中 X ( s ) X(s) X(s)是输入信号的拉普拉斯变换, H ( s ) H(s) H(s)为冲激响应的拉普拉斯变换。 -
系统函数和信号的拉普拉斯变换
我们把拉普拉斯变换定义为
F ( s ) = L [ f ( t ) ] = ∫ − ∞ ∞ f ( t ) e − s 0 t d t , s ∈ R O C F(s)=\mathscr{L}[f(t)]=\int_{-\infty}^{\infty}f(t)e^{-s_{0}t}dt,s\in{ROC} F(s)=L[f(t)]=∫−∞∞f(t)e−s0tdt,s∈ROC
R O C ROC ROC称为收敛域,是指此无限积分有意义的区域。
拉普拉斯逆变换定义为
f ( t ) = 1 2 π j ∫ σ − j ∞ σ + j ∞ F ( s ) e s t d s f(t)=\frac{1}{2\pi{j}}\int_{\sigma-j\infty}^{\sigma+j\infty}F(s)e^{st}ds f(t)=2πj1∫σ−j∞σ+j∞F(s)estds
如果这里的 f ( t ) f(t) f(t)是信号的话,那我们这个拉普拉斯变换称为信号的拉普拉斯变换,如果这个 f ( t ) f(t) f(t)是系统的冲激响应,也就是 h ( t ) h(t) h(t),那么我们把这个拉普拉斯变换称为系统函数或转移函数。 -
拉普拉斯变换的 M A T L A B MATLAB MATLAB计算
在计算信号的拉普拉斯变换时,我们一般使用符号计算,我们来看一段代码。
syms t
x=exp(-t);
y=x*cos(10*t)
X=laplace(x)
Y=laplace(y)
这段代码计算了 x 和 y x和y x和y的拉普拉斯变换,结果为
X =1/(s + 1)
Y =(s + 1)/((s + 1)^2 + 100)
这里需要注意一点,matlab里的拉普拉斯变换是单边拉普拉斯变换,所以我们在做相应的变换时需要考虑到这一点。
如果对于一个非因果信号进行拉氏变换。例如对
e
2
t
u
(
−
t
)
e^{2t}u(-t)
e2tu(−t)进行拉普拉斯变换。那么我们便不能使用
m
a
t
l
a
b
matlab
matlab里的
l
a
p
l
a
c
e
laplace
laplace函数,我们应该换一种方法,即
L
[
e
2
t
u
(
−
t
)
]
=
L
[
e
−
2
t
u
(
t
)
]
∣
−
s
\mathscr{L}[e^{2t}u(-t)]={\mathscr{L}[e^{-2t}u(t)]}|_{-s}
L[e2tu(−t)]=L[e−2tu(t)]∣−s