momftfrt.m

来自「matlab时频工具箱」· M 代码 · 共 55 行

%function momftfrt
%MOMFTFRT Unit test for the function MOMFTFR.

%	O. Lemoine - March 1996.

N=128;

% For a perfect line
[sig,ifl]=fmlin(N);
tfr=tfrideal(ifl);
[tm,D2]=momftfr(tfr); 
tmth=(1:N);
if any(abs(tm-tmth')>sqrt(eps)),
  error('momftfr test 1 failed');
elseif any(abs(D2)>sqrt(eps)),
  error('momftfr test 2 failed');
end

% For a sinusoid
[sig,ifl]=fmsin(N,0,.5,N*2,1,0);
tfr=tfrideal(ifl);
[tm,D2]=momftfr(tfr'); 
tmth=round(ifl*2*(N-1))+1;
if any(abs(tm-tmth)>sqrt(1/N)),
  error('momftfr test 3 failed');
elseif any(abs(D2)>sqrt(eps)),
  error('momftfr test 4 failed');
end


N=117;

% For a perfect line
[sig,ifl]=fmlin(N);
tfr=tfrideal(ifl);
[tm,D2]=momftfr(tfr); 
tmth=(1:N);
if any(abs(tm-tmth')>sqrt(eps)),
  error('momftfr test 5 failed');
elseif any(abs(D2)>sqrt(eps)),
  error('momftfr test 6 failed');
end

% For a sinusoid
[sig,ifl]=fmsin(N,0,.5,N*2,1,0);
tfr=tfrideal(ifl);
[tm,D2]=momftfr(tfr'); 
tmth=round(ifl*2*(N-1))+1;
if any(abs(tm-tmth)>sqrt(1/N)),
  error('momftfr test 7 failed');
elseif any(abs(D2)>sqrt(eps)),
  error('momftfr test 8 failed');
end