动态系统稳定性分析【2】-信号的正交特性和复数表示

1. 信号的正交特性

1.1 信号的正交特性

函数$u(t)$和$v(t)$为实空间上以$2\pi$为周期的函数，两函数在周期$[-\pi ,\pi ]$上正交，可表示为如下形式：

$${\int }_{-\pi }^{\pi }u(t)v(t)dt=0$$

1.2 应用1（无线通信）：sinx和cosx的正交性

在无线通信等领域，sinx和cosx基底常常被用做最基础的信号单元，下述是一个简单的demo展示了两者的正交特性。
$${\int }_{-\pi }^{\pi }\sin (t)\cos \left(t\right)dt=0$$

% 定义积分区间
a = -pi;
b = pi;

% 计算sin(x)和cos(x)函数乘积的积分结果
integral_result = integral(@(x) sin(x).*cos(x), a, b);
disp(integral_result);

结果如下：

  -1.8215e-17

1.3 应用2（信号成分检验）

积分的和等于和的积分的条件： 由于积分区间相同，并且信号$f(t)$具有全正交成分，因此具有性质：
$${\int }_{-\pi }^{\pi }\left(f_1\left(t\right)f_k\left(t\right)+f_2\left(t\right)f_k\left(t\right)+...+f_j\left(t\right)\right)dt={\int }_{-\pi }^{\pi }{f}_1\left(t\right)f_k\left(t\right)dt+...+{\int }_{-\pi }^{\pi }{f}_j\left(t\right)dt$$
下面是一个简单例子：

% Time vector
t = linspace(-pi, pi, 1000);

% Define the signals
f1 = sin(t);
f2 = cos(t);
f3 = sin(2*t);
fk = sin(3*t);

% Compute the left-hand side of the equation
lhs = trapz(t, (f1.*fk + f2.*fk + f3.*fk));

% Compute the right-hand side of the equation
rhs = trapz(t, f1.*fk) + trapz(t, f2.*fk) + trapz(t, f3.*fk);

% Check if lhs and rhs are almost equal (upto some tolerance)
tol = 1e-10;
if abs(lhs - rhs) < tol
    disp('The signals are orthogonal.');
else
    disp('The signals are not orthogonal.');
end

结果如下：

The signals are orthogonal.

1.4 复数底 正交信号

1.4.1 复数形式

常用的复数表示形式有以下几种：

• 矩形形式：
  $$c=a+jb$$
• 三角形式：
  $$c=m\left(\cos \left(\varphi \right)+j\sin \left(\varphi \right)\right)$$
• 基于欧拉公式，进而可以获得其极坐标形式：
  $$C=M{e}^{j\varphi }$$

PS: 值得一提的是，欧拉虽然证明了实正余弦函数与复数的等价，但人们对于这个等价关系的理解如同现代的物理学家们不知道什么是电子，但能理解它的特性一样。（这句话非常有哲理，体现了目前我们人类对于物质认识的程度仅限于只知道现象，而并不知道其本质！） 因为我们人类对某些物理现象还没有一个合理的解释，因此我们不用担心j是什么，而只需要理解它们的作用就够了。对于我们来说，j算子表示把一个复数逆时针旋转90度。

1.4.2 为什么信号需要用复数表示？

首先，实际现实的信号是实信号的，但是在数学建模过程中通常处理的是复数形式的信号。

实信号，例如正弦波信号，是正负频率的共轭对称。例如cos(wt)=cos[(-w)t]，对cos(wt)进行傅里叶变换就会得到正负频率点上的信息，回忆一下cos(wt)的傅里叶变换，实部相等、虚部相反。

正负频率点都含有信息，信息就会有冗余，如果只保留正频率点的信息依然能恢复实信号，以达到简化的目的。可以考虑通过滤波器将负频率点给滤掉。这个滤波器就是希尔伯特变换。

实际上，不管是信号处理中的傅里叶变换还是希尔波特变换，等都涉及到了复实数的转换。

1.4.3 复数信号Python测试

代码如下：

% Define the complex signal representing a spiral
t = linspace(0, 4*pi, 1000); % Time vector
r = exp(0.1*t); % Magnitude of the spiral
theta = 5*t; % Phase of the spiral
complex_signal = r .* exp(1i*theta); % Complex signal with spiral behavior

% Create 3D plot to visualize the spiral in the complex plane
figure;
plot3(real(complex_signal), imag(complex_signal), t, 'b', 'LineWidth', 1.5);
xlabel('Real Part');
ylabel('Imaginary Part');
zlabel('Time');
title('Spiral in the Complex Plane with Time as the Third Dimension');
grid on;

结果如下：

参考资料

【1】知乎-信号处理专业必读：是复数但并不复杂的正交信号（万字长文）
【2】为什么信号要用复数表示？