Cubic Phase Impulse Response

From Optical Imaging and Spectroscopy (Wiley - OSA 2009)
By David J. Brady
www.opticalimaging.org
Figure 10.7 and 10.8

ap=10^-5;
hr[thetax_,thetaz_]:=AiryAi[(-2 Pi thetax +Pi^2 thetaz^2/(3 ap))/Power[3 ap,1/3]];

GraphicsColumn[{Plot[hr[tx,0]^2,{tx,-.05,.5},PlotRange->All,PlotPoints->300,Axes->False, Frame->True, FrameLabel->{"(Subscript[θ, x])","(Subscript[h, r])"},PlotStyle->{Black},PlotLabel->{(a)}],
DensityPlot[hr[tx,0]^2 hr[ty,0]^2,{tx,-.05,.1},{ty,-.05,.1},PlotRange->All,PlotPoints->100,ColorFunction->GrayLevel,PlotLabel->{(b)}]}]

DensityPlot[hr[tx,0]^2 hr[ty,0]^2,{tx,-.05,.5},{ty,-.05,.5},PlotRange->All,PlotPoints->100]

ap=10^-5;
hrl[thetax_]:=NIntegrate[AiryAi[-2 Pi thetax /Power[3 ap l^2,1/3]]^2,{l,.9,1.1}];

abd=Plot[hr[tx,0]^2,{tx,-.005,.1},PlotRange->All,PlotPoints->300,Axes->False, Frame->True, FrameLabel->{"(Subscript[θ, x])","(Subscript[h, r])"},PlotStyle->{Black},PlotLabel->{(a)}]

GraphicsRow[{abd,dog}]

dog=Plot[hrl[tx],{tx,-.005,.1},PlotRange->All,PlotPoints->300,Axes->False, Frame->True, FrameLabel->{"(Subscript[θ, x])","(Subscript[h, r])"},PlotStyle->{Black}]