把添加的特判修改到用MATLAB对心电信号进行小波变换，检测R波峰并计算平均RR间期的示例代码

好的，你可以将原来的代码修改为：

% 读取心电信号
ecg = load('100m.mat');
ecg = ecg.val(1:216000);

% 小波变换
wname = 'db6'; % 小波类型
level = 6; % 分解层数
[C, L] = wavedec(ecg, level, wname);
thr = wthrmngr('sw1ddenoLVL',C,L); % 计算阈值
sorh = 's'; % 软阈值去噪
xd = wden(C, L, thr, sorh, level, wname); % 去噪后的信号

% 检测R波峰
fs = 360; % 采样率
d1 = diff(xd); % 一阶差分
d2 = diff(d1); % 二阶差分
d3 = diff(d2); % 三阶差分
thr = 0.5 * max(d3); % 设置阈值
if length(d3) < 3
    % 处理样本数不足3个的情况
    disp('样本数不足3个');
else
    [pks, locs] = findpeaks(d3, 'MinPeakHeight', thr, 'MinPeakDistance', 0.2 * fs); % 检测模极大值
end

% 计算平均RR间期
rr = diff(locs) / fs; % 计算RR间期
avg_rr = mean(rr); % 计算平均RR间期
disp(['平均RR间期为：', num2str(avg_rr), '秒']);

这样，当样本数不足3个时，就会输出一条提示信息，而不会调用`findpeaks`函数，避免了错误的发生。