.001:0.5 Sx0 M1000 for k1:M-1 SxSxRkexp-jwkRkexpjwk由自相关函数的傅里叶变换求 功率谱P89页4.1.2 end Sxlog10absSx plotwabsSx title间接法谱估计 grid on 通信工程 孙大江 通信工程 孙大江 2Levinson法程序 clc clear all fs100 设采样频率为100 N150 数据长度改变数据长度会导致分辨率的变化 n0:0.01:2 变化范围 xsqrt20sin2pi0.2nsqrt10sin2pi0.3nrandnsizen for m1:N 接下来的循环计算自相关函数的无偏估计仍使用P89也方法 Rm0 end for m1:N s0 for n1:N-m-1 ssxnxnm-1 end Rms/N end for M1:N-5 接下来的循环定义初值----阶数M要小于数据长度N aMM0 FPEM0 PM0 end 接下来初始化由Levinson公式 a11-R2/R1 计算一阶模型参数a11见P66页公式3.3.15 P1R1a11R2 计算一阶模型的误差功率P1p66页底部误差功率PG2 FPE1P1N11/N-1-1 计算一阶模型的最终预测误差准则函数FPE1 sum0 接下来进行递推 for M2:N-5 for k1:M-1 sumsumaM-1kRM-k1 end aMM-RM1sum/PM-1 P67页计算反射系数公式3.3.16a for k1:M-1 aMkaM-1kaMMaM-1M-k 由反射系数通过levinson算法递推求系统函数 系数a3.3.16b end PM1-absakk2PM-1 由反射系数求预测误差功率即G23.3.16c FPEMPMNM1/N-M-1 计算模型的最终预测误差准则函数FPEM sum0 sum归零 end 确定阶数:求出使FPE最小时的阶数M1 minFPE1 确定比较对象为FPE1 for M2:N-5 依次比较 if FPEM
上一篇:
浅析朝鲜式社会主义 政治学国际关系学等专业毕业论文 原创论文
下一篇:
汝康皮肤护理液医治烧伤的可行性和效果分析