rtestfun.m

Matlab时频分析工具箱,希望能对大家有所帮助啊

function sig2=rtestfun(L)
%RTESTFUN  Real 1-D test function.
%   Usage  sig=rtestfun(L);
%
%   RTESTFUN(L) return a synthetic test function of length L.
% 
%   The function is a superposition of chirp and an indicator function.

% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
% 
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
% GNU General Public License for more details.
% 
% You should have received a copy of the GNU General Public License
% along with this program.  If not, see <http://www.gnu.org/licenses/>.

tt=(0:floor(L/4)-1).';
f1=sin(2*pi*tt/ceil(L/12)); % sinusoid

f2=zeros(floor(L/8),1);
f2(round(L/16))=1; % delta-like

f3=triang(floor(L/8));

f4 = fftshift(pgauss(floor(L/4)));
f4 = f4/max(f4);

sig = 0.9*[f1; f2; f3; -f3; f3; f4];

sig2=zeros(L,1);
sig2(1:length(sig))=sig;

%ftest=zeros(L,1);

%sp=round(L/4);
%lchirp=round(L/2);
%ftest(sp+1:sp+lchirp)=sin(2*linspace(0,2*pi*sqrt(lchirp)/10,lchirp).^2)';

%s=round(L*7/16);
%l=round(L/16);
%ftest(s:s+l)=ftest(s:s+l)+ones(l+1,1);