请写出心电信号获得希尔伯特边际谱的方法和代码

心电信号获得希尔伯特边际谱的方法：

1. 首先，对心电信号进行希尔伯特变换，得到其解析信号。
2. 对解析信号进行快速傅里叶变换，得到频域信号。
3. 计算频域信号的振幅谱，即幅度平方，得到边际谱。

代码：

import numpy as np
from scipy.signal import hilbert

# 生成心电信号
fs = 1000  # 采样率
T = 10  # 时长
t = np.linspace(0, T, T*fs, endpoint=False)
ecg = np.sin(2*np.pi*60*t) + 0.5*np.sin(2*np.pi*120*t) + 0.2*np.sin(2*np.pi*240*t)

# 希尔伯特变换
analytic_signal = hilbert(ecg)

# 快速傅里叶变换
fft_signal = np.fft.fft(analytic_signal)

# 计算边际谱
amplitude_spectrum = np.abs(fft_signal)**2

# 绘制结果
import matplotlib.pyplot as plt

freq = np.fft.fftfreq(len(t), 1/fs)
plt.plot(freq[:len(freq)//2], amplitude_spectrum[:len(freq)//2])
plt.xlabel('Frequency (Hz)')
plt.ylabel('Amplitude Spectrum')
plt.show()

上述代码中，我们首先生成了一个包含三个频率成分的心电信号，然后使用 `hilbert` 函数进行希尔伯特变换，得到解析信号。接着使用 `np.fft.fft` 函数进行快速傅里叶变换，得到频域信号。最后计算频域信号的振幅谱，即幅度平方，得到边际谱，并绘制了其频谱图。