標(biāo)題: 關(guān)于matlab頻域心率計(jì)算 [打印本頁(yè)]
作者: permanent    時(shí)間: 2022-12-22 16:21
標(biāo)題: 關(guān)于matlab頻域心率計(jì)算
本帖最后由 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)+"次");

幅頻特性.png (12 KB, 下載次數(shù): 55)

振幅圖

振幅圖

代碼截圖.png (22.69 KB, 下載次數(shù): 59)

代碼截圖

代碼截圖

DataFIR.zip

588 Bytes, 下載次數(shù): 3

心電數(shù)據(jù)
作者: glinfei    時(shí)間: 2022-12-22 21:37
先看看采樣個(gè)數(shù)和頻率，再直接看看分解圖，再處理。另外，直流分量應(yīng)該是0頻率的值處理來(lái)的，不用管它啊
作者: permanent    時(shí)間: 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é)果么
作者: glinfei    時(shí)間: 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)看看，就知道怎么處理了
作者: permanent    時(shí)間: 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)題了么還是我的精度不夠呀？

幅頻特性曲線.png (12 KB, 下載次數(shù): 61)

您好，這個(gè)是快速傅里葉變換后畫(huà)的圖

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

它是周期對(duì)稱的，而且你只要0到3Hz區(qū)間，放大，應(yīng)該有四個(gè)峰值。我覺(jué)得是你頻率區(qū)間取的太寬造成的誤差。
作者: Y97    時(shí)間: 2022-12-23 17:29
你把快速傅立葉分解振幅圖譜打出來(lái)看一下，就知道怎么去處理了
作者: permanent    時(shí)間: 2022-12-23 18:39
glinfei 發(fā)表于 2022-12-23 16:57
它是周期對(duì)稱的，而且你只要0到3Hz區(qū)間，放大，應(yīng)該有四個(gè)峰值。我覺(jué)得是你頻率區(qū)間取的太寬造成的誤差。

您好，我這個(gè)程序是先做心電數(shù)據(jù)的傅里葉變換，然后因?yàn)樽鐾旰笫亲笥覍?duì)稱的，所以我把幅度譜取了一半的長(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)題，謝謝您了
作者: permanent    時(shí)間: 2022-12-23 18:42
glinfei 發(fā)表于 2022-12-23 16:57
它是周期對(duì)稱的，而且你只要0到3Hz區(qū)間，放大，應(yīng)該有四個(gè)峰值。我覺(jué)得是你頻率區(qū)間取的太寬造成的誤差。

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

幅頻特性.png (12 KB, 下載次數(shù): 60)

振幅圖

振幅圖




歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1