f7_8.m

来自「《时频分析技术》这本书中的MATLAB CODE」· M 代码 · 共 46 行

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');