清华大学自动化系 张贤达 zxd-dau@tsinghua.edu.cn 电话:62794875 第三章 现代谱估计 清华大学自动化系 张贤达 zxd-dau@tsinghua.edu.cn 电话:62794875.

Slides:



Advertisements
Similar presentations
一、 一阶线性微分方程及其解法 二、 一阶线性微分方程的简单应用 三、 小结及作业 §6.2 一阶线性微分方程.
Advertisements

第五节 函数的微分 一、微分的定义 二、微分的几何意义 三、基本初等函数的微分公式与微分运算 法则 四、微分形式不变性 五、微分在近似计算中的应用 六、小结.
第九章 常微分方程数值解法 §1 、引言. 微分方程的数值解:设方程问题的解 y(x) 的存在区间是 [a,b] ,令 a= x 0 < x 1
第八章 第四节 机动 目录 上页 下页 返回 结束 一个方程所确定的隐函数 及其导数 隐函数的微分法.
2.6 隐函数微分法 第二章 第二章 二、高阶导数 一、隐式定义的函数 三、可微函数的有理幂. 一、隐函数的导数 若由方程 可确定 y 是 x 的函数, 由 表示的函数, 称为显函数. 例如, 可确定显函数 可确定 y 是 x 的函数, 但此隐函数不能显化. 函数为隐函数. 则称此 隐函数求导方法.
2012江苏历史高考 重点与热点考点分析与复习.
第六章 动态模型.
信号与系统 第三章 傅里叶变换 东北大学 2017/2/27.
«地学建模» 之 “随机时间序列分析模型”.
非线性时间序列模型 一般非线性时间序列模型介绍 条件异方差模型 上海财经大学 统计与管理学院.
西方行政学说史 导论:西方行政学的产生与发展历程.
岩層中的奧秘與寶藏.
第五章 二次型. 第五章 二次型 知识点1---二次型及其矩阵表示 二次型的基本概念 1. 线性变换与合同矩阵 2.
第3节 二次型与二次型的化简 一、二次型的定义 二、二次型的化简(矩阵的合同) 下页.
一、能线性化的多元非线性回归 二、多元多项式回归(线性化)
一、二阶行列式的引入 用消元法解二元线性方程组. 一、二阶行列式的引入 用消元法解二元线性方程组.
第三章 函数逼近 — 最佳平方逼近.
第7章 数字信号处理中的有效字长效应.
第2章 时域离散信号和系统的频域分析 教学内容包括: 序列的傅立叶变换定义及性质 Z变换的定义与收敛域 利用z变换分析信号和系统的频域特性.
第二章 平稳时间序列模型(单变量) 选择单变量时间序列的原因
恰当方程(全微分方程) 一、概念 二、全微分方程的解法.
不确定度的传递与合成 间接测量结果不确定度的评估
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
2-7、函数的微分 教学要求 教学要点.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
ARIMA模型在电量预测中的应用 蔡跳
Signals and Systems Lecture 28
第2章 Z变换 Z变换的定义与收敛域 Z反变换 系统的稳定性和H(z) 系统函数.
Matlab 中IIR数字滤波器设计相关函数
7 平稳时间序列预测法 7.1 概述 7.2 时间序列的自相关分析 7.3 单位根检验和协整检验 7.4 ARMA模型的建模 回总目录.
全国高校数学微课程教学设计竞赛 知识点名称: 导数的定义.
计算机数学基础 主讲老师: 邓辉文.
现代信号处理 何继爱.
第七章 参数估计 7.3 参数的区间估计.
第4章 非线性规划 4.5 约束最优化方法 2019/4/6 山东大学 软件学院.
第四章 功率谱估计.
抽样和抽样分布 基本计算 Sampling & Sampling distribution
Partial Differential Equations §2 Separation of variables
第二十二章 曲面积分 §1 第一型曲面积分 §2 第二型曲面积分 §3 高斯公式与斯托克斯公式.
Module_4_Unit_11_ppt Unit11:系统动态特性和闭环频率特性的关系 东北大学《自动控制原理》课程组.
長虹虹頂新建工程 中鹿營造/ 宏林營造廠- 聯合承攬
§2.4 零输入响应和零状态响应 零输入响应 零状态响应 对系统线性的进一步认识.
第十章 双线性型 Bilinear Form 厦门大学数学科学学院 网址: gdjpkc.xmu.edu.cn
第4章 Excel电子表格制作软件 4.4 函数(一).
实验一 熟悉MATLAB环境 常用离散时间信号的仿真.
定理21.9(可满足性定理)设A是P(Y)的协调子集,则存在P(Y)的解释域U和项解释,使得赋值函数v(A){1}。
第16讲 相似矩阵与方阵的对角化 主要内容: 1.相似矩阵 2. 方阵的对角化.
§8.3 不变因子 一、行列式因子 二、不变因子.
§6.7 子空间的直和 一、直和的定义 二、直和的判定 三、多个子空间的直和.
相关与回归 非确定关系 在宏观上存在关系,但并未精确到可以用函数关系来表达。青少年身高与年龄,体重与体表面积 非确定关系:
2019/5/11 实验三 线性相位FIR滤波器的特性 05:31:30.
第三章 函数的微分学 第二节 导数的四则运算法则 一、导数的四则运算 二、偏导数的求法.
第一部分:概率 产生随机样本:对分布采样 均匀分布 其他分布 伪随机数 很多统计软件包中都有此工具 如在Matlab中:rand
建模常见问题MATLAB求解  .
第三章 随机信号的功率谱估计 郑宝玉.
第15讲 特征值与特征向量的性质 主要内容:特征值与特征向量的性质.
线 性 代 数 厦门大学线性代数教学组 2019年5月12日4时19分 / 45.
概率论与数理统计B.
§2 方阵的特征值与特征向量.
§7.3 离散时间系统的数学 模型—差分方程 线性时不变离散系统 由微分方程导出差分方程 由系统框图写差分方程 差分方程的特点.
现代通信原理 吉林大学远程教育学院.
第三节 函数的微分 3.1 微分的概念 3.2 微分的计算 3.3 微分的应用.
第三章 从概率分布函数的抽样 (Sampling from Probability Distribution Functions)
第三章 平稳时间序列分析.
教学大纲(甲型,54学时 ) 教学大纲(乙型, 36学时 )
三角 三角 三角 函数 余弦函数的图象和性质.
数学模型实验课(二) 最小二乘法与直线拟合.
第三章 线性方程组 §4 n维向量及其线性相关性(续7)
§2 自由代数 定义19.7:设X是集合,G是一个T-代数,为X到G的函数,若对每个T-代数A和X到A的函数,都存在唯一的G到A的同态映射,使得=,则称G(更严格的说是(G,))是生成集X上的自由T-代数。X中的元素称为生成元。 A变, 变 变, 也变 对给定的 和A,是唯一的.
第三章 高阶谱估计 3.1 累积量及高阶谱 3.2 高阶谱估计 3.3 有色噪声背景下的频率估计 3.4 高阶谱的应用 2019/12/16.
Presentation transcript:

