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

数字信号处理 实验报告 实验一：相关运算和FFT的Matlab实现 实验操作记录 实验步骤 完成情况 签字确认 心电图数据的读入与显示 工频噪声信号的生成与显示 退化数据的生成、显示与保存 干扰信号频率的检测 扩展实验 实验要求 实验内容 在心电图数据中加入干扰工频噪声，然后利用相关运算和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、实验结果 三、退化数据的生成、显示与保存 1、实验代码及说明 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×a?? %?¨á￠input.dat???t′?′￠í??ˉêy?Y1?dspê1ó? N=1024; n=0:N-1; t=n/fs; plot(t,x1); xlabel(t/s) ylabel(U/mv)