% power delay angular profile Ex 6.8
p=[.1 .31 .1 .032];
d=0:.2:6;
ang=-90:5:90;
[x y]=meshgrid(d,ang);
Z1=(abs(x-0)+abs(y-(-80))<=.01)*p(1);
Z2=(abs(x-1)+abs(y-(-45))<=.01 )*p(2);
Z3=(abs(x-3)+abs(y-(40))<=.01)*p(3);
Z4=(abs(x-5)+abs(y-(60))<=.01)*p(4);
Z=Z1+Z2+Z3+Z4;
figure;
mesh(y,x,Z)
colormap([0 0 0])
set(gca,'ydir','reverse')
axis([-90 90 0 6 0 .4])
ylabel('delay (\mus)','rotation',-30)
xlabel('arrival angle','rotation',20)
Zlabel('PDAP')