matlab写个离散时间序列希尔博特变换程序

好的，下面是一个简单的 MATLAB 程序，用于计算离散时间序列的 Hilbert 变换：

function [ht, ht_abs, ht_phase] = hilbert_transform(x)
% 输入参数：
% x: 输入的离散时间序列，行向量或列向量。
% 输出参数：
% ht: Hilbert 变换后的复数序列。
% ht_abs: Hilbert 变换后的幅度序列。
% ht_phase: Hilbert 变换后的相位序列。

% 离散时间序列的长度
N = length(x);

% 计算输入序列的傅里叶变换
X = fft(x);

% 构造 Hilbert 变换的频域滤波器
H = zeros(1, N);
if mod(N, 2) == 0  % N 是偶数
    H(1) = 1;
    H(2:N/2) = 2*ones(1, N/2-1);
    H(N/2+1) = 1;
    H(N/2+2:N) = zeros(1, N/2-1);
else  % N 是奇数
    H(1) = 1;
    H(2:(N+1)/2) = 2*ones(1, (N-1)/2);
    H((N+3)/2:N) = zeros(1, (N-1)/2);
end

% 对输入序列的傅里叶变换进行频域滤波
HX = H .* X;

% 计算 Hilbert 变换后的复数序列
ht = ifft(complex(real(HX), imag(HX)));

% 计算 Hilbert 变换后的幅度序列和相位序列
ht_abs = abs(ht);
ht_phase = angle(ht);
end

使用时，只需要传入一个离散时间序列即可，例如：

% 生成一个长度为 100 的随机离散时间序列
x = randn(1, 100);

% 计算 Hilbert 变换
[ht, ht_abs, ht_phase] = hilbert_transform(x);

其中 `ht` 是 Hilbert 变换后的复数序列，`ht_abs` 是 Hilbert 变换后的幅度序列，`ht_phase` 是 Hilbert 变换后的相位序列。