方法一：原始方法

T = 10e-6;               %脈寬
B = 4e6;                %信號(hào)帶寬
K = B/T;                 %調(diào)頻斜率
Fs = 80e6;                %采樣頻率
Ts = 1/Fs;               %采樣間隔
N = T/Ts;
t = linspace(-T/2,T/2,N);
lfm = exp(j*pi*K*t.^2);  %生成線性調(diào)頻信號(hào)
%以上可以分解為exp(j*pi*(K*t)*t),即K*t=f，f隨著時(shí)間線性變化；

方法二：進(jìn)階方法，利用matlab封裝好的對(duì)象

clc;
clear;
fs = 80e6;
waveform = phased.LinearFMWaveform('SampleRate',fs,...
    'SweepBandwidth',4e6,...
    'PulseWidth',1e-5,'PRF',4e3);
% 'SweepBandwidth'為掃頻帶寬 
% 'PulseWidth'為脈沖寬度 
% 'PRF'為脈沖重復(fù)周期
waveform.CoefficientsOutputPort=true;
% plot(waveform,'PlotType','real');
% plot(waveform,'PlotType','imag');
plot(waveform,'PlotType','complex');

下圖為上圖線性調(diào)頻信號(hào)對(duì)應(yīng)額的匹配濾波系數(shù)，該系數(shù)與原始波形的關(guān)系為共軛（虛部反相）鏡像（時(shí)間上翻轉(zhuǎn)）。匹配系數(shù)在相關(guān)接收時(shí)需要用到，如線性調(diào)頻脈沖雷達(dá)的脈沖壓縮（即匹配濾波）系數(shù)；

方法三：高階方法，利用matlab工具箱（雷達(dá)波形分析工具箱）

工具箱產(chǎn)生波形后可導(dǎo)出到matlab工作區(qū)間做仿真使用。