- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 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)