我在用matlab求心电信号的功率谱和频谱图
来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/15 19:54:17
我在用matlab求心电信号的功率谱和频谱图
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;
txt中存放的是心电图中R_R的一系列数据 但是出来的结果不对 麻烦看下问题出在哪里
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;
txt中存放的是心电图中R_R的一系列数据 但是出来的结果不对 麻烦看下问题出在哪里
matlab里面有直接求出功率谱的函数呀~用周期图法或者间接法计算就行.
cx=xcorr(x,'unbiased');
CXk=fft(cx,N);
Pxx=abs(CXk);
index=0:round(N/2-1);
k=index*Fs/nfft;
plot_Pxx=10*log10(Pxx(index+1));
figure(3);
plot(k,plot_Pxx);
这只是计算功率谱的一种方法,不知道对你有没有帮助
或者你可以直接用periodogram函数计算,函数的格式你在matlab的帮助可以看到
cx=xcorr(x,'unbiased');
CXk=fft(cx,N);
Pxx=abs(CXk);
index=0:round(N/2-1);
k=index*Fs/nfft;
plot_Pxx=10*log10(Pxx(index+1));
figure(3);
plot(k,plot_Pxx);
这只是计算功率谱的一种方法,不知道对你有没有帮助
或者你可以直接用periodogram函数计算,函数的格式你在matlab的帮助可以看到
我在用matlab求心电信号的功率谱和频谱图
用matlab求信号的频谱
在MATLAB中信号分析的频谱与功率谱密度相似是怎么回事?
用MATLAB实现函数的频谱图
matlab sinx的频谱图
在matlab中,如何画出一个信号的频谱图?
为什么我用matlab做的这个1hz信号的频谱图不对啊,
MATLAB里用傅里叶变换求频谱
频谱和功率谱密度是什么关系?
matlab 画频谱MATLAB画时域离散信号频谱时,时域坐标和频谱坐标怎么对应.求程序
跪求:matlab频谱分析,如何得到功率谱图? 离散数据信号为p,长度为n.请高手指教!
求matlab周期三角波信号频谱分析的代码,能画出三角波信号、幅度谱和相位谱.