0.3*randn(1,n));

yf = fftshift(fft(y));
ynf = fftshift(fft(yn));
xf = fdaxis(dt,n);

subplot(2,1,1)
plot(time,imag(y),time,imag(yn)); xlabel('time (ns)')
subplot(2,1,2)
plot(xf,real(yf),xf,real(ynf)); xlabel('frequency (MHz)'); xlim([-100 100])
</pre>

<p>
Perform lpsvd on the noisy data
</p>

<pre class="matlab">
yp = lpsvd(yn,time,'ss','aic');
ypf = fftshift(fft(yp));

subplot(2,1,1)
plot(time,imag(yn),time,imag(yp)); xlabel('time (ns)')
subplot(2,1,2)
plot(xf,real(ynf),xf,real(ypf)); xlabel('frequency (