出错 untitled (第 9 行) y = y + phaseretard(x, d(i)*phi, fs); % 添加相位延迟

这个错误可能是由于`phaseretard`函数没有被识别，如果你没有添加信号处理工具箱，那么这个函数可能不可用。你可以尝试手动编写一个相位延迟函数，以下是一个简单的相位延迟函数的示例代码：

function y = phase_delay(x, phi, fs)
% x: 输入信号
% phi: 相位延迟（弧度）
% fs: 采样率
d = phi/(2*pi*fs); % 延迟时间（单位：秒）
n = 1:length(x); % 时间向量
delay = round(d*fs); % 延迟样本数
y = zeros(size(x)); % 初始化输出信号
if delay >= 0
    y(delay+1:end) = x(1:end-delay);
else
    y(1:end+delay) = x(-delay+1:end);
end
end

你可以将这个函数添加到你的代码中，然后使用以下代码替换原来的代码行：

y = y + phase_delay(x, d(i)*phi, fs); % 添加相位延迟

这个函数的原理是将输入信号按照一定的延迟时间进行移位，然后用零填充未移位的位置。如果你仍然遇到问题，请尝试查看MATLAB文档或联系MATLAB技术支持。