q713.asv

Cycle slip probability of the phase unrapping algorithm

clear;close all
n=0:255;N=4096;
x=cos(0.04*pi*n)+cos(0.08*pi*n)+cos(0.4*pi*n);
w=randn(size(x));
x=x+0.3*w;
b=[1,2,1];
B=0.0003738*conv(conv(b,b),b);
a1=[1,-1.2686,0.7051];
a2=[1,-1.0106,0.3583];
a3=[1,-0.9044,0.2155];
A=conv(conv(a1,a2),a3);
y=filter(B,A,x);
X=fft(x,N);
Y=fft(y,N);
[H,f]=freqz(B,A,N,'whole');
Ym=H'.*X;
k=0:N-1;f=2*k/N;
stem(x,'.');
subplot(3,2,1)
stem(y,'.')
s