清华大学自动化系 张贤达 zxd-dau@tsinghua.edu.cn 电话:62794875 第三章 现代谱估计 清华大学自动化系 张贤达 zxd-dau@tsinghua.edu.cn 电话:62794875

经典谱估计 样本 直接法 间接法 假设已零均值化, 周期函数

周期图法 有偏估计,平滑性差 加窗函数 功率谱曲线平滑,但分辨率下降 要提高分辨率,使用参数化的谱估计! 经典谱估计:使用FFT的谱估计 数据窗 谱窗 要提高分辨率,使用参数化的谱估计! 经典谱估计:使用FFT的谱估计 现代谱估计:参数化谱估计

3.1 ARMA谱估计与系统辨识 平稳ARMA过程 为离散白噪声,则称 为ARMA过程。 离散随机过程 服从线性差分方程: 为离散白噪声,则称 为ARMA过程。 自回归(autoregressive)—滑动平均(moving average)过程 AR阶数 MA阶数 AR参数 MA参数

ARMA模型描述的线性时不变(LTI)系统 传递函数:

(1)冲激响应系数必须绝对可求和: (系统稳定) (2)A(z)和B(z)无公共因子(p,q唯一) (3)系统是物理可实现的(因果系统) 满足ARMA模型的条件: (1)冲激响应系数必须绝对可求和: (系统稳定) (2)A(z)和B(z)无公共因子(p,q唯一) (3)系统是物理可实现的(因果系统) 极点的作用:决定系统的稳定性和因果性 即极点不在单位圆上 因果性:称x(n)是e(n)的因果函数,若 即因果系统要求极点在单位圆以内,A(z)的根|z|<1 零点部分 极点部分

