Yule-Walker方程-实验报告.docVIP

评分 大理大学实验报告 课程名称 生物医学信号处理 实验名称 Yule-Walker方程 专业班级 姓 名 羽卒兰cl 学 号 实验日期 2016年5月27日 实验地点 2015—2016学年度第 3 学期 实验目的 学习求解Yule-Walker方程，建立随机信号的AR模型。 二、实验环境 1、硬件配置：Intel（R） Core(TM) i5-4210U CPU @1.7GHz 1.7GHz 安装内存（RAM）：4.00GB 系统类型：64位操作系统 2、软件环境：MATLAB R2013b软件 三、实验内容 (包括本实验要完成的实验问题及需要的相关知识简单概述) 编写求解Yule-Walker方程的程序，并对实际生理信号（例如心电、脑电）建立AR模型。 对同一数据，使用Matlab信号处理工具箱自带函数aryule计算相同阶数AR模型系数，检验程序是否正确。 用伪随机序列（白噪声）驱动AR模型，观察输出是否和真实心电、脑电信号相似，对比真实信号和仿真信号的功率谱。 四、实验结果和分析 （包括实验原理、数据的准备、运行过程分析、源程序（代码）、图形图象界面等） 实验原理 随机信号可以看作是由当前激励白噪声w(n)以及若干次以往信号x(n-k)的线性组合产生，即所谓自回归模型（AR模型） 模型参数满足Yule-Walker方程 矩阵形式 求解Yule-Walker方程，就可以得到AR模型系数 当模型阶次较大时，直接用矩阵运算求解的计算量大，不利于实时运算。利用系数矩阵的特性，人们提出了如L-D算法等快速算法。 源程序： clear; clc; M=1024; load ecgdata;x = ecgdata(1:1024); %导入心电信号的数据 %load eegdata;x = eegdata(1:M); %导入脑电信号的数据 %load icpdata;x = icpdata(1:1024); %导入颅内压信号的数据 %load respdata;x = respdata (1:1024);%导入个呼吸信号的数据 p=1:60; %P的取值范围 Sw=zeros(1,length(p)); %创建一个1行列长为length(p)的Sw零矩阵或数组 E=zeros(1,length(p)); %创建一个1行列长为length(p)的E零矩阵或数组 FPE=zeros(1,length(p)); %创建一个1行列长为length(p)的FPE零矩阵或数组 for i=1:60 %尝试改变模型阶数，观察效果 Rxx = xcorr(x,biased);%估计随机过程中的互相关序列，biased为有偏的互相关函数估计； Rtemp = zeros(1,i); %创建一个i行1列的Rtemp零矩阵或数组 Rl = zeros(i,1); %创建一个i行1列的Rl零矩阵或数组 for k = 1:length(Rtemp) %for循环从1到length(Rtemp) Rtemp(k) = Rxx(length(x)-1+k); %取Rxx（0）到Rxx（p-1）赋值 Rl(k) = Rxx(length(x)+k); %取Rxx（1）到Rxx（p）赋值 end Rs = toeplitz(Rtemp); %生成自相关系数矩阵(Toeplitz型) A = -inv(Rs)*Rl; %AR模型系数估计 Sw(i)= [Rtemp(1),Rl]*[1;A]; %白噪声方差估计 % 采用mala