工作频率1GHz、目标1（距离75km、速度100m/s）、目标2（距离50km、速度80m/s）发射信号时宽为1μs，带宽为1GHz，采样频率为2GHz，脉冲重复周期为10μs。1.线性调频波形产生、②回波信号产生、③脉冲压缩、④MTI/MTD、⑤目标探测、⑥参数估计，最终获得目标的距离、速度信息。MATLAB仿真

这是一个比较复杂的雷达信号处理问题，需要分几步来实现。以下是一种可能的实现方案，仅供参考。

1. 线性调频波形产生

首先需要产生一个带宽为1GHz、持续时间为1μs的线性调频信号。可以使用MATLAB中的`chirp`函数来实现：

fc = 1e9; % 工作频率
T = 1e-6; % 信号持续时间
B = 1e9; % 带宽
t = linspace(0, T, 2*B*T); % 时间向量
signal_tx = chirp(t, 0, T, B, 'linear', 90); % 产生线性调频信号

2. 回波信号产生

将信号发射到目标上，得到回波信号。假设有两个目标，分别位于距离75km和50km处，速度分别为100m/s和80m/s，可以使用以下代码产生回波信号：

lambda = physconst('LightSpeed')/fc; % 波长
R = [75e3, 50e3]; % 两个目标距离
v = [100, 80]; % 两个目标速度
tau = 2*R/c; % 两个目标回波时间延迟
signal_rx = zeros(1, length(signal_tx)); % 初始化回波信号
for i = 1:length(R)
    signal_tx_delayed = [zeros(1, floor(tau(i)*fs)), signal_tx(1:end-floor(tau(i)*fs))]; % 延迟发射信号
    signal_rx = signal_rx + exp(-1j*2*pi*fc*tau(i))*signal_tx_delayed.*exp(1j*2*pi*fc*2*v(i)/lambda*(t-tau(i))); % 加入多普勒频移
end

3. 脉冲压缩

接下来需要对回波信号进行脉冲压缩，将1μs宽的信号压缩到更短的时间内。可以使用线性调频信号的相关性质来实现，以下是一种可能的实现方案：

T_p = 10e-6; % 脉冲重复周期
signal_tx = chirp(linspace(0, T_p, 2*B*T_p), 0, T, B, 'linear', 90); % 重新生成线性调频信号
signal_tx_delayed = [zeros(1, floor(tau(1)*fs)), signal_tx(1:end-floor(tau(1)*fs))]; % 延迟发射信号
signal_tx_comp = conj(fliplr(signal_tx)); % 产生压缩脉冲
signal_rx_comp = conv(signal_tx_delayed, signal_tx_comp, 'same'); % 对回波信号进行脉冲压缩

4. MTI/MTD

对于目标探测，可以使用MTI（运动目标指示）或MTD（运动目标检测）算法来实现。以下是一种可能的MTI算法实现方案：

M = 10; % MTI中的通道数
window = ones(1, M)/M; % 移动平均窗口
signal_mti = zeros(size(signal_rx_comp)); % 初始化MTI信号
for i = M+1:length(signal_rx_comp)
    signal_mti(i) = abs(sum(signal_rx_comp(i-M:i).*window)); % 对每个时刻的信号进行移动平均
end

5. 目标探测

对于目标探测，可以使用一些常见的方法，例如判决门限法、恒虚警率法等。以下是一种基于判决门限的实现方案：

SNR_threshold = 10; % 判决门限（信噪比阈值）
signal_detected = signal_mti > SNR_threshold*std(signal_mti); % 目标探测

6. 参数估计

对于探测到的每个目标，可以通过一些信号处理算法来估计其距离和速度。以下是一种可能的实现方案：

[N, p] = findpeaks(signal_mti, 'MinPeakHeight', SNR_threshold*std(signal_mti)); % 找到所有峰值
for i = 1:length(N)
    R_est = c*(p(i)/fs + tau(1))/2; % 估计距离
    v_est = lambda/4*(angle(signal_rx(p(i)+1))-angle(signal_rx(p(i)-1)))/(2*pi/fs); % 估计速度
    disp(['目标 ', num2str(i), ' 估计距离为 ', num2str(R_est/1000), ' km，速度为 ', num2str(v_est), ' m/s']);
end

以上是一个可能的实现方案，但由于雷达信号处理涉及到许多细节和参数，实际的实现可能需要根据具体情况进行调整和优化。