随机信号的数字特征分析教程.docx

评分大理大学实验报告课程名称 生物医学信号处理 实验名称随机信号的数字特征分析 专业班级姓名 学号实验日期 实验地点2015—2016学年度第3学期实验目的1、了解随机信号的特征。2、掌握随机信号的数字特征分析算法。实验环境1、硬件配置：处理器：AMD A10-5750M APU with Radeon(tm) Graphics 2.50GHz 安装内存：（RAM）4.00GB 系统类型：64位操作系统，基于x64位处理器 2、软件环境：Matlab R2012b实验原理对于平稳各态遍历随机过程，可以用单一样本函数的时间平均代替集总平均，即通过测量过程的单一样本来估计信号的统计特征量。样本均值：样本均方值：样本方差：实验内容用Matlab编制程序，分析信号的数字特征，包括均值、方差、均方值、协方差。可以使用Matlab自带函数。观察信号的直方图，粗略估计其概率分布。信号1：利用Matlab中的伪随机序列产生函数randn()产生的长1024点的序列；信号2：实际采集的生物医学信号（脑电，心电，颅内压，心电等）。五、实验结果与分析实验程序：disp(请选择信号); %字符串原样显示disp(1 ---- 伪随机序列randn());disp(2 ---- 实际测量的心电信号);disp(3 ---- 实际测量的脑电信号);disp(4---- 实际测量的颅内压信号);disp(5 ---- 实际测量的呼吸信号);b = input(信号：); %用 input函数输入对应信号switch b % 输入序号，用switch 函数有条件地执行语句的选项。每个选项用case 语句表示，产生相应信号。case 1 %随机信号的case 语句 L = input(每段数据长度 L \n); %用 input函数输入每段数据长度 L N = input(数据共多少段 N \n); %用 input函数输入数据段 数N x = randn(1, L*N); %随机产生 L*N个数据case 2 %实际测量的心电信号的case 语句 load ecgdata; %载入心电信号的数据 display([数据总长度,num2str(length(ecgdata)),点]); L = input(每段数据长度 L \n); %用 input函数输入每段数据长度 L N = input(数据共多少段 N \n); %用 input函数输入数据段 数N x = ecgdata (1:(N*L)); %产生 L*N个心电数据case 3 %实际测量的脑电信号的case 语句 load eegdata; %载入脑电信号的数据 display([数据总长度,num2str(length(eegdata)),点]); L = input(每段数据长度 L \n); %用 input函数输入每段数据长度 L N = input(数据共多少段 N \n); %用 input函数输入数据段 数N x = eegdata (1:(N*L)); %产生 L*N个脑电数据case 4 %实际测量的颅内压信号的case 语句 load icpdata; %载入颅内压信号的数据 display([数据总长度,num2str(length(icpdata)),点]); L = input(每段数据长度 L \n); %用 input函数输入每段数据长度 L N = input(数据共多少段 N \n); %用 input函数输入数据段 数N x =icpdata (1:(N*L)