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

以下是短时分数阶傅里叶变换提取线性调频信号的脉冲宽度的MATLAB代码：

clear all;
close all;

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

% 对信号进行STFT
wlen = 128; % 窗长
nfft = 256; % FFT点数
hop = wlen/4; % 帧移
win = hamming(wlen); % 汉明窗
[S,F,T] = spectrogram(x,win,hop,nfft,fs); % 短时傅里叶变换

% 计算分数阶阶数
alpha = 0.8; % 阶数
N = nfft; % FFT点数
M = wlen; % 窗长
K = (N/M)*(M-1); % 帧数
beta = 2*pi*alpha/K;

% 计算分数阶傅里叶变换
Y = zeros(N,K);
for k = 1:K
    for n = 1:N
        for m = 1:M
            Y(n,k) = Y(n,k) + win(m)*x((k-1)*hop+m)*exp(-1i*beta*(k-1)*(n-1))*((k-1)*hop+m-1)^alpha;
        end
    end
end

% 提取脉冲宽度
pw = zeros(1,K);
for k = 1:K
    [~,ind] = max(abs(Y(:,k))); % 找到最大值的索引
    pw(k) = (ind-1)*fs/N; % 计算脉冲宽度
end

% 画图
figure;
subplot(211);
plot(t,x);
title('线性调频信号');
xlabel('时间(s)');
ylabel('幅度');
subplot(212);
mesh(T,F,abs(S));
title('短时傅里叶变换');
xlabel('时间(s)');
ylabel('频率(Hz)');
zlabel('幅度');
view(0,90);
figure;
plot(pw);
title('脉冲宽度');
xlabel('帧数');
ylabel('宽度(Hz)');

该代码首先生成一个线性调频信号，然后对该信号进行短时傅里叶变换，并计算分数阶傅里叶变换。接着，提取每个帧的脉冲宽度，并绘制出来。