tfrideat.m

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

%function tfrideat
%TFRIDEAT Unit test for the function TFRIDEAL.

%	O. Lemoine - March 1996.

N=128;

[sig1,ifl1]=fmlin(N); [sig2,ifl2]=fmsin(N);
[sig3,ifl3]=fmhyp(N,'f',[1 0.5],[N/2 0.1]); 
tfr=tfrideal([ifl1;ifl2;ifl3]); 
if any(any(tfr~=1 & tfr~=0)),
 error('tfrideal test 1 failed');
elseif any(sum(tfr)~=1),
 error('tfrideal test 2 failed');
end

iflaws1=[ifl1;ifl2;ifl3]';
for t=1:3*N,
  [Max iflaws2(t)]=max(tfr(:,t));
end
iflaws2=(iflaws2-1)*0.5/(3*N-1);

if any(abs(iflaws1-iflaws2)>1/(3*N)), 
 error('tfrideal test 3 failed');
end


clear; N=111;

[sig1,ifl1]=fmlin(N); [sig2,ifl2]=fmsin(N);
[sig3,ifl3]=fmhyp(N,'f',[1 0.5],[round(N/2) 0.1]); 
tfr=tfrideal([ifl1;ifl2;ifl3]); 
if any(any(tfr~=1 & tfr~=0)),
 error('tfrideal test 4 failed');
elseif any(sum(tfr)~=1),
 error('tfrideal test 5 failed');
end

iflaws1=[ifl1;ifl2;ifl3]';
for t=1:3*N,
  [Max iflaws2(t)]=max(tfr(:,t));
end
iflaws2=(iflaws2-1)*0.5/(3*N-1);

if any(abs(iflaws1-iflaws2)>1/(3*N)), 
 error('tfrideal test 6 failed');
end