杜芬振子提取线性调频信号的脉冲宽度matlab

杜芬振子可以用于提取线性调频（LFM）信号中的脉冲宽度。下面是用MATLAB实现的示例代码：

% 生成一个LFM信号
fs = 1000; % 采样率
t = 0:1/fs:1; % 时间范围
f0 = 10; % 起始频率
f1 = 100; % 终止频率
T = 0.5; % 脉冲宽度
s = chirp(t,f0,1,f1,'linear',[],T);

% 杜芬振子提取脉冲宽度
N = length(s); % 信号长度
k = (0:N-1) - (N-1)/2; % 构造k序列
w = exp(-1i*2*pi*k/N); % 杜芬振子
y = fft(s.*w); % 信号乘上杜芬振子后进行FFT
p = abs(y); % 取振幅谱
[maxval,maxidx] = max(p); % 找到最大值及其下标
pw = abs(k(maxidx))/N*fs; % 计算脉冲宽度

% 绘图
figure;
plot(t,s);
title('LFM信号');
xlabel('时间（秒）');
ylabel('幅度');
figure;
plot(k,abs(y));
title('LFM信号频谱');
xlabel('频率');
ylabel('幅度');
figure;
plot(p);
title('LFM信号杜芬振子谱');
xlabel('杜芬振子频率');
ylabel('幅度');
fprintf('脉冲宽度：%.2f 秒\n',pw);

上述代码中，首先生成一个LFM信号，然后构造杜芬振子并将其应用于信号上，进行FFT得到频谱，最后在频谱中找到最大值并计算脉冲宽度。