關(guān)于matlab頻域心率計(jì)算

ID:1058749 · 發(fā)表于 2022-12-22 16:21

本帖最后由 permanent 于 2022-12-23 19:11 編輯

現(xiàn)在要在頻域用matlab做一個(gè)心率計(jì)算，我.m文件里是我在arduino濾波完以后的數(shù)據(jù),我給這個(gè)數(shù)組存的變量名是DataFIR,然后我基本思路是先給這個(gè)數(shù)據(jù)去掉直流分量，然后用fft()這個(gè)函數(shù)給他做了一個(gè)傅里葉變換，然后找到最大值處對(duì)應(yīng)的頻率，乘以60就是我要的每分鐘的心跳次數(shù)，但是運(yùn)行結(jié)果差的太多了，竟然是327！！！求助大佬看看我的程序是哪里有問(wèn)題么！！！改了好久了嗚嗚嗚！！附件是心電數(shù)據(jù)，圖片是代碼截圖和運(yùn)行結(jié)果

clear;
fs=500;
L=fs;
load('DataFIR.mat')
x=DataFIR;
n=1:1:length(x);
y0=x(n)-mean(x);%去掉直流分量
y1=abs(fft(y0));
m = y1(1:length(y1)/2);
N= length(x);
% for i=1:length(m)
%    if m(i)==max(m)%找峰值
%       f=i*(L/N);
%    end
% end
f=find(m==max(m))*L/N;
xinlv=round(f*60)
m = y1(1:length(y1));
x= (0:N-1)*(L/N);
figure;
plot(y0)
figure;
plot(x,m)
xlabel('Frequency (Hz)');ylabel('Magnitude (dB)');title('Filter output');
text(200,2000,"心率"+num2str(xinlv)+"次");

ID:844772 · 發(fā)表于 2022-12-22 21:37

先看看采樣個(gè)數(shù)和頻率，再直接看看分解圖，再處理。另外，直流分量應(yīng)該是0頻率的值處理來(lái)的，不用管它啊

ID:1058749 · 發(fā)表于 2022-12-22 22:42

glinfei 發(fā)表于 2022-12-22 21:37
先看看采樣個(gè)數(shù)和頻率，再直接看看分解圖，再處理。另外，直流分量應(yīng)該是0頻率的值處理來(lái)的，不用管它啊

您好，這個(gè)采樣頻率是題目要求固定的不能變，請(qǐng)問(wèn)這個(gè)采樣個(gè)數(shù)和采樣頻率會(huì)影響我的心率計(jì)算結(jié)果么

ID:844772 · 發(fā)表于 2022-12-23 07:22

permanent 發(fā)表于 2022-12-22 22:42
您好，這個(gè)采樣頻率是題目要求固定的不能變，請(qǐng)問(wèn)這個(gè)采樣個(gè)數(shù)和采樣頻率會(huì)影響我的心率計(jì)算結(jié)果么

采樣頻率大于一倍就行了，采樣個(gè)數(shù)只是影響精度而已，你先把快速傅立葉分解振幅圖譜打出來(lái)看看，就知道怎么處理了

ID:1058749 · 發(fā)表于 2022-12-23 09:13

glinfei 發(fā)表于 2022-12-23 07:22
采樣頻率大于一倍就行了，采樣個(gè)數(shù)只是影響精度而已，你先把快速傅立葉分解振幅圖譜打出來(lái)看看，就知道怎 ...

您好，這是快速傅里葉變換后的幅頻特性曲線，我把最高的頻率分量點(diǎn)找到乘以60是不是就是每分鐘的心跳次數(shù)，可是我這個(gè)找到算出來(lái)是327，應(yīng)該是60才對(duì)，是我的程序哪里出現(xiàn)問(wèn)題了么還是我的精度不夠呀？

ID:844772 · 發(fā)表于 2022-12-23 16:57

permanent 發(fā)表于 2022-12-23 09:13
您好，這是快速傅里葉變換后的幅頻特性曲線，我把最高的頻率分量點(diǎn)找到乘以60是不是就是每分鐘的心跳次數(shù) ...

它是周期對(duì)稱(chēng)的，而且你只要0到3Hz區(qū)間，放大，應(yīng)該有四個(gè)峰值。我覺(jué)得是你頻率區(qū)間取的太寬造成的誤差。

ID:1058494 · 發(fā)表于 2022-12-23 17:29

你把快速傅立葉分解振幅圖譜打出來(lái)看一下，就知道怎么去處理了

ID:1058749 · 發(fā)表于 2022-12-23 18:39

glinfei 發(fā)表于 2022-12-23 16:57
它是周期對(duì)稱(chēng)的，而且你只要0到3Hz區(qū)間，放大，應(yīng)該有四個(gè)峰值。我覺(jué)得是你頻率區(qū)間取的太寬造成的誤差。

您好，我這個(gè)程序是先做心電數(shù)據(jù)的傅里葉變換，然后因?yàn)樽鐾旰笫亲笥覍?duì)稱(chēng)的，所以我把幅度譜取了一半的長(zhǎng)度，然后通過(guò)max函數(shù)找到一個(gè)最大值，通過(guò)find函數(shù)找對(duì)應(yīng)頻率，并不是看圖自己找的，所以振幅圖的區(qū)間應(yīng)該并不影響計(jì)算結(jié)果，因?yàn)閒ind找到的應(yīng)該就是那個(gè)精確的最大值，而振幅圖就是我的一個(gè)參考，那能不能麻煩您幫我看一下這到底是哪里的問(wèn)題，謝謝您了

ID:1058749 · 發(fā)表于 2022-12-23 18:42

glinfei 發(fā)表于 2022-12-23 16:57
它是周期對(duì)稱(chēng)的，而且你只要0到3Hz區(qū)間，放大，應(yīng)該有四個(gè)峰值。我覺(jué)得是你頻率區(qū)間取的太寬造成的誤差。

您好，我這個(gè)程序是先做心電數(shù)據(jù)的傅里葉變換，然后因?yàn)樽鐾旰笫亲笥覍?duì)稱(chēng)的，所以我把幅度譜取了一半的長(zhǎng)度，然后通過(guò)max函數(shù)找到那個(gè)最大值，find找對(duì)應(yīng)頻率，并不是看圖自己找的，所以振幅圖的區(qū)間應(yīng)該并不影響計(jì)算結(jié)果，因?yàn)閙ax找到的應(yīng)該就是那個(gè)精確的最大值，而振幅圖就是我的一個(gè)參考，那能不能麻煩您幫我看一下這到底是哪里的問(wèn)題，謝謝您了

ID:1058749 · 發(fā)表于 2022-12-23 18:54

Y97 發(fā)表于 2022-12-23 17:29
你把快速傅立葉分解振幅圖譜打出來(lái)看一下，就知道怎么去處理了

您好，這是我的振幅圖，但是顯然最高點(diǎn)對(duì)應(yīng)的頻率乘以60并不是正確的計(jì)算結(jié)果，可不可以麻煩您幫我看一下程序是哪里出現(xiàn)了問(wèn)題么，謝謝您了。

帳號(hào)		自動(dòng)登錄	找回密碼
密碼			立即注冊(cè)

久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

關(guān)于matlab頻域心率計(jì)算

相關(guān)帖子