信息与通信的数学基础——第三次小班

1. 卷积

计算$x(t)=\{\begin{array}{l}1,0<t<1\\ 0,else\end{array},h(t)=\{\begin{array}{l}t,0<t<1\\ 0,else\end{array}$的卷积

（1）定义函数

x[t_] := If[0 <= t && t < 1, 1, 0];
h[t_] := If[0 < t && t < 2, t, 0];
Plot[{x[t], h[t]}, {t, -2, 2}, PlotStyle -> {Thick, Dashed}]

（2）卷积计算

1. 定义计算

${\int }_{-\infty }^{+\infty }x(\tau )h(t-\tau )d\tau ={\int }_{-\infty }^{+\infty }x(t-\tau )h(\tau )d\tau$

f0 = Integrate[x[w]*h[t - w], {w, -Infinity, +Infinity}]
f1 = Integrate[x[t - w]*h[w], {w, -Infinity, +Infinity}]
Plot[{f0, f1}, {t, -3, 3}, PlotStyle -> {Thick, Dashed}]

2. 系统函数计算

f3 = Convolve[x[w], h[w], w, t]

3. 傅里叶变换

$h(t)\ast x(t)=F^{-1}(H(w)\cdot X(w))$

fx = FourierTransform[x[t], t, w]
fh = FourierTransform[h[t], t, w]
f = Sqrt[2*Pi]*InverseFourierTransform[fx*fh, w, t]
Plot[{f, f1}, {t, -3, 3}, PlotStyle -> {Thick, Dashed}]

问题：如何将两个图形打印在一副图形中？
使用Plot函数，设定两个函数图像，并设置PlotStyle分别为实现与虚线。

Plot[{f0, f1}, {t, -3, 3}, PlotStyle -> {Thick, Dashed}]

2. Laplace变换

求$f(t)=t^4,g(t)=e^t\sin t$得Laplace变换

（1）定义函数

f[t_] := t^4;
g[t_] := E^t*Sin[t];

（2）函数变换

2.1 函数正变换

l1 = LaplaceTransform[f[t], t, s]
l2 = LaplaceTransform[g[t], t, s]

2.2 函数逆变换

InverseLaplaceTransform[l1, s, t]
InverseLaplaceTransform[l2, s, t]

问题：结果的三角和指数表达

（1）三角转指数表达

TrigToExp[g[t]]

（2）全化简
全化简用于结果语句的化简。

InverseLaplaceTransform[l2, s, t]
FullSimplify[%]

3. Laplace变换解非齐次微分方程

$$\begin{gathered} x^{\prime \prime \prime }+3x^{\prime \prime }+3x^{\prime }+x=1 \\ {x}^{\prime \prime }(0)=x^{\prime }(0)=x(0)=0 \end{gathered}$$

（1）Laplace变换

$3\ast s^{2}F(s)+3\ast sF(s)+F(s)=\frac{1}{s}$

f1 = LaplaceTransform[x'''[t] + 3 x''[t] + 3 x'[t] + x[t], t, s]
f2 = LaplaceTransform[1, t, s]

（2）给出初始条件

X[0] = 0;
X''[0] = 0;
X'[0] = 0;

（3）解方程获取Laplace变换函数

Solve[f1 == f2, LaplaceTransform[X[t], t, s]]
LaplaceTransform[X[t], t, 
  s] /. {{LaplaceTransform[X[t], t, s] -> 1/(s (1 + s)^3)}}

（4）Laplace反变换求原函数

x = InverseLaplaceTransform[sol, s, t]
FullSimplify[x]

所有代码

ClearAll;
X[0] = 0;
X''[0] = 0;
X'[0] = 0;
f1 = LaplaceTransform[X'''[t] + 3*X''[t] + 3*X'[t] + X[t], t, s]
f2 = LaplaceTransform[1, t, s]
Solve[f1 == f2, LaplaceTransform[X[t], t, s]]
sol = LaplaceTransform[X[t], t, 
   s] /. {{LaplaceTransform[X[t], t, s] -> 1/(s (1 + s)^3)}}
x = InverseLaplaceTransform[sol, s, t]
FullSimplify[x]

例题

$$\{\begin{array}{l}{x}^{\prime \prime }(t)-2x^{\prime }(t)-y^{\prime }(t)+2y(t)=0\\ x^{\prime }(t)+y^{\prime }(t)-2x(t)=-2e^{-t}\end{array}$$

ClearAll;
f1 = LaplaceTransform[{X''[t] - 2 X'[t] - Y'[t] + 2 Y[t], 
   X'[t] + Y'[t] - 2 X[t]}, t, s]
f2 = LaplaceTransform[{0, -2*E^{-t}}, t, s]
X[0] = 2;
X'[0] = 3;
Y[0] = 1;
Solve[f1 == f2, {LaplaceTransform[X[t], t, s], 
  LaplaceTransform[Y[t], t, s]}]
sol = {LaplaceTransform[X[t], t, s], LaplaceTransform[Y[t], t, s]} /. %
x = InverseLaplaceTransform[sol, s, t]
FullSimplify[x]

4. Laplace变换求解积分方程

$a{t}^2-{\int }_0^t\cos (x-t)f(x)dx=f(t)$