%cross-terms;
t=[1:128];
x1=exp(j*0.2*t);
x2=exp(j*0.8*t);
A1=ambigu(x1,128);
A2=ambigu(x2,128);
x=x1+x2;
A=ambigu(x,128);
subplot(221);
imagesc(abs(A1));
subplot(222);
imagesc(abs(A2));
subplot(223);
imagesc(abs(A));
An=A.*[zeros(128,63) ones(128,1) zeros(128,64)];
tfn=dochange(An,'amb','tfd');
subplot(224);
imagesc(tfn);
pause;
%example 2;
x1=logon(32,0.4,128);
x2=logon(64,0.8,128);
x=x1+x2;
A1=ambigu(x1,128);
x2=ambigu(x2,128);
A=ambigu(x,128);
figure
subplot(221);
imagesc(abs(A1));
subplot(222);
imagesc(abs(A));
w=dochange(A,'amb','tfd');
subplot(223);
imagesc(w);
An=A.*[zeros(128,50) ones(128,30) zeros(128,48)];
tfn=dochange(An,'amb','tfd');
subplot(224);
imagesc(tfn);
pause;
%three signals;
figure;
x=freqhops(32,0.2,64,0.8,32,-0.6);
A=ambigu(x,128);
W=wigner(x,128);
subplot(221);
imagesc(W);
subplot(222);
imagesc(abs(A));
subplot(223);
An=A.*[zeros(128,40) ones(128,40) zeros(128,48)];
imagesc(abs(An));
subplot(224);
tf=dochange(An,'amb','tfd');
imagesc(tf);