基于积分公式手动实现正弦函数Hilbert变换的Matlab问题咨询

代码的核心错误 #

你当前的代码存在几个致命问题，导致无法运行：

• 未定义omega和t变量，integral调用缺少闭合括号
• 积分变量与时间轴变量同名（都用t），导致逻辑混乱
• 直接积分得到的是标量，无法与时间数组t绘图
• 未处理积分的奇点（当t=τ时分母为0），需要用柯西主值积分

关于你困惑的相位偏移相关符号 #

Hilbert变换的本质是对实信号的频谱进行相位调制：正频率分量加90°（π/2）相位，负频率分量减90°相位。对于正弦信号sin(ωt)，其Hilbert变换的结果是cos(ωt)——这是相位偏移π/2的结果，和周期2π无关，你之前代入2π是对符号含义的误解。

正确的手动实现代码 #

以下是可运行的手动实现，包含柯西主值积分处理、时间轴定义和绘图：

% 1. 定义基础参数
omega = 2*pi;          % 正弦信号的角频率
t = linspace(-5, 5, 1000);  % 时间轴，覆盖5个周期
hilbert_imag = zeros(size(t));  % 存储Hilbert变换结果（解析信号虚部）

% 2. 逐点计算Hilbert变换（柯西主值积分）
for idx = 1:length(t)
    t_current = t(idx);
    % 被积函数：积分变量为tau，避免与时间t冲突
    integrand = @(tau) sin(omega*tau) ./ (t_current - tau);
    % 计算柯西主值积分，处理tau=t_current的奇点
    hilbert_imag(idx) = (1/pi) * integral(integrand, -Inf, Inf, 'PrincipalValue', true);
end

% 3. 绘制解析信号的实部（原信号）和虚部（Hilbert变换）
figure;
plot(t, sin(omega*t), 'b-', t, hilbert_imag, 'r-', 'LineWidth', 1.5);
legend('解析信号实部：sin(ωt)', '解析信号虚部：Hilbert变换结果');
xlabel('时间 t');
ylabel('幅值');
grid on;

关键原理说明 #

Hilbert变换的积分公式是柯西主值积分：
H[x(t)] = \frac{1}{\pi} \text{P.V.} \int_{-\infty}^{\infty} \frac{x(τ)}{t - τ} dτ
其中P.V.表示取主值，用来绕过t=τ处的奇点。对于正弦信号，积分结果恰好是相位偏移π/2的余弦信号，和Matlab自带hilbert函数的虚部结果一致。

内容的提问来源于stack exchange，提问作者newbie-mat