?? f7_7.m
字號:
%正常人心律信號
%fid=fopen('normalecg.txt','r');
%充血性心臟心律信號
fid=fopen('heartfailure.txt','r');
%心室纖維顫動心律信號
%fid=fopen('atrialfib.txt','r');
x=fscanf(fid,'%f %f %f',[300 3]);
fclose(fid);
%采樣頻率
fs=1/0.008;% 正常信號下
%fs=1/0.004;%充血性心臟
%fs=1/0.004;%心室纖維顫動
N=length(x);
for i=0:N/3-1
t(i+1)=x(3*i+1);
sig1(i+1)=x(3*i+2);
sig2(i+1)=x(3*i+3);
end
%時域波形
figure(1);
subplot(2,2,1);
plot(t,sig1,'LineWidth',2);
xlabel('時間 t/s');
ylabel('電壓 V/mV');
%估計瞬時頻率
sig=hilbert(sig1);
sig=sig-mean(sig);
[ifr,tt]=instfreq(sig',2:length(sig)-1);
%顯示結果
subplot(2,2,2);
plot(tt/fs,(1-ifr)*fs/2, 'LineWidth',2);
xlabel('時間 t/s');
ylabel('瞬時頻率 f/Hz');
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -