资源描述：
《工程软件大作业-matlab在信号处理中的应用》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、嘉兴学院机电学院《工程软件》大作业报告第9页共9页MATLAB应用实验报告模板题目3基于BMD101心电信号Soc的QRS波群定位的仿真实现作业类型MATLAB应用完成人201233285128项方康联系方式xiangfangkang@163.com自己设计并实现的作业本作业是基于BMD101心电信号前端芯片的QRS波群定位的仿真实现。主要功能为将BMD101Soc发回的心电信号预处理数据包作为系统输入，做其进行滤波及特征提取。作业涉及到MATLAB的基本操作、数字滤波器的设计、QRS波群检测算法的MATLAB实现。M文件中的指令代
2、码%[可选作业3]%[可选作业3中的指令代码]clc;clearall;closeall;%%%%%%%%%-------文件读取-------%%%%%%%%%%%%%%%%%%-------补码转换-------%%%%%%%%%PATH='E:国创国创论文matlab原始数据';DATAFILE='dataLog-2015-3-15-13-49-31.txt';DATAQUANTITY=15000;%绘制坐标XXXX=zeros(DATAQUANTITY,1);XX=1:DATAQUANTITY;XX=XX';XX=X
3、X/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;whilei<=DATAQUANTITY;char(1,:)=fread(msg,10,'uint8=>char');%%文件读取ifchar(1,1)==':'&&char(1,2)==''&&
4、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);%%十进制数转为二进制字符串iftemp2(1,1)=='0'%%补码首位为-1ecg_data(i,1)=bin2dec(temp2);else%%补码首位为1forj=1:16iftemp2(1,j)=='1'temp2(1,j
5、)='0';elsetemp2(1,j)='1';endendecg_data(i,1)=-(bin2dec(temp2)+1);endi=i+1;end嘉兴学院机电学院《工程软件》大作业报告第9页共9页fseek(msg,-9,'cof');%%文件指针移动endfclose(msg);%%%%%%%%----------FFT----------%%%%%%%%%%%%%%%%%------巴特沃斯低通------%%%%%%%%%%%%%%%%%--------滤波后FFT-------%%%%%%%%%Fs=512;%采样频
6、率512HZ采样频率len=512;%fft点数%IIR滤波器设计Nf=0;%阶数Wn=0;%3dB截止频率Fp=36.5;%通带截止频率37HzFc=50;%阻带截止频率500HzRp=3;%通带波纹最大衰减为1dBRs=60;%阻带衰减为60dBF0=18;BW=33;%----------低通滤波器----------%巴特沃兹figure(1);%画图subplot(2,1,1);plot(XX(1:DISPLAYMAX-DISPLAYMIN+1,1),ecg_data(DISPLAYMIN:DISPLAYMAX,1),'b
7、lue');grid;legend('巴特沃斯低通滤波前');axis([0XX(DISPLAYMAX-DISPLAYMIN+500,1)-500018000]);xlabel('时间/秒');ylabel('幅值');[Nf,Wn]=buttord(Fp*2/Fs,Fc*2/Fs,Rp,Rs);[BbBa]=butter(Nf,Wn,'low');%调用MATLABbutter函数快速设计滤波器Hf=freqz(Bb,Ba,len/2);%绘制频率响应曲线Bf=filter(Bb,Ba,ecg_data);%进行低通滤波%By=f
8、ft(Bf,len);%对信号f1做len点FFT变换subplot(2,1,2);plot(XX(1:DISPLAYMAX-DISPLAYMIN+1,1),Bf(DISPLAYMIN:DISPLAYMAX,1),'blue');gr
