五步相移相位提取算法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_diff = diff(phi); % 相位差
phi_diff = [phi_diff(1), phi_diff]; % 补充首项
omega = 2*pi*f; % 角频率
tau = phi_diff./(omega(2)-omega(1)); % 时延谱
tau = tau - min(tau); % 时延归零
t = 0:1/fs:(N-1)/fs; % 时延向量

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

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