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

以下是一个利用杜芬振子提取线性调频信号脉冲宽度的MATLAB代码示例：

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

% 利用杜芬振子提取脉冲宽度
fc = 150; % 杜芬振子中心频率
bw = 50; % 杜芬振子带宽
[b,a] = butter(4,2*fc/Fs); % 设计4阶低通滤波器
y = filter(b,a,x); % 信号经过低通滤波器
z = hilbert(y); % 信号解析成复信号
w = abs(z); % 取模得到幅值
bw_index = find(w > max(w)/2); % 找到幅值大于一半的部分
pw = (t(bw_index(end)) - t(bw_index(1))) * 1000; % 计算脉冲宽度，单位：ms

% 绘制信号和杜芬振子中心频率、带宽
figure;
subplot(2,1,1);
plot(t,x);
xlabel('时间（秒）');
ylabel('幅值');
title('线性调频信号');
subplot(2,1,2);
plot(t,w);
hold on;
line([t(1),t(end)],[max(w)/2,max(w)/2],'Color','r');
line([t(bw_index(1)),t(bw_index(end))],[max(w)/2,max(w)/2],'Color','g');
xlabel('时间（秒）');
ylabel('幅值');
title(['杜芬振子中心频率：',num2str(fc),'Hz，带宽：',num2str(bw),'Hz']);

这段代码首先生成一个1秒钟的线性调频信号，起始频率为100Hz，终止频率为200Hz。然后设计一个中心频率为150Hz，带宽为50Hz的杜芬振子，将信号经过低通滤波器后解析成复信号，取模后得到幅值，找到幅值大于一半的部分，计算脉冲宽度并输出。最后绘制信号和杜芬振子的幅值随时间变化的曲线，以及杜芬振子的中心频率和带宽。