零点的作用:决定系统的可逆性,即 是否存在。 可逆性:称e(n)是x(n)的可逆函数,若 (1)存在序列 ,并满足 ——可逆系统的稳定性 (2) ——可逆性条件

ARMA过程的功率谱密度 则功率谱 其中

ARMA功率谱估计的两种线性方法 Cadzow谱估计子 又 其中

则 两边同乘 ,比较系数得 所以,Cadzow谱估计子的关键:估计AR阶数p和AR参数

Kaveh谱估计子 非线性方程,MA参数辨识 (Newton-Raphson迭代)

协方差函数的Fourier变换 Kaveh谱估计子:

ARMA功率谱密度的特例 特例一:MA过程 有限冲激响应(FIR)系统

特例二:AR过程 中含有 的无数多项 无限冲激响应(IIR)系统 白噪声中的AR过程: ARMA(p,p)过程

特例三:完全可预测过程 线 谱 加性白噪声中的可预测过程: 特殊的ARMA

所以: 白噪声中的AR过程 = ARMA过程 白噪声中的可预测过程 = 特殊的ARMA过程

修正Yule-Walker方程 等价 高斯白噪

BBR公式:

修正Yule-Walker方程(MYW方程)

定理(AR参数的可辨识性): 若A(z)和B(z)无可对消公共因子,且 ,则AR参数 可由p个修正Yule-Walker方程唯一确定或辨识。

若构造: 使得 ,则

主奇异值:p个大的奇异值(p个信号分量的能量) AR阶数确定的奇异值分解方法 奇异值分解(SVD): 酉矩阵: 主奇异值:p个大的奇异值(p个信号分量的能量) 次奇异值:其它小奇异值(扰动或误差的能量)

信号与噪声的分离: 准则一:归一化比值 准则二:使用归一化奇异值 若阈值=0.995,v(k)>阈值的最小整数k定为矩阵A的“有效秩”。 准则二:使用归一化奇异值 <某个很小的阈值(0.05)的最小整数k定为有效秩。

AR阶数确定的信息量准则法 最终预报误差方法(FPE, Finite Prediction Error): FPE准则选择使FPE(p,q)最小,作为AR模型的阶数。 AIC(Akaike’s Information Criterion) “信息量准则”遵循“吝啬原则”

扩展阶MYW方程 若 ,则 选

总体最小二乘TLS: Total Least Squares 扰动矩阵 总体最小二乘TLS: Total Least Squares 思想:寻求一个解z,使得

定义代价函数

方法2:只包含p个参数(主要因素) :用秩为p的矩阵对B的最佳逼近 令 ,则

构造代价函数

由于存在误差

AR定阶与AR参数估计的SVD-TLS算法: 步骤2:确定 的有效秩p,给出AR阶数估计值 步骤3:计算 步骤4:计算