五步移相法matlab程序

以下是五步移相法的MATLAB程序示例：

% 输入信号
fs = 1000; % 采样率
t = 0:1/fs:1-1/fs; % 时间向量
f = 50; % 信号频率
x = cos(2*pi*f*t); % 输入信号

% 五步移相法
N = length(x); % 信号长度
k = 0:N-1; % 频率索引
f = k*fs/N; % 频率向量
X = fft(x); % 傅里叶变换
phi = angle(X); % 相位谱
phi = unwrap(phi); % 相位展开
df = f(2) - f(1); % 频率分辨率
dt = 1/(N*df); % 时间分辨率
tau = phi/(2*pi*f); % 时延谱
tau = tau - min(tau); % 时延归零
t = 0:dt:(N-1)*dt; % 时延向量

% 显示结果
subplot(2,1,1);
plot(t,x); xlabel('时间/s'); ylabel('幅值');
title('输入信号');
subplot(2,1,2);
plot(tau,abs(X)); xlabel('时延/s'); ylabel('幅值');
title('时延谱');

在此示例中，我们首先生成了一个50Hz的正弦波，并使用五步移相法计算了其时延谱。结果显示在两个子图中：第一个子图显示了输入信号的时域波形，第二个子图显示了信号的时延谱。