短时分数阶傅里叶变换提取线性调频信号的脉冲宽度matlab

短时分数阶傅里叶变换（S-fracFT）是一种用于分析非平稳信号的工具，可以提取信号中的瞬态特征。下面是使用S-fracFT提取线性调频信号脉冲宽度的matlab代码示例：

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

% 设置S-fracFT参数
alpha = 0.5; % 分数阶阶数
wlen = 100; % 窗口长度
hop = 50; % 帧移
nfft = 1024; % FFT点数

% 计算S-fracFT
[sfracft,f,t] = sfracft(x,alpha,wlen,hop,nfft,fs);

% 提取脉冲宽度
[~,idx] = max(abs(sfracft),[],1); % 找到每帧最大值的位置
pw = zeros(1,length(t)); % 初始化脉冲宽度向量
for i = 1:length(t)
    k = idx(i); % 最大值位置
    while (k > 1) && (abs(sfracft(k)) > 0.1*abs(sfracft(idx(i)))) % 找到左边界
        k = k - 1;
    end
    l = k;
    k = idx(i);
    while (k < length(f)) && (abs(sfracft(k)) > 0.1*abs(sfracft(idx(i)))) % 找到右边界
        k = k + 1;
    end
    r = k;
    pw(i) = (f(r)-f(l))/(f2-f1)*length(x)/fs; % 计算脉冲宽度
end

% 绘制结果
figure;
subplot(2,1,1);
spectrogram(x,wlen,hop,nfft,fs,'yaxis');
title('原始信号');
subplot(2,1,2);
plot(t,pw);
xlabel('时间（s）');
ylabel('脉冲宽度（s）');
title('脉冲宽度随时间变化');

该代码首先生成一个线性调频信号，然后使用S-fracFT计算信号的时频分布，最后根据S-fracFT结果提取脉冲宽度并绘制结果。