clc
clear all
close all

w=0.5e-3;
x=(-512:511)/1024*w;
D=D.*(D>0)/0.5;
D=fftshift(D);
plot(x,D)