anaaskt.m

matlab时频工具箱

%function anaaskt
%ANAASKT Unit test for the function ANAASK.

%	O. Lemoine - February 1996.

N=256;

% Output frequency law
f0=0.02;
[signal,am]=anaask(N,32,f0);
iflaw=instfreq(signal);
if any(abs(iflaw-f0*ones(N-2,1))>sqrt(eps))~=0,
  error('anaask test 1 failed');
end

f0=0.25;
[signal,am]=anaask(N,41,f0);
iflaw=instfreq(signal);
if any(abs(iflaw-f0*ones(N-2,1))>sqrt(eps))~=0,
  error('anaask test 2 failed');
end

f0=0.49;
[signal,am]=anaask(N,57,f0);
iflaw=instfreq(signal);
if any(abs(iflaw-f0*ones(N-2,1))>sqrt(eps))~=0,
  error('anaask test 3 failed');
end

% Output amplitude
f0=0.01;
[signal,am]=anaask(N,26,f0);
if any(am>1)~=0|any(am<0)~=0,
  error('anaask test 4 failed');
end

f0=0.17;
[signal,am]=anaask(N,37,f0);
if any(am>1)~=0|any(am<0)~=0,
  error('anaask test 5 failed');
end

f0=0.43;
[signal,am]=anaask(N,12,f0);
if any(am>1)~=0|any(am<0)~=0,
  error('anaask test 6 failed');
end

N=221; 

% Output frequency law
f0=0.02;
[signal,am]=anaask(N,32,f0);
iflaw=instfreq(signal);
if any(abs(iflaw-f0*ones(N-2,1))>sqrt(eps))~=0,
  error('anaask test 7 failed');
end

f0=0.25;
[signal,am]=anaask(N,41,f0);
iflaw=instfreq(signal);
if any(abs(iflaw-f0*ones(N-2,1))>sqrt(eps))~=0,
  error('anaask test 8 failed');
end

f0=0.49;
[signal,am]=anaask(N,57,f0);
iflaw=instfreq(signal);
if any(abs(iflaw-f0*ones(N-2,1))>sqrt(eps))~=0,
  error('anaask test 9 failed');
end

% Output amplitude
f0=0.01;
[signal,am]=anaask(N,26,f0);
if any(am>1)~=0|any(am<0)~=0,
  error('anaask test 10 failed');
end

f0=0.17;
[signal,am]=anaask(N,37,f0);
if any(am>1)~=0|any(am<0)~=0,
  error('anaask test 11 failed');
end

f0=0.43;
[signal,am]=anaask(N,12,f0);
if any(am>1)~=0|any(am<0)~=0,
  error('anaask test 12 failed');
end

N=221;