后附有源码可直接运行，x,y,z坐标尚未添加

等高线图为上图横切面，横轴为延迟时间/s，纵轴为多普勒频移/Hz

clc;
close all;
clear all;
taup = 1;
b = 10;
up_down = 1;
eps = 0.000001;
i = 0;
mu = up_down * b / 2./ taup;
delt = 2.2*taup/250;
delf = 2*b /250;
for tau = -1.1*taup:.05:1.1*taup
   i = i + 1;
   j = 0;
   for fd = -b:.05:b
      j = j + 1;
      val1 = 1. - abs(tau) / taup;
      val2 = pi * taup * (1.0 - abs(tau) / taup);
      val3 = (fd + mu * tau);
      val = val2 * val3;
      x(j,i) = abs( val1 * (sin(val+eps)/(val+eps))).^2;
   end
end
figure(1)
mesh(x);%三维网格图形绘制
figure(2)
contour(x)%等高线图形绘制