随机信号的数字特征实验资料.doc

评分 大理大学实验报告 课程名称 实验名称 随机信号的数字特征分析 专业班级 姓 名 羽卒兰cl 学 号 实验日期 实验地点 2015—2016学年度第 学期 实验目的 二、实验环境 1、硬件配置：IntelR） Core(TM) i5-4210U CPU @1.7GHz 1.7GHz 安装内存（RAM）：4.00GB 系统类型：64位操作系统 2、软件环境： 四、实验结果与分析 实验原理： 1.对于平稳各态遍历随机过程，可以用单一样本函数的时间平均代替集总平均，即通过测量过程的单一样本来估计信号的统计特征量。 样本均值： 样本均方值： 样本方差： 样本自相关函数： 样本自协方差： 2.平稳随机过程:如果随机信号的概率特性不随时间变化而变化，则称为平稳随机过程，否则称为非平稳随机过程。 （1）.一阶平稳过程：信号的平均值与t无关； （2）.二阶平稳过程： ①信号的平均值与t无关； ②信号的均方值与t无关； ③信号的协方差只是时间间隔的函数，而与时间的原点无关； （3）.我们至少把二阶平稳过程叫做准平稳过程或广义平稳过程。 源代码： % 选择信号类型并设定参数，产生信号x(n) clear; clc;%clear清空内存变量,clc只清空工作区内的内容 disp(请选择信号);%显示：请选择信号 disp(1 ---- 伪随机序列randn()); %显示：1 ---- 伪随机序列randn() disp(2 ---- 实际测量的心电信号); %显示：2 ---- 实际测量的心电信号 disp(3 ---- 实际测量的脑电信号); %显示：3 ---- 实际测量的脑电信号 disp(4 ---- 实际测量的呼吸信号); %显示：4 ---- 实际测量的呼吸信号 disp(5 ---- 实际测量的颅内压信号);%显示：5 ---- 实际测量的颅内压信号 b = input(信号：); %在屏幕上显示：信号：，等待输入内容 switch b % 输入序号，产生相应信号switch条件语句 case 1 %情况1 L = input(每段数据长度 L \n);%每段数据长度 L，作为一个字符（串），输入后回车赋值 N = input(数据共多少段 N \n);%数据共多少段 N，作为一个字符（串），输入后回车赋值 x = randn(1, L*N);%randn是生成随机数的函数，randn(1,L*N)表示生成L*N的矩阵，矩阵 的每个元素都是随机数。 case 2 %情况2 load ecgdata;%导入ECG心电信号 display([数据总长度,num2str(length(ecgdata)),点]);%把ECG心电信号长度数值转换 成字符串，转换后可以使用fprintf或disp函数进行输出。 L = input(每段数据长度 L \n); %输入每段数据长度L N = input(数据共多少段 N \n); %输入段数N x = ecgdata (1:(N*L)); %表示生成N*L的矩阵的ECG心电信号 case 3%情况3 load eegdata; %导入EEG脑电信号 display([数据总长度,num2str(length(eegdata)),点]);%把EEG脑电信号长度数值 换成字符串，转换后可以使用fprintf或disp函数进行输出。 L = input(每段数据长度 L \n); %输入每段数据长度L N = input(数据共多少段 N \n); %输入段数N x = eegdata (1:(N*L));%表示生成N*L的矩阵的EEG脑电信号 case 4 %情况4 load icpdata; %导入实际测量的呼吸信号 display([数据总长度,num2str(length(icpdata)),点]);%