function [Xi,fsi]=FecgInterp(X,fs,interpFact,cName,graph)
% ----------------------------------------------------------------------------------------------
%   ECG  interpolation
%
% [Xi,fsi]=FecgInterp(X,fs,interpFact,cName,graph)
%
% X          : input signal matrix (one signal per column)
% fs         : sampling frequency
% interpFact : interpolation factor
% cName      : record name
% graph      : flag enabling figure drawing
%
% Xi         : interpolated signal matrix
% fsi        : sampling frequency of interpolated signals
%
% Author: Maurizio Varanini, Clinical Physiology Institute, CNR, Pisa, Italy
% For any comment or bug report, please send e-mail to: maurizio.varanini@ifc.cnr.it
% ----------------------------------------------------------------------------------------------

if(nargin<3), cName=''; end
if(nargin<4), graph=0; end

[progpath, progname] = fileparts(which(mfilename));
fprintf('\n --------------------------------------------------------- \n');
fprintf('Program: %s,  record name: %s\n', progname, cName);

%-------------------------------------------------------------
% recording duration
[ndt, ns]=size(X);

% interpolation by Fourier Transform
Xi=interpft(X,round(interpFact*ndt));
fsi=interpFact*fs;

if(graph)
    figure, set(gcf,'Color','white');
    vtime= [1:ndt]/fs;
    vtimei= [1:size(Xi,1)]/(fsi);
    for is=1:ns,
        subplot(ns,1,is), hold on, plot(vtime,X(:,is),'r'), plot(vtimei,Xi(:,is),'b');
        wgmimaV=mimaxscG(Xi(:,is),0,0,.1);
        ylim(wgmimaV);
        if(is~=ns), set(gca,'XTickLabel',''); end
        if(is==1), title([cName,': original & interpolated ECG']); end
        shg
    end
end
end %== function ================================================================
%