工程软件大作业_Matlab在信号处理中的应用.docx

MATLAB应用实验报告模板题目3 基于BMD101心电信号Soc的QRS波群定位的仿真实现作业类型MATLAB应用完成人201233285128 项方康联系方式xiangfangkang@163.com自己设计并实现的作业本作业是基于BMD101心电信号前端芯片的QRS波群定位的仿真实现。主要功能为将BMD101Soc发回的心电信号预处理数据包作为系统输入，做其进行滤波及特征提取。作业涉及到MATLAB 的基本操作、数字滤波器的设计、QRS波群检测算法的MATLAB实现。M文件中的指令代码%[可选作业3]%[可选作业3中的指令代码]clc;clear all;close all;%%%%%%%%%-------文件读取-------%%%%%%%%%%%%%%%%%%-------补码转换-------%%%%%%%%%PATH=E:\国创\国创论文\matlab\原始数据;DATAFILE=dataLog-2015-3-15-13-49-31.txt;DATAQUANTITY=15000; %绘制坐标 XXXX=zeros(DATAQUANTITY,1);XX=1:DATAQUANTITY;XX=XX;XX=XX/512;%显示范围DISPLAYMIN=3001;DISPLAYMAX=6500;FULLPATH=fullfile(PATH,DATAFILE);msg=fopen(FULLPATH,r);ecg_char=zeros(DATAQUANTITY,4);ecg_data=zeros(DATAQUANTITY,1);%读入文件i=1;while i=DATAQUANTITY; char(1,:)=fread(msg,10,uint8=char);%%文件读取 if char(1,1)==: char(1,2)== char(1,3)==8 char(1,4)==0 ecg_char(i,:)=[char(1,6) char(1,7) char(1,9) char(1,10)]; temp1=hex2dec(ecg_char(i,:));%%十六进制字符串转换为十进制数 temp2=dec2bin(temp1,16); %%十进制数转为二进制字符串 if temp2(1,1)==0 %%补码首位为-1 ecg_data(i,1)=bin2dec(temp2); else %%补码首位为1 for j=1:16 if temp2(1,j)==1 temp2(1,j)=0; else temp2(1,j)=1; end end ecg_data(i,1)=-(bin2dec(temp2)+1); end i=i+1; end fseek(msg,-9,cof); %%文件指针移动endfclose(msg);%%%%%%%%----------FFT----------%%%%%%%%%%%%%%%%%------巴特沃斯低通------%%%%%%%%%%%%%%%%%--------滤波后FFT-------%%%%%%%%%Fs=512;%采样频率 512HZ 采样频率len=512;%fft点数%IIR滤波器设计 Nf=0; % 阶数 Wn=0; % 3dB截止频率Fp=36.5; % 通带截止频率37Hz Fc=50; % 阻带截止频率500Hz Rp=3; % 通带波纹最大衰减为1dB Rs=60; % 阻带衰减为60dBF0=18;BW=33;%----------低通滤波器----------%巴特沃兹figure(1); % 画图subplot(2,1,1);plot(XX(1:DISPLAYMAX-DISPLAYMIN+1,1),ecg_data(DISPLAYMIN:DISPLAYMAX,1),blue);grid; legend(巴特沃斯低通滤波前); axis([0 XX(DISPLAYMAX-DISPLAYMIN+500,1) -5000 18000]);xlabel(时间/秒);ylabel(幅值); [Nf,Wn]=buttord(Fp*2/Fs,Fc*2/Fs,Rp,Rs);[Bb Ba]=butter(Nf,Wn,low); % 调用MATLAB butter函数快速设计滤波器 Hf=freqz(Bb,Ba,len/2); % 绘制频率响应曲线Bf=filter(Bb,Ba,ecg_data); % 进行低通滤波 % By=fft(Bf,len); % 对信号f1做len点FFT变换 subplot(2,1,2);plot(XX(1:DISPLAYMAX-DISPLAYMIN+1,1),Bf(DISPLAYM