心电信号的频率普和频谱分析

N=128

fs=100%设定采样频率

N=128;

n=0:N-1;

%进行FFT 变换并做频谱图

x=load('F:\ZB.txt');

y=fft(x,N);%进行fft 变换

mag=abs(y);%求幅值

f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换 figure(1);

subplot(2,2,1);

plot(f,mag);%做频谱图

axis([0,100,0,50]);

xlabel('频率(Hz)');

ylabel('幅值');

title('心电信号幅频谱图');

grid;

%求功率谱

x=load('F:\ZB.txt');

c=xcorr(x);

p=fft(c,N);

figure(1);

subplot(2,2,2);

plot(f,p);

xlabel('频率(Hz)');

ylabel('功率谱');

title('心电信号功率谱');

grid;

N=256

N=128

fs=100%设定采样频率

N=128;

n=0:N-1;

%进行FFT 变换并做频谱图

x=load('F:\ZB.txt');

y=fft(x,N);%进行fft 变换

mag=abs(y);%求幅值

f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换 figure(1);

subplot(2,2,1);

plot(f,mag);%做频谱图

axis([0,100,0,50]);

xlabel('频率(Hz)');

ylabel('幅值');

title('心电信号幅频谱图');

grid;

%求功率谱

x=load('F:\ZB.txt');

c=xcorr(x);

p=fft(c,N);

figure(1);

subplot(2,2,2);

plot(f,p);

xlabel('频率(Hz)');

ylabel('功率谱');

title('心电信号功率谱');

grid;

N=256