在该代码的基础上继续编写，找出t波位置并画图

在该代码的基础上，可以通过以下步骤来找出T波位置并画图：

1. 在R波位置向量rvalue的基础上，计算R波到R波之间的距离（即RR间期），并计算平均RR间期。

2. 在平均RR间期的基础上，计算T波检测窗口的长度。例如，可以将T波检测窗口设置为平均RR间期的一半。

3. 在R波位置向量rvalue的基础上，遍历每个R波位置，计算对应的T波位置。可以通过以下方法来计算T波位置：

a. 在当前R波位置的基础上，向后搜索T波检测窗口的长度。

b. 在搜索窗口内，搜索信号的最小值，以确定T波位置。可以使用MATLAB中的min函数来实现。

c. 如果在搜索窗口内未发现明显的最小值，则可以继续向后搜索，直到找到为止。

4. 将T波位置存储为向量或矩阵，以供后续分析和处理使用。

5. 在原始信号上标记出R波和T波位置，并进行画图展示。

下面是代码的具体实现：

% 计算平均RR间期
RR_interval = diff(rvalue);
mean_RR_interval = mean(RR_interval);

% 计算T波检测窗口长度
T_window_length = round(mean_RR_interval / 2);

% 初始化T波位置向量
tvalue = zeros(size(rvalue));

% 搜索T波位置
for i = 1:length(rvalue)
    % 计算搜索范围
    search_range = round(T_window_length * 1.5); % 向后搜索1.5倍的T波检测窗口长度
    
    % 截取信号
    signal_segment = x(rvalue(i):rvalue(i)+search_range);
    
    % 搜索最小值
    [~, tvalue_index] = min(signal_segment);
    
    % 计算T波位置
    tvalue(i) = rvalue(i) + tvalue_index - 1;
end

% 画图展示R波和T波位置
figure
plot((1:length(x))/fs, x); hold on;
plot(rvalue/fs, x(rvalue), 'r^', 'MarkerFaceColor', 'r');
plot(tvalue/fs, x(tvalue), 'bv', 'MarkerFaceColor', 'b');
xlabel('Time (s)'); ylabel('Amplitude');
legend('ECG signal', 'R wave', 'T wave');

需要注意的是，这里假设心电信号数据已经读取并存储在变量x中，R波位置存储在向量rvalue中，采样率为1000Hz。在实际应用中，需要根据具体情况进行调整和优化。