目标距离为100km，场景中心为99km，目标速度为10m/s，脉冲数为64，阵元数为16，目标角度为15°，雷达载频为5GHz，带宽为10MHz，采样率为20MHz，PRF为100Hz，利用matlab雷达测速、测距、测角

下面是一个MATLAB示例代码，用于雷达测速、测距和测角：

% 雷达参数设置
speedOfLight = 3e8; % 光速 (m/s)
carrierFrequency = 5e9; % 载波频率 (Hz)
bandwidth = 10e6; % 带宽 (Hz)
pulseWidth = 1/bandwidth; % 脉冲宽度 (s)
chirpSlope = bandwidth/pulseWidth; % 扫频斜率 (Hz/s)
samplingRate = 20e6; % 采样率 (Hz)
prf = 100; % PRF (Hz)
numPulses = 64; % 脉冲数
numElements = 16; % 阵元数

% 目标参数设置
targetDistance = 100000; % 目标距离 (m)
sceneCenter = 99000; % 场景中心距离 (m)
targetSpeed = 10; % 目标速度 (m/s)
targetAngle = 15; % 目标角度 (度)

% 计算回波信号
timeDelay = (2 * (targetDistance - sceneCenter)) / speedOfLight; % 信号传播时间延迟 (s)
frequencyShift = 2 * targetSpeed * carrierFrequency / speedOfLight; % 频率偏移 (Hz)
t = 0:1/samplingRate:pulseWidth-1/samplingRate; % 时间轴
echoSignal = exp(1i * (2*pi*(carrierFrequency*t + chirpSlope*t.^2/2) + frequencyShift * t)); % 回波信号

% 生成阵列波束
theta = deg2rad(targetAngle); % 目标角度 (弧度)
arrayResponse = exp(-1i * 2 * pi * (0:numElements-1)' * sin(theta)); % 阵列波束

% 重复脉冲和波束
pulseRepetition = samplingRate/prf; % 脉冲重复周期 (样本数)
pulseTrain = repmat(echoSignal, [1, numPulses]); % 重复脉冲
arrayBeam = repmat(arrayResponse, [1, numPulses]); % 重复波束

% 接收信号
receivedSignal = pulseTrain .* arrayBeam; % 接收信号

% 快时间傅里叶变换
rangeProfile = fft(receivedSignal, [], 1); % 范围剖面
rangeProfile = rangeProfile(1:length(rangeProfile)/2, :); % 保留一半的范围剖面

% 显示范围剖面
figure;
imagesc(1:numPulses, linspace(0, targetDistance - sceneCenter, size(rangeProfile, 1)), abs(rangeProfile));
xlabel('脉冲序号');
ylabel('距离 (m)');
title('范围剖面');

% 提取速度信息
dopplerShift = fftshift(fft(rangeProfile, [], 2), 2); % 多普勒频移
dopplerShift = dopplerShift(:, numPulses/2 + 1:end); % 保留正频率部分
dopplerAxis = linspace(-prf/2, prf/2, size(dopplerShift, 2)); % 多普勒频率轴

% 显示速度频谱
figure;
imagesc(dopplerAxis, linspace(0, targetDistance - sceneCenter, size(rangeProfile, 1)), abs(dopplerShift));
xlabel('速度 (m/s)');
ylabel('距离 (m)');
title('速度频谱');

% 计算目标距离
maxRangeIndex = round((targetDistance - sceneCenter) / speedOfLight * samplingRate); % 目标距离索引
peakRangeProfile = rangeProfile(maxRangeIndex, :); % 目标距离处的范围剖面
targetRange = (sceneCenter + peakRangeProfile * speedOfLight / (2 * chirpSlope)); % 目标距离 (m)
disp(['目标距离为：', num2str(targetRange), ' 米']);

% 计算目标速度
[maxValue, maxIndex] = max(abs(dopplerShift(maxRangeIndex, :))); % 多普勒频移处的最大值和索引
targetVelocity = dopplerAxis(maxIndex); % 目标速度 (m/s)
disp(['目标速度为：', num2str(targetVelocity), ' m/s']);

% 计算目标角度
beamPattern = abs(fftshift(fft(arrayResponse))); % 阵列波束图案
angleAxis = linspace(-90, 90, numElements); % 角度轴
[maxValue, maxIndex] = max(abs(beamPattern)); % 阵列波束图案的最大值和索引
targetAngle = angleAxis(maxIndex); % 目标角度 (度)
disp(['目标角度为：', num2str(targetAngle), ' 度']);

这段代码首先设置了雷达的参数，包括光速、载波频率、带宽、脉冲宽度、扫频斜率、采样率和PRF等。然后设置了目标车辆的距离、场景中心距离、速度和角度。通过计算信号传播时间延迟和频率偏移，生成了回波信号。使用阵列波束对回波信号进行波束形成。然后进行脉冲重复和波束重复，得到接收信号。接下来进行快时间傅里叶变换，得到范围剖面。通过多普勒频移分析速度信息，并显示在速度频谱图中。计算目标距离、速度和角度，并在命令窗口中显示出来。

请注意，这只是一个简单的示例代码，实际的雷达测速系统可能会涉及更多的参数和算法。此外，对于测角部分的计算，需要根据具体的阵列形状和波束形成算法进行调整。