实验一相关运算和FFT的Matlab实现与应用.doc

实验一相关运算和FFT的Matlab实现与应用.doc

  1. 1、本文档共8页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验一相关运算和FFT的Matlab实现与应用

数字信号处理实验报告 实验一 相关运算和FFT的Matlab实现与应用 PAGE 7 MACROBUTTON MTEditEquationSection2 SEQ MTEqn \r \h \* MERGEFORMAT SEQ MTSec \r 1 \h \* MERGEFORMAT SEQ MTChap \r 1 \h \* MERGEFORMAT 数字信号处理 实验报告 实验一:相关运算和FFT的Matlab实现与应用 实验人: ++++ 学号: 090802076 指导教师: 谢晓春 实验日期: 2012年 月 日 实验地点: 6-602 实验预习 实验操作 实验报告 实验成绩 实验操作记录 实验步骤 完成情况 签字确认 心电图数据的读入与显示 工频噪声信号的生成与显示 退化数据的生成、显示与保存 干扰信号频率的检测 扩展实验 实验要求 实验内容 在心电图数据中加入干扰工频噪声,然后利用相关运算和FFT从中检测出干扰周期信号的频率。 实验目的 理解并掌握相关运算在周期信号检测中的应用; 掌握在matlab环境中相关运算和FFT的实现; 掌握在matlab环境中的数据读取、写入及显示的方法。 三、实验步骤 读入心电图数据myECG.mat,并将其显示; clc;close all; load(ecgdata.mat); x=mydata; plot(x); (提示:将ecgdata.mat拷到工程目录下,直接查看变量x中数据) 生成幅度为0.1mV的50Hz工频噪声信号; clc;close all; fs=360; load(ecgdata.mat); x=mydata; noise=0.1*sin(2*pi*50*(0:length(x)-1)/fs); plot(noise); 将噪声信号与心电图数据相加,得到退化的心电图数据,并将其保存为DSP能够处理的数据格式,为后面的实验做准备; clc;close all; clear; fs=360; load(ecgdata.mat); x=mydata; noise=0.1*sin(2*pi*50*(0:length(x)-1)/fs); x1=x+noise; x1=x1;%x1转置 %建立input.dat文件存储退化数据供dsp使用 fid=fopen(input.dat,w); %打开文件,w是将此文件定义为可写的,fid是此文件的整数标示 fprintf(fid,1651 4 0 1 0\n); fprintf(fid,%f\n,x1); fclose(fid); 利用相关运算和FFT检测干扰信号的频率; clc;clear;close all; fs=360; N=1024; n=0:N-1;t=n*fs/N; load(ecgdata.mat); x=mydata; noise=0.1*sin(2*pi*50*(0:length(x)-1)/fs); x1=x+noise; x2=xcorr(x1,noise);%相关运算 xk=fft(x2,N); plot(t,abs(xk)) 扩展实验:编程直接获得干扰信号频率。 四、相关说明 1、心电图数据的采样频率是360Hz,共有1024点数据,数据的单位是mV; 2、在绘图时要正确标注横坐标与纵坐标; 3、DSP的dat文件格式是由文件头(1651 4 0 1 0)和数据构成,其中1651为供DSP识别的固定值,4表示文件中的数据为十进制浮点型。 4、实验中涉及到的matlab指令有 数据读写指令:load,write,fopen, fprintf, fclose; 绘图指令:figure, plot, xlabel, ylabel, title,axis; FFT运算: fft, fftshift; 相关运算: xcorr。 实验结果 一、心电图数据的读入与显示 1、实验代码及说明 clc;close all; N=1024; fs=360; n=0:N-1; t=n/fs; load(ecgdata.mat); x=mydata; plot(t,x); xlabel(t/s); ylabel(U/mV); title(D?μ?í?); 2、实验结果 二、工频噪声信号的生成与显示 1、实验代码及说明 clc;close all; fs=360; N=1024; n=0:N-1; t=n/fs; load(ecgdata.mat); x=mydata; noise=0.1*sin(2*pi*50*(0:length(x)-1)/fs); plot(t,noise); xlabel(t/s) ylabel(U/mv) title(1¤?μ??éùD?o?); 2、实验结果 三、退

文档评论(0)

shenlan118 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档