写出ＭＡＴＬＡＢ实现MTI动目标显示的函数

function [output] = MTI_Doppler_Filter(input,PRF,fc,PRI,tau)

% input: 输入的信号矩阵，每列为一帧数据
% PRF: 脉冲重复频率，单位为Hz
% fc: 雷达中心频率，单位为Hz
% PRI: 脉冲重复间隔，单位为秒
% tau: MTI滤波延迟时间，单位为秒

% output: 经过MTI滤波后的信号矩阵

% 计算滤波器系数
N = size(input,1); % 每帧数据的点数
fs = 1/PRI; % 采样频率，单位为Hz
f_step = fs/N; % 频率分辨率，单位为Hz
f = zeros(N,1);
for i = 1:N
f(i) = (i-N/2-1)*f_step; % 构造频率轴
end
H = zeros(N,1);
for i = 1:N
if abs(f(i)) < PRF/2
H(i) = 1; % 带通滤波器
end
end
H = H.*exp(-1j*2*pi*f*tau); % 带通滤波器乘以相位因子
h = ifft(H); % 滤波器时域系数

% MTI滤波
output = zeros(size(input));
for i = 1:size(input,2)
input_frame = input(:,i); % 取出一帧数据
output_frame = conv(input_frame,h,'same'); % 卷积滤波
output(:,i) = output_frame;
end

% 显示结果
figure;
imagesc(abs(output));
xlabel('Range (m)');
ylabel('Time (s)');
title('MTI Doppler Filtered Signal');
colorbar;
end