注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

highspeedlogic

FPGA/MATLAB/Simulink

 
 
 

日志

 
 

matlab代做|matlab专业代做|ECG信号仿真源码程序  

2017-04-24 04:29:30|  分类: 默认分类 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
matlab代做|matlab专业代做|ECG信号仿真源码程序
QQ: 1224848052
网址1:www.hslogic.com
网址2:www.mat7lab.com
ECG其全称为electrocardiogram,即心电图。所谓ECG(心电图)就是指心脏在每个心动周期中,由起搏点、心房、心室相继兴奋,伴随着心电图生物电的变化,通过心电描记器从体表引出多种形式的电位变化的图形(简称ECG)。心电图是心脏兴奋的发生、传播及恢复过程的客观指标。

ECG的检查意义在于:用于对各种心律失常、心室心房肥大、心肌梗死、心率异常、心肌缺血、电解质紊乱(对血钾不正常变化有快速直视的临床参考意义)、心衰等病症检查,可用于床边24小时监视病人心脏功能。

  心电图是反映心脏兴奋的电活动过程,它对心脏基本功能及其病理研究方面,具有重要的参考价值。心电图可以分析与鉴别各种心律失常;也可以反映心肌受损的程度和发展过程和心房、心室的功能结构情况。在指导心脏手术进行及指示必要的药物处理上有参考价值。然而,心电图并非检查心脏功能状态必不可少的指标。

%-----------------------------------------------------------------------------
%    experiment 1:ECG filter
%    ecgfilter.m 
%    written by Ren Yan,no. 20070130406
%    last updated: 01/01/2008
%-----------------------------------------------------------------------------
clear all;
load ecg.mat; 
wp = 0.2*pi; ws = 0.3*pi;
rp = 1; rs = 15;
% design
[n,wn] = buttord(wp/pi,ws/pi,rp,rs);
[b,a] = butter(n,wn);
[h1,w1] = freqz(b,a,512,1);
H1 = 20*log10(abs(h1));
figure(1);
plot(w1*2,H1);
xlabel('\omega/\pi');
ylabel('|H(e^j^\omega)|/dB');
title('Magnitude and frequency response of butterworth filter(N=6)');
grid on;

%butterworth filter of different orders in time-domain 
figure(2);
for m=1:4
    switch m
        case 1
            n = 3;
        case 2
            n = 6;
        case 3
            n = 20;
        case 4
            n = 40;     
  end 
  t=0:1:55;
  [bm,am] = butter(n,wn);
  ym = filter(bm,am,x);
  hold on;
  subplot(2,2,m)
  plot(t,x,'b',t,ym,'r');
  xlabel('n');
  ylabel('x(n)');
  title(['N= ',num2str(n)]);
  grid on;
  t=0:1:55;
end

%butterworth filter of different orders in frequency-domain 
figure(3);
for k=1:4
    switch k
        case 1
            n = 3;
        case 2
            n = 6;
        case 3
            n = 20;
        case 4
            n = 40; 
    end 
    t=0:1:55;
    [bn,an] = butter(n,wn);
    yn = filter(bn,an,x);
    X_Fn = fftshift(abs(fft(x)));
    Y_Fn = fftshift(abs(fft(yn)));
 
    hold on;
    subplot(2,2,k)
    plot(t,X_Fn,'b',t,Y_Fn,'r');
    xlabel('\omega/2\pi');
    ylabel('|X(k)|');
    title(['N= ',num2str(n)]);
    grid on;
end

matlab代做|matlab专业代做|ECG信号仿真源码程序 - highspeedlogicblog - highspeedlogic
 
  评论这张
 
阅读(5483)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017