基于matlab心电信号的功率谱和频谱分析课程设计

N=128fs=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