f7_7.m

matlab时频分析配套源码第6章的程序

%正常人心律信号
%fid=fopen('normalecg.txt','r');
%充血性心脏心律信号
fid=fopen('heartfailure.txt','r');
%心室纤维颤动心律信号
%fid=fopen('atrialfib.txt','r');
x=fscanf(fid,'%f %f %f',[300 3]);
fclose(fid);

%采样频率
fs=1/0.008;% 正常信号下
%fs=1/0.004;%充血性心脏
%fs=1/0.004;%心室纤维颤动

N=length(x);
for i=0:N/3-1
    t(i+1)=x(3*i+1);
    sig1(i+1)=x(3*i+2);
    sig2(i+1)=x(3*i+3);
end

%时域波形
figure(1);
subplot(2,2,1);
plot(t,sig1,'LineWidth',2);
xlabel('时间 t/s');
ylabel('电压 V/mV');

%估计瞬时频率
sig=hilbert(sig1);
sig=sig-mean(sig);
[ifr,tt]=instfreq(sig',2:length(sig)-1);

%显示结果
subplot(2,2,2);
plot(tt/fs,(1-ifr)*fs/2, 'LineWidth',2);
xlabel('时间 t/s');
ylabel('瞬时频率 f/Hz');