医学信号处理综合实验

医学信号处理 实验指导

河北工业大学电气工程学院生物医学工程专业专用

实验一、随机信号的产生及其数字特征分析

实验项目的目的和任务：了解随机信号的产生以及随机信号的特征。 实验内容：

1. 随机信号序列的产生 2. 随机信号的数字特征分析

实验步骤：

一、随机信号序列的产生

1、 利用matlab函数rand和randn产生随机信号序列 2、 利用同余算法生成随机信号序列 同余算法：

k(i+1)=（c*k(i)）(mod M) i=0 1 2 ....m

k(0) 给定，后面的数由k与c的乘积与M取余求得，c与

M都可以任意指定，一般情况下M=2^b，k(0)为一奇数 c=8t(+/-)3 t取整数。

具体b和t取值，可由试验得出。

随机信号不重复的长度：2^(b-2)次，b取值越大，随机信号

序列长度越长

建议 b=12 即2^10次不会重复，此时c=2^9+3 即 t=2^6。

k的取值范围，0~M-1之间。k在0~M-1上是均匀分布的。

生成随机信号序列要求：

1 用同余算法产生均匀分布的随机序列，长度

一、 语音数字滤波系统设计

1. 原理

实验中所需理论包括采样定理、时域信号的FFT 分析及数字滤波器设计原理和方法 采样定理：在进行模拟/数字信号的转换过程中，当采样频率s f 大于等于信号最高截止频率h f 的2倍时，即：2s h f f ，则采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高截止频率的5～10倍；采样定理又称奈奎斯特定理。

2. 任务

试编制一MATLAB 或LabVIEW 程序，实现语音数字滤波系统的用户界面设计，在所设计的系统界面上完成的功能包括：

（1）对于任意的语音信号进行采样并对加入加性噪声的信号作频谱分析；

（2）通过频谱分析选择合适的滤波器性能指标，设计合适的数字滤波器，并对含噪音的语音信号进行数字滤波，得出滤波器的时域波形和频谱；

（3）对滤波前后的信号进行对比，分析信号的变化；回放语音信号，对滤波前后的声音进行比较。

3. 思考题

（1）试分析各种不同类型滤波器的性能。

（2）试分析语音信号中的低频、中频和高频成分的特点及适合采用的滤波器类型。

二、 音频信号处理系统设计

1. 原理

实验中所需理论包括采样定理（可参照语音数字滤波系统设计）、时域信号的FFT 分析、AM 调制解调原理及IIR

生物医学信号处理

实验指导书

2009年8月

目录

实验一 随机信号的数字特征分析............................................................................ 1 实验二 数字相关和数字卷积.................................................................................... 4 实验三 维纳－霍夫方程............................................................................................ 7 实验四

Yule-Walker方程 ........................................................................................ 11 实验一 随机信号的数字特征分析

（一）实验目的

了解随机信号的特征。

掌握随机信号的数字特征分析算法。

（二）实验原理

对于平稳各态遍历随机过程，可以用单一样本函数的时间平

数字信号处理综合实验

班级： 姓名： 学号：

一、实验目的

1．掌握MATLAB的程序设计方法；

2．掌握数字信号处理的基本理论和基本方法； 3．掌握语音信号的采集与处理方法；

4．掌握用MATLAB设计FIR和IIR数字滤波器的方法； 5．掌握用MATLAB对信号进行分析和处理的方法；

二、实验原理

y=exp(x): 以e为底的指数。

conj(x): 取x的共轭，即改变x的虚部符号。 real(x): 取复数x的实部。

rand(1,N): 生成在0和1之间均匀分布的随机序列，长度为N。 randn(1,N): 生成正态分布(高斯分布)的随机序列，长度为N。 sound(f,fs): 输入参量是音频数据向量、采样频率和转换位数。

X=fft(x,N): 计算序列x的N点FFT。如果x的长度小于N，则在x后面补零；如果x的长度大于N，则对x进行截取；如果不指定参数N，则以x的实际长度作为FFT的点数。

x=ifft(X,N): 计算序列X的N点IFFT。

Y=fftshift(X): 将序列X分成左右两部分并交换位置。

[N, Wn] = buttord(Wp, Ws, Rp, Rs, 's'): 求巴特沃

课 程 设 计

系统分析与设计

信息系统就是为了有效的传输和交换信息而对载荷信息的信号进行加工、处理和变换的设备的总称。系统根据输入、输出信号的不同分为连续系统和离散系统。其分析方法有卷积法，零输入、零状态法，变换域法等。

3.1 连续系统分析

3.1.1 卷积积分法求解系统的响应

对LTI连续系统，输入信号为x(t)，输出信号为y(t)，系统的冲激响应为h(t)，则输入、输出之间的关系为：

时域：y(t)=x(t)*h(t) 频域：Y(jω)=X(jω) H(jω)

其中线性卷积的计算引入了DFT分析连续信号的频谱。

例1：图（A）和（B）分别表示一个LTI连续系统的输入x(t)和单位冲激响应h(t)，试利用DFT的快速算法计算该系统的零状态响应。

x(t) h(t) （A） （B） 1 1

3 4

实验一 MATLAB仿真软件的基本操作命令和使用方法

实验内容

1、帮助命令

使用 help 命令，查找 sqrt（开方）函数的使用方法；

2、MATLAB命令窗口

（1）在MATLAB命令窗口直接输入命令行计算y1?2sin(0.5?)1?3的值；

（2）求多项式 p(x) = x3 + 2x+ 4的根；

3、矩阵运算 （1）矩阵的乘法

已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B

（2）矩阵的行列式

已知 A=[1 2 3;4 5 6;7 8 9]， 求A

（3）矩阵的转置及共轭转置 已知A=[1 2 3;4 5 6;7 8 9]，求A'

已知 B=[5+i,2-i,1;6*i,4,9-i], 求 B.' , B'

（4）特征值、特征向量、特征多项式

已知 A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵 A的特征值、特征向量、特征多项式；

（5）使用冒号选出指定元素

已知： A=[1 2 3;4 5 6;7 8 9]；求 A 中第 3 列前 2 个元素；A 中所有列第 2，3 行的元素；

4、Matlab 基本编程方法

（1） 编写命令

实验一 MATLAB仿真软件的基本操作命令和使用方法

实验内容

1、帮助命令

使用 help 命令，查找 sqrt（开方）函数的使用方法；

2、MATLAB命令窗口

（1）在MATLAB命令窗口直接输入命令行计算y1?2sin(0.5?)1?3的值；

（2）求多项式 p(x) = x3 + 2x+ 4的根；

3、矩阵运算 （1）矩阵的乘法

已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B

（2）矩阵的行列式

已知 A=[1 2 3;4 5 6;7 8 9]， 求A

（3）矩阵的转置及共轭转置 已知A=[1 2 3;4 5 6;7 8 9]，求A'

已知 B=[5+i,2-i,1;6*i,4,9-i], 求 B.' , B'

（4）特征值、特征向量、特征多项式

已知 A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵 A的特征值、特征向量、特征多项式；

（5）使用冒号选出指定元素

已知： A=[1 2 3;4 5 6;7 8 9]；求 A 中第 3 列前 2 个元素；A 中所有列第 2，3 行的元素；

4、Matlab 基本编程方法

（1） 编写命令

太原理工大学

数字信号处理课程 实验报告

专业班级

学 号2013000000 姓 名 XXX 指导教师XXX

实验一: 系统响应及系统稳定性

1.实验目的

（1）掌握 求系统响应的方法。 （2）掌握时域离散系统的时域特性。 （3）分析、观察及检验系统的稳定性。

2.实验原理与方法

在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。在计算机上适合用递推法求差分方程的解，最简单的方法是采用MATLAB语言的工具箱函数filter函数。也可以用MATLAB语言的工具箱函数conv函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。

系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。

系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。或者系统的单位脉冲响应满足绝对可和的条件。系统的稳定性由其差分方程的系数决定。

实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否都是

信号专业资料

一、测试系统基本组成及类型1 基本组成传感器：拾取信号， 将振动信号变为电荷信号 电荷放大器：放大信号， 将电荷信号变为电压信号

A/D卡：转换信号， 将电压信号变为数字信号 计算机：分析信号 进行诊断

信号专业资料

2.测量量

– 振动位移– 振动速度 – 振动加速度

3、振动的位移、速度、加速度指标 位移： x(t)=A*sin(2πft+φ) A-振幅，大小 f-频率，快慢 速度： v(t)=dx(t)/dt=fAcos(2πft+φ) 加速度：a(t)=dv(t)/dt=-f2Asin(2πft+φ) 三者间频率 f 不变，最大幅值呈 f 倍递增 位移xmax=A 速度 vmax=fA 加速度 amax=f2A

信号专业资料

安装方法

压电加速度计的安装方法 1.钢螺栓 2.绝缘螺栓和云母垫 3.磁铁吸附 4.胶合 5.蜡和橡胶泥粘附 6手持探头

信号专业资料

二、采样定理为保证采样后信号能真实地保留原始模拟信号信 息，信号采样频率必须至少为原信号中最高频率成 分的2倍，亦称仙农(香农)定理。

f s 2 f max注意：满足采样定理时，只保证不发生混叠，而不 能保证采样信号能真实地反映原信号 x(t) 。工程实际 中采样频率通常大

数字信号处理实验

报告

实验一 信号、系统及系统响应

一．实验目的

（1） 熟悉连续信号理想采样前后的频谱变化关系，加深对时域采样定理的理解； （2） 熟悉时域离散系统的时域特性；

（3） 利用卷积方法观察分析系统的时域特性；

（4） 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离

散信号及系统响应进行频域分析。

二．实验原理与方法

采样时连续信号数字处理的第一个关键环节。对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生变化以及信号信息不丢失的条件，而且可以加深对傅里叶变换、Z变换和序列傅里叶变换之间关系式的理解。

对一个连续信号

xa(t)

进行理想采样的过程可用下式表示：

其中

a(t) xa(t)p(t)x

a(t)x

为

xa(t)

的理想采样，p(t)为周期脉冲，即

p(t)

m

(t nT)

jwk

N 1n 0

a(t)x

X(e

的傅里叶变换为

N 1n 0

) x(m)e jwkn

X(e

其中，

jwk

) x(m)e jwknwk

2

k

M ，k=0,1, M-1

时域离散线性非时变系统的输入输出关系为

y(n) x(n)*h(n)

m

x(m)h(n m)

jwjwjwY(e) X(e)H(e) 卷积运算也可在频域实现

三．实