test2.m

来自「用matlab设计FIR和IIR滤波器的设计函数与例程」· M 代码 · 共 22 行

clear;
fs=400;
gs=50;

M=fs/gs;
b=[1 zeros(1,95) -49 zeros(1,3) -49 zeros(1,95) 1];

a=[1 zeros(1,3) 1];
figure(1);
freqz(b,a,1024,fs);
figure(2);
zplane(b,a);

filename='D:\数据\蔡刿.mat';
load(filename);
ecgdata=data(:,2);
ecg=ecgdata(1:2*fs);


x=filter(b,a,ecg);
figure(3);
plot(x);