اعمال فیلتر بر روی ارایه ecg

[hd_uni_pro]

با سلام
دوستان من می خوام فیلتر طراحی شده با ابزار fdatool رو که یک فیلتر چبیشف نوع 1 مرتبه 4 هستش رو روی داده ای سی جی که از فیزیو نت گرفتم اعمال کنم با چه دستوری می تونم این کارو بکنم ؟ من داده رو در متلب لود کردم و همچنین کد فیلتر دارم حالا چه طور فیلتر رو اعمال کنم :
کد چبیشف که ابزار متلب ساخته :
function Hd = chebyshef 1 order 4
%CHEBYSHEF 1 ORDER 4 Returns a discrete-time filter object.

%
% MATLAB Code
% Generated by MATLAB(R) 7.11 and the Signal Processing Toolbox 6.14.
%
% Generated on: 13-Mar-2012 1005
%

% Chebyshev Type I Bandpass filter designed using FDESIGN.BANDPASS.

% All frequency values are in Hz.
Fs = 55; % Sampling Frequency

N = 4; % Order
Fpass1 = 6; % First Passband Frequency
Fpass2 = 16; % Second Passband Frequency
Apass = 1; % Passband Ripple (dB)

% Construct an FDESIGN object and call its CHEBY1 method.
h = fdesign.bandpass('N,Fp1,Fp2,Ap', N, Fpass1, Fpass2, Apass, Fs);
Hd = design(h, 'cheby1'

% [EOF]

[BlueAbble]

پاسخ : اعمال فیلتر بر روی ارایه ecg

سلام
شما اصل کار را انجام داده اید. تنها یک دستور مانده است که باید اجرا شود و آن دستور filter است.
اگر فرض کنیم که داده های ابتدایی در آرایه x ریخته شده و میخواهید داده های فیلتر شده را در آرایه y بریزید، باید از دستور زیر استفاده کنید:

y = filter(Hd,x);

به عنوان مثال کد زیر را اجرا کنید:

t = 0:1/Fs:5;
x = 5*sin(2*pi*30*t);
y = filter(Hd,x);
plot(t,x,'blue',t,y,'red');

فرکانس سینوسی را در خط دوم تغییر دهید تا اثر فیلتر را مشاهده کنید

[matlab]

پاسخ : اعمال فیلتر بر روی ارایه ecg

این نکته نباید فراموش شود که فرکانس نمونه برداری در طراحی فیلتر باید با فرکانس نمونه برداری سیگنال(ٍECG) یکی باشه، به همین جهت بهتره تابعی که ضرایب فیلترHd را ایجاد می کند به جای آنکه Fs را ثابت تعریف بکنیم، به عنوان ورودی در تابع تعریف بشود.