Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 其它 > 资料存档
资料存档 资料存档
回复
 
主题工具 显示模式
旧 2019-11-27, 05:26   #1
poster
高级会员
 
注册日期: 2019-11-21
帖子: 3,006
声望力: 66
poster 正向着好的方向发展
默认 Extracting QT from EKG signal

Basically, I want to extract the QT from an EKG signal, for this I am following some steps from an article. Firstly, I am filtering the signal with a pass-band filter, then I have to apply the coiflet-1 CWT transform(never used it before) due to QT Pattern = (CWT Output)^2 and then choose
QRS Pattern(i) = {0 , QT Pattern(i) < A
QT Pattern(i), otherwise.
This is what I did until now, but i found myself in a moment when i don't know what to do next:



**clear all
clc
close all
[filename, pathname] = uigetfile('semnal.dat');
if isequal(filename, 0) || isequal(pathname, 0)
disp('File input canceled.');
ECG_Data = [];
else
fid=fopen(filename,'r');
end;
time=10;
f=fread(fid,2*360*time,'ubit12');
semnal=f(1:2:length(f));
figure(1)
plot(semnal)

xlabel('samples');
ylabel('time')


fs = 250; % Sampling Frequency (Hz)



N=length(semnal);
fs=250;
t=[0:N-1]/fs;
fcutlow=0.5; %low cut frequency in Hz
fcuthigh=40; %high cut frequency in Hz
[b,a]=fir1(330,[fcutlow,fcuthigh]/(fs/2),'bandpass');
semnal_filtrat=filter(b,a,semnal)
figure(2)
plot(t,semnal_filtrat)

%CTW


x = [1, 2, 3, 4];
scales = [9];
wavelet_name = 'coif2';
coefs = cwt(x,scales, wavelet_name)**


I am waiting for any idea, thank you!!





More answer...
poster 当前离线   回复时引用此帖
回复


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛禁用 表情符号
论坛启用 [IMG] 代码
论坛启用 HTML 代码



所有时间均为北京时间。现在的时间是 23:13


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.