clear;
%正常人心律信號
fid=fopen('normalecg.txt','r');
%x=fscanf(fid,' %f %f %f',[300 3]);
%充血性心臟心律信號
%fid=fopen('heartfailure.txt','r');
x=fscanf(fid,' %f %f %f',[300 3]);
%心室纖維顫動心律信號
%fid=fopen('atrialfib.txt','r');
%x=fscanf(fid,' %f %f %f',[300 3]);
fclose(fid);

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

subplot(2,1,1);
plot(t,sig1,'LineWidth',2);
xlabel('時間 t/s');
ylabel('電壓 V/mV');

subplot(2,1,2);
plot(t,sig2,'LineWidth',2);
xlabel('時間 t/s');
ylabel('電壓 V/mV');

%采樣頻率
fs=1/0.008;% 正常信號下
%fs=1/0.004;%充血性心臟
%fs=1/0.004;%心室纖維顫動

x1=sig1;
x2=sig2;

%計算互Wigner-Ville分布
[tfr,t,f]=tfrwv([x1 x2]',1:length(x1),256);

%顯示結果
figure(2)
contour(t/fs,f(1:128)*fs,tfr(1:128,:));
xlabel('時間 t/s');
ylabel('頻率 f/Hz');