在上面使用线性调频信号

好的，如果使用线性调频信号，我们需要修改生成脉冲信号的代码。具体来说，我们需要先定义线性调频信号的参数，如起始频率、终止频率、调频时间等：

% 线性调频信号参数
f0 = 9e9; % 起始频率
f1 = 11e9; % 终止频率
t_chirp = 10e-6; % 调频时间

然后，我们生成线性调频信号并进行相参积累：

% 生成线性调频信号
txsig = zeros(8, ceil(prf/pulselen)*pulselen*fs);
for i = 1:3
    for j = 1:8
        t = (0:1/fs:pulselen-1/fs)';
        k = (f1-f0)/t_chirp;
        txsig(j,:) = txsig(j,:) + exp(1j*2*pi*(f0*t + 0.5*k*t.^2 + 2*d(i)/lambda)) .* rectpuls(t-pulselen/2, pulselen);
    end
end

% 相参积累
n = 8;
sumsig = sum(txsig, 1);
for i = 2:n
    sumsig(i:end) = sumsig(i:end) .* conj(txsig(1:end-i+1));
end

接下来，我们进行 MTI 检测和 CFAR 检测，与之前的方法相同：

% MTI 检测
tau = 1/prf;
mti = abs(filter([1 -1], 1, sumsig));
mti(1) = 0; % 第一个点要设为0
mti = mti .* exp(-1j*2*pi*fc*tau*(0:length(mti)-1));

% CFAR 检测
winlen = 100; % 窗口长度
guardlen = 10; % 前后保护间隔
threshold = 10; % 阈值

cfar = zeros(size(mti));
for i = 1+guardlen : length(mti)-guardlen-winlen+1
    noise = [mti(i-guardlen:i-1), mti(i+winlen:i+guardlen)];
    threshold_value = threshold * mean(abs(noise));
    if abs(mti(i)) > threshold_value
        cfar(i) = mti(i);
    end
end

这样，我们就完成了使用线性调频信号进行相参积累、MTI 和 CFAR 检测的代码实现。