24小时热门版块排行榜    

查看: 1631  |  回复: 2

q3264585

新虫 (初入文坛)

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

matlab中有官方提供时频分析工具tftb,请问有大神会使用其中的ecgpuwave这个函数么?
回复此楼

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

程序

» 猜你喜欢

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

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +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
2楼2016-08-04 06:58:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

3楼2016-08-04 07:08:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 学员bG1DZw 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 今天基金会出结果吗?20260819 +13 kkkl_v 2026-08-19 13/650 2026-08-19 22:37 by 宝贝虫子
[基金申请] 时间戳变了,能看出什么问题? +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息,中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见