解多普勒模糊matlab代码

以下是一个用MATLAB实现多普勒模糊处理的示例代码：

% 读入图像
img = imread('image.bmp');
% 显示原图
figure,imshow(img);
% 计算频率域运动补偿函数
[M, N] = size(img);
dx = 1.0;
fx = linspace(-1, 1, N) / dx;
fy = linspace(-1, 1, M) / dx;
[FX, FY] = meshgrid(fx, fy);
velocity = [10, 20]; % 目标速度
H = exp(-2j * pi * (FX * velocity(1) + FY * velocity(2)));
% 频率域滤波
img_fft = fft2(img);
img_fft_mc = img_fft .* H;
% 逆傅里叶变换
img_mc = ifft2(img_fft_mc);
% 显示运动补偿后的图像
figure,imshow(abs(img_mc),[]);

以上代码实现了频域运动补偿法处理多普勒模糊。需要注意的是，由于MATLAB中的FFT函数默认将原点放在左上角，而不是在中心位置，因此在计算频率域运动补偿函数时需要将频率坐标系做相应的调整。