fsmooth.m

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

function y=fsmooth(x)

%FSMOOTH Taajuustasoon lis刾isteiden aproksimointi ja silitys t-f-spektrille
%      y=fsmooth(x)     

% PK 4.8.-93  PK 22.9.1993


x=x.';
[m,n]=size(x);
nz=10;

z=zeros((m+nz)*n,1);
z(:)=[x;zeros(nz,n)];
zz=zeros(2*size(z,1),1);
zz(:)=[z.';zeros(size(z.'))];

[b,a]=butter(2,.5);

y=zeros(2*(m+nz),n);
y(:)=2*filtfilt(b,a,zz);

y=y(1:(2*m-1),:).';
y(y<0)=0*y(y<0);

% VANHA JA HIDAS
%  [m,n]=size(x);y=zeros(m,4*n-3); 
%  zer0=zeros(10,1);
%  zer1=zeros(3,n);
%  zer2=zeros(20,1);
%  x1=zeros(4*n,1);
%  [b,a]=butter(3,.25);
%  for ii=1:m;
%      x1(:)=[x(ii,:);zer1];
%      z=filtfilt(b,a,[zer0;x1;zer2]);
%      y(ii,:)=4*z(11:4*n-3+10)';
%  end