%正常人心律信號
%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');