xd actd.ilsb1.x_ax Excursion
x1 (1/T0)*timeint((NT-1)*T0,NT*T0,xd*exp(-i*2*pi*f0*t)*2,'nointerp') "Fundamental, displacement"
p1 (1/T0)*timeint((NT-1)*T0,NT*T0,p*exp(-i*2*pi*f0*t)*2,'nointerp') "Fundamental, Pressure"
p2 (1/T0)*timeint((NT-1)*T0,NT*T0,p*exp(-i*2*pi*2*f0*t)*2,'nointerp') "Second harmonic, Pressure"
p3 (1/T0)*timeint((NT-1)*T0,NT*T0,p*exp(-i*2*pi*3*f0*t)*2,'nointerp') "Third harmonic, Pressure"
p1dB 20*log10((abs(p1)/sqrt(2))/actd.pref_SPL)[dB] "Fundamental, dB"
p2dB 20*log10((abs(p2)/sqrt(2))/actd.pref_SPL)[dB] "Second harmonic, dB"
p3dB 20*log10((abs(p3)/sqrt(2))/actd.pref_SPL)[dB] "Third harmonic, dB"