资源说明:求助Matlab关于心电信号中基线纠漂程序出错-ECG.txt
要对心电心电信号中基线纠漂,使用截止频率为0.7Hz的高通滤波器,一下程序有什么问题,怎么不出来图阿
%载入信号
clear all
clc
ECG=load;
N=1001;
Fs=250;
n=ECG;
y=ECG;
% 基于频率抽样方法1实现Ⅰ型FIR数字高通滤波器
T1=0.1095;
T2=0.598;
alpha= /2;
l=0:N-1;
wl= *l;
Hrs=[zeros,T1,T2,ones,T2,T1,zeros]; %理想振幅采样响应
Hdr=[0,0,1,1];
wdl=[0,0.75,0.75,1];
k1=0:floor/2);
k2=floor/2) 1:N-1;
angH=[-alpha*/N*k1,alpha*/N*]; %相位约束条件
Hdk=Hrs.*exp; %构成Hd
h2=real); %实际单位冲激响应
[db,mag,pha,w]=freqz_m2;
[Hr,ww,a,L]=hr_type1; %实际振幅响应
y1=filter;
plot(y1);
function [db, mag, pha, w]=freqz_m2;
%Modified version of freqz subroutine
[H,w]=freqz;
H=)'; w=)';
mag=abs;
db=20*log10/max);
pha=angle;
%grd=grpdelay;
function [Hr, w, a,L]=hr_type1;
%Computes Amplitude response of Type-1 LP FIR filter
M=length;
L=/2;
a=[h 2*h];
n=[0:L];
w=[0:500]'*pi/500;
Hr=cos*a';
本源码包内暂不包含可直接显示的源代码文件,请下载源码包。