24小时热门版块排行榜

返回列表

q3264585

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 4.5
帖子: 9
在线: 3.3小时
虫号: 4734168
注册: 2016-05-29
性别: GG
专业: 人工智能与知识工程

[求助] 请教matlab中时频分析工具已有1人参与

matlab中有官方提供时频分析工具tftb，请问有大神会使用其中的ecgpuwave这个函数么？

回复此楼

» 收录本帖的淘帖专辑推荐

程序

» 猜你喜欢

1楼 2016-08-03 17:35:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

信息EPI: 1
应助: 591 (博士)
贵宾: 0.03
金币: 18173.2
红花: 97
帖子: 8891
在线: 1487.4小时
虫号: 2400059
注册: 2013-04-04
专业: 计算机应用技术

【答案】应助回帖

感谢参与，应助指数 +1

%% Write an ECG header in MIT format
% This function receives a "header" struct with ECG signal propeties and
% writes [header.name '.hea'] file into "header_path".
%
% Example
%
%       MIT_filename  = 'your_filename';
%       MIT_path = ['.' filesep];
%       ECG_header.recname = MIT_filename;
%       ECG_header.nsig = 8;
%       ECG_header.nsamp = 9978;
%       ECG_header.freq = 1000; %Hz
%       ECG_header.desc = char( {'I','II','V1','V2','V3','V4','V5','V6'} );
%       ECG_header.adczero = zeros( ECG_header.nsig,1 );
%       ECG_header.gain = ones( ECG_header.nsig,1 );
%       ECG_header.units = repmat( 'uV',ECG_header.nsig,1 );
%       ECG_header.btime = '00:00:00';
%       ECG_header.bdate = '01/01/2000';
%
%       writeheader( MIT_path, ECG_header );
%
% See also read_ECG, ECGwrapper
%
% Author: Salvador Olmos, Juan Pablo Mart穩nez and/or Juan Bolea
% adapted to ECGkit by: Mariano Llamedo Soria
% <matlab:web( 'mailto:llamedom@electron.frba.utn.edu.ar','-browser' ) ( email )>
% Version: 0.1 beta
% Birthdate: 5/01/2014
% Last update: 19/11/2014
%
function writeheader( header_path, header )
%
% WRITHEAD function writes header file for signal data struct in directory header_path
% Input parameters.
%    header_path: directory of work
%    header: struct of header to write

if( header_path( end ) ~= filesep  )
header_path = [header_path filesep];
end

% Opening ASCII header file
if ~isempty( header.recname )
fid=fopen( [header_path  header.recname '.hea'], 'w' );
header.fname = repmat( [ header.recname '.dat'], header.nsig,1 );
else
fid = fopen( [header_path  header.fname( 1,: ) '.hea'], 'w' );
end

if ( fid<0 )
error( ['Can not write ' [header_path  header.recname '.hea'] '\n'] );
end

try

% writing first line of record_name, # signals and so on
fprintf( fid,'%s %d %d %d\n',header.recname,header.nsig );
if isfield( header,'freq' )
      fprintf( fid,'%d ',header.freq );
end
if isfield( header,'nsamp' )
      fprintf( fid,'%d ',header.nsamp );
end
if isfield( header,'btime' )
      fprintf( fid,'%s ',header.btime );
end
if isfield( header,'bdate' );
      fprintf( fid,'%s ',header.bdate );
end

fmt = 16;
% writing the rest of lines
for ii=1:header.nsig

      if isfield( header,'spf' ) || isfield( header,'skew' ) || isfield( header,'offset' )
         fprintf( fid,'\n%s %d',header.fname( ii,: ),fmt );
         if isfield( header,'spf' ) && ~isnan( header.spf( ii ) )
            fprintf( fid,'x%d',header.spf( ii ) );
         end
         if isfield( header,'skew' ) && ~isnan( header.skew( ii ) )
            fprintf( fid,':%d',header.skew( ii ) );
         end
         if isfield( header,'offset' ) && ~isnan( header.offset( ii ) )
            fprintf( fid,'+%d',header.offset( ii ) );
         end
         fprintf( fid,' ' );
      else
         fprintf( fid,'\n%s %d ',header.fname( ii,: ),fmt );
      end

      if isfield( header,'gain' )
         fprintf( fid,'%f',header.gain( ii )  );
      else
         fprintf( fid,'%f',200  );
      end

      if( isfield( header,'baseline' ) && ~isnan( header.baseline( ii ) )  )
         fprintf( fid,'( %d )',header.baseline( ii ) );
      end

      if ( isfield( header,'units' )  )
         fprintf( fid,'/%s ',header.units( ii,: ) );
      else
         fprintf( fid,' '  );
      end

      if isfield( header,'adcres' )
         fprintf( fid,'%d ',header.adcres( ii ) );
      else
         fprintf( fid,'%d ', 12 );
      end

      if isfield( header,'adczero' )
         fprintf( fid,'%d ',header.adczero( ii ) );
      else
         fprintf( fid,'%d ', 0 );
      end

      if isfield( header,'initval' )
         fprintf( fid,'%d ',header.adczero( ii ) );
      else
         fprintf( fid,'%d ', 0 );
      end

      if isfield( header,'cksum' )
         fprintf( fid,'%d ',header.cksum( ii ) );
      else
         fprintf( fid,'%d ', 0 );
      end

      if isfield( header,'bsize' )
         fprintf( fid,'%d ',header.bsize( ii ) );
      else
         fprintf( fid,'%d ', 0 );
      end

      if isfield( header,'desc' )
         fprintf( fid,'%s',header.desc( ii,: ) );
      end

end

fprintf( fid,'\n' );

fclose( fid );

catch MEE
fclose( fid );
rethrow( MEE );
end

參考:

https://groups.google.com/forum/ ... U0xsMY/CMXs2ZOvNoEJ