limnums.m

来自「A MATLAB tool for analysis of Transient 」· M 代码 · 共 47 行

function [mme, smme] = limnums(Pr, ttf, n_lev, resplim);

%  limnums    finds figures of time-dependent spectrum
%
%   [mme, smme] = limnums(Pr, ttf, n_lev, resplim);
%
%   in:   Pr        T-D Spectrum
%         ttf       time in milliseconds
%         n_lev     noise level dB
%         resplim   limit of start and end points
%
%  out:   mme       Matrix of figures 
%                   [ Maximum values dB
%                     time of the maxima
%                     start of response
%                     end of response
%                     maximum derivative
%                     minimum derivative ]
%
%         smme      Same figures as strings in rows
%                   separeted with setstr(124)

% (c) Pekka Kumpulainen 21.2.1995
%     PK 26.11.1996 maksimit ja muut korjattu

[m, n] = size(Pr);
mme = zeros(6,n);

% max dB
[mp, mi] = max(Pr);
limind = mp-n_lev<0; % maksimin oltava suurempi kuin kohinataso
mp(limind) = NaN*mp(limind);
smaxim = strrep(sprintf('%3.1f|',mp),'NaN',' ');
%smaxim = smaxim(1:(length(smaxim)-1));

% max ms
maxtime = ttf(mi)';
maxtime(limind) = NaN*maxtime(limind);
stmaxim = strrep(sprintf('%3.1f|',maxtime),'NaN',' ');
%stmaxim = stmaxim(1:(length(stmaxim)-1));

% start and end points
sttime = zeros(1,n); % 24.9.97
entime = zeros(1,n); % 24.9.97
ssttime = []; 
sentime = [];
limit = n_lev+resplim; % miss