生物医学信号处理实验指导山东大学控制学院生医专业2011-2-5实验目录实验一 离散时间信号与系统分析3实验二 离散时间信号与系统的Z变换分析11实验三IIR滤波器的设计与信号滤波17实验四 用窗函数法设计FIR数字滤波器21实验五 用FFT作谱分析24实验六综合实验27实验一 离散时间信号与系统分析一、实验目的1掌握离散时间信号与系统的时域分析方法。2掌握序列傅氏变换的计算机实现方法，利用序列的傅氏变换对离散信号、系统及系统响应进行频域分析。 3熟悉理想采样的性质，了解信号采样前后的频谱变化，加深对采样定理的理解。 二、实验原理1离散时间系统一个离散时间系统是将输入序列变换成输出序列的一种运算。若以来表示这种运算，则一个离散时间系统可由下图来表示： 图 离散时间系统即输出与输入之间关系用下式表示 离散时间系统中最重要、最常用的是线性时不变系统。2离散时间系统的单位脉冲响应设系统输入，系统输出的初始状态为零，这是系统输出用表示，即，则称为系统的单位脉冲响应。可得到： 该式说明线性时不变系统的响应等于输入序列与单位脉冲序列的卷积。3连续时间信号的采样采样是从连续信号到离散时间信号的过渡桥梁，对采样过程的研究不仅可以了解采样前后信号时域何频域特性发生的变化以及信号内容不丢失的条件，而且有助于加深对拉氏变换、傅氏变换、Z变换和序列傅氏变换之间关系的理解。对一个连续时间信号进行理想采样的过程可以表示为信号与一个周期冲激脉冲的乘积，即：其中，是连续信号的理想采样，是周期冲激脉冲设模拟信号，冲激函数序列以及抽样信号的傅立叶变换分别为、和，即根据连续时间信号与系统中的频域卷积定理，式（2.59）表示的时域相乘，变换到频域为卷积运算，即 其中 由此可以推导出 由上式可知，信号理想采样后的频谱是原来信号频谱的周期延拓，其延拓周期等于采样频率。根据香农定理，如果原信号是带限信号，且采样频率高于原信号最高频率的2倍，则采样后的离散序列不会发生频谱混叠现象。4有限长序列的分析对于长度为N的有限长序列，我们只观察、分析在某些频率点上的值。 一般只需要在之间均匀的取M个频率点，计算这些点上的序列傅立叶变换： 其中，。是一个复函数，它的模就是幅频特性曲线。三、主要实验仪器及材料微型计算机、Matlab软件（或TC编程环境）。四、实验内容1知识准备认真复习离散信号与系统、单位脉冲响应、抽样定理等有关内容，阅读本实验原理与方法。2编制信号产生子程序，用于产生实验中要用到的信号序列（1）单位脉冲序列单位脉冲序列（2）系统单位脉冲响应序列（3）理想采样信号序列对信号进行理想采样，可以得到一个理想的采样信号序列，。其中为幅度因子，a是衰减因子，W是频率，为采样周期。这几个参数要在实验过程中输入，以产生不同的。首先产生理想采样信号序列，使444.128, =50,=50。然后改变参数A1，=0.4,=2.0734，产生理想采样信号序列。3离散信号、系统和系统响应的分析观察信号xb(n)和系统hb(n)的时域和频域特性；利用线性卷积求信号通过系统以后的响应。比较系统响应和信号的时域和幅频特性。注意它们之间有无差异，绘出图形。4分析理想采样信号序列的特性产生理想采样信号序列，使：（1）首先选用采样频率为1000Hz，T=1/1000，观察所得理想采样信号的幅频特性，在折叠频率以内和给定的理想幅频特性无明显差异，并作记录。（2）改变采样频率为300Hz，T=1/300，观察所得理想采样信号的幅频特性曲线的变化，并作记录。（3）进一步减小采样频率为200Hz，T=1/200，观察频谱混叠现象是否明显存在，说明原因，并记录此时的幅频特性曲线。5. 卷积定律的验证。采用参数444.128, =50,=50， T=1/1000，将和系统的傅氏变换相乘，直接求得，将得到的幅频特性曲线和先求后再求得的幅频特性曲线进行比较，观察二者有无差异。验证卷积定律。五、思考题1线性时不变系统的输出的长度与输入和系统的单位冲激响应的长度有什么关系？2. 对信号进行理想抽样时，抽样频率不同，相应理想采样序列傅立叶变换频谱的数字频率度量是否都相同7它们所对应的模拟频率是否相同?为什么？六、实验报告要求 1简述实验原理及目的。2. 总结在上机实验内容中要求比较时域、幅频曲线差异部分内容的结果，定性分析它们正确与否，并简要说明这些结果的含义。3.总结实验所得主要结论。4.简要回答思考题。附：实验一 离散时间系统分析参考程序部分参考程序四、2.（1）,（2）的参考程序：%单位脉冲序列的时域和幅频特性%在MatLab中，这一函数可以用zeros函数实现：n=1:50; x=zeros(1,50); %MatLab中数组的下标从1开始x(1)=1;close all;subplot(3,1,1);stem(x);title(单位冲击信号序列x(n);k=-24:25;X=x*(exp(-j*pi/25).(n-1)*k);magX=abs(X);subplot(3,1,2);stem(magX);title(单位冲击信号序列的幅度谱);angX=angle(X);subplot(3,1,3);stem(angX);title(单位冲击信号序列的相位谱);%以下是的时域和幅频特性n=1:50;x=zeros(1,50);x(1)=1;x(2)=2.5;x(3)=2.5;x(4)=1;close all;subplot(3,1,1);stem(x);title(系统单位脉冲响应信号序列);k=-24:25;X=x*(exp(-j*pi/25).(n-1)*k);magX=abs(X);subplot(3,1,2);stem(magX);title(系统频率响应的幅度谱);angX=angle(X);subplot(3,1,3);stem(angX);title(系统频率响应的相位谱)%以下矩形脉冲序列的时域和幅频特性(没有要求做，仅作参考)n=1:50; x=sign(sign(10-n)+1);close all;subplot(3,1,1);stem(x);title(矩形脉冲序列);k=-24:25;X=x*(exp(-j*pi/25).(n-1)*k);magX=abs(X);subplot(3,1,2);stem(magX);title(矩形脉冲序列傅立叶变换的幅度谱);angX=angle(X);subplot(3,1,3);stem(angX);title(矩形脉冲序列傅立叶变换的相位谱);四、2.（3）的参考程序：n=1:50; %定义序列的长度是50A=444.128; %设置信号有关的参数a=50*sqrt(2.0)*pi;T=0.001; %采样率 w0=50*sqrt(2.0)*pi; %符号在MatLab中不能输入，用w代替x=A*exp(-a*(n-1)*T).*cos(w0*(n-1)*T); %pi是MatLab中定义的%close allfiguresubplot(3,1,1);stem(n-1,x); %绘制x（n）的图形title(理想采样信号序列x(n); %设置结果图形的标题k=-24:25;%W=(pi/25)*k;X=x*(exp(-j*pi/25).(n-1)*k);magX=abs(X); %绘制x（n）的幅度谱subplot(3,1,2);stem(n-1,magX);title(理想采样信号序列的幅度谱);angX=angle(X); %绘制x（n）的相位谱subplot(3,1,3);stem(n-1,angX);title(理想采样信号序列的相位谱);四、3. 的参考程序（要用代替）：%卷积计算%在MatLab中提供了卷积函数conv，即y=conv（n，h），调用十分方便。%信号和系统单位脉冲响应的卷积n=1:50;hb=zeros(1,50);hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1;close all;subplot(3,1,1);stem(hb);title(系统单位脉冲响应hbn);m=1:50;T=0.001;A=444.128;a=50*sqrt(2.0)*pi;w0=50*sqrt(2.0)*pi;x=A*exp(-a*(m-1)*T).*cos(w0*(m-1)*T);subplot(3,1,2);stem(x);title(输入信号xn);y=conv(x,hb);subplot(3,1,3);stem(y);title(输出信号yn);四、5. 的参考程序%参考程序clearhb=zeros(1,50);hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1;m=1:50;T=0.001;A=444.128;a=50*sqrt(2.0)*pi;w0=50*sqrt(2.0)*pi;x=A*exp(-a*(m-1)*T).*cos(w0*(m-1)*T);n=1:50;k=-24:25;X=x*(exp(-j*pi/25).(n-1)*k);magX=abs(X);%绘制信号x(n)傅立叶变换的幅度谱subplot(3,2,1);stem(magX);title(输入信号傅立叶变换的幅度谱);angX=angle(X); %绘制x(n)的相位谱subplot(3,2,2);stem(angX);title(输入信号傅立叶变换的相位谱);Hb=hb*(exp(-j*pi/25).(n-1)*k);magHb=abs(Hb); %绘制hb(n)的幅度谱subplot(3,2,3);stem(magHb);