生成一段心电信号，并用matlab dwt 提取心电信号的P波

生成一段心电信号可以使用 Matlab 中的 `ecgsyn` 函数，提取心电信号的 P 波可以先进行小波变换，然后根据 P 波的特征选择相应的细节系数进行重构。

以下是一个示例代码，演示了如何生成一段心电信号并提取其中的 P 波：

% 生成心电信号
fs = 200; % 采样频率
t = 0:1/fs:5; % 时间序列
ecg = ecgsyn(fs*(length(t)-1), 1); % 生成心电信号

% 提取 P 波
[c, l] = wavedec(ecg, 5, 'db4'); % 进行小波变换
d4 = detcoef(c, l, 4); % 取出第四层细节系数
[pks,locs] = findpeaks(d4, 'MinPeakHeight', 0.1); % 寻找峰值
p_wave = zeros(size(ecg)); % 初始化 P 波
for i = 1:length(locs)
    idx = wcodemat(c(l(4)+1:l(3)), length(ecg)); % 将细节系数重构为长度与原信号相同的向量
    p_wave(locs(i):locs(i)+length(idx)-1) = idx; % 将第 i 个 P 波插入到信号中
end

% 绘制结果
figure;
subplot(2,1,1);
plot(t, ecg);
title('Original ECG');
subplot(2,1,2);
plot(t, p_wave);
title('P Wave');

在这个示例中，我们首先使用 `ecgsyn` 函数生成了一段模拟心电信号，然后使用小波变换进行信号分析，并通过寻找细节系数的峰值来确定 P 波的位置。接着，我们将每个 P 波的细节系数重构为一个与原信号长度相同的向量，并将它们插入到一个全零的向量中，最终得到了一个只包含 P 波的信号。最后，我们将原始信号和提取后的 P 波绘制在同一张图上进行比较。