短时分数阶傅里叶变换提取线性调频信号的脉冲周期matlab代码

下面是一个简单的 MATLAB 代码来提取线性调频信号的脉冲周期，使用短时分数阶傅里叶变换（S-FRFT）。

% 生成线性调频信号
fs = 1000; % 采样率
t = 0:1/fs:1; % 时间轴
f0 = 10; % 起始频率
f1 = 100; % 终止频率
x = chirp(t, f0, 1, f1);

% 生成脉冲周期
N = 256; % 每个窗口的长度
M = floor(length(x)/N); % 总窗口数
tau = zeros(M, 1); % 脉冲周期
for i = 1:M
    xi = x((i-1)*N+1:i*N); % 获取当前窗口信号
    alpha = 0.5; % 分数阶系数
    sf = frft(xi, alpha); % 短时分数阶傅里叶变换
    [pks, locs] = findpeaks(abs(sf)); % 找到峰值点
    [pks_sort, idx_sort] = sort(pks, 'descend'); % 按降序排序
    if length(pks_sort) >= 2 % 至少有两个峰值点
        tau(i) = abs(locs(idx_sort(1)) - locs(idx_sort(2))); % 计算脉冲周期
    end
end

% 绘制脉冲周期图像
t = 0:N/fs:(M-1)*N/fs;
figure;
plot(t, tau);
xlabel('时间 (s)');
ylabel('脉冲周期 (样本点)');

该代码中使用了 `chirp` 函数生成线性调频信号，使用了 `frft` 函数进行短时分数阶傅里叶变换，使用了 `findpeaks` 函数找到峰值点，并使用了简单的算法计算脉冲周期。