随机信号的数字特征分析实验报告.docx

《随机信号的数字特征分析实验报告.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、评分大理大学实验报告课程名称生物医学信号处理实验名称随机信号的数字特征分析专业班级姓名羽卒兰cl学号实验日期实验地点2015—2016学年度第3学期.一、实验目的了解随机信号的特征。掌握随机信号的数字特征分析算法。二、实验环境1、硬件配置：Intel（R）Core(TM)i5-4210UCPU@1.7GHz1.7GHz安装内存（RAM）：4.00GB系统类型：64位操作系统2、软件环境：MATLABR2013b软件三、实验内容用Matlab编制程序，分析信号的数字特征，包括均值、方差、均方值、协方差。可以使用Matlab自带函数。观察信号的直方图，粗略估计其概率分布。信号1：利用Mat

2、lab中的伪随机序列产生函数randn()产生的长1024点的序列；信号2：实际采集的生物医学信号（脑电，心电等）。四、实验结果与分析实验原理：1.对于平稳各态遍历随机过程，可以用单一样本函数的时间平均代替集总平均，即通过测量过程的单一样本来估计信号的统计特征量。样本均值：1n?ximxni1样本均方值：Exn21nnxi2i1样本方差：?x21nnxim?x2i1样本自相关函数：1nRxx(m)xiximni11n样本自协方差：Cxx(m)(ximx)(ximmx)ni12.平稳随机过程:如果随机信号的概率特性不随时间变化而变化，则称为平稳随机过程，否则称为非平稳随机过程。（1）

3、.一阶平稳过程：信号的平均值与t无关；（2）.二阶平稳过程：①信号的平均值与t无关；②信号的均方值与t无关；③信号的协方差只是时间间隔的函数，而与时间的原点无关；（3）.我们至少把二阶平稳过程叫做准平稳过程或广义平稳过程。...源代码：%选择信号类型并设定参数，产生信号x(n)clear;clc;%clear清空内存变量,clc只清空工作区内的内容disp('请选择信号');%显示：请选择信号disp('1----伪随机序列randn()');%显示：1----伪随机序列randn()disp('2----实际测量的心电信号');%显示：2----实际测量的心电信号disp('3---

4、-实际测量的脑电信号');%显示：3----实际测量的脑电信号disp('4----实际测量的呼吸信号');%显示：4----实际测量的呼吸信号disp('5----实际测量的颅内压信号');%显示：5----实际测量的颅内压信号b=input('信号：');%在屏幕上显示：信号：，等待输入内容switchb%输入序号，产生相应信号switch条件语句case1%情况1L=input('每段数据长度L');%每段数据长度L，作为一个字符（串），输入后回车赋值N=input('数据共多少段N');%数据共多少段N，作为一个字符（串），输入后回车赋值x=randn(1,L*N);%

5、randn是生成随机数的函数，randn(1,L*N)表示生成L*N的矩阵，矩阵的每个元素都是随机数。case2%情况2loadecgdata;%导入ECG心电信号display(['数据总长度',num2str(length(ecgdata)),'点']);%把ECG心电信号长度数值转换成字符串，转换后可以使用fprintf或disp函数进行输出。L=input('每段数据长度L');%输入每段数据长度LN=input('数据共多少段N');%输入段数Nx=ecgdata(1:(N*L));%表示生成N*L的矩阵的ECG心电信号case3%情况3loadeegdata;%导入

6、EEG脑电信号display(['数据总长度',num2str(length(eegdata)),'点']);%把EEG脑电信号长度数值换成字符串，转换后可以使用fprintf或disp函数进行输出。L=input('每段数据长度L');%输入每段数据长度LN=input('数据共多少段N');%输入段数Nx=eegdata(1:(N*L));%表示生成N*L的矩阵的EEG脑电信号case4%情况4loadicpdata;%导入实际测量的呼吸信号display(['数据总长度',num2str(length(icpdata)),'点']);%把ECG心电信号长度数值转换成字符串

7、（转换后可以使用fprintf或disp函数进行输出）L=input('每段数据长度L');%输入每段数据长度LN=input('数据共多少段N');%输入段数Nx=icpdata(1:(N*L));%b=5时，产生N*L的矩阵的呼吸信号case5%情况5loadrespdata;%导入实际测量的颅内压信号display(['数据总长度',num2str(length(respdata)),'点']);%把EEG脑电信号长度数值转...