实验一医学信号.docxVIP

西安交通大学实验报告课 程 医学信号处理 第页 共 页 系 别 生物医学工程 实 验 日 期 2016年3 月 21日 专业班级 医电31 班 组别________ 交 报 告 日 期 2016年4 月 18日 姓 名 麻聃 学号 2131201005 报 告退 发 ( 订正 、 重做 ) 同 组 人 韩瑞芳 教师审批签字 实验一 用FFT 对信号作频谱分析一、实验目的应用离散傅里叶变换DFT 分析模拟信号x(t)的频谱，深刻理解利用DFT 分析模拟信号频谱的原理、分析过程中出现的现象及解决方法。二、实验原理与内容连续周期信号相对于离散周期信号，连续非周期信号相对于离散非周期信号，都可以通过时域抽样定理建立相互关系。因此，在离散信号的DFT 分析方法基础上，增加时域抽样的步骤，就可以实现连续信号的DFT 分析。1.利用 FFT 计算连续周期信号的频谱已知：f1=100Hz f2=110Hz1）设fs≧4fh=4×110=440Hz 取 fs=500Hz2）①取 N1=36，对~x(t)抽样得 x1(n);②用N1 点hamming 窗乘以x1(n)得x2(n)；③取 N2=50，对~x(t)抽样得 x3(n)。3）分别对x1(n)、x2(n)、x3(n)作FFT 得z1(k)、z2(k)、z3(k)数字频谱。4）分别用stem 语句显示各自的幅度谱。横坐标为f。程序及结果：n1=0:35;f1=100;f2=110;fs=500;n2=0:49;x=cos(2*pi*f1/fs*n1)+0.75*cos(2*pi*f2/fs*n1);x1=cos(2*pi*f1/fs*n2)+0.75*cos(2*pi*f2/fs*n2);x2=x.*hamming(36);z1=fft(x1);z=fft(x);z2=fft(x2);subplot(1,3,1);stem(n1(1:18),abs(z(1:18)))subplot(1,3,2);stem(n1(1:18),abs(z2(1:18)))subplot(1,3,3);stem(n2(1:25),abs(z1(1:25)))5）分析讨论：①z1(k)、z2(k)、z3(k)数字频谱有无频谱泄漏？为什么？答：有频谱泄漏，因为只截取了信号的有限点数，相当于在时域信号上加窗，频域上做卷积，使得在以原来频率为中心的频带范围内都有谱线出现，这样能量就泄露了。②对x1(n)加hamming 窗后，频谱z2 相对z1 有何变化？泄漏有无减少？答：z2相对于z1来说频谱泄漏减少，但是主瓣变宽了。③ z1(k)、z2(k)、z3(k)哪一个最接近真实的频谱？请说明理由。答：z3最接近真实的频谱，因为其采样点数最多，泄漏较少，频率分辨率好。2.利用FFT 计算保存在数据文件shiyanyi.mat 中的信号的频谱。该数据抽样频率fs=250Hz1） 使用load 语句将数据文件shiyanyi 中的变量ecg1,ecg2,ecg3 调入matlab内存空间。其中ecg2 与ecg3 分别为从ecg1 中截取的一段。2） ①对ecg2 作FFT 得频谱z1；②对ecg2 中的下标从1 到100 共100 个点作FFT 得频谱z2；③对ecg2 中的下标从101 到150 共50 个点作FFT 得频谱z3；④用plot 语句分别显示z1 、z2、z3 的幅度谱。横坐标为f 。3） ①对ecg3 作DFT 得频谱z4；②将ecg3 乘以hamming 窗，然后作FFT 得z5；③用plot 语句分别显示z1、z4、z5 的幅度谱。横坐标为f 。程序及结果：load ecgfs=250;z1=fft(ecg2);z2=fft(ecg2(1:100));z3=fft(ecg2(101:150));subplot(3,2,2);plot((0:125)*fs/251,abs(z1(1:126)));title(z1)subplot(3,2,4);plot((0:49)*fs/100,abs(z2(1:50)));title(z2)subplot(3,2,6);plot((0:24)*fs/50,abs(z3(1:25)));title(z3)subplot(3,2,1);plot(ecg2);title(ecg2)subplot(3,2,3);plot(ecg2(1:100));title(ecg2(1:100))subplot(3,2,5);plot(ecg2(101:150));title(ecg2(101:150))z4=fft(ecg2);z5=fft(ecg3);z6=fft(ecg3.*hamming(196));figuresubplot(3,2,2);plot((0:125)*fs